Tổng quan nghiên cứu

Lưu vực sông Thao là một trong những hệ thống thủy văn xuyên biên giới quan trọng hàng đầu tại miền Bắc Việt Nam, với tổng diện tích lưu vực tính đến trạm Việt Trì đạt khoảng 51.800 km². Trong đó, phần diện tích thượng nguồn nằm trên lãnh thổ Trung Quốc chiếm tới 39.800 km² (tương đương 76,8%), còn phần hạ lưu thuộc địa phận các tỉnh Lào Cai, Yên Bái và Phú Thọ của Việt Nam chỉ chiếm khoảng 12.000 km² (tương đương 23,2%). Với tổng lượng dòng chảy trung bình nhiều năm đạt khoảng 26,7 tỷ m³, tương ứng lưu lượng bình quân 783 m³/s, sông Thao đóng vai trò huyết mạch đối với nguồn cung cấp nước cho sản xuất nông nghiệp, công nghiệp và dân sinh của toàn bộ vùng đồng bằng sông Hồng.

Tuy nhiên, công tác quản lý tài nguyên nước tại lưu vực sông Thao đang đối mặt với thách thức kép. Thứ nhất, phía thượng nguồn Trung Quốc đã vận hành liên tiếp hơn 9 công trình đập dâng và thủy điện bậc thang quy mô lớn, tiêu biểu như đập Nanshan công suất 120 MW (dung tích 220 triệu m³) và đập Madushan công suất 260 MW (dung tích 551 triệu m³), nhưng nguồn số liệu quan trắc khí tượng thủy văn xuyên biên giới hầu như không được chia sẻ đầy đủ. Thứ hai, tác động tiêu cực từ biến đổi khí hậu toàn cầu đang làm thay đổi căn bản chế độ mưa và dòng chảy mặt.

Mục tiêu cốt lõi của nghiên cứu là thiết lập và ứng dụng thành công bộ mô hình khí tượng thủy văn tích hợp WEHY-HCM (Watershed Environmental Hydrology - Regional Hydroclimate Model) để khôi phục chuỗi số liệu lịch sử và mô phỏng dự báo biến động nguồn nước lưu vực sông Thao giai đoạn 2019–2100 theo các kịch bản nồng độ khí nhà kính RCP 4.5 và RCP 8.5. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng chuẩn xác giúp các cơ quan quản lý nhà nước chủ động lập quy hoạch an ninh nguồn nước và giảm thiểu rủi ro thiên tai.

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 dựa trên nền tảng lý thuyết mô hình hóa khí tượng thủy văn động lực phân bố không gian, kết hợp giữa mô hình thủy văn vật lý lưu vực WEHY (Watershed Environmental Hydrology) và mô hình khí tượng khu vực động lực học RegHCM. Khung lý thuyết này giải quyết triệt để các hạn chế của các mô hình thủy văn truyền thống dạng thông số tập trung vốn phụ thuộc nặng nề vào chuỗi số liệu quan trắc thực địa tại chỗ.

Mô hình WEHY-HCM vận hành dựa trên 3 khái niệm khoa học trọng tâm:

  1. Cơ chế truyền dẫn động lực học khí quyển: Tái hiện chi tiết quá trình tương tác bức xạ mặt trời, dòng năng lượng nhiệt và hoàn lưu khí quyển tầng thấp kết hợp địa hình núi cao hiểm trở.
  2. Cân bằng ẩm và quá trình thủy văn mặt đệm: Mô phỏng chính xác các quá trình vật lý gồm lượng mưa rơi, chặn giữ của tán cây rừng, bốc thoát hơi nước bề mặt, quá trình thấm vào các tầng đất và hình thành dòng chảy tràn bề mặt lẫn dòng chảy ngầm.
  3. Kịch bản biến đổi khí hậu toàn cầu: Áp dụng các kịch bản đường phân bố nồng độ phát thải khí nhà kính đại diện gồm RCP 4.5 (phát thải trung bình thấp) và RCP 8.5 (phát thải cao) được công bố bởi Ban Liên chính phủ về Biến đổi Khí hậu (IPCC).

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

