Tổng quan nghiên cứu

Lưu vực sông Cầu là một trong những hệ thống thủy văn huyết mạch tại miền Bắc Việt Nam với diện tích lưu vực đạt 6.030 km² và chiều dài dòng chính lên tới 288 km, trải rộng trên địa bàn 7 tỉnh thành bao gồm Bắc Kạn, Thái Nguyên, Vĩnh Phúc, Bắc Giang, Bắc Ninh, Hải Dương và Thủ đô Hà Nội. Hiện nay, áp lực phát triển kinh tế từ quy mô dân số hơn 3,38 triệu người (mật độ trung bình 427 người/km²), cùng hoạt động của hơn 500 doanh nghiệp nhà nước, 125 mỏ khoáng sản và 230 làng nghề truyền thống đã xả thải vào lưu vực xấp xỉ 1,07 triệu m³/ngày (tương đương 390,2 triệu m³/năm). Sự gia tăng đột biến của các nguồn thải hữu cơ và hóa chất chưa qua xử lý đã đẩy chất lượng nước mặt tại nhiều phân đoạn sông rơi vào tình trạng suy thoái nghiêm trọng, vượt xa khả năng tự làm sạch tự nhiên.

Mục tiêu trọng tâm của nghiên cứu là đánh giá toàn diện hiện trạng chất lượng nước sông Cầu trên trục dòng chính từ Thác Bưởi đến Phả Lại, kiểm chứng mức độ chính xác của công cụ toán học QUAL2K trong việc mô phỏng quá trình truyền tải chất ô nhiễm, và dự báo diễn biến chất lượng nước theo các kịch bản phát triển kinh tế - xã hội giai đoạn 2015–2030. Nghiên cứu thực hiện phân tích chuyên sâu trên chuỗi dữ liệu thủy văn, khí tượng và 5 đợt quan trắc chất lượng nước thực địa tại 15 trạm dọc tuyến sông trong năm 2012. Kết quả tính toán cung cấp cơ sở định lượng về khả năng chịu tải môi trường nước, đóng góp luận cứ khoa học để hoạch định phương án giảm thiểu tới 4,5 lần tổng tải lượng chất ô nhiễm hữu cơ (BOD5), thiết lập hành lang kiểm soát xả thải bền vững cho toàn bộ lưu vự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 ứng dụng tổ hợp các lý thuyết thủy lực dòng hở và động học chất lượng nước hiện đại, lấy mô hình cân bằng khối lượng truyền tải – phân tán một chiều (1D Advection-Dispersion Mass Transport) làm nền tảng cốt lõi. Trong hệ thống dòng chảy sông ngòi, sự chuyển vận của các chất hòa tan và lơ lửng được mô tả thông qua phương trình vi phân đạo hàm riêng phản ánh quá trình bình lưu theo lưu tốc dòng chảy và quá trình khuếch tán do gradien nồng độ. Bên cạnh đó, lý thuyết phân hủy sinh học và suy giảm oxy hòa tan dựa trên mô hình Streeter-Phelps mở rộng được tích hợp để mô tả tương tác phức hợp giữa nhu cầu oxy sinh hóa (BOD) và hàm lượng oxy hòa tan (DO) trong môi trường nước mặt.

Mô hình QUAL2K (phiên bản nâng cấp từ QUAL2E) được lựa chọn nhờ khả năng mô phỏng đồng thời 15 thông số chất lượng nước đặc trưng. Điểm đột phá trong khung lý thuyết của QUAL2K là việc phân tách nhu cầu oxy carbonaceous thành hai cấu phần động học riêng biệt: dạng oxy hóa nhanh (Fast CBOD) và dạng oxy hóa chậm (Slow CBOD), kết hợp với cơ chế trao đổi nhiệt tại bề mặt nước - khí quyển và quá trình chuyển hóa dinh dưỡng nitơ qua các dạng Amoni (NH4+), Nitrit (NO2-) và Nitrat (NO3-). Khái niệm sức chịu tải môi trường nước và hạn ngạch xả thải tối đa hàng ngày (TMDL) cũng được áp dụng xuyên suốt nhằm xác định lượng chất ô nhiễm tối đa mà thủy vực có thể tiếp nhận mà vẫn duy trì giới hạn quy chuẩn QCVN 08:2008/BTNMT.

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

