Tổng quan nghiên cứu

Hồ chứa thủy điện Hòa Bình trên dòng chính sông Đà là công trình thủy điện chiến lược quốc gia với dung tích toàn bộ 9,45 tỷ m³ và dung tích hữu ích 5,65 tỷ m³, giữ vai trò then chốt trong việc cắt lũ cho đồng bằng Bắc Bộ, phát điện với sản lượng 8,16 tỷ kWh/năm và cung cấp nước tưới tiêu. Kể từ khi tích nước đến mực nước dâng bình thường ở cao trình 115–117 m vào năm 1990, hiện tượng sa bồi lòng hồ đã diễn ra với tốc độ đáng báo động. Trong giai đoạn 1990–2013, tổng lượng bùn cát bồi lắng tích lũy đã vượt mức 1.387,1 triệu m³, trung bình mỗi năm giữ lại 57,8 triệu m³ bùn cát, làm bồi lấp tới 37% dung tích chết và làm suy giảm trực tiếp tuổi thọ công trình.

Vấn đề nghiên cứu trọng tâm là đánh giá định lượng hiện trạng bồi lắng lòng dẫn và biến động bùn cát trong bối cảnh thủy văn thay đổi khi có thêm các hồ chứa bậc thang thượng lưu. Mục tiêu cụ thể của đề tài là phân tích nguyên nhân, cơ chế bồi xói lòng hồ Hòa Bình trên đoạn sông dài hơn 200 km từ trạm thủy văn Tạ Bú đến đập thủy điện, ứng dụng mô hình một chiều HEC-RAS và phương pháp so sánh thể tích địa hình để dự báo diễn biến bồi lắng đến năm 2020. Nghiên cứu mang ý nghĩa thực tiễn sâu sắc, cung cấp cơ sở dữ liệu khoa học giúp duy trì khả năng trữ nước hữu ích trên 5,5 tỷ m³, nâng cao hiệu quả vận hành liên hồ chứa theo Quyết định 1622/QĐ-TTg và kiểm soát rủi ro sạt lở bờ hồ.

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 thủy lực sông ngòi và động lực học dòng bùn cát, kết hợp hai khung lý thuyết cốt lõi:

  • Lý thuyết dòng chảy không ổn định một chiều: Dựa trên hệ phương trình Saint-Venant gồm phương trình liên tục bảo toàn khối lượng và phương trình động lượng biến thiên, mô tả sự thay đổi lưu lượng và mực nước dọc theo lòng dẫn tự nhiên.
  • Lý thuyết chuyển tải và bảo toàn khối lượng bùn cát: Dựa trên phương trình Exner mô tả tương quan giữa độ dốc dòng chảy, biến thiên thể tích bùn cát đáy và lưu lượng bùn cát qua từng tiết diện khống chế để xác định trạng thái bồi tụ hay xói lở.

Các khái niệm chuyên ngành then chốt bao gồm: dung tích chết (khoang chứa bùn cát dưới cao trình 80 m), dung tích hữu ích (khoang điều tiết từ 80 m đến 117 m), hệ số nhám thủy lực Manning, cấp phối hạt trầm tích và quá trình sa bồi vùng hồi thủy. Mô hình tích hợp quy trình vận hành liên hồ chứa trên lưu vực sông Hồng nhằm kiểm soát chế độ dòng chảy mùa lũ.

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

  • Nguồn dữ liệu và cỡ mẫu: Đề tài thu thập chuỗi dữ liệu địa hình mặt cắt ngang lòng hồ gồm 59 mặt cắt chuẩn trải dài 213 km từ trạm Tạ Bú về tuyến đập. Chuỗi dữ liệu khí tượng thủy văn, lưu lượng nước và bùn cát lơ lửng được đo đạc liên tục hàng ngày từ năm 1990 đến năm 2013, tập trung sâu vào giai đoạn hiệu chỉnh và kiểm định từ ngày 01/01/2007 đến ngày 31/12/2013.
  • Phương pháp chọn mẫu: Khảo sát địa hình đáy hồ cố định vào tháng 12 hàng năm, thời điểm dòng chảy và nồng độ bùn cát ở mức thấp nhất, lòng dẫn ổn định nhất để đảm bảo tính đồng nhất của số liệu đo đạc trắc đạc lòng hồ.
  • Phương pháp phân tích: Kết hợp phương pháp so sánh thể tích hình học không gian giữa các kỳ đo đạc với mô hình số trị HEC-RAS 4.1. Lý do lựa chọn mô hình 1D HEC-RAS là vì hồ Hòa Bình có hình thái dạng thung lũng sông kéo dài hẹp và sâu (tỷ lệ chiều dài trên chiều rộng lớn), giúp tối ưu hóa tốc độ tính toán cho mô phỏng dài hạn nhiều năm mà vẫn đảm bảo độ chính xác cao hơn so với mô hình 2D/3D vốn đòi hỏi khối lượng dữ liệu đầu vào cực lớn. Quá trình mô phỏng được hiệu chỉnh trong giai đoạn 2007–2010 và kiểm định trong giai đoạn 2011–2013 qua hệ số tương quan R² đạt trên 0,85.

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

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

  • Bồi lắng làm suy giảm nghiêm trọng dung tích thiết kế: Sau 24 năm vận hành (1990–2013), tổng khối lượng bùn cát tích tụ trong lòng hồ đạt 1.387,1 triệu m³, tương đương tốc độ bồi lắng bình quân 57,8 triệu m³/năm, chiếm trọn 37% dung tích chết ban đầu. Lớp bồi tích nâng cao trình đáy hồ lên trung bình 0,24 m/năm trên toàn tuyến.
  • Phân hóa không gian rõ rệt: Khu vực trung lưu từ Bản Khộc (mặt cắt 37) đến Suối Lúa – Nà Giang (mặt cắt 19) dài 56,1 km là điểm hội tụ sa bồi nặng nề nhất, chiếm tới 77,9% tổng lượng bùn cát toàn hồ (khoảng 1.080,5 triệu m³). Bề dày tầng phù sa lắng đọng tại khu vực này lên tới 40 m, hình thành bãi bồi lớn dịch chuyển về phía hạ lưu với tốc độ trung bình 3,9 km/năm. Trong khi đó, khu vực thượng lưu chỉ chiếm 5,78% (80,12 triệu m³) và khu vực hạ lưu trước đập chiếm 16,3% (226,57 triệu m³).
  • Tác động chặn bùn cát rõ rệt của thủy điện Sơn La: Từ khi thủy điện Sơn La tích nước và phát điện (giai đoạn 2010–2013), tổng lượng bùn cát về hồ Hòa Bình giảm mạnh chỉ còn 54,6 triệu m³ trong 4 năm, bình quân đạt 13,7 triệu m³/năm, giảm khoảng 76,3% so với mức trung bình nhiều năm trước đó.

Thảo luận kết quả

Hiện tượng bồi lắng tập trung ở khu vực giữa hồ bắt nguồn từ quy luật thủy lực khi dòng chảy xiết từ sông Đà đổ vào vùng nước dềnh có diện tích mặt cắt ướt mở rộng đột ngột, làm vận tốc dòng chảy suy giảm mạnh khiến các cấp hạt cát sỏi thô và phù sa hạt mịn lắng đọng nhanh chóng. Về mặt tự nhiên, lưu vực sông Đà có tới 81,7% diện tích là đất xám trên địa hình đồi dốc bị chia cắt mạnh, kết hợp các đợt mưa lũ cực đoan như trận lũ lịch sử năm 1996 với đỉnh lũ 33.671 m³/s đã cuốn trôi 87,5 triệu m³ bùn cát xuống lòng hồ. Về mặt nhân sinh, tập quán canh tác nông nghiệp trên đất dốc, nạn phá rừng đầu nguồn và 13.000 ha đất bán ngập bị khai thác thiếu quy hoạch đã gia tăng xói mòn bề mặt.

Khi so sánh với các hồ chứa khác, hồ Hòa Bình có mức độ bồi lắng tuyệt đối lớn hơn rất nhiều so với hồ Thác Bà (163–175 triệu m³ sau 30 năm) hay hồ Trị An (bồi lắng trung bình 12 cm sau 23 năm). Kết quả nghiên cứu có thể được trực quan hóa tối ưu thông qua:

  • Biểu đồ mặt cắt dọc địa hình đáy hồ qua các thời kỳ (1990, 1996, 2009, 2013) để làm nổi bật sự hình thành và dịch chuyển của đỉnh bãi bồi.
  • Bảng tương quan lưu lượng nước thực đo và tính toán mô phỏng tại trạm Hòa Bình với chỉ số R² tiệm cận 0,9.
  • Bản đồ số hóa GIS phân vùng nguy cơ bồi xói 12 mặt cắt trọng điểm, phục vụ công tác cảnh báo trực quan cho đơn vị quản lý vận hành.

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

  • Đẩy mạnh phục hồi thảm phủ thực vật lưu vực thượng nguồn: Ủy ban nhân dân các tỉnh Hòa Bình, Sơn La, Lai Châu và Điện Biên cần phối hợp thực hiện chính sách chi trả dịch vụ môi trường rừng, nâng cao độ che phủ rừng phòng hộ đầu nguồn lên trên 65% giai đoạn 2018–2025 nhằm giảm thiểu 25–30% lượng đất xói mòn trôi vào các nhánh suối nhập lưu.
  • Xây dựng công trình bậc thang kiểm soát bùn cát nhánh: Bộ Nông nghiệp và Phát triển Nông thôn cùng Tập đoàn Điện lực Việt Nam (EVN) cần quy hoạch và đẩy nhanh tiến độ các hồ thủy điện nhỏ và đập giữ cát trên các nhánh nhập lưu chính như Suối Sập, Nậm Pia, Nậm Chim, hoàn thành trước năm 2025 để giữ lại khoảng 1,0–1,9 triệu tấn bùn cát khu giữa mỗi năm.
  • Áp dụng quy chế xả cát và điều tiết mực nước linh hoạt: Công ty Thủy điện Hòa Bình cần phối hợp với Trung tâm Điều độ Hệ thống điện Quốc gia tối ưu hóa quy trình vận hành liên hồ chứa theo Quyết định 1622/QĐ-TTg, thực hiện hạ mực nước hồ đón lũ sớm kết hợp mở các cửa xả đáy để tạo dòng chảy động lực đẩy bùn cát mịn qua thân đập trong các đợt lũ chính vụ.
  • Ban hành bản đồ số hóa nguy cơ bồi xói và kiểm soát khai thác khoáng sản: Sở Tài nguyên và Môi trường tỉnh Hòa Bình cần ứng dụng bản đồ nguy cơ bồi lắng tỷ lệ 1:25.000 để khoanh vùng cấm khai thác cát sỏi trái phép tại các chân mái dốc từ năm 2020, đồng thời xây dựng kè bảo vệ tại các mặt cắt có nguy cơ sạt lở từ 6–10 m ven lòng hồ.

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

  • Đơn vị quản lý và vận hành nhà máy thủy điện: Cán bộ kỹ thuật của Công ty Thủy điện Hòa Bình, Thủy điện Sơn La và EVN có thể ứng dụng số liệu bồi lắng và quy luật dịch chuyển bãi bồi để lập kế hoạch duy tu cửa lấy nước, nạo vét cục bộ và vận hành cửa xả đáy an toàn.
  • Cơ quan quản lý nhà nước về tài nguyên nước và môi trường: Các chuyên viên thuộc Bộ Tài nguyên và Môi trường, Tổng cục Thủy lợi tham khảo phương pháp tính toán để xây dựng chính sách quản lý lưu vực sông liên tỉnh, quy hoạch vùng đất bán ngập 13.000 ha và thẩm định quy trình vận hành liên hồ.
  • Các nhà khoa học, giảng viên và học viên cao học: Giảng viên và học viên chuyên ngành Khoa học Môi trường, Kỹ thuật Tài nguyên nước, Thủy văn tại các trường đại học khối kỹ thuật sử dụng luận văn làm tài liệu tham khảo chuẩn mực về ứng dụng HEC-RAS trong mô phỏng bùn cát 1D.
  • Đơn vị tư vấn quy hoạch phòng chống thiên tai: Các viện nghiên cứu quy hoạch thủy lợi, đơn vị tư vấn xây dựng công trình ngầm và thủy điện sử dụng kết quả dự báo bồi lắng đến năm 2020 để tính toán dung tích bồi lắng thiết kế cho các công trình hạ tầng hạ lưu sông Đà và sông Hồng.

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

  • Bồi lắng lòng hồ Hòa Bình đã làm suy giảm bao nhiêu dung tích chết của hồ? Theo số liệu quan trắc và tính toán địa hình giai đoạn 1990–2013, tổng lượng bùn cát lắng đọng đạt 1.387,1 triệu m³, làm mất đi 37% dung tích chết thiết kế ban đầu (3,85 tỷ m³), trực tiếp nâng cao trình đáy hồ trung bình 0,24 m mỗi năm.

  • Sự xuất hiện của thủy điện Sơn La ảnh hưởng như thế nào đến tốc độ bồi lắng hồ Hòa Bình? Kể từ khi hồ Sơn La tích nước và đi vào vận hành từ năm 2010, phần lớn lượng bùn cát thô ở thượng nguồn sông Đà đã bị giữ lại, làm tốc độ bồi lắng hồ Hòa Bình giảm từ mức bình quân 57,8 triệu m³/năm xuống còn 13,7 triệu m³/năm trong giai đoạn 2010–2013, giảm khoảng 76,3%.

  • Vị trí nào trong lòng hồ Hòa Bình xảy ra hiện tượng bồi tích nghiêm trọng nhất? Đoạn lòng hồ từ Bản Khộc đến Suối Lúa – Nà Giang (dài 56,1 km ở khu vực trung lưu) là nơi bồi lắng nặng nề nhất, chiếm 77,9% tổng lượng phù sa toàn hồ, hình thành lớp bồi tích dày tới 40 m và tạo nên bãi bồi dịch chuyển về phía hạ lưu 3,9 km/năm.

  • Tại sao tác giả lại lựa chọn mô hình 1D HEC-RAS thay vì các mô hình 2D hay 3D? Hồ Hòa Bình có hình thái thung lũng sông điển hình với chiều dài hơn 200 km nhưng chiều ngang hẹp. Mô hình HEC-RAS 1D giải tốt hệ phương trình Saint-Venant và bảo toàn khối lượng bùn cát, tiết kiệm tài nguyên tính toán cho chuỗi thời gian dài nhiều năm với độ tin cậy R² trên 0,85.

  • Những giải pháp phi công trình nào đóng vai trò cốt lõi để hạn chế bồi lắng? Hai giải pháp then chốt là tăng cường giao đất giao rừng kết hợp chi trả dịch vụ môi trường rừng để nâng độ che phủ lưu vực lên trên 65%, cùng với việc quy hoạch chuyển đổi cơ cấu cây trồng lâu năm trên 13.000 ha đất bán ngập nhằm ngăn chặn xói mòn mặt đất.

Kết luận

  • Luận văn đã xác định chính xác tổng lượng bùn cát bồi lắng tích lũy trong hồ Hòa Bình giai đoạn 1990–2013 là 1.387,1 triệu m³, làm suy giảm 37% dung tích chết và nâng cao trình đáy hồ trung bình 0,24 m/năm.
  • Nghiên cứu đã làm rõ quy luật phân hóa không gian với 77,9% lượng bồi lắng tập trung tại trung lưu (mặt cắt 19 đến 37), tạo nên dải bãi bồi dài 56,1 km có tốc độ dịch chuyển 3,9 km/năm về phía đập.
  • Đề tài chứng minh định lượng vai trò của công trình thủy điện Sơn La giúp giảm tới 76,3% lượng sa bồi hàng năm đổ về hồ Hòa Bình trong giai đoạn 2010–2013.
  • Ứng dụng thành công chuỗi dữ liệu 59 mặt cắt trên mô hình thủy lực bùn cát HEC-RAS 4.1 kết hợp GIS để xây dựng bản đồ nguy cơ bồi xói phục vụ công tác quản lý an toàn hồ đập.
  • Hệ thống giải pháp kết hợp giữa trồng rừng phòng hộ, xây dựng đập phụ giữ cát trên các nhánh suối và điều tiết xả lũ hợp lý theo Quyết định 1622/QĐ-TTg cần được các cơ quan quản lý triển khai đồng bộ ngay trong giai đoạn tới để bảo vệ bền vững nguồn nước sông Đà.