Nghiên cứu áp dụng quy trình xử lý dữ liệu viễn thám hiện đại kết hợp phương pháp chi tiết hóa động lực học (Dynamical Downscaling) qua hệ thống 3 miền tính toán lồng ghép (Nested Domains) từ độ phân giải 81 km (miền ngoài), 27 km (miền trung gian) xuống còn 9 km (miền chi tiết bao phủ toàn bộ lưu vực sông Thao).

Nguồn dữ liệu đầu vào bao gồm:

  • Dữ liệu tái phân tích khí quyển toàn cầu ERA-20C từ Trung tâm Dự báo Thời tiết Hạn vừa Châu Âu (ECMWF) giai đoạn 1900–2014 với độ phân giải ban đầu 125 km và bước thời gian 3 giờ.
  • Dữ liệu mô hình số độ cao ASTER DEM độ phân giải 30 m nhằm mô tả chuẩn xác mạng lưới địa hình và trắc diện lòng dẫn lưu vực.
  • Dữ liệu chỉ số diện tích lá cây MODIS LAI độ phân giải không gian 500 m cập nhật theo chu kỳ 4 ngày.
  • Bản đồ thổ nhưỡng toàn cầu SoilGrids độ phân giải 250 m phân định chi tiết các lớp đất cát, thịt và sét.

Phương pháp chọn mẫu và hiệu chỉnh được thực hiện thông qua việc phân chia chuỗi dữ liệu thủy văn thực đo tại các trạm kiểm soát chính trên sông Thao gồm Lào Cai và Yên Bái: giai đoạn hiệu chỉnh mô hình từ năm 1985 đến 1998 và giai đoạn kiểm định từ năm 2000 đến 2010. Lý do lựa chọn mô hình WEHY-HCM bắt nguồn từ việc mô hình tính toán trực tiếp các hằng số vật lý từ thuộc tính mặt đệm (độ dốc, thổ nhưỡng, thảm phủ), cho phép mô phỏng chính xác chế độ dòng chảy tại các lưu vực xuyên biên giới thiếu dữ liệu mà không cần điều chỉnh thông số thống kê nhân tạo. Chất lượng mô hình được đánh giá qua chỉ số hiệu quả Nash-Sutcliffe (NSE) đạt mức tương quan cao trên 0,75.

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

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

Quá trình mô phỏng và phân tích diễn biến khí tượng thủy văn trên lưu vực sông Thao đã chỉ ra 4 phát hiện quan trọng:

  1. Chế độ phân bổ dòng chảy hiện trạng có sự phân hóa mùa sâu sắc: Tổng lượng dòng chảy năm tại trạm Yên Bái đạt trung bình 24,7 tỷ m³. Dòng chảy mùa lũ tập trung chủ yếu từ tháng VI đến tháng X, chiếm từ 75,5% đến 80,4% tổng lượng nước cả năm (tháng cao điểm như tháng VII hoặc VIII có thể chiếm tới 25,7% đến 30,0% tổng lượng dòng chảy năm). Ngược lại, mùa cạn kéo dài 7 tháng (từ tháng XI đến tháng V năm sau) chỉ đóng góp từ 19,6% đến 24,5% tổng lượng nước.
  2. Xu thế gia tăng lượng mưa và dòng chảy lũ cực đoan dưới tác động của biến đổi khí hậu: Dưới kịch bản RCP 4.5 và RCP 8.5 trong giai đoạn 2019–2100, tổng lượng mưa mùa mưa có xu hướng tăng từ 6,2% đến 15,8% so với thời kỳ nền. Điều này dẫn tới sự gia tăng tương ứng của lưu lượng đỉnh lũ lớn nhất từ 10,5% đến 22,3% tại trạm Yên Bái, làm tăng tần suất xuất hiện của các trận lũ cực lớn chu kỳ 10 năm và 20 năm.
  3. Nguy cơ suy giảm nghiêm trọng dòng chảy kiệt: Lưu lượng dòng chảy trung bình mùa cạn trong các thập niên tương lai giảm từ 8,5% đến 18,3% dưới kịch bản RCP 4.5 và giảm mạnh tới 24,6% dưới kịch bản RCP 8.5. Tình trạng cạn kiệt nguồn nước trở nên gay gắt nhất vào các tháng II và III, đe dọa trực tiếp đến việc lấy nước của các hệ thống thủy nông hạ du.
  4. Gia tăng biên độ biến động dòng chảy chu kỳ ngắn: Tổng lượng dòng chảy lũ tích lũy trong các khoảng thời gian 1 ngày, 3 ngày và 7 ngày lớn nhất đều có xu thế tăng mạnh, rút ngắn thời gian tập trung dòng chảy trên sườn dốc do tác động cộng hưởng của mưa cường suất lớn và đặc điểm địa hình chia cắt mạnh.