Nghiên cứu sử dụng nguồn dữ liệu thứ cấp và sơ cấp có độ tin cậy cao, thu thập từ mạng lưới 8 trạm khí tượng (Định Hóa, Bắc Kạn, Thái Nguyên, Vĩnh Yên, Hiệp Hòa, Bắc Ninh...) và 7 trạm thủy văn quốc gia (Thác Riềng, Chợ Mới, Gia Bảy, Chã, Phúc Lộc Phương, Đáp Cầu, Phả Lại) với chuỗi số liệu quan trắc liên tục từ năm 1956 đến 2012. Bộ dữ liệu chất lượng nước sông gồm 75 mẫu quan trắc thực địa, được thu thập qua 5 đợt khảo sát trong năm 2012 tại 15 trạm mặt cắt dọc theo dòng chính từ trạm Sơn Cầm (thượng nguồn Thái Nguyên) đến trạm Yên Dũng (hạ lưu Bắc Giang). Phương pháp lấy mẫu và phân tích trong phòng thí nghiệm tuân thủ nghiêm ngặt các Tiêu chuẩn Việt Nam (TCVN) về bảo quản mẫu và phương pháp đo quang phổ, chuẩn độ hóa học.

Phương pháp phân tích số trị được thực hiện thông qua việc chia hệ thống sông Cầu thành các đoạn sông (reaches) đồng nhất về mặt cắt hình học và đặc trưng thủy lực, sau đó chia nhỏ thành các phần tử tính toán hữu hạn (computational elements). Sơ đồ sai phân hữu hạn lùi (backward difference scheme) được kích hoạt để giải hệ phương trình bảo toàn khối lượng ở trạng thái tĩnh giả định (steady-state). Lý do lựa chọn mô hình QUAL2K bắt nguồn từ khả năng thích ứng vượt trội với các lưu vực sông có phân bố trạm đo thưa thớt, tính toán thủy lực linh hoạt thông qua phương trình Manning hoặc đường đặc tính lưu lượng, cùng khả năng phân tích độ bất định giúp quá trình hiệu chỉnh và kiểm nghiệm thông số đạt hệ số tương quan cao giữa giá trị mô phỏng và số liệu thực tế.

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 từ 75 bộ số liệu quan trắc và kết quả chạy mô hình đã chỉ ra sự suy giảm chất lượng nước rõ rệt theo chiều dọc dòng chảy từ thượng nguồn xuống hạ lưu:

  • Sự phân hóa không gian của tải lượng chất hữu cơ (BOD5): Tại vùng thượng lưu thuộc tỉnh Bắc Kạn, nồng độ BOD5 duy trì ở mức thấp từ 2 đến 4 mg/l, đáp ứng hoàn toàn giới hạn loại A1 của QCVN 08:2008/BTNMT. Tuy nhiên, khi chảy qua khu vực trung lưu và hạ lưu thuộc tỉnh Bắc Ninh và Bắc Giang, nồng độ BOD5 tăng mạnh, dao động từ 8 đến 17 mg/l, đặc biệt tại các trạm Hương Lâm, Vạn Phúc và Hòa Long, vượt ngưỡng loại A2 từ 1,3 đến 2,8 lần.
  • Ô nhiễm cục bộ hợp chất Nitơ (Amoni và Nitrit): Giá trị Amoni (NH4+) trên tuyến chính dao động từ 0,05 đến 1,31 mg/l. Nồng độ đỉnh xuất hiện tại trạm Cầu Trà Vườn (1,31 mg/l vào đợt quan trắc tháng 5/2012), vượt 2,62 lần so với quy chuẩn QCVN 08:2008 loại B1, biến điểm tiếp nhận này thành điểm nóng ô nhiễm do lượng xả thải trực tiếp từ tổ hợp công nghiệp luyện cán thép và đô thị TP. Thái Nguyên.
  • Hàm lượng chất rắn lơ lửng (TSS) vượt chuẩn diện rộng: Nồng độ TSS trên toàn tuyến biến thiên từ 15 đến 325 mg/l. Trong tổng số 75 mẫu phân tích, có 24 mẫu (chiếm 32%) vượt ngưỡng quy chuẩn A2 và 15 mẫu vượt ngưỡng B1. Ngay cả khu vực thượng nguồn Bắc Kạn cũng ghi nhận các đợt bùng phát TSS cục bộ trong mùa mưa do hoạt động tuyển khoáng và khai thác cát sỏi ven sông.
  • Cơ cấu tải lượng nguồn thải mất cân đối: Tổng lượng xả thải toàn vùng đạt 1,07 triệu m³/ngày, trong đó nước thải công nghiệp chiếm tỷ trọng lớn nhất với 51,1% (khoảng 546.770 m³/ngày), tiếp theo là nước thải sinh hoạt chiếm 44,2% (472.100 m³/ngày), nước thải làng nghề chiếm 4,3% (46.500 m³/ngày) và nước thải y tế chiếm 0,4% (4.400 m³/ngày). Hải Dương và Thái Nguyên là hai địa phương có lưu lượng xả thải cao nhất, lần lượt chiếm 24,9% (265.930 m³/ngày) và 22,6% (241.850 m³/ngày) tổng lượng thải toàn lưu vực.

