Tổng quan nghiên cứu

Lưu vực sông Mê Kông đóng vai trò huyết mạch đối với sự phát triển kinh tế và an ninh nguồn nước của khu vực Đồng bằng sông Cửu Long với diện tích tự nhiên khoảng 40.548 km² cùng mạng lưới sông ngòi chằng chịt. Trung bình mỗi năm, dòng sông này chuyển tải hơn 500 tỷ m³ nước và nguồn phù sa màu mỡ bồi đắp cho toàn bộ châu thổ. Tuy nhiên, chất lượng nước mặt tại khu vực đang đối mặt với những biến động gay gắt do tác động kép từ biến đổi khí hậu và việc phát triển thủy điện thượng nguồn. Tính đến giai đoạn nghiên cứu, trên toàn lưu vực đã có 94 công trình thủy điện vận hành, khiến tổng lượng phù sa mịn và chất rắn lơ lửng chảy về hạ lưu sụt giảm nghiêm trọng từ 160 triệu tấn vào năm 1990 xuống chỉ còn khoảng 75 triệu tấn vào năm 2015 (mức giảm hơn 53%). Sự suy giảm lưu lượng mùa khô (chỉ chiếm 15% đến 30% tổng lượng dòng chảy năm) đã thúc đẩy hiện tượng xâm nhập mặn tiến sâu vào nội đồng, gia tăng xói lở bờ sông và làm suy thoái chất lượng nước mặt tại hai nhánh sông chính là sông Tiền và sông Hậu.

Trước áp lực đó, đề tài nghiên cứu tập trung ứng dụng công nghệ viễn thám đa phổ kết hợp số liệu quan trắc mặt đất nhằm phân tích biến động không gian và thời gian của hai chỉ tiêu then chốt: độ đục và tổng chất rắn lơ lửng (TSS). Nghiên cứu được triển khai thực nghiệm trên hệ thống sông Tiền, sông Hậu và phân lưu sông Cổ Chiên trong khung thời gian từ tháng 3 năm 2014 đến tháng 9 năm 2016. Mục tiêu cụ thể là xây dựng phương trình hồi quy tương quan giữa giá trị phản xạ phổ từ ảnh vệ tinh Landsat 7, Landsat 8 với dữ liệu thực nghiệm, từ đó mô phỏng định lượng phân bố chất lượng nước trên diện rộng. Kết quả nghiên cứu cung cấp cơ sở khoa học chính xác với độ tin cậy thống kê cao, giúp các cơ quan quản lý nâng cao năng lực cảnh báo sớm, tối ưu hóa chi phí giám sát môi trường nước mặt và hỗ trợ quy hoạch sử dụng nguồn nước bền vững cho hơn 17,3 triệu người dân trong vùng.

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 ứng dụng nền tảng lý thuyết quang học thủy văn của Mobley kết hợp định luật truyền xạ bức xạ mặt trời qua môi trường nước. Theo nguyên lý viễn thám quang học, bức xạ rời khỏi mặt nước mang thông tin phản ánh các thành phần vật chất lơ lửng và độ đục thông qua tỷ số phản xạ viễn thám $R_{rs} = L_w / E_d$ (trong đó $L_w$ là bức xạ rời khỏi mặt nước và $E_d$ là bức xạ tới mặt nước). Nước sông chứa nhiều hạt vật chất vô cơ và hữu cơ có xu hướng phản xạ mạnh trong dải sóng nhìn thấy (0,4 - 0,7 µm), đặc biệt là kênh lam (0,4 - 0,5 µm) và kênh lục (0,5 - 0,6 µm), đồng thời bị hấp thụ gần như hoàn toàn ở dải cận hồng ngoại (0,8 - 1,1 µm). Độ thấu quang của cột nước thay đổi rõ rệt theo bước sóng: đạt tới 10m ở dải sóng 0,5 - 0,6 µm nhưng giảm xuống dưới 10cm ở dải sóng 0,8 - 1,1 µm.

Bên cạnh mô hình truyền xạ quang học, nghiên cứu sử dụng lý thuyết phân tích thống kê suy diễn và mô hình hồi quy đa biến nhằm thiết lập mối liên hệ định lượng giữa biến phụ thuộc (nồng độ TSS tính bằng mg/l và độ đục tính bằng NTU) với biến độc lập là giá trị phản xạ phổ của các pixel ảnh vệ tinh. Các khái niệm cốt lõi được chuẩn hóa trong nghiên cứu bao gồm:

  • Tổng chất rắn lơ lửng (TSS): tập hợp các hạt có kích thước lớn hơn 2 micron lơ lửng trong nước, xác định theo khối lượng khô trên một đơn vị thể tích mẫu.
  • Độ đục: thước đo mức độ tán xạ ánh sáng của các hạt keo và trầm tích trong nước, đo đạc theo tiêu chuẩn quốc tế ISO 7027:1999 và TCVN 6184:2008 với đơn vị NTU.
  • Trường nhìn tức thì (IFOV) và độ phân giải không gian: tham số đặc trưng cho kích thước điểm ảnh bề mặt (pixel) của đầu thu vệ tinh.

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

Nghiên cứu khai thác nguồn dữ liệu ảnh vệ tinh đa thời gian Landsat 7 ETM+ và Landsat 8 OLI/TIRS có độ phân giải không gian 30m, được thu thập từ Cục Khảo sát Địa chất Hoa Kỳ (USGS) qua cổng EarthExplorer trong giai đoạn 2014 - 2016. Cỡ mẫu nghiên cứu thực địa bao gồm 19 vị trí giám sát chiến lược (ký hiệu từ M1 đến M19) phân bố trên dòng chính sông Hậu, sông Tiền, sông Cổ Chiên và 2 trạm khí tượng thủy văn Cần Thơ, Mỹ Thuận. Phương pháp chọn mẫu là chọn mẫu có chủ đích kết hợp phân tầng không gian tại các điểm xung yếu: cửa lấy nước nhà máy nước sinh hoạt (Trà Ôn, Vĩnh Long), ranh giới hành chính liên tỉnh, ngã ba sông giao thoa dòng chảy và khu vực lồng bè nuôi trồng thủy sản tập trung.

Phương pháp lấy mẫu nước tại hiện trường được thực hiện đồng bộ với thời gian vệ tinh chụp qua khu vực nghiên cứu, chia thành 3 đợt quan trắc định kỳ hàng năm vào tháng 3 (cao điểm mùa khô), tháng 6 (đầu mùa mưa) và tháng 9 (đỉnh lũ). Toàn bộ mẫu nước được bảo quản và phân tích trong phòng thí nghiệm theo đúng quy chuẩn kỹ thuật quốc gia. Lý do lựa chọn phương pháp phân tích hồi quy thực nghiệm trên ảnh Landsat là nhờ tính khả thi cao, chi phí thấp, chu kỳ lặp lại 16 ngày và khả năng nội suy số liệu chính xác tại những khu vực sông rộng lớn, địa hình hiểm trở mà các trạm quan trắc mặt đất truyền thống khó tiếp cận thường xuyên.

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

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

Phân tích thống kê đối sánh giữa giá trị phổ phản xạ của ảnh Landsat và số liệu thực nghiệm đã mang lại những phát hiện quan trọng về mối tương quan và quy luật phân bố ô nhiễm:

  • Mô hình toán học phản ánh độ đục đạt độ chính xác rất cao với hệ số tương quan $R = 0,923$, khẳng định mối liên hệ chặt chẽ giữa dải phổ bức xạ rời mặt nước với mật độ hạt tán xạ. Phương trình hồi quy đối với chỉ tiêu TSS cũng ghi nhận hệ số tương quan thuận lợi, cho phép nội suy nồng độ chất lơ lửng từ ảnh vệ tinh.
  • Phân bố không gian của TSS và độ đục thể hiện rõ xu hướng tăng dần từ thượng nguồn xuống hạ lưu và khu vực cửa sông ven biển (như cửa Định An, Trần Đề, Cung Hầu), sau đó có xu hướng giảm nhanh khi ra xa bờ do quá trình hòa loãng với khối nước biển đại dương.
  • Kết quả tính toán nồng độ và đối chiếu với Quy chuẩn kỹ thuật quốc gia về chất lượng nước mặt (QCVN 08-MT:2015/BTNMT) cho thấy: trong năm 2014 và năm 2015 chỉ ghi nhận 2 điểm quan trắc vượt ngưỡng giới hạn cho phép; tuy nhiên đến năm 2016, số lượng điểm vượt quy chuẩn đã tăng gấp đôi lên 4 điểm quan trắc.
  • Nồng độ bùn cát và TSS biến động mạnh mẽ theo mùa: trong mùa lũ (tháng 7 đến tháng 11), lượng dòng chảy chiếm 70% đến 85% tổng lượng cả năm khiến độ đục trung bình trên sông chính tăng cao từ 500 đến 1.660 g/m³, trong khi vào mùa kiệt (tháng 12 đến tháng 6 năm sau) dòng chảy suy giảm chỉ còn 15% đến 30% làm độ đục tại các kênh rạch nội đồng hạ xuống dưới 50 g/m³.

Thảo luận kết quả

Sự gia tăng đột biến của số điểm vượt quy chuẩn môi trường trong năm 2016 bắt nguồn trực tiếp từ hiện tượng El Nino lịch sử, gây hạn hán khốc liệt và khiến lưu lượng dòng chảy mùa kiệt của sông Mê Kông xuống mức thấp nhất trong nhiều thập kỷ (dòng chảy trung bình mùa kiệt tại Tân Châu và Châu Đốc chỉ đạt khoảng 1.700 m³/s đến 2.600 m³/s). Lưu lượng nước ngọt từ thượng nguồn giảm mạnh đã tạo điều kiện cho thủy triều biển Đông (biên độ triều lớn từ 2,5m đến 3,5m) đẩy mặn và xâm nhập sâu vào nội địa, khuấy động trầm tích đáy và làm gia tăng độ đục tầng mặt tại các đoạn sông hạ lưu.

Dữ liệu mô phỏng trong luận văn được thể hiện trực quan qua hệ thống bản đồ phân bố không gian GIS đa thời gian trên 9 cảnh ảnh Landsat (ngày 21/3/2014, 22/6/2014, 18/9/2014; ngày 13/3/2015, 17/6/2015, 22/9/2015; ngày 23/3/2016, 16/6/2016, 19/9/2016) kết hợp các bảng thống kê phân cấp chất lượng nước theo chuẩn màu quốc tế. Khi so sánh với các công trình quốc tế như nghiên cứu của Doxaran trên cửa sông Macrotidal tại Pháp (độ không đảm bảo dao động 18% đến 22%) hay nghiên cứu của Nguyễn Quốc Phi tại khu vực ven biển Cửa Đáy (hệ số hồi quy $R = 0,955$), mô hình của luận văn đạt độ tương thích hoàn toàn tương đương, chứng minh thuật toán viễn thám có khả năng phản ánh khách quan hiện trạng thủy văn châu thổ.

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

Dựa trên kết quả định lượng không gian và thời gian của độ đục và TSS, các giải pháp quản trị môi trường nước mặt cần được triển khai đồng bộ:

  1. Xây dựng hệ thống giám sát chất lượng nước mặt tự động bằng viễn thám: Bộ Tài nguyên và Môi trường cùng Sở Tài nguyên và Môi trường 13 tỉnh thành ĐBSCL cần tích hợp công nghệ ảnh vệ tinh Landsat 8/9 và Sentinel-2 vào quy trình quan trắc định kỳ từ giai đoạn 2024 - 2028, hướng tới mục tiêu phủ sóng giám sát 100% diện tích mặt nước lưu vực sông Tiền và sông Hậu với chu kỳ cập nhật 5 đến 10 ngày/lần.
  2. Kiểm soát chặt chẽ hoạt động khai thác cát và nạo vét luồng lạch: Ủy ban nhân dân các tỉnh An Giang, Đồng Tháp và Vĩnh Long cần thiết lập hành lang bảo vệ bờ sông tại 2 đoạn uốn khúc có nguy cơ xói lở nghiêm trọng nhất (Hồng Ngự - Sa Đéc và Châu Đốc - Long Xuyên), đặt mục tiêu cắt giảm ít nhất 40% các điểm nạo vét không phép trước năm 2026 nhằm hạn chế xáo trộn trầm tích đáy làm gia tăng nồng độ TSS.
  3. Điều chỉnh cơ cấu mùa vụ nông nghiệp và công trình tích trữ nước ngọt: Bộ Nông nghiệp và Phát triển nông thôn cùng các địa phương cần điều chỉnh lịch thời vụ cho 1,2 triệu ha đất phù sa ngọt và 2,5 triệu ha đất phèn, đất mặn; chủ động vận hành hệ thống cống ngăn mặn và hồ trữ nước ngọt trong mùa khô (từ tháng 12 đến tháng 4) để giảm thiểu 25% nhu cầu khai thác nước thô khi sông có độ đục và độ mặn cao.
  4. Nâng cấp mạng lưới 19 trạm kiểm chứng thực địa: Đầu tư bổ sung thiết bị đo độ đục quang học trực tuyến tại các trạm khí tượng thủy văn Cần Thơ, Mỹ Thuận và các trạm cửa sông ven biển đến năm 2027, cung cấp chuỗi số liệu nền liên tục để hiệu chỉnh độ chính xác của các thuật toán giải đoán viễn thám.

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

Nội dung và kết quả thực nghiệm của luận văn mang lại giá trị ứng dụng thiết thực cho nhiều nhóm chuyên môn:

  • Cơ quan quản lý nhà nước về môi trường và tài nguyên nước: Lãnh đạo và chuyên viên tại Cục Quản lý Tài nguyên Nước, Chi cục Bảo vệ Môi trường các tỉnh ĐBSCL có thể sử dụng quy trình xử lý ảnh và phương trình hồi quy để xây dựng bản đồ cảnh báo ô nhiễm nguồn nước mặt với chi phí vận hành thấp hơn 60% so với phương pháp đo thủ công.
  • Nhà nghiên cứu, giảng viên và nghiên cứu sinh chuyên ngành GIS - Viễn thám: Tài liệu cung cấp cơ sở phương pháp luận chi tiết về chuyển đổi bức xạ phổ, xử lý nhiễu khí quyển và thuật toán tích hợp số liệu thực địa trong đánh giá môi trường nước nội địa.
  • Kỹ sư vận hành nhà máy nước và quản lý vùng nuôi trồng thủy sản: Cung cấp dữ liệu dự báo biến động độ đục theo mùa tại các điểm lấy nước thô (như Trà Ôn, Vĩnh Long), giúp chủ động tính toán liều lượng hóa chất keo tụ (phèn nhôm, PAC) và điều tiết nước cấp an toàn cho ao nuôi tôm, cá.
  • Các tổ chức phi chính phủ và viện nghiên cứu biến đổi khí hậu: Nguồn tài liệu tham khảo giá trị về tác động tích hợp của thủy điện thượng nguồn Mê Kông và hiện tượng El Nino đối với cán cân trầm tích châu thổ.

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

Tại sao ảnh vệ tinh Landsat được lựa chọn để đánh giá độ đục và TSS thay vì các phương pháp viễn thám khác?
Ảnh Landsat 7 và Landsat 8 cung cấp độ phân giải không gian 30m hoàn toàn phù hợp với độ rộng lòng sông Tiền và sông Hậu (thường rộng từ 500m đến hơn 2.000m). Dữ liệu được cung cấp miễn phí với chuỗi ảnh lịch sử liên tục từ USGS, giúp tối ưu hóa chi phí nghiên cứu và cho phép theo dõi diễn biến môi trường trong nhiều năm.

Mô hình hồi quy trong luận văn đạt độ chính xác như thế nào khi đối soát với thực tế?
Phương trình tương quan độ đục đạt hệ số $R = 0,923$, phản ánh sự trùng khớp cao giữa giá trị số của kênh phổ và phép đo thực địa tại 19 trạm quan trắc. Sai số tính toán nằm trong ngưỡng cho phép của kỹ thuật viễn thám quang học quốc tế (dưới 20%), đảm bảo độ tin cậy để lập bản đồ phân vùng ô nhiễm.

Đâu là trở ngại lớn nhất khi ứng dụng viễn thám quang học quan trắc nước sông tại ĐBSCL?
Hạn chế lớn nhất là độ che phủ của mây trong mùa mưa lũ (tháng 5 đến tháng 11) và độ thấu quang của nước bị suy giảm nhanh ở dải sóng cận hồng ngoại (dưới 10cm đối với bước sóng 0,8 - 1,1 µm), đòi hỏi phải kết hợp lọc mây kỹ thuật số và đồng bộ hóa thời gian lấy mẫu mặt đất chuẩn xác.

Các đập thủy điện thượng nguồn ảnh hưởng cụ thể thế nào đến độ đục nước sông tại ĐBSCL?
94 đập thủy điện giữ lại phần lớn bùn cát hạt thô, làm tổng lượng phù sa về hạ lưu giảm hơn 53% (từ 160 triệu tấn xuống 75 triệu tấn). Dòng chảy bị chặn làm biến đổi quy luật bồi lắng tự nhiên, gia tăng hiện tượng xói lở lòng dẫn và làm suy giảm độ đục tự nhiên có lợi cho nông nghiệp.

Kết quả nghiên cứu có thể chuyển giao trực tiếp cho các địa phương trong vùng không?
Toàn bộ quy trình giải đoán ảnh, tọa độ 19 trạm quan trắc chuẩn VN-2000 và các phương trình hồi quy có thể tích hợp trực tiếp vào phần mềm GIS của các Sở Tài nguyên và Môi trường để tự động hóa công tác lập báo cáo hiện trạng môi trường nước mặt định kỳ.

Kết luận

Luận văn đã giải quyết trọn vẹn mục tiêu khoa học và thực tiễn thông qua các kết quả then chốt:

  • Thiết lập thành công mô hình hồi quy viễn thám xác định độ đục và tổng chất rắn lơ lửng (TSS) trên hệ thống sông ĐBSCL với hệ số tương quan nổi bật đạt $R = 0,923$.
  • Xây dựng chuỗi bản đồ phân bố không gian - thời gian chất lượng nước chi tiết qua 9 thời điểm trong giai đoạn 2014 - 2016, phản ánh chân thực ảnh hưởng của hạn hán El Nino năm 2016 làm gia tăng gấp đôi số điểm vượt chuẩn QCVN 08-MT:2015/BTNMT.
  • Chứng minh tính ưu việt của công nghệ viễn thám trong việc quan trắc diện rộng, giảm thiểu chi phí và cho phép nội suy dữ liệu tại các vị trí khó tiếp cận dọc sông Tiền, sông Hậu và sông Cổ Chiên.
  • Đề xuất hệ thống giải pháp liên ngành đồng bộ từ quy hoạch giám sát tự động giai đoạn 2024 - 2028, quản lý khai thác cát đến thích ứng sản xuất nông nghiệp mùa khô.
  • Mở ra hướng nghiên cứu tiếp theo về việc tích hợp ảnh vệ tinh siêu phổ và dữ liệu radar (SAR) để quan trắc chất lượng nước xuyên suốt cả trong điều kiện thời tiết nhiều mây mùa lũ.

Các cơ quan quản lý, viện nghiên cứu và doanh nghiệp quan tâm có thể khai thác phương pháp luận và cơ sở dữ liệu của luận văn để nâng cao hiệu quả quản trị tài nguyên nước và phát triển bền vững lưu vực Đồng bằng sông Cửu Long.