Thảo luận kết quả

Nguyên nhân căn bản của sự phân bố dòng chảy không đều bắt nguồn từ cấu trúc địa hình lòng máng kẹp giữa hai dãy núi cao: dãy Hoàng Liên Sơn ở phía Tây (đỉnh Fansipan cao 3.143 m) và dãy Con Voi ở phía Đông. Cấu trúc này tạo thành hành lang hút ẩm cực mạnh từ gió mùa mùa hè Tây Nam và Đông Nam, gây ra các tâm mưa lớn cục bộ có lượng mưa năm vượt trên 2.200 mm tại khu vực Sa Pa và Thượng Than Uyên, trong khi lượng mưa phía thượng nguồn Trung Quốc chỉ dao động từ 1.000 mm đến 1.400 mm.

Dữ liệu kết quả nghiên cứu có thể được trực quan hóa tối ưu thông qua các biểu đồ phân bố đường tần suất dòng chảy (Flow Duration Curves), đồ thị thủy văn mô phỏng dạng Hydrograph so sánh giữa lưu lượng thực đo và lưu lượng tính toán tại trạm Lào Cai, Yên Bái, cũng như các bản đồ phân vùng lượng mưa đẳng trị không gian trên toàn lưu vực.

So sánh với các nghiên cứu quốc tế sử dụng mô hình khí tượng thủy văn động lực như nghiên cứu tại lưu vực sông Kunhar của Pakistan (tăng dòng chảy năm nhưng giảm dòng chảy mùa kiệt) hay lưu vực sông Kabul tại Afghanistan (dòng chảy biến động từ -14% đến +95% tùy kịch bản), kết quả tại lưu vực sông Thao thể hiện sự tương đồng rõ rệt về quy luật phân cực nguồn nước dưới tác động của biến đổi khí hậu toàn cầu. Sự phối hợp giữa mô hình động lực học khí quyển và mô hình thủy văn vật lý trong luận văn đã khắc phục hoàn toàn điểm mù thông tin tại vùng thượng nguồn lãnh thổ Trung Quốc.

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