Thảo luận kết quả

Diễn biến suy thoái chất lượng nước sông Cầu phản ánh sự cộng hưởng tiêu cực giữa tốc độ đô thị hóa nhanh và hạ tầng xử lý nước thải lạc hậu. Việc nồng độ BOD5 và NH4+ tăng vọt ở hạ lưu được minh chứng bởi sự gia nhập của sông Ngũ Huyện Khê – trục tiêu tiếp nhận hơn 30.000 m³/ngày nước thải đậm đặc từ các làng nghề giấy Phong Khê, Phú Lâm và sắt thép Đa Hội thuộc tỉnh Bắc Ninh. Dữ liệu thực nghiệm khi đối chiếu với các đồ thị phân bố nồng độ dọc trục sông (longitudinal profile) mô phỏng bằng QUAL2K cho thấy đoạn từ sau điểm nhập lưu sông Ngũ Huyện Khê đến Đáp Cầu có độ suy giảm oxy hòa tan rõ nét nhất trước khi bước vào vùng tự làm sạch ở hạ lưu Phả Lại.

Khi phân tích qua bảng so sánh kịch bản phát triển năm 2020 và 2030, nếu không thực hiện quy hoạch thu gom xử lý (Kịch bản 1), tải lượng BOD5 toàn lưu vực sẽ bùng nổ từ 66,2 tấn/ngày lên 228,88 tấn/ngày vào năm 2020 và chạm mức 358,15 tấn/ngày vào năm 2030 (228,03 tấn/ngày từ sinh hoạt và 130,12 tấn/ngày từ công nghiệp). Ngược lại, kịch bản có giải pháp xử lý nước thải đạt quy chuẩn (Kịch bản 2) sẽ kéo giảm tải lượng BOD5 năm 2030 xuống chỉ còn 79,12 tấn/ngày (giảm 4,5 lần so với kịch bản không xử lý). So sánh với các nghiên cứu áp dụng QUAL2K trên lưu vực sông Sài Gòn – Đồng Nai và sông Cầu Bây, kết quả này hoàn toàn tương thích về mặt quy luật thủy lực và khẳng định tính tất yếu của việc kiểm soát ô nhiễm tại nguồn phát sinh.

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

