Tổng quan nghiên cứu

Hệ thống sông ngòi Việt Nam sở hữu hơn 3.450 con sông, suối có chiều dài từ 10 km trở lên, phân bố trên 108 lưu vực với tổng lượng nước mặt trung bình hàng năm đạt khoảng 830 tỉ m³. Tuy nhiên, nguồn tài nguyên này đang đứng trước thách thức lớn khi khoảng 63% tổng lượng dòng chảy bắt nguồn từ ngoài biên giới quốc gia, khiến nguồn nước nội địa chịu áp lực nặng nề từ các hoạt động nhân sinh. Sông Gianh là hệ thống sông lớn nhất tỉnh Quảng Bình với diện tích lưu vực 4.680 km² và chiều dài dòng chính đạt 158 km, đóng vai trò biểu trưng địa lý và là nguồn cấp nước sống còn cho nông nghiệp, thủy sản và sinh hoạt của hàng trăm nghìn người dân.

Trong bối cảnh công nghiệp hóa, đô thị hóa gia tăng cùng những biến cố môi trường biển phức tạp, áp lực nước thải chưa qua xử lý đổ vào lưu vực sông Gianh ngày càng nghiêm trọng. Vấn đề nghiên cứu trọng tâm của luận văn là giải quyết khoảng trống dữ liệu quan trắc toàn diện, đánh giá đúng thực trạng ô nhiễm và các yếu tố chi phối động thái chất lượng nước mặt.

Mục tiêu cụ thể của công trình bao gồm việc xác định 13 thông số lý hóa, vi sinh quan trọng (nhiệt độ, pH, EC, TDS, DO, COD, BOD5, amoni, nitrat, photphat, tổng sắt tan, độ cứng, coliform), so sánh đối chiếu với Quy chuẩn Kỹ thuật Quốc gia QCVN 08-MT:2015/BTNMT, QCVN 39:2011/BTNMT và QCVN 01:2009/BYT. Nghiên cứu được triển khai đồng bộ qua 4 đợt quan trắc thực địa tại 5 vị trí trọng điểm trải dài từ thượng nguồn xã Thanh Hóa đến hạ lưu Cửa Gianh. Kết quả nghiên cứu cung cấp cơ sở khoa học tin cậy, giúp nâng cao độ chính xác trong công tác dự báo chất lượng nước lên trên 95% và hỗ trợ các nhà quản lý đưa ra giải pháp bảo vệ tài nguyên nước bền vững cho địa phươ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 dựa trên nền tảng lý thuyết về cân bằng vật chất thủy vực và động lực học phân hủy sinh - hóa học của các chất ô nhiễm trong môi trường nước mặt. Về mô hình tải lượng ô nhiễm, luận văn phân tích theo hai nhóm nguồn phát thải chính: nguồn điểm (Point Sources) bao gồm nước thải sinh hoạt đô thị với tỷ lệ chuyển hóa từ 70% đến 90% lượng nước cấp và nước thải công nghiệp từ các cơ sở chế biến; nguồn phân tán (Non-Point Sources) phát sinh do hiện tượng nước mưa chảy tràn qua các vùng sản xuất nông nghiệp và khai khoáng kéo theo chất dinh dưỡng dư thừa.

Các khái niệm then chốt được chuẩn hóa bao gồm: Nhu cầu oxy hóa học (COD) phản ánh lượng chất hữu cơ bị oxy hóa bởi chất oxy hóa mạnh; Nhu cầu oxy sinh hóa sau 5 ngày (BOD5) đo lường mức tiêu thụ oxy của vi sinh vật hiếu khí ở 20°C; Tổng chất rắn hòa tan (TDS) và độ dẫn điện (EC) biểu thị độ khoáng hóa và mức độ ion hóa trong nước; cùng hàm lượng các ion dưỡng chất như amoni, nitrat và photphat đại diện cho nguy cơ phú dưỡng hóa.

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

Nguồn dữ liệu thực nghiệm được thu thập từ cỡ mẫu gồm 20 mẫu nước tổ hợp, phân bố đều tại 5 vị trí quan trắc mang tính đại diện cao dọc lưu vực sông Gianh (ký hiệu từ SG1 đến SG5) trong 4 đợt khảo sát định kỳ 2 tháng một lần vào các tháng 2, 4, 6 và 8. Phương pháp chọn mẫu áp dụng kỹ thuật lấy mẫu tổ hợp hai bên bờ sông theo tỷ lệ thể tích 1:1 tại tầng mặt ở độ sâu 50 cm, sau đó bảo quản nghiêm ngặt ở nhiệt độ âm 20°C hoặc axit hóa mẫu theo Tiêu chuẩn Quốc gia TCVN nhằm ngăn chặn hiện tượng biến tính sinh hóa trước khi phân tích.

Hệ thống phương pháp phân tích phòng thí nghiệm được lựa chọn dựa trên độ nhạy cao và tính chuẩn mực quốc tế:

  • Phương pháp trắc quang bicromat xác định COD tại bước sóng 420 nm sau khi đun hồi lưu 2 giờ ở 150°C.
  • Phương pháp cấy ủ 5 ngày ở 20°C kết hợp đo điện cực màng quang học để xác định BOD5.
  • Phương pháp so màu phenat tạo phức indophenol xanh xác định amoni ở bước sóng 640 nm.
  • Phương pháp natri salixylat xác định nitrat ở bước sóng 415 nm.
  • Phương pháp xanh molypden xác định photphat ở bước sóng 880 nm.
  • Phương pháp tạo phức với 1,10-phenanthroline ở bước sóng 510 nm xác định tổng sắt tan.
  • Phương pháp chuẩn độ complexon III với chỉ thị murexit xác định độ cứng canxi ở pH từ 12 đến 13.
  • Phương pháp lên men nhiều ống nghiệm (MPN) trên môi trường chọn lọc để định lượng tổng coliform.

Toàn bộ quy trình phân tích được kiểm soát chất lượng (QA/QC) nghiêm ngặt với hệ số tương quan đường chuẩn R đạt từ 0,9996 đến 0,9999; giới hạn định lượng LOQ của COD đạt 4,74 mg/L và nitrat đạt 0,06 mg/L; độ lệch chuẩn tương đối RSD đều nhỏ hơn giá trị giới hạn tính theo hàm Horwitz (từ 1,13% đến 7,14%); độ thu hồi phép đo đạt mức tối ưu từ 95% đến 102%. Để đánh giá tác động đồng thời của vị trí không gian và thời gian quan trắc, tác giả sử dụng mô hình phân tích phương sai hai yếu tố (ANOVA 2 chiều) ở mức ý nghĩa thống kê p < 0,05 kết hợp tính toán độ lệch nhỏ nhất LSD.

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

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

Thứ nhất, mức độ khoáng hóa và độ dẫn điện có sự phân hóa không gian cực kỳ sâu sắc giữa các vùng nước. Độ dẫn điện EC dao động trong khoảng 300,09 đến 625,01 µS/cm và tổng chất rắn hòa tan TDS nằm trong khoảng 156,02 đến 331,10 mg/L. Tại trạm hạ lưu cửa sông Cảng Gianh (SG1), giá trị trung bình của TDS đạt 311,59 mg/L và EC đạt 589,21 µS/cm, cao hơn khoảng 88% về TDS và 85% về EC so với trạm thượng lưu Đồng Lê (SG5 có TDS trung bình 165,80 mg/L và EC 313,56 µS/cm).

Thứ hai, độ cứng của nước sông Gianh biến thiên từ 123,2 đến 218,6 mg/L, phân loại thuộc nhóm nước cứng vừa đến nước cứng. Điểm đặc biệt là nồng độ độ cứng trung bình tại thượng nguồn (SG5) đạt mức cao nhất lên tới 203,80 mg/L, cao hơn khoảng 39,6% so với trạm cửa sông SG1 (trung bình 145,90 mg/L).

Thứ ba, nồng độ các chất ô nhiễm hữu cơ và vi sinh biểu hiện xu hướng gia tăng cục bộ tại các khu vực tập trung dân cư và hoạt động thương mại. Hàm lượng COD dao động từ 4,80 đến 30,70 mg/L, BOD5 từ 2,50 đến 17,00 mg/L. Trạm cầu Gianh (SG2) ghi nhận giá trị COD cao nhất vào đợt 4 đạt 30,70 mg/L (vượt ngưỡng A1 của QCVN 08-MT:2015/BTNMT). Mật độ coliform dao động rộng từ 470 đến 7.900 MPN/100mL, với đỉnh điểm xuất hiện tại trạm SG3 và SG2 do nguồn xả thải sinh hoạt chưa qua xử lý.

Thứ tư, chỉ số oxy hòa tan DO duy trì ở mức tương đối tốt, dao động từ 3,50 đến 8,37 mg/L, trong khi pH dao động ổn định trong khoảng 7,10 đến 7,67, hoàn toàn thỏa mãn giới hạn cho phép của nguồn nước mặt loại A1 phục vụ mục đích bảo tồn thủy sinh và sinh hoạt.

Thảo luận kết quả

Dữ liệu thực nghiệm khi biểu diễn qua biểu đồ biến thiên không gian - thời gian và đồ thị phân tích phương sai đã phản ánh rõ nét cơ chế tác động của môi trường tự nhiên lẫn hoạt động nhân sinh. Kết quả kiểm định ANOVA 2 chiều cho thấy yếu tố vị trí địa lý có ý nghĩa thống kê vượt trội đối với EC và TDS với giá trị F tính toán đạt 108,70 (vượt xa giá trị F lý thuyết là 3,26 ở mức ý nghĩa 0,05). Sự gia tăng nồng độ khoáng hóa tại SG1 và SG2 được lý giải bởi tác động lan truyền của chế độ bán nhật triều không đều với biên độ triều cực đại 1,66 m đưa nước biển xâm nhập sâu vào lòng sông tới 40 - 60 km.

Ngược lại, hiện tượng độ cứng tăng vọt tại trạm thượng nguồn SG5 (F tính toán = 19,51 > F lý thuyết = 3,26) phản ánh đặc thù địa chất của vùng núi đá vôi kỷ Devonian và Carbon-Pecmi chiếm hơn 33% diện tích lưu vực sông Gianh. Quá trình phong hóa và hòa tan khoáng chất canxi, magie từ các khối núi đá vôi dạng địa hình karst đã làm giàu ion kim loại kiềm thổ vào dòng chảy thượng lưu.

Đối với các thông số ô nhiễm hữu cơ (COD, BOD5, Coliform), biến động thời gian giữa các đợt quan trắc cho thấy tải lượng ô nhiễm tăng cao vào các tháng mùa khô nắng nóng (tháng 6 và tháng 8) do lưu lượng dòng chảy suy giảm và nhiệt độ nước đạt mức 29,2 - 32,1°C thúc đẩy quá trình phân hủy kỵ khí cục bộ. Các kết quả này hoàn toàn tương thích với các nghiên cứu thủy văn hiện đại, khẳng định tính tất yếu của việc kiểm soát ô nhiễm nước thải đô thị thị xã Ba Đồn và các khu dân cư ven sông.

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

Nhằm kiểm soát ô nhiễm và nâng cao chất lượng môi trường nước lưu vực sông Gianh, bốn nhóm giải pháp hành động cụ thể được đề xuất:

  • Thiết lập mạng lưới trạm quan trắc tự động liên tục: Lắp đặt 3 trạm quan trắc tự động tại cửa sông Gianh (SG1), khu vực thị xã Ba Đồn (SG2) và thượng lưu Đồng Lê (SG5) để giám sát trực tuyến 24/7 các thông số pH, DO, TDS, EC và COD. Mục tiêu nâng cao độ bao phủ giám sát lên 100% các điểm xung yếu trước năm 2027, do Sở Tài nguyên và Môi trường tỉnh Quảng Bình chủ trì thực hiện.
  • Đầu tư và hoàn thiện hệ thống xử lý nước thải tập trung: Xây dựng mới 2 nhà máy xử lý nước thải sinh hoạt đô thị tại thị xã Ba Đồn và huyện Quảng Trạch với công suất thiết kế 10.000 m³/ngày-đêm, ứng dụng công nghệ lọc sinh học hiếu khí kết hợp hồ sinh học nhân tạo nhằm cắt giảm 85% tải lượng COD và BOD5 xả vào lòng sông trước năm 2028, do Ủy ban nhân dân các huyện thị phối hợp nhà đầu tư triển khai.
  • Quy hoạch hành lang bảo vệ nguồn nước và phân vùng xả thải: Ban hành quy định phân vùng tiếp nhận nước thải dựa trên sức chịu tải của từng đoạn sông, bảo vệ nghiêm ngặt vùng đệm karst Vườn Quốc gia Phong Nha - Kẻ Bàng với diện tích lưu vực 4.680 km². Tăng cường thanh tra, nâng mức xử phạt vi phạm hành chính đối với các cơ sở chế biến xả thải chưa đạt chuẩn cột A theo quy chuẩn môi trường từ năm 2026.
  • Đẩy mạnh tuyên truyền giáo dục và chuyển đổi nông nghiệp sinh thái: Tổ chức 30 lớp tập huấn hàng năm cho cộng đồng dân cư ven sông về quản lý chất thải rắn và nước thải chăn nuôi, giảm thiểu 30% lượng phân bón vô cơ và thuốc bảo vệ thực vật nhằm hạn chế ô nhiễm diện trước năm 2027, do Chi cục Trồng trọt và Bảo vệ Thực vật cùng Hội đồng nhân dân các cấp chỉ đạo.

Đố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 và học thuật cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về môi trường và tài nguyên nước: Điển hình là Sở Tài nguyên và Môi trường, Chi cục Thủy lợi và Quản lý Môi trường tỉnh Quảng Bình. Luận văn cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác để xây dựng quy hoạch tài nguyên nước mặt giai đoạn 2026-2035, cấp phép xả thải công bằng và lập kế hoạch ứng phó biến đổi khí hậu, xâm nhập mặn.
  • Giảng viên, học viên cao học và sinh viên ngành Hóa học, Môi trường: Tài liệu là cẩm nang hữu ích về quy trình thẩm định phương pháp phân tích (QA/QC), cách xác định đường chuẩn, LOD, LOQ theo chuẩn AOAC và kỹ thuật ứng dụng mô hình ANOVA hai chiều để xử lý số liệu môi trường phức tạp.
  • Các doanh nghiệp và hợp tác xã nuôi trồng thủy sản, sản xuất nước sạch: Cung cấp số liệu chi tiết về biến động độ cứng, TDS, DO theo không gian và mùa vụ để doanh nghiệp chủ động điều chỉnh quy trình lấy nước thô, giảm thiểu rủi ro sốc mặn hoặc phú dưỡng trong nuôi thủy sản lồng bè.
  • Các viện nghiên cứu và tổ chức phi chính phủ về bảo tồn thiên nhiên: Kế thừa dữ liệu phân tích để xây dựng các dự án bảo tồn đa dạng sinh học lưu vực sông, bảo vệ hệ sinh thái rừng ngập mặn và hệ thống hang động karst đặc sắc của di sản thế giới.

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

Tại sao độ cứng của nước sông Gianh lại có xu hướng cao nhất ở vùng thượng lưu?

Hàm lượng độ cứng tại trạm thượng nguồn SG5 đạt mức trung bình 203,80 mg/L do đặc trưng địa hình lưu vực có hơn 33% diện tích tự nhiên là núi đá vôi kỷ Devonian và Carbon-Pecmi. Nước mưa và dòng ngầm chảy qua các khối karst hòa tan một lượng lớn khoáng chất canxi và magie trước khi đổ vào dòng chính, tạo nên đặc tính nước cứng tự nhiên cho khu vực thượng nguồn.

Sự gia tăng tổng chất rắn hòa tan (TDS) và độ dẫn điện (EC) ở hạ lưu phản ánh điều gì?

Giá trị TDS đạt trung bình 311,59 mg/L và EC đạt 589,21 µS/cm tại trạm SG1 Cảng Gianh phản ánh quá trình xâm nhập mặn rõ rệt từ biển Đông. Do sông Gianh chịu chế độ bán nhật triều với biên độ triều lên tới 1,66 m, nước biển chứa nồng độ ion hòa tan cao khuếch tán ngược vào lòng sông khoảng 40 - 60 km, làm tăng mạnh độ khoáng hóa ở hạ lưu.

Mô hình phân tích ANOVA 2 chiều đóng góp vai trò gì trong việc đánh giá số liệu?

Mô hình phân tích phương sai ANOVA hai chiều với mức ý nghĩa p = 0,05 giúp phân tách chính xác ảnh hưởng của không gian (5 trạm quan trắc) và thời gian (4 đợt khảo sát). Kết quả kiểm định thống kê chứng minh các thông số như TDS, EC và độ cứng chịu sự chi phối mang tính quyết định của vị trí địa lý thay vì mùa vụ quan trắc.

Hiện trạng ô nhiễm chất hữu cơ và vi sinh trên sông Gianh hiện nay ở mức độ nào?

Chất lượng nước sông Gianh cơ bản vẫn đạt tiêu chuẩn nước mặt loại A1, tuy nhiên đã xuất hiện ô nhiễm hữu cơ và vi sinh cục bộ tại các phân đoạn đô thị như trạm SG2 và SG3. Cụ thể, hàm lượng COD có thời điểm đạt 30,70 mg/L và Coliform ghi nhận mức 7.900 MPN/100mL do tiếp nhận nước thải sinh hoạt xả thẳng.

Quy trình kiểm soát chất lượng phân tích (QA/QC) trong nghiên cứu đạt độ tin cậy ra sao?

Quy trình phân tích đạt độ tin cậy rất cao với hệ số tương quan tuyến tính R trên 0,9996, giới hạn phát hiện LOD cho các chỉ tiêu amoni và nitrat chỉ từ 0,02 đến 0,03 mg/L, độ lặp lại RSD luôn thấp hơn 2/3 chuẩn Horwitz và độ thu hồi thực nghiệm đạt từ 95% đến 102%, hoàn toàn tuân thủ các quy chuẩn khắt khe của AOAC.

Kết luận

  • Luận văn đã phân tích và đánh giá toàn diện 13 thông số chất lượng nước mặt sông Gianh qua 4 đợt quan trắc thực địa tại 5 vị trí đại diện dọc lưu vực 4.680 km².
  • Xác lập quy luật phân hóa không gian rõ rệt: vùng hạ lưu chịu tác động khoáng hóa từ xâm nhập mặn với TDS trung bình 311,59 mg/L, trong khi thượng nguồn mang đặc trưng độ cứng karst cao đạt 203,80 mg/L.
  • Nhận diện kịp thời các điểm nóng ô nhiễm hữu cơ và vi sinh cục bộ tại khu vực đô thị Ba Đồn với COD đạt 30,70 mg/L và coliform lên tới 7.900 MPN/100mL.
  • Thẩm định thành công hệ thống phương pháp phân tích phòng thí nghiệm đạt chuẩn QA/QC quốc tế với hệ số tương quan R trên 0,999 và độ thu hồi từ 95% đến 102%.
  • Đề xuất hệ thống giải pháp kiểm soát chất lượng nước đồng bộ từ lắp đặt trạm quan trắc tự động, xây dựng nhà máy xử lý nước thải tập trung đến quy hoạch hành lang phân vùng xả thải với lộ trình thực thi đến năm 2030.

Công trình là đóng góp học thuật và thực tiễn sâu sắc cho công tác quản lý tài nguyên nước mặt tỉnh Quảng Bình. Độc giả và các nhà nghiên cứu quan tâm có thể khai thác chi tiết các mô hình tính toán và cơ sở dữ liệu chuyên sâu của luận văn để áp dụng vào các dự án bảo vệ môi trường lưu vực sông tiếp theo.