Nhằm chủ động thích ứng với tình trạng biến động nguồn nước phức tạp trên lưu vực sông Thao, luận văn đưa ra 4 nhóm khuyến nghị mang tính thực thi cao:

  1. Nâng cấp và số hóa hệ thống quan trắc cảnh báo sớm lũ lụt: Bộ Tài nguyên và Môi trường cùng Ủy ban nhân dân các tỉnh Lào Cai, Yên Bái, Phú Thọ cần ưu tiên đầu tư lắp đặt mạng lưới radar thời tiết và trạm đo mưa tự động mật độ cao dọc hành lang sông Thao. Mục tiêu đến năm 2025 nâng cao độ chính xác dự báo lũ trước 24–48 giờ, giúp giảm thiểu ít nhất 25% thiệt hại về người và tài sản do lũ quét và sạt lở đất.
  2. Thúc đẩy cơ chế hợp tác chia sẻ thông tin thủy văn xuyên biên giới: Bộ Ngoại giao phối hợp cùng Bộ Nông nghiệp và Phát triển Nông thôn đẩy mạnh đàm phán song phương với các cơ quan quản lý nước tỉnh Vân Nam (Trung Quốc). Đích ngắm giai đoạn 2024–2030 là đạt được thỏa thuận chia sẻ dữ liệu xả lũ theo thời gian thực từ các đập Nanshan và Madushan, đảm bảo duy trì dòng chảy môi trường tối thiểu không dưới 150 m³/s về hạ lưu trong mùa cạn.
  3. Tái cơ cấu quy hoạch tưới tiêu và sản xuất nông nghiệp mùa khô: Sở Nông nghiệp và Phát triển Nông thôn các tỉnh vùng hạ lưu cần chủ động chuyển đổi cơ cấu cây trồng tiêu thụ ít nước và nâng cấp công trình trạm bơm dã chiến lấy nước trên sông Thao. Giai đoạn 2024–2028 phấn đấu tiết kiệm 20% lượng nước tưới mặt ruộng và hạn chế tối đa nguy cơ treo trạm bơm do mực nước sông hạ thấp dưới chuẩn thiết kế.
  4. Mở rộng ứng dụng công nghệ mô hình WEHY-HCM cho các lưu vực sông liên quốc gia lân cận: Viện Khoa học Thủy lợi Việt Nam và Trường Đại học Thủy lợi cần tiếp tục chuẩn hóa và chuyển giao bộ công cụ mô hình WEHY-HCM cho các hệ thống sông Đà và sông Lô. Lộ trình đến năm 2030 hoàn thiện hệ thống hỗ trợ ra quyết định (DSS) vận hành liên hồ chứa theo thời gian thực trên toàn lưu vực sông Hồng - Thái Bình.

Đố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ị thực tiễn sâu rộng cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước về tài nguyên nước và khí tượng thủy văn: Sử dụng khung tính toán và nguồn số liệu dự báo của luận văn làm căn cứ khoa học xây dựng Quy hoạch lưu vực sông Hồng - Thái Bình tầm nhìn 2050, đồng thời ban hành các phương án ứng phó hạn hán và xâm nhập mặn vùng hạ du.
  • Đơn vị quản lý, vận hành hệ thống hồ chứa và công trình thủy lợi: Ứng dụng số liệu dự báo dòng chảy đến theo chu kỳ 10 ngày, tháng và mùa để tối ưu hóa quy trình điều tiết nước, gia tăng hiệu suất tích nước phát điện và đảm bảo an toàn đập mùa mưa lũ.
  • Các nhà khoa học, giảng viên và học viên cao học chuyên ngành Thủy văn, Biến đổi khí hậu, Địa lý tài nguyên: Tham khảo quy trình thiết lập mô hình kết hợp khí quyển - thủy văn động lực WEHY-HCM, phương pháp xử lý dữ liệu tái phân tích ERA-20C và kỹ thuật chi tiết hóa động lực học không gian.
  • Các tổ chức quốc tế và chuyên gia quản lý lưu vực sông xuyên biên giới: Sử dụng luận văn như một trường hợp điển hình (case study) mẫu mực về việc ứng dụng công nghệ viễn thám và mô hình toán để giải quyết bài toán thiếu dữ liệu quan trắc tại các lưu vực sông quốc tế.

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

1. Mô hình khí tượng thủy văn WEHY-HCM có ưu điểm gì vượt trội so với các mô hình thủy văn truyền thống?
Mô hình WEHY-HCM kết hợp trực tiếp mô hình khí quyển động lực học khu vực với mô hình thủy văn dựa trên bản chất vật lý phân bố không gian. Khác với các mô hình truyền thống (như HEC-HMS, NAM hay SWAT) phụ thuộc lớn vào việc hiệu chỉnh chuỗi số liệu quan trắc mưa mặt đất, WEHY-HCM trích xuất trực tiếp các tham số vật lý từ dữ liệu vệ tinh, cho phép mô phỏng chính xác cả những lưu vực hoàn toàn không có trạm đo.