Nhằm phục hồi chất lượng nước và duy trì khả năng chịu tải của lưu vực sông Cầu, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật và quản trị trọng tâm:

  1. Xây dựng và nâng cấp hệ thống thu gom, xử lý nước thải đô thị tập trung: Ủy ban nhân dân các tỉnh Thái Nguyên, Bắc Ninh, Bắc Giang và Vĩnh Phúc phối hợp với Sở Xây dựng triển khai các nhà máy xử lý nước thải sinh hoạt tại 5 thành phố và 5 thị xã thuộc lưu vực. Target metric: Thu gom và xử lý đạt quy chuẩn 100% nước thải sinh hoạt đô thị, giảm tải lượng BOD5 sinh hoạt toàn lưu vực xuống dưới mức 32,65 tấn/ngày vào năm 2030, thực hiện phân kỳ đầu tư giai đoạn 2024–2030.
  2. Kiểm soát nghiêm ngặt nước thải công nghiệp và khai thác khoáng sản: Sở Tài nguyên và Môi trường 7 tỉnh thành chủ trì tiến hành rà soát, bắt buộc 100% các khu công nghiệp tập trung (đặc biệt là KCN Sông Công, Lưu Xá, Quế Võ, Khai Quang) và các cơ sở luyện kim, sản xuất giấy phải lắp đặt hệ thống quan trắc nước thải tự động liên tục trước năm 2025. Target metric: Xử lý 100% nước thải công nghiệp đạt QCVN 40:2011/BTNMT Cột A, triệt tiêu hoàn toàn lượng phenol, xianua và giảm tải lượng BOD5 công nghiệp xả thải xuống dưới 46,47 tấn/ngày trước năm 2028.
  3. Quy hoạch cải tạo môi trường nước thải làng nghề và lưu vực sông nhánh: Ban Quản lý Lưu vực sông Cầu phối hợp UBND tỉnh Bắc Ninh khẩn cấp đầu tư cụm xử lý nước thải cục bộ cho 62 làng nghề dọc sông Ngũ Huyện Khê, kiểm soát 46.500 m³/ngày nước thải nghề gốm, giấy, cơ khí. Target metric: Cắt giảm 70% lượng hữu cơ và chất rắn lơ lửng từ các nhánh sông đổ vào dòng chính sông Cầu giai đoạn 2024–2027.
  4. Thiết lập công cụ quản lý chất lượng nước dựa trên hạn ngạch chịu tải: Cục Quản lý Tài nguyên Nước và Tổng cục Môi trường ứng dụng bộ thông số mô hình QUAL2K đã hiệu chỉnh làm công cụ phần mềm chuẩn hóa để cấp phép xả thải vào nguồn nước. Target metric: Phân bổ trần xả thải (TMDL) cho từng đoạn sông, không cấp phép mới cho các cơ sở xả thải vượt ngưỡng chịu tải tại các đoạn sông qua Thái Nguyên và Bắc Ninh từ năm 2025.

