Tổng quan nghiên cứu

Lưu vực sông Hồng là một trong 10 hệ thống sông có tải lượng dòng chảy và phù sa lớn nhất thế giới, cung cấp nguồn nước và bồi đắp dinh dưỡng thiết yếu cho hơn 14 triệu người tại vùng Đồng bằng Bắc Bộ. Tổng lượng phù sa trung bình nhiều năm đo tại trạm thủy văn Sơn Tây đạt khoảng 114 triệu tấn/năm. Tuy nhiên, hơn 50% nguồn nước sông Hồng bắt nguồn từ ngoài lãnh thổ Việt Nam, khiến chất lượng nước và chế độ dòng chảy vô cùng nhạy cảm với các hoạt động khai thác thượng nguồn. Sự phát triển mạnh mẽ của hơn 29 công trình đập thủy điện bậc thang phía Trung Quốc, cùng mạng lưới 105 điểm khai thác cát quy mô lớn và hàng chục nghìn mét khối nước thải công nghiệp chưa qua xử lý xả trực tiếp mỗi ngày đã làm biến đổi sâu sắc động lực học bùn cát.

Phương pháp quan trắc truyền thống bằng lấy mẫu thủ công định kỳ tồn tại nhiều hạn chế lớn về mặt chi phí, thời gian và phạm vi không gian đại diện. Nhằm khắc phục khoảng trống dữ liệu này, nghiên cứu đặt mục tiêu xây dựng thuật toán định lượng hàm lượng phù sa lơ lửng (SSC) dựa trên tư liệu ảnh vệ tinh Sentinel-2A. Phạm vi nghiên cứu tập trung dọc theo 270 km dòng chính sông Thao từ Hà Khẩu (Lào Cai) đến Việt Trì (Phú Thọ) trong giai đoạn 2017 - 2018.

Ý nghĩa cốt lõi của nghiên cứu là xây dựng quy trình giám sát chất lượng nước liên tục với độ phân giải không gian 10 m và chu kỳ chụp lặp từ 5 đến 10 ngày. Ứng dụng mô hình viễn thám quang học này giúp giảm hơn 60% chi phí vận hành so với trạm đo cố định, đồng thời mang lại cơ sở khoa học định lượng phục vụ công tác quản lý tài nguyên nước xuyên biên giới và cảnh báo sớm biến động môi trườ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 được xây dựng trên nền tảng lý thuyết quang học thủy văn (Hydrological Optics) kết hợp với các nguyên lý thủy - thạch động lực học dòng chảy sông ngòi. Theo lý thuyết quang học thủy văn, phổ phản xạ rời mặt nước là một hàm số phụ thuộc vào các hợp phần quang hoạt tính trong cột nước, bao gồm vật chất lơ lửng vô cơ (ISS), vật chất lơ lửng hữu cơ (OSS), độ thấu quang và sắc tố diệp lục. Các hạt phù sa vô cơ trong nước sông tán xạ mạnh bức xạ mặt trời ở dải sóng ánh sáng nhìn thấy và cận hồng ngoại.

Bên cạnh đó, lý thuyết cân bằng trầm tích sông ngòi giải thích cơ chế bồi xói, lắng đọng và vận chuyển vật chất rắn dọc theo độ dốc địa hình thung lũng đới đứt gãy sông Hồng. Ba khái niệm chuyên ngành then chốt được chuẩn hóa trong nghiên cứu gồm:

  • Hàm lượng phù sa lơ lửng (SSC): Khối lượng chất rắn lơ lửng trên một đơn vị thể tích nước, tính bằng đơn vị mg/L.
  • Độ thấu quang đĩa Secchi (SD): Chỉ số biểu thị độ trong quang học của cột nước, đo bằng chiều sâu biến mất của đĩa Secchi chuẩn 20 cm.
  • Tỷ số kênh phổ phản xạ (Band Ratio): Tỷ lệ phản xạ giữa dải sóng đỏ, rìa đỏ (Red Edge) và xanh lục nhằm chuẩn hóa các tác động nhiễu của khí quyển và bề mặt gợn sóng.

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

Nguồn dữ liệu thực nghiệm được thu thập từ 3 chiến dịch khảo sát thực địa với tổng cỡ mẫu là 77 điểm đo phân bố trên toàn tuyến sông:

  • Đợt 1 (ngày 11/11/2017 - đại diện mùa khô) gồm 13 điểm đo tại khu vực hợp lưu ngã ba sông Lô - Đà - Thao.
  • Đợt 2 (ngày 23/03/2018 và 01/04/2018 - thời kỳ chuyển mùa khô sang mưa) gồm 40 điểm đo từ cửa khẩu Hà Khẩu đến cầu Giang Đông (Lào Cai).
  • Đợt 3 (ngày 29/07/2018 và 10/08/2018 - cao điểm mùa mưa) gồm 24 điểm đo từ cửa khẩu Lào Cai đến cầu Phố Mới.

Phương pháp chọn mẫu áp dụng là chọn mẫu ngẫu nhiên phân tầng theo trắc diện ngang và dọc dòng chảy, ưu tiên các đoạn sông hợp lưu, khu vực khai thác khoáng sản và cửa xả nước thải lớn nhằm đảm bảo tính đại diện quang học. Tại mỗi điểm đo, phổ phản xạ mặt nước được thu nhận bằng máy đo bức xạ hiện trường GER 1500 (dải phổ 350 - 1050 nm, độ phân giải 1,5 nm) theo góc nghiêng 40 - 45 độ chuẩn Mobley. Mẫu nước mặt được lấy ở độ sâu biểu kiến, bảo quản lạnh ở 4°C và phân tích trong phòng thí nghiệm theo tiêu chuẩn APHA 2540D bằng màng lọc Whatman GF/F 47 mm, sấy 104°C để xác định SSC và nung 500°C trong 4 giờ để xác định hàm lượng chất vô cơ lơ lửng (ISS).

Dữ liệu ảnh vệ tinh Sentinel-2A đa phổ được xử lý hiệu chỉnh khí quyển bằng phương pháp trừ vật tối (DOS) và trích xuất mặt nước qua tỷ số kênh hồng ngoại sóng ngắn B11/B3. Phân tích hồi quy phi tuyến hàm mũ giữa dữ liệu phổ mặt đất và tỷ số kênh ảnh được lựa chọn vì tính ưu việt trong việc khắc phục hiện tượng bão hòa phổ khi nồng độ phù sa vượt ngưỡng 200 mg/L.

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

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

Mô hình thực nghiệm thiết lập mối tương quan chặt chẽ giữa tỷ số phổ phản xạ kết hợp kênh 4 (665 nm), kênh 5 (705 nm) và kênh 3 (560 nm) của ảnh Sentinel-2A với hàm lượng SSC thực tế. Phương trình hồi quy dạng hàm mũ đạt hệ số xác định R² > 0,85 và sai số toàn phương trung bình RMSE dưới 15% trên tập dữ liệu 45 điểm kiểm chứng, cho phép định lượng chính xác dải SSC từ 15 mg/L đến trên 600 mg/L.

Sự biến động hàm lượng phù sa lơ lửng thể hiện tính phân hóa sâu sắc theo mùa:

  • Trong mùa mưa (tháng 6 đến tháng 10), hàm lượng SSC trung bình toàn tuyến dao động từ 250 mg/L đến 550 mg/L, cục bộ tại các đỉnh lũ có thể vượt 650 mg/L.
  • Trong mùa khô (tháng 11 đến tháng 5 năm sau), hàm lượng SSC giảm mạnh xuống ngưỡng 40 mg/L đến 120 mg/L. Mức chênh lệch tải lượng phù sa giữa mùa mưa và mùa khô đạt từ 2,5 đến 3,2 lần.

Phân bố không gian của SSC có xu hướng giảm dần dọc theo hành trình dòng chảy từ Tây Bắc xuống Đông Nam:

  • Đoạn sông Thao chảy qua Lào Cai ghi nhận mức SSC cao nhất trên toàn lưu vực, đóng góp tới 78% tổng lượng vật chất lơ lửng từ vùng thượng nguồn bị xói mòn mạnh.
  • Khi dòng chảy qua địa phận Yên Bái, hàm lượng SSC giảm trung bình khoảng 25% đến 35% do địa hình mở rộng và các nhánh sông phụ lưu pha loãng.
  • Đến đoạn qua tỉnh Phú Thọ, hàm lượng SSC tiếp tục giảm thêm 15% đến 20% do quá trình bồi lắng lòng dẫn trước khi hòa nhập vào khu vực hợp lưu Việt Trì.

Hoạt động nhân sinh tạo ra các dị thường cục bộ rõ nét:

  • Tại các mỏ khai thác cát ở huyện Bảo Thắng (công suất 80.000 m³/năm) và thành phố Lào Cai (công suất 60.000 m³/năm), hàm lượng SSC tức thời tại hạ lưu mỏ tăng từ 40% đến 60% so với đoạn thượng lưu liền kề.
  • Tại các điểm xả thải công nghiệp ven sông ở Lâm Thao và Việt Trì với lưu lượng xả từ 3.800 đến 4.500 m³/ngày, nồng độ chất rắn lơ lửng hữu cơ và vô cơ tăng đột biến, làm giảm độ thấu quang đĩa Secchi từ 45 cm xuống dưới 15 cm.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự tập trung phù sa cực cao tại Lào Cai xuất phát từ đặc điểm địa hình núi cao chia cắt mạnh, độ dốc lớn kết hợp với lớp vỏ phong hóa nhiệt đới dày trên nền đá biến chất yếu. Khi xuất hiện các đợt mưa lớn với lượng mưa từ 1.400 mm đến 2.200 mm/năm, hiện tượng xói mòn bề mặt và trượt lở đất diễn ra dữ dội, cuốn trôi một khối lượng đất đá khổng lồ vào lòng sông chính.

Các dữ liệu phân tích được mô hình hóa trực quan qua hệ thống bản đồ chuyên đề đa thời gian trên nền ảnh Sentinel-2A tổ hợp màu thực và màu giả. Biểu đồ phân tán hồi quy minh chứng rõ ràng sự tương thích cao giữa số liệu tính toán từ viễn thám và số liệu đo đạc thực tế tại trạm quan trắc.

So sánh với các nghiên cứu quốc tế trên lưu vực sông Amazon hay sông Rhone, thuật toán tổ hợp đa kênh (Đỏ + Rìa đỏ / Xanh lục) trong luận văn này đã giải quyết triệt để sự sai lệch do ảnh hưởng của màu nước nền hữu cơ, mang lại độ tin cậy vượt trội so với các mô hình đơn kênh truyền thống. Kết quả này khẳng định sông Thao giữ vai trò chủ đạo tuyệt đối trong việc cấp phù sa cho toàn bộ vùng hạ du sông Hồng, ngay cả khi các hồ chứa thủy điện lớn đã giữ lại tới 50% đến 70% lượng trầm tích trên lưu vực sông Đà.

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

  1. Triển khai hệ thống tự động hóa giám sát chất lượng nước mặt bằng ảnh vệ tinh Sentinel-2: Bộ Tài nguyên và Môi trường cùng Sở Tài nguyên và Môi trường 3 tỉnh Lào Cai, Yên Bái, Phú Thọ cần tích hợp thuật toán xử lý ảnh viễn thám vào quy trình quan trắc định kỳ 10 ngày/lần. Mục tiêu hướng tới nâng cao 80% độ bao phủ không gian giám sát và tiết kiệm ít nhất 40% ngân sách đo kiểm thực địa định kỳ trong giai đoạn 2026 - 2028.

  2. Quy hoạch và siết chặt cấp phép khai thác cát lòng sông: Ủy ban nhân dân các tỉnh Lào Cai và Phú Thọ cần tiến hành rà soát toàn diện 105 điểm khai thác cát sỏi hiện hữu, kiên quyết thu hồi giấy phép đối với các mỏ vi phạm hành lang bảo vệ đê điều. Yêu cầu lắp đặt thiết bị định vị và camera giám sát tải trọng tại 100% tàu hút cát, cắt giảm 30% công suất khai thác vào mùa kiệt trước năm 2027 để hạn chế xói lở bờ sông.

  3. Nâng cấp quy chuẩn xử lý nước thải công nghiệp tại các lưu vực trọng điểm: Ban quản lý các khu công nghiệp tại Phú Thọ và Lào Cai cần buộc các nhà máy sản xuất giấy, hóa chất và dệt may có lưu lượng xả thải trên 1.000 m³/ngày phải lắp đặt trạm quan trắc tự động liên tục truyền dữ liệu về cơ quan quản lý. Đảm bảo 100% nước thải xả ra sông Hồng đạt quy chuẩn QCVN 40:2011/BTNMT trước quý IV năm 2027.

  4. Thiết lập cơ chế hợp tác chia sẻ dữ liệu thủy văn xuyên biên giới: Ủy ban sông Mê Kông - sông Hồng phối hợp với các cơ quan ngoại giao thúc đẩy ký kết hiệp định chia sẻ số liệu vận hành xả lũ và điều tiết bùn cát của 29 đập thủy điện thượng nguồn với các quốc gia láng giềng. Mục tiêu duy trì lưu lượng dòng chảy mùa cạn tối thiểu đạt 745 m³/s tại trạm thủy văn khống chế nhằm bảo tồn hệ sinh thái hạ du trước năm 2028.

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

  • Cơ quan quản lý nhà nước về môi trường và tài nguyên nước: Cán bộ thuộc Cục Quản lý Tài nguyên nước, Tổng cục Môi trường và các Sở Tài nguyên và Môi trường địa phương có thể ứng dụng trực tiếp bản đồ phân vùng SSC và cơ sở dữ liệu viễn thám để thẩm định báo cáo đánh giá tác động môi trường (ĐTM), cấp phép khai thác tài nguyên và kiểm soát các nguồn xả thải.

  • Kỹ sư và chuyên viên phân tích dữ liệu không gian (GIS & Remote Sensing): Các chuyên viên địa tin học có thể khai thác thuật toán tỷ số kênh phổ, quy trình hiệu chỉnh bức xạ khí quyển DOS và phương pháp trích xuất diện tích mặt nước để phát triển các ứng dụng webGIS theo dõi tài nguyên nước theo thời gian thực.

  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Môi trường, Địa lý: Tài liệu cung cấp khung phương pháp luận chuẩn xác về đo phổ phản xạ hiện trường bằng máy GER 1500, kỹ thuật phân tích mẫu phòng thí nghiệm theo chuẩn APHA và phương pháp xử lý thống kê đa biến phục vụ các đề tài nghiên cứu liên quan.

  • Ban quản lý các nhà máy nước sạch và công trình thủy lợi, thủy điện: Đội ngũ kỹ thuật có thể sử dụng các chỉ số dự báo độ đục và tải lượng phù sa theo mùa để lập kế hoạch vận hành lắng lọc, tối ưu hóa hóa chất xử lý nước và tính toán chu kỳ nạo vét bồi lắng lòng hồ chứa.

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

  1. Dữ liệu vệ tinh Sentinel-2A có ưu điểm gì vượt trội trong quan trắc phù sa lơ lửng sông ngòi? Vệ tinh Sentinel-2A sở hữu độ phân giải không gian cao 10 m ở các kênh phổ ánh sáng nhìn thấy và dải rìa đỏ (Red Edge), chu kỳ lặp 5 ngày khi kết hợp Sentinel-2B và hoàn toàn miễn phí từ ESA. Độ phân giải này rất thích hợp với các lòng sông hẹp từ 100 m đến 500 m như sông Thao.

  2. Phương trình tính toán SSC trong luận văn sử dụng tổ hợp kênh phổ nào? Nghiên cứu sử dụng tổ hợp hàm mũ dựa trên tỷ số giữa tổng phổ kênh Đỏ (B4 - 665 nm) cộng kênh Rìa đỏ (B5 - 705 nm) chia cho kênh Xanh lục (B3 - 560 nm). Tổ hợp này triệt tiêu hiện tượng phản xạ nền và đạt tương quan R² > 0,85 với sai số toàn phương trung bình RMSE dưới 15%.

  3. Tại sao hàm lượng phù sa sông Thao lại suy giảm đáng kể khi chảy qua Phú Thọ? Khi chảy về hạ lưu tỉnh Phú Thọ, độ dốc dòng sông giảm mạnh, lòng dẫn mở rộng và vận tốc dòng chảy suy giảm khiến các hạt phù sa có kích thước lớn nhanh chóng lắng đọng xuống đáy. Đồng thời, sự gia nhập của các nhánh sông có độ đục thấp hơn góp phần làm loãng nồng độ lơ lửng.

  4. Mối quan hệ giữa độ thấu quang đĩa Secchi và hàm lượng SSC được thể hiện thế nào? Độ thấu quang đĩa Secchi tỷ lệ nghịch phi tuyến với hàm lượng SSC. Trong mùa lũ, khi SSC vượt 400 mg/L, độ thấu quang ghi nhận dưới 15 cm; ngược lại trong mùa kiệt, khi SSC giảm xuống dưới 50 mg/L, độ thấu quang của nước sông tăng lên từ 40 cm đến 65 cm.

  5. Quy trình xử lý viễn thám trong nghiên cứu có thể áp dụng cho các sông khác ở Việt Nam không? Hoàn toàn có thể áp dụng. Tuy nhiên, do đặc tính thổ nhưỡng và thành phần quang hoạt tính của mỗi lưu vực khác nhau, cần tiến hành đo đạc hiệu chỉnh phổ mặt đất với tối thiểu 30 điểm mẫu thực địa để chuẩn hóa lại các hệ số trong phương trình hồi quy.

Kết luận

  • Xây dựng thành công thuật toán định lượng hàm lượng phù sa lơ lửng (SSC) cho sông Thao bằng tư liệu ảnh Sentinel-2A với độ chính xác cao (R² > 0,85, RMSE < 15%).
  • Xác định rõ quy luật phân hóa phù sa theo mùa với nồng độ mùa mưa (250 - 550 mg/L) cao gấp 2,5 đến 3,2 lần mùa khô (40 - 120 mg/L).
  • Chứng minh sông Thao đoạn qua Lào Cai là nguồn cung cấp phù sa chủ đạo cho hạ lưu sông Hồng, chiếm tới 78% tổng tải lượng xói mòn thượng nguồn.
  • Định vị chính xác các điểm nóng gia tăng độ đục cục bộ do hoạt động khai thác cát tại Bảo Thắng và xả thải công nghiệp tại Lâm Thao, Việt Trì.
  • Khẳng định tính khả thi của công nghệ viễn thám quang học trong việc thay thế và bổ trợ các phương pháp quan trắc thủy văn truyền thống.

Đóng góp lớn nhất của luận văn là hoàn thiện khung phương pháp luận viễn thám kết hợp thực địa để giám sát liên tục chất lượng môi trường nước trên các lưu vực sông xuyên biên giới. Trong giai đoạn 2026 - 2030, hướng phát triển tiếp theo cần mở rộng tích hợp ảnh vệ tinh Landsat 9 và dữ liệu siêu phổ của Sentinel-3 nhằm tự động hóa quy trình cảnh báo ô nhiễm nguồn nước theo thời gian thực. Các cơ quan quản lý và nhà khoa học hãy chủ động ứng dụng bộ công cụ viễn thám này vào công tác quản trị lưu vực sông để bảo vệ bền vững an ninh tài nguyên nước quốc gia.