2. Tại sao nghiên cứu có thể đánh giá được dòng chảy sông Thao dù thiếu số liệu quan trắc từ phía Trung Quốc?
Luận văn đã khai thác bộ dữ liệu tái phân tích khí quyển toàn cầu ERA-20C với chuỗi thời gian 115 năm kết hợp ảnh viễn thám độ phân giải cao (ASTER DEM 30 m, SoilGrids 250 m, MODIS LAI). Thông qua phương pháp chi tiết hóa động lực học 3 miền tính toán lồng ghép (từ 81 km xuống 9 km), mô hình đã tự động tái tạo hoàn chỉnh trường mưa và các yếu tố khí tượng trên toàn bộ 39.800 km² lãnh thổ thượng nguồn.

3. Biến đổi khí hậu theo kịch bản RCP 8.5 sẽ tác động thế nào đến tài nguyên nước sông Thao?
Theo kịch bản phát thải cao RCP 8.5 giai đoạn 2019–2100, tài nguyên nước sông Thao sẽ chuyển biến theo hướng cực đoan hóa mạnh mẽ: lượng mưa mùa lũ tăng tới 15,8%, lưu lượng đỉnh lũ lớn nhất tăng 22,3% làm gia tăng nguy cơ ngập lụt nghiêm trọng; ngược lại, dòng chảy mùa cạn suy giảm sâu tới 24,6%, đẩy hạ lưu vào tình trạng thiếu nước sản xuất gay gắt.

4. Dữ liệu viễn thám mặt đệm nào giữ vai trò quan trọng nhất trong việc thiết lập mô hình WEHY?
Mô hình độ cao số ASTER DEM 30 m và bản đồ thổ nhưỡng SoilGrids 250 m giữ vai trò then chốt nhất. ASTER DEM quyết định việc phân chia cấu trúc lưu vực, hướng dốc và vận tốc dòng chảy tràn, trong khi SoilGrids cung cấp chính xác tỷ lệ các thành phần hạt đất, kiểm soát toàn bộ cơ chế thấm và hình thành tầng trữ nước ngầm của lưu vực.

5. Kết quả và quy trình mô hình hóa trong luận văn có thể chuyển giao cho các lưu vực khác không?
Hoàn toàn có thể. Quy trình tích hợp bộ công cụ WEHY-HCM kết hợp dữ liệu tái phân tích khí quyển và dữ liệu viễn thám toàn cầu là một giải pháp công nghệ mở, có khả năng áp dụng hiệu quả cho bất kỳ lưu vực sông xuyên biên giới nào tại Việt Nam như sông Đà, sông Lô, sông Mã hay lưu vực sông Mê Công.

Kết luận

  • Luận văn đã thiết lập thành công mô hình khí tượng thủy văn động lực tích hợp WEHY-HCM cho toàn bộ lưu vực sông Thao (51.800 km²), giải quyết triệt để bài toán thiếu hụt số liệu đo đạc xuyên biên giới từ phía thượng nguồn Trung Quốc.
  • Khôi phục chuẩn xác chuỗi dòng chảy lịch sử và làm rõ đặc điểm phân bố nguồn nước hiện trạng với 75,5% – 80,4% lượng nước tập trung trong mùa mưa lũ và chỉ 19,6% – 24,5% trong mùa cạn.
  • Định lượng chi tiết tác động của biến đổi khí hậu theo các kịch bản RCP 4.5 và RCP 8.5 đến năm 2100, chỉ ra nguy cơ gia tăng 10,5% – 22,3% lưu lượng đỉnh lũ và sự suy giảm tới 24,6% dòng chảy kiệt mùa khô.
  • Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, đặt mục tiêu đến năm 2030 số hóa hoàn chỉnh hệ thống dự báo cảnh báo lũ và thiết lập cơ chế giám sát nguồn nước xuyên biên giới hiệu quả.
  • Nghiên cứu là tài liệu tham khảo khoa học có giá trị cao, thúc đẩy việc mở rộng ứng dụng công nghệ viễn thám và mô hình toán động lực học trong quản lý bền vững tài nguyên nước lưu vực sông tại Việt Nam.