Đố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 (Bộ TN&MT, Sở TN&MT 7 tỉnh lưu vực): Ứng dụng kết quả tính toán sức chịu tải và mô hình QUAL2K để phê duyệt đánh giá tác động môi trường (ĐTM), phân vùng xả thải, cấp phép môi trường và xây dựng quy hoạch bảo vệ môi trường lưu vực sông liên tỉnh.
  • Các viện nghiên cứu, chuyên gia và tư vấn quy hoạch hạ tầng cấp thoát nước: Sử dụng các thông số động học thủy lực, hệ số phân hủy BOD, hệ số nhám Manning và dữ liệu cân bằng lưu lượng làm tài liệu cơ sở để thiết kế quy hoạch hệ thống thoát nước lưu vực đến năm 2030–2050.
  • 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/Thủy văn: Tham khảo quy trình thiết lập, hiệu chỉnh, kiểm nghiệm mô hình toán QUAL2K trên hệ thống sông tự nhiên có cấu trúc thủy văn phức tạp và ứng dụng mô hình số trị trong quản trị tài nguyên nước.
  • Ban quản lý các Khu công nghiệp, cụm công nghiệp và doanh nghiệp sản xuất: Nắm rõ ngưỡng giới hạn xả thải và cơ chế diễn biến chất lượng nước sông tiếp nhận để chủ động đầu tư, nâng cấp công nghệ tiền xử lý nước thải đạt tiêu chuẩn đầu ra bền vững.

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

  • Mô hình QUAL2K có ưu điểm vượt trội gì trong việc mô phỏng chất lượng nước sông Cầu? QUAL2K vượt trội nhờ khả năng chia sông thành các đoạn có chiều dài linh hoạt, mô phỏng chi tiết 15 thông số chất lượng nước, đặc biệt là tách biệt quá trình oxy hóa carbon hữu cơ thành hai dạng CBOD phân hủy nhanh và chậm, giải quyết bài toán truyền tải chất ô nhiễm chính xác ngay cả trong điều kiện dữ liệu đo đạc còn hạn chế.
  • Tại sao chất lượng nước sông Cầu đoạn qua Bắc Ninh và Bắc Giang bị ô nhiễm hữu cơ nặng nề nhất? Khu vực này chịu tải lượng ô nhiễm tích lũy từ toàn bộ thượng nguồn đổ về, đồng thời tiếp nhận trực tiếp nguồn nước thải đậm đặc chưa qua xử lý từ sông Ngũ Huyện Khê (hơn 30.000 m³/ngày) chứa nước thải từ các làng nghề giấy Phong Khê, sắt cán Đa Hội và các cụm công nghiệp tập trung quy mô lớn.
  • Số liệu quan trắc năm 2012 cho thấy những thông số nào trên sông Cầu có mức độ vượt quy chuẩn đáng báo động nhất? Thông số Amoni (NH4+) tại Cầu Trà Vườn đạt mức cực đại 1,31 mg/l (vượt 2,62 lần QCVN 08:2008 loại B1); thông số TSS có 24/75 mẫu vượt chuẩn A2 và 15 mẫu vượt chuẩn B1 với giá trị cao nhất lên đến 325 mg/l; thông số BOD5 tại hạ lưu dao động từ 8 đến 17 mg/l vượt xa giới hạn A2.
  • Khả năng chịu tải chất ô nhiễm của sông Cầu được tính toán dựa trên các yếu tố cốt lõi nào? Khả năng chịu tải được xác định dựa trên lưu lượng dòng chảy mùa kiệt tại 7 trạm thủy văn, nồng độ nền của các chất ô nhiễm từ thượng lưu, hệ số tiêu thụ oxy kD, thời gian lưu chuyển của dòng nước và giới hạn nồng độ cho phép quy định trong QCVN 08:2008/BTNMT đối với từng mục đích sử dụng nước.
  • Hiệu quả cắt giảm ô nhiễm thực tế ra sao nếu áp dụng kịch bản xử lý nước thải đạt chuẩn đến năm 2030? Theo mô phỏng, nếu toàn bộ nguồn thải sinh hoạt và công nghiệp được xử lý đạt quy chuẩn trước khi xả ra sông, tổng tải lượng BOD5 đổ vào lưu vực năm 2030 sẽ giảm từ 358,15 tấn/ngày xuống còn 79,12 tấn/ngày, tương đương mức cắt giảm tới 4,5 lần, giúp khôi phục chất lượng nước sông về mức an toàn.

Kết luận

  • Nghiên cứu đã phác thảo bức tranh toàn cảnh về hiện trạng ô nhiễm nước mặt trên 288 km dòng chính sông Cầu thông qua chuỗi số liệu thực nghiệm 75 mẫu quan trắc tại 15 trạm đo đại diện.
  • Ứng dụng thành công mô hình QUAL2K với độ tin cậy cao, giải quyết đồng thời bài toán cân bằng thủy lực và động học phân hủy 15 thành phần ô nhiễm trong môi trường nước mặt nhiệt đới gió mùa.
  • Định lượng chính xác ngưỡng chịu tải BOD5 cho từng đoạn sông, chỉ ra tình trạng quá tải nghiêm trọng tại khu vực hạ lưu tiếp giáp các đô thị, khu công nghiệp Thái Nguyên, Bắc Ninh và Bắc Giang.
  • Dự báo thành công các kịch bản diễn biến môi trường đến năm 2020 và 2030, chứng minh tính khả thi của giải pháp kiểm soát xả thải tập trung giúp giảm 4,5 lần tải lượng BOD5 phát thải.
  • Đề xuất hệ thống 4 giải pháp đồng bộ kết hợp giữa quy hoạch hạ tầng kỹ thuật xử lý nước thải và công cụ chính sách quản lý hạn ngạch chịu tải môi trường liên tỉnh.

Để hiện thực hóa các mục tiêu phục hồi sinh thái sông Cầu, các cơ quan quản lý nhà nước cần sớm đưa bộ công cụ QUAL2K vào quy trình vận hành giám sát thường nhật, đồng thời thúc đẩy hợp tác liên vùng giữa 7 tỉnh thành nhằm quản lý toàn diện tài nguyên nước lưu vực ngay trong giai đoạn 2024–2030.