Tổng quan nghiên cứu

Lưu vực sông Vu Gia – Thu Bồn sở hữu diện tích tự nhiên 10.350 km², trải rộng trên địa bàn 17 huyện, thành phố thuộc thành phố Đà Nẵng, tỉnh Quảng Nam và một phần tỉnh Kon Tum. Đây là hệ thống thủy văn huyết mạch, cung cấp nguồn nước ngọt phục vụ đời sống cho hơn 2 triệu dân cư, chiếm 73,7% tổng dân số vùng dự án, đồng thời duy trì hoạt động sản xuất nông nghiệp, công nghiệp và hệ sinh thái đa dạng. Mặc dù tổng lượng dòng chảy toàn lưu vực hàng năm đạt xấp xỉ 19,45 tỷ m³ (trong đó sông Thu Bồn tại Giao Thủy đạt 9,25 tỷ m³ và sông Vu Gia tại Ái Nghĩa đạt 8,55 tỷ m³), tài nguyên nước lại phân bố cực kỳ bất đối xứng theo thời gian. Mùa lũ kéo dài từ tháng X đến tháng XII chiếm tới 62,5% - 69,2% tổng lượng nước năm, trong khi mùa cạn kéo dài từ tháng I đến tháng VIII chỉ chiếm 20% - 25%, đặc biệt dòng chảy kiệt tháng IV chỉ đạt 2,1% - 2,6% tổng lượng dòng chảy năm.

Vấn đề nghiên cứu then chốt xuất phát từ sự bùng nổ của các bậc thang thủy điện thượng nguồn, kết hợp hiện tượng biến đổi lòng dẫn và xói lở nghiêm trọng. Điển hình là hiện tượng cướp dòng tại thôn 9, xã Đại Cường (Đại Lộc), khiến lưu lượng từ sông Vu Gia đổ về sông Yên cấp cho Đà Nẵng giảm xuống dưới 20 m³/s trong mùa kiệt, đẩy ranh mặn 2 g/l xâm nhập sâu 13 - 16 km. Mục tiêu cụ thể của luận văn là ứng dụng mô hình MIKE 11 mô phỏng chế độ thủy lực mùa cạn, xác định ngưỡng dòng chảy tối thiểu cần thiết trên dòng chính nhằm duy trì sự sống của dòng sông, bảo vệ 134 loài cá bản địa và đáp ứng nhu cầu sử dụng nước hạ du. Nghiên cứu sử dụng chuỗi dữ liệu thủy văn 1976 - 2006, tạo cơ sở khoa học để hiện thực hóa các quy định tại Nghị định 120/2008/NĐ-CP và Nghị định 112/2008/NĐ-CP về quản lý tài nguyên nước.

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 dòng chảy môi trường (Environmental Flows) và nguyên lý dòng chảy tối thiểu được chuẩn hóa theo khung tiếp cận của Liên minh Quốc tế Bảo tồn Thiên nhiên và Tài nguyên Thiên nhiên (IUCN) cùng Viện Quản lý Nước Quốc tế (IWMI). Dòng chảy tối thiểu được xác định là lưu lượng thấp nhất cần thiết nhằm duy trì các chức năng sinh thái cơ bản của hệ thống sông ngòi, bảo tồn quần thể thủy sinh vật và đảm bảo mức khai thác tối thiểu cho các ngành kinh tế - xã hội.

Về mặt thủy lực, nghiên cứu vận dụng hệ phương trình vi phân Saint-Venant một chiều mô tả dòng chảy không ổn định trong lòng dẫn hở, được tích hợp trong phần mềm MIKE 11 do Viện Thủy lực Đan Mạch (DHI) phát triển. Khung lý thuyết bao quát 4 khái niệm cốt lõi:

  • Khái niệm modul dòng chảy kiệt ($M_0$ tính bằng l/s.km²);
  • Khái niệm chu vi ướt (Wetted Perimeter) làm ngấn sinh thái đáy;
  • Hệ thống chỉ thị chất lượng nước bằng sinh vật đáy không xương sống BMWP (Biological Monitoring Working Party);
  • Khái niệm dung tích hữu ích và dung tích điều tiết phòng lũ của các hồ chứa thượng lưu.

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

Nguồn dữ liệu đầu vào bao gồm chuỗi số liệu quan trắc khí tượng liên tục từ năm 1977 đến 2006 tại 3 trạm đại diện (Đà Nẵng, Trà My, Tam Kỳ) và dữ liệu thủy văn thực đo tại 11 trạm phân bố trên toàn lưu vực (Nông Sơn, Thành Mỹ, Ái Nghĩa, Hội Khách, Giao Thủy, Câu Lâu, Cẩm Lệ...). Dữ liệu địa hình bao gồm hệ thống mặt cắt ngang lòng sông thực đo dọc theo các nhánh sông chính Thu Bồn, Vu Gia, sông Yên, sông Vĩnh Điện và sông Quảng Huế.

Phương pháp mô hình toán được triển khai thông qua module thủy động lực học MIKE 11 HD. Quá trình thiết lập mạng sông và điều kiện biên sử dụng lưu lượng thực đo tại 2 trạm kiểm soát thượng lưu là Thành Mỹ ($F = 1.850\text{ km}^2$) và Nông Sơn ($F = 3.150\text{ km}^2$), biên dưới là diễn biến mực nước triều tại các cửa biển Cửa Đại và Cửa Hàn. Phương pháp chọn mẫu áp dụng toàn bộ chuỗi số liệu quan trắc thủy văn liên tục 30 năm (cỡ mẫu 30 năm) nhằm phân tích tần suất dòng chảy kiệt thiết kế $P = 75%$ và $P = 90%$. Quá trình hiệu chỉnh và kiểm định mô hình được thực hiện dựa trên hệ số nhám Manning ($n$) dao động từ 0,025 đến 0,045, đánh giá tính phù hợp qua chỉ số tương quan pha và đỉnh mực nước tại 5 trạm thủy văn hạ lưu. Nghiên cứu kết hợp phương pháp điều tra thực địa 595 công trình thủy lợi và khảo sát hệ sinh thái thủy sinh vật tại 50 điểm lấy mẫu sinh học.

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

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

Mô phỏng thủy lực và phân tích chuỗi dữ liệu thực tế đã mang lại 4 phát hiện khoa học quan trọng:

  • Sự bất cân đối nghiêm trọng trong phân phối dòng chảy mùa cạn: Dòng chảy mùa kiệt phân hóa sâu sắc qua 2 giai đoạn: dòng chảy ổn định (tháng I - IV) và không ổn định (tháng V - VIII). Tại trạm Nông Sơn, lưu lượng bình quân tháng IV xuống mức thấp kỷ lục 72,5 m³/s (chiếm 2,2% lượng nước năm), chênh lệch gấp 13,8 lần so với tháng đỉnh lũ XI (999,7 m³/s). Tại trạm Thành Mỹ, lưu lượng tháng IV chỉ đạt 38,0 m³/s (chiếm 2,63% lượng nước năm) so với tháng XI đạt 385,7 m³/s. Hệ số biến động dòng chảy năm đạt $Cv = 0,32 - 0,35$, cao hơn hẳn các lưu vực sông khu vực Tây Nguyên ($Cv = 0,15 - 0,25$).

  • Suy giảm lưu lượng cục bộ do thay đổi hình thái lòng sông: Hiện tượng cướp dòng tại xã Đại Cường đã chuyển hướng phần lớn nước từ Vu Gia sang Thu Bồn qua sông Quảng Huế, khiến lưu lượng chuyển tiếp về sông Yên chỉ còn dưới 20 m³/s trong các tháng kiệt đỉnh điểm. Sự sụt giảm này làm tê liệt 150 trạm bơm nông nghiệp hạ lưu, khiến diện tích thực tế tưới chỉ đạt 14.613 ha (đạt 60% công suất thiết kế 24.228 ha).

  • Xâm nhập mặn gia tăng về cường độ và phạm vi: Khi dòng chảy ngọt từ thượng lưu sụt giảm dưới ngưỡng kiểm soát, mặn 2 g/l đã xâm nhập qua cầu Câu Lâu 2 km, chiều dài xâm nhập đạt 13 - 16 km từ cửa biển, ảnh hưởng trên phạm vi 24 km. Các trạm bơm sinh hoạt và nông nghiệp dọc sông Cẩm Lệ, Vĩnh Điện buộc phải ngừng vận hành từ 5 đến 7 ngày trong mỗi chu kỳ triều cường.

  • Ngưỡng dòng chảy tối thiểu xác định qua mô hình: Kết quả tính toán chỉ ra rằng, để ngăn mặn và bảo vệ hệ sinh thái thủy sinh (với 50 loài tảo và 21 loài cá hồ Đồng Nghệ, 134 loài cá sông chính), dòng chảy tối thiểu cần duy trì thường xuyên tại trạm Ái Nghĩa phải đạt từ 35 đến 42 m³/s, và tại trạm Giao Thủy phải đạt từ 65 đến 75 m³/s.

Thảo luận kết quả

Nguyên nhân trực tiếp dẫn tới tình trạng kiệt nước hạ du là sự kết hợp giữa địa hình đồi núi dốc đứng (độ dốc trung bình 21,3% - 25,5%) làm nước lũ thoát nhanh ra biển, và sự gia tăng của 11 bậc thang thủy điện lớn trên dòng chính (A Vương, Đak Mi 4, Sông Bung 4, Sông Tranh 2...). Đáng chú ý, công trình thủy điện Đak Mi 4 chuyển nước từ sông Vu Gia sang sông Thu Bồn đã làm suy thoái nghiêm trọng chế độ dòng chảy tự nhiên của phân lưu Vu Gia.

So sánh với các nghiên cứu dòng chảy môi trường trên sông Hương của IUCN (năm 2004) hay sông Ba và Trà Khúc của Trường Đại học Thủy lợi (năm 2006), nghiên cứu này khẳng định việc ứng dụng mô hình MIKE 11 giúp phản ánh chính xác tương quan phi tuyến giữa lưu lượng xả và độ sâu ngập của chu vi ướt lòng dẫn.

Các dữ liệu nghiên cứu này được trực quan hóa tối ưu thông qua:

  • Biểu đồ đường quá trình mực nước thực đo và tính toán tại các trạm Hội Khách, Ái Nghĩa, Câu Lâu;
  • Bảng ma trận phân phối dòng chảy thiết kế ứng với các cấp tần suất $P = 50%, 75%, 90%$;
  • Bản đồ không gian hóa phạm vi xâm nhập mặn 1 g/l và 2 g/l trên toàn dải hạ lưu ven biển.

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

Nhằm bảo đảm duy trì dòng chảy tối thiểu trên dòng chính sông Vu Gia – Thu Bồn, luận văn đề xuất 4 nhóm giải pháp mang tính hành động cao:

  • Nâng cấp và xây dựng công trình chỉnh trị phân lưu: Khẩn trương thi công kè chỉnh trị và đập dâng điều tiết tại cửa sông Quảng Huế nhằm tái lập tỷ lệ phân chia dòng chảy hợp lý (50:50) giữa Vu Gia và Thu Bồn. Cải tạo nâng cấp cụm đập ngăn mặn An Trạch trên sông Yên, đập Thanh Quýt và đập Hà Thanh, hướng tới target duy trì lưu lượng xả về hạ lưu sông Yên luôn đạt trên 30 m³/s trước năm 2025. Chủ thể thực hiện: Bộ Nông nghiệp và Phát triển nông thôn chủ trì, phối hợp cùng Ủy ban nhân dân thành phố Đà Nẵng và tỉnh Quảng Nam.

  • Thiết lập quy trình vận hành liên hồ chứa mùa kiệt: Ban hành quy chế điều tiết bắt buộc đối với các nhà máy thủy điện thượng nguồn (A Vương, Đak Mi 4, Sông Bung 2, 4, 5, Sông Tranh 1, 2). Buộc các hồ chứa xả dòng chảy sau đập không thấp hơn lưu lượng kiệt nhỏ nhất trung bình nhiều năm ($M_{kiet} \ge 8,76\text{ l/s.km}^2$), thực hiện các đợt xả nước tập trung theo chu kỳ mùa vụ gieo cấy lúa Đông Xuân và Hè Thu. Timeline hoàn thành: Giai đoạn 2013 - 2015. Chủ thể thực hiện: Bộ Công Thương, Bộ Tài nguyên và Môi trường, Tập đoàn Điện lực Việt Nam (EVN).

  • Thành lập cơ quan điều phối quản lý tổng hợp lưu vực sông: Thành lập Ban Quản lý lưu vực sông Vu Gia – Thu Bồn thống nhất quyền hạn giữa hai địa phương, lắp đặt hệ thống 11 trạm quan trắc tự động SCADA truyền tin thời gian thực để kiểm soát việc khai thác của 595 công trình thủy lợi, cắt giảm tỷ lệ thất thoát nước tưới xuống dưới 20% vào năm 2020. Chủ thể thực hiện: Ủy ban nhân dân thành phố Đà Nẵng và Ủy ban nhân dân tỉnh Quảng Nam.

  • Kiểm soát ô nhiễm và phục hồi rừng đầu nguồn: Đẩy mạnh chương trình đóng cửa và hoàn nguyên môi trường đối với 100% tụ điểm khai thác vàng sa khoáng trái phép dùng hóa chất Xyanua tại các huyện Phước Sơn, Trà My, Nam Giang; nâng độ che phủ rừng tự nhiên đầu nguồn từ 38% lên trên 45% nhằm nâng cao khả năng giữ nước mặt và bổ cập nước ngầm. Chủ thể thực hiện: Sở Tài nguyên và Môi trường cùng Chi cục Kiểm lâm tỉnh Quảng Nam.

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

  • Cơ quan hoạch định chính sách và quản lý nhà nước: Bộ Tài nguyên và Môi trường, Cục Quản lý Tài nguyên nước, Sở Tài nguyên và Môi trường Đà Nẵng và Quảng Nam. Luận văn cung cấp luận cứ thực tiễn để thẩm định hồ sơ cấp phép khai thác nước mặt, xây dựng quy định dòng chảy môi trường và khung pháp lý vận hành liên hồ chứa.

  • Đơn vị quản lý, vận hành thủy điện và công trình thủy nông: Các công ty cổ phần thủy điện (A Vương, Sông Bung, Đak Mi, Sông Tranh) và các Công ty Trách nhiệm hữu hạn Một thành viên Khai thác Công trình Thủy lợi. Tài liệu giúp xác định chuẩn xác lưu lượng xả đáy tối thiểu, điều hòa hài hòa lợi ích kinh tế phát điện và an toàn cấp nước tưới tiêu.

  • Viện nghiên cứu và cơ sở đào tạo chuyên ngành: Các trường đại học khối kỹ thuật, môi trường, thủy lợi và các viện nghiên cứu khí tượng thủy văn. Đây là tài liệu tham khảo chất lượng về phương pháp luận ứng dụng mô hình MIKE 11 HD trong mô phỏng thủy lực mạng sông phức tạp miền Trung.

  • Các tổ chức phi chính phủ và chuyên gia bảo tồn sinh thái: Các nhóm nghiên cứu bảo tồn nguồn lợi thủy sản, đa dạng sinh học nguồn nước. Luận văn cung cấp số liệu toàn diện về danh mục 134 loài cá, 50 loài thực vật phiêu sinh và chỉ số chất lượng nước sinh thái BMWP làm cơ sở quan trắc môi trường lâu dài.

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

1. Dòng chảy tối thiểu trên sông Vu Gia – Thu Bồn được định lượng như thế nào?
Dòng chảy tối thiểu là lưu lượng nước thấp nhất ở hạ lưu cần được duy trì để bảo tồn hệ sinh thái thủy sinh và cấp nước dân sinh theo Nghị định 112/2008/NĐ-CP. Qua tính toán mô hình thủy lực MIKE 11, lưu lượng kiệt kiểm soát tại trạm Ái Nghĩa phải đạt tối thiểu 35 - 42 m³/s và tại trạm Giao Thủy đạt 65 - 75 m³/s để duy trì chu vi ướt an toàn cho sinh thái đáy.

2. Tại sao mô hình MIKE 11 được chọn làm công cụ tính toán chính trong đề tài?
MIKE 11 do DHI phát triển có khả năng mô phỏng thủy động lực học một chiều (module HD) cực kỳ chính xác cho mạng sông phức tạp có nhiều phân lưu, hợp lưu, công trình cống, đập dâng và chịu tác động triều mạnh. Bộ thông số mô hình được hiệu chỉnh với hệ số nhám Manning $n = 0,025 - 0,045$, đạt độ trùng khớp cao giữa đường mực nước thực đo và tính toán tại 5 trạm thủy văn.

3. Hiện tượng cướp dòng tại xã Đại Cường gây ra hậu quả gì đối với hạ lưu Đà Nẵng?
Hiện tượng sông cướp dòng tự nhiên tại thôn 9, xã Đại Cường đã chuyển hướng dòng chảy chính từ Vu Gia sang Thu Bồn, làm lưu lượng nước cấp về sông Yên chảy ra Cửa Hàn tụt sâu xuống dưới 20 m³/s trong mùa khô. Hệ quả là trạm bơm Cầu Đỏ bị nhiễm mặn nghiêm trọng, gây thiếu hụt nguồn nước sạch sinh hoạt trên diện rộng cho thành phố Đà Nẵng.

4. Khai thác vàng sa khoáng thượng nguồn tác động tiêu cực thế nào đến nguồn nước?
Hoạt động khai khoáng trái phép tại Phước Sơn, Trà My, Nam Giang làm sạt lở hàng ngàn hecta rừng phòng hộ, gây bồi lấp lòng sông và giải phóng các kim loại nặng. Đặc biệt, việc lạm dụng hóa chất cực độc Xyanua trong tuyển quặng đã làm suy giảm nghiêm trọng quần thể 134 loài cá và làm gia tăng nguy cơ suy thoái tầng nước ngầm chứa 4,8 triệu m³/ngày trên lưu vực.

5. Các công trình thủy lợi hiện hữu đáp ứng năng lực cấp nước mùa cạn ra sao?
Toàn lưu vực có 595 công trình thủy lợi (78 hồ chứa, 367 đập dâng, 150 trạm bơm), nhưng hiệu suất phục vụ thực tế chỉ đạt 50,2% năng lực thiết kế (thực tưới 22.675 ha trên thiết kế 44.632 ha). Vào mùa kiệt tháng IV và V, nhiều trạm bơm dọc sông La Thọ, Cổ Cò, Vĩnh Điện phải dừng bơm từ 5 đến 7 ngày do mực nước kiệt sâu và nhiễm mặn.

Kết luận

  • Luận văn đã thiết lập thành công mô hình thủy lực MIKE 11 HD mô phỏng mạng sông Vu Gia – Thu Bồn diện tích 10.350 km², giải quyết bài toán thủy lực mùa kiệt phức tạp.
  • Đã lượng hóa chi tiết mức độ mất cân bằng dòng chảy với mùa cạn chỉ chiếm 20% - 25% tổng lượng nước năm, đồng thời làm rõ tác động tiêu cực của hiện tượng cướp dòng tại Đại Cường và chuyển nước của thủy điện Đak Mi 4.
  • Đã xác định căn cứ khoa học cho ngưỡng dòng chảy tối thiểu tại Ái Nghĩa (35 - 42 m³/s) và Giao Thủy (65 - 75 m³/s) nhằm đẩy mặn và bảo vệ 134 loài cá bản địa.
  • Đã phân tích toàn diện hiện trạng suy thoái chất lượng nước do khai thác vàng sa khoáng bằng Xyanua và hiệu suất khai thác thực tế của 595 công trình thủy lợi (đạt 50,2%).
  • Đã xây dựng bộ giải pháp tích hợp giữa công trình chỉnh trị phân lưu, quy chế vận hành liên hồ chứa và thành lập cơ quan quản lý lưu vực sông thống nhất.

Về định hướng tiếp theo, các cơ quan chức năng cần tích hợp mô hình dòng chảy môi trường hai chiều (2D) cho vùng cửa sông ven biển và hoàn thiện hệ thống quan trắc tự động toàn lưu vực trong giai đoạn đến năm 2025. Các nhà quản lý, doanh nghiệp năng lượng và cộng đồng nghiên cứu cần chung tay áp dụng ngay các giải pháp duy trì dòng chảy tối thiểu để đảm bảo an ninh nguồn nước và phát triển bền vững lưu vực Vu Gia – Thu Bồn.