Tổng quan nghiên cứu

Lưu vực sông Cầu là trục thoát nước và cung cấp nguồn nước mặt thiết yếu thuộc hệ thống sông Thái Bình với tổng diện tích lưu vực khoảng 6.030 km², chiều dài dòng chính đạt 288 km chảy qua 6 tỉnh thành phía Bắc. Tổng lượng dòng chảy mặt hàng năm của lưu vực đạt khoảng 4,5 tỷ m³, giữ vai trò chiến lược đối với sự phát triển kinh tế và an sinh xã hội toàn vùng. Trong đó, đoạn sông Cầu chảy qua địa bàn tỉnh Thái Nguyên giữ vị trí đặc biệt nhạy cảm về môi trường. Với tốc độ tăng trưởng tổng sản phẩm trên địa bàn (GDP) đạt 15,2% năm 2016 cùng dòng vốn FDI lũy kế vượt 7,2 tỷ USD, hoạt động sản xuất công nghiệp luyện kim, khai khoáng và mở rộng đô thị đang tạo áp lực gay gắt lên hệ sinh thái nước mặt.

Vấn đề cốt lõi được đặt ra là sự suy giảm chất lượng nước nghiêm trọng vào mùa kiệt do tiếp nhận lượng nước thải sinh hoạt trung bình khoảng 75 lít/người/ngày cùng nguồn thải từ các khu công nghiệp chưa qua xử lý đồng bộ. Mục tiêu của nghiên cứu là đánh giá hiện trạng ô nhiễm, thiết lập công cụ mô hình hóa chuyên sâu để mô phỏng và dự báo diễn biến chất lượng nước sông Cầu đoạn từ Trạm thủy văn Gia Bảy đến Cầu Tân Phú. Nghiên cứu tập trung vào 06 thông số chỉ thị chính gồm DO, BOD5, NH4+, NO3-, PO43- và Coliform trong giai đoạn 2016 - 2017, đồng thời dự phóng các kịch bản chất lượng nước đến năm 2020 và 2030. Kết quả nghiên cứu mang ý nghĩa thực tiễn lớn khi cung cấp bộ công cụ định lượng hỗ trợ các nhà quản lý kiểm soát ô nhiễm, duy trì sai số tính toán dưới 10% và đảm bảo cấp nước an toàn cho hơn 1,2 triệu người dân đị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 trên nền tảng lý thuyết thủy động lực học dòng chảy một chiều không ổn định trong kênh hở tự nhiên và lý thuyết truyền tải khuếch tán chất ô nhiễm. Hệ phương trình Saint-Venant đầy đủ (gồm phương trình bảo toàn khối lượng liên tục và phương trình bảo toàn động lượng) được sử dụng làm cốt lõi để mô tả trường vận tốc và mực nước. Quá trình dịch chuyển và phân tán chất hòa tan trong nước sông được mô tả bằng phương trình truyền tải - khuếch tán Advection-Dispersion kết hợp mô-đun sinh thái Ecolab cấp độ 6 mở rộng.

Khung mô hình Ecolab cấp độ 6 cho phép tích hợp đồng thời 13 quá trình phản ứng sinh - hóa học phức tạp, bao gồm chu trình tiêu thụ và tái tạo oxy hòa tan (DO), quá trình phân hủy các hợp phần nhu cầu oxy sinh hóa (BOD5 dạng hòa tan, lơ lửng và tầng đáy), chu trình chuyển hóa các hợp chất Nitơ (Amoni hóa, Nitrat hóa) và Phốt pho (PO43-), cùng quy luật suy giảm vi sinh vật Coliform theo thời gian. Các khái niệm và thông số then chốt được chuẩn hóa gồm hệ số nhám lòng dẫn Manning (giá trị khởi tạo n = 0,04), hệ số phân tán dọc trục D, chỉ số hiệu quả mô hình NASH với ngưỡng tin cậy bắt buộc đạt từ 90% trở lên và quy chuẩn so sánh QCVN 08-MT:2015/BTNMT cột A2, B1.

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

Nghiên cứu áp dụng quy trình tiếp cận hệ thống, kết hợp chặt chẽ giữa quan trắc thực địa, phân tích phòng thí nghiệm và mô phỏng số trị:

  • Cỡ mẫu và mạng lưới thu thập: Thiết lập mạng lưới gồm 14 vị trí quan trắc đại diện trải dài trên đoạn sông nghiên cứu từ Gia Bảy đến Tân Phú. Tần suất quan trắc gồm 5 đợt phân bổ trong năm 2017 (vào các tháng 3, 5, 7, 9 và 11), đảm bảo bao quát cả mùa mưa và mùa khô.
  • Phương pháp lấy mẫu và phân tích: Quy trình lấy mẫu, bảo quản và thử nghiệm tuân thủ nghiêm ngặt các tiêu chuẩn quốc gia và quốc tế như TCVN 6663-1:2011, TCVN 6663-6:2008, SMEWW 5210B:2012 cho 11 chỉ tiêu lý hóa (Nhiệt độ, pH, DO, Độ đục, TDS, EC, BOD5, NH4+, NO3-, PO43-, Coliform).
  • Phương pháp mô hình hóa: Thu thập và số hóa dữ liệu địa hình từ 30 mặt cắt ngang lòng sông với khoảng cách trung bình 1,5 km/mặt cắt. Biên trên của mô hình là chuỗi dữ liệu lưu lượng thực đo tại Trạm thủy văn Gia Bảy.
  • Lý do lựa chọn phương pháp: Mô hình MIKE 11 (do DHI phát triển) được lựa chọn vì khả năng ghép nối đồng bộ 3 mô-đun HD, AD và Ecolab. Công cụ này xử lý hiệu quả bài toán thủy lực phi tuyến phức tạp của sông tự nhiên, cho phép kiểm soát tốt hiện tượng khuếch tán số và dự báo chính xác các kịch bản tải lượng trong tương lai.

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

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

Kết quả phân tích chuỗi dữ liệu quan trắc thực tế và mô phỏng trên phần mềm MIKE 11 đã chỉ ra những đặc trưng biến đổi rõ nét của chất lượng nước sông Cầu:

  • Oxy hòa tan (DO) và ô nhiễm hữu cơ (BOD5): Giá trị DO tại trạm đầu nguồn duy trì ở mức khá tốt từ 4,5 đến 7,8 mg/L, nhưng giảm mạnh khi đi qua các khu vực đông dân cư. Nồng độ BOD5 dao động từ 2,0 đến 31,0 mg/L, trong đó khoảng 9% số mẫu quan trắc vào các tháng mùa khô và đầu mùa mưa vượt ngưỡng quy chuẩn QCVN 08-MT:2015/BTNMT loại B1 (15,0 mg/L).
  • Ô nhiễm dinh dưỡng Nitơ và Phốt pho: Nồng độ Amoni (N-NH4+) diễn biến phức tạp trong khoảng 0,10 đến 1,92 mg/L. Tại một số điểm quan trắc hạ lưu khu dân cư tập trung, giá trị N-NH4+ cao gấp 3,8 lần giới hạn cho phép của chuẩn B1 (0,5 mg/L). Tỷ lệ mẫu vượt ngưỡng chuẩn A2 đối với Nitrit (N-NO2-) lên tới 59% vào đợt quan trắc tháng 3.
  • Ô nhiễm vi sinh vật: Chỉ tiêu Coliform có xu hướng tăng vọt tại các đoạn sông tiếp nhận nước thải sinh hoạt đô thị và làng nghề, nhiều thời điểm vượt ngưỡng cho phép phục vụ mục đích cấp nước sinh hoạt.
  • Độ chính xác mô hình: Quá trình hiệu chỉnh và kiểm định bộ thông số mô hình thủy lực HD cho hệ số hiệu quả NASH đạt trên 90%. Sai số tương đối giữa kết quả mô phỏng và giá trị thực đo của 06 thông số chất lượng nước trong mô-đun Ecolab đều duy trì ở mức dưới 10%, khẳng định độ tin cậy khoa học của mô hình.

Thảo luận kết quả

Sự suy giảm chất lượng nước mặt bắt nguồn từ hai nhóm nguyên nhân chính: áp lực xả thải gia tăng và đặc điểm thủy văn phân mùa rõ rệt. Về nguồn thải, việc tiếp nhận nước thải từ các đô thị dọc lưu vực (với định mức phát thải khoảng 75 lít/người/ngày) cùng hơn 20 triệu m³/năm nước thải công nghiệp từ các cơ sở luyện kim lớn đã vượt quá khả năng tự làm sạch tự nhiên của đoạn sông. Về chế độ dòng chảy, lưu lượng dòng chảy mùa kiệt (chiếm chưa đầy 20% đến 30% tổng lượng nước cả năm) khiến khả năng hòa loãng các chất ô nhiễm giảm xuống mức thấp nhất, dẫn đến hiện tượng tích tụ nồng độ Amoni và BOD5 cục bộ.

Khi so sánh với các nghiên cứu mô phỏng chất lượng nước trên thế giới như lưu vực sông Dender (Bỉ) hay sông Bertam (Malaysia), mô hình MIKE 11 ứng dụng tại sông Cầu thể hiện sự thích ứng vượt trội trong việc tái hiện động lực lan truyền chất trong điều kiện thủy văn biến đổi mạnh. Trong các báo cáo chuyên sâu, kết quả này được minh họa tối ưu thông qua hệ thống biểu đồ chuỗi thời gian so sánh giữa đường cong tính toán và điểm thực đo của 6 thông số, kết hợp sơ đồ mặt cắt dọc phân bố nồng độ chất ô nhiễm theo khoảng cách từ km số 0 (Gia Bảy) đến km 45 (Tân Phú). Cách trình bày này giúp khoanh vùng chính xác các điểm nóng tiếp nhận tải lượng ô nhiễm cao sau điểm hợp lưu sông Đu và vùng ven thành phố Thái Nguyên.

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

Dựa trên kết quả mô phỏng các kịch bản biến đổi chất lượng nước đến năm 2020 và 2030, luận văn đưa ra 4 nhóm giải pháp quản lý đồng bộ:

  • Đầu tư và nâng cấp hạ tầng xử lý nước thải đô thị: Xây dựng mới và hoàn thiện hệ thống gom, xử lý nước thải tập trung cho thành phố Thái Nguyên và thị xã Sông Công với tổng công suất thiết kế đạt tối thiểu 30.000 m³/ngày đêm; đảm bảo 100% lượng nước thải sinh hoạt đô thị đạt quy chuẩn trước khi thải ra lưu vực sông Cầu trong giai đoạn 2020 - 2025 (Chủ thể thực hiện: UBND tỉnh Thái Nguyên, Sở Xây dựng).
  • Tăng cường kiểm soát nguồn thải công nghiệp tự động: Yêu cầu 100% các nhà máy và cơ sở sản xuất có lưu lượng xả thải từ 1.000 m³/ngày đêm trở lên (đặc biệt là các tổ hợp gang thép, luyện kim, khai khoáng) phải lắp đặt trạm quan trắc tự động liên tục truyền dữ liệu trực tiếp về Sở Tài nguyên và Môi trường trước năm 2022 (Chủ thể thực hiện: Sở Tài nguyên và Môi trường, Ban Quản lý các Khu công nghiệp).
  • Ứng dụng công cụ mô hình hóa trong cấp phép xả thải: Sử dụng định kỳ bộ công cụ MIKE 11 với mạng lưới 30 mặt cắt đã hiệu chỉnh để tính toán sức chịu tải môi trường nước 2 năm/lần; tuyệt đối không cấp phép xả thải mới cho các dự án phát sinh tải lượng vượt quá khả năng tiếp nhận của từng đoạn sông (Chủ thể thực hiện: Chi cục Bảo vệ Môi trường Thái Nguyên).
  • Thiết lập cơ chế phối hợp kiểm soát ô nhiễm liên tỉnh: Xây dựng quy chế chia sẻ dữ liệu quan trắc và điều tiết dòng chảy mùa kiệt giữa các tỉnh Thái Nguyên, Bắc Kạn, Bắc Giang và Bắc Ninh, hướng tới mục tiêu giảm 25% đến 30% tổng lượng phát thải các chất ô nhiễm hữu cơ trên toàn lưu vực sông Cầu đến năm 2030 (Chủ thể thực hiện: Ủy ban Bảo vệ Môi trường Lưu vực sông Cầu).

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

Luận văn cung cấp nguồn tài liệu khoa học giá trị và có khả năng ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:

  • Cơ quan quản lý nhà nước về tài nguyên và môi trường: Cán bộ thuộc Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường có thể ứng dụng trực tiếp bộ thông số thủy lực và kết quả phân vùng tải lượng ô nhiễm 6 chỉ tiêu để thẩm định báo cáo đánh giá tác động môi trường (ĐTM), cấp giấy phép môi trường và xây dựng quy hoạch phân vùng chất lượng nước mặt.
  • Ban quản lý khu công nghiệp và doanh nghiệp sản xuất: Đội ngũ kỹ sư môi trường tại các khu, cụm công nghiệp trên địa bàn khai thác dữ liệu dự báo tải lượng BOD5 và NH4+ để thiết kế công nghệ xử lý nước thải nội bộ tối ưu, tuân thủ nghiêm quy chuẩn QCVN 08-MT:2015/BTNMT.
  • Giảng viên, nhà nghiên cứu và học viên sau đại học: Nhóm đối tượng học thuật chuyên ngành Khoa học Môi trường, Kỹ thuật Tài nguyên Nước kế thừa phương pháp luận thiết lập mô hình MIKE 11 (HD, AD, Ecolab), quy trình hiệu chỉnh hệ số Manning và chỉ số NASH trên 90% cho các nghiên cứu mô phỏng lưu vực sông khác.
  • Các đơn vị tư vấn và tổ chức phi chính phủ: Khai thác chuỗi dữ liệu 14 điểm quan trắc qua 5 đợt khảo sát thực tế để xây dựng các đề án bảo vệ nguồn nước và phát triển sinh kế bền vững cho cộng đồng dân cư ven sông.

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

  • Vì sao mô hình MIKE 11 được lựa chọn để mô phỏng chất lượng nước sông Cầu? Mô hình MIKE 11 tích hợp toàn diện 3 mô-đun HD, AD và Ecolab, có khả năng giải bài toán thủy lực dòng không ổn định và mô phỏng 13 quá trình phản ứng sinh - hóa học. Kết quả kiểm định tại sông Cầu đạt hệ số NASH trên 90% và sai số dưới 10%, đáp ứng hoàn hảo yêu cầu dự báo khoa học.

  • Những thông số chất lượng nước nào giữ vai trò chỉ thị chính trong nghiên cứu? Nghiên cứu tập trung vào 06 thông số cốt lõi gồm DO, BOD5, NH4+, NO3-, PO43- và Coliform. Trong đó, DO và BOD5 phản ánh mức độ suy giảm oxy và tải lượng hữu cơ, còn NH4+ và PO43- là các chỉ thị nhạy cảm về nguy cơ phú dưỡng và ô nhiễm nước thải sinh hoạt, công nghiệp.

  • Diễn biến chất lượng nước sông Cầu thay đổi như thế nào giữa mùa mưa và mùa khô? Chất lượng nước suy giảm mạnh vào mùa khô (tháng 11 đến tháng 4) khi dòng chảy chỉ chiếm khoảng 20% đến 30% tổng lượng cả năm. Điển hình, nồng độ NH4+ trong đợt quan trắc mùa kiệt ghi nhận giá trị lên tới 1,92 mg/L, vượt 3,8 lần quy chuẩn loại B1, trong khi mùa mưa có lưu lượng dòng chảy lớn giúp tăng cường khả năng tự làm sạch.

  • Dữ liệu đầu vào phục vụ mô hình hóa chất lượng nước gồm những thành phần nào? Dữ liệu đầu vào gồm 30 mặt cắt ngang địa hình (khoảng cách trung bình 1,5 km), chuỗi lưu lượng biên trên tại Trạm thủy văn Gia Bảy, số liệu quan trắc 14 điểm thực địa qua 5 đợt năm 2017 và dữ liệu thống kê tải lượng từ các nguồn thải sinh hoạt, làng nghề, công nghiệp luyện kim.

  • Kết quả dự báo chất lượng nước đến năm 2020 và 2030 mang lại giá trị gì cho địa phương? Kết quả mô phỏng các kịch bản tương lai giúp chính quyền định lượng trước mức độ gia tăng chất ô nhiễm khi GDP tăng trưởng trên 14%/năm. Đây là cơ sở khoa học then chốt để quy hoạch công suất nhà máy xử lý nước thải, kiểm soát điểm xả thải và đảm bảo an ninh nguồn nước bền vững.

Kết luận

  • Đã đánh giá toàn diện hiện trạng chất lượng nước mặt sông Cầu đoạn qua Thái Nguyên qua 14 điểm quan trắc với 11 thông số lý hóa và sinh học.
  • Thiết lập và kiểm định thành công mô hình MIKE 11 trên 30 mặt cắt địa hình, đạt độ tin cậy rất cao với chỉ số NASH trên 90% và sai số dưới 10%.
  • Nhận diện rõ các điểm nóng ô nhiễm hữu cơ và dinh dưỡng tập trung vào mùa khô, với nồng độ BOD5 lên tới 31,0 mg/L và NH4+ đạt 1,92 mg/L.
  • Dự báo định lượng diễn biến 06 thông số chất lượng nước then chốt theo các kịch bản phát triển kinh tế - xã hội giai đoạn 2020 - 2030.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật, kiểm soát quan trắc và quản lý liên vùng mang tính khả thi cao.

Công trình luận văn đóng góp cơ sở khoa học và công cụ mô hình hóa định lượng chuẩn xác, phục vụ đắc lực cho công tác quản trị môi trường nước lưu vực sông Cầu trong giai đoạn 2020 - 2030. Quý cơ quan quản lý, nhà nghiên cứu và doanh nghiệp quan tâm có thể khai thác, ứng dụng bộ thông số mô hình và phương pháp luận của luận văn vào thực tiễn quản lý nguồn nước tại địa phương.