Tổng quan nghiên cứu

Hệ thống sông Thái Bình là mạng lưới thủy văn lớn thứ hai tại miền Bắc Việt Nam, giữ vai trò huyết mạch trong việc cung cấp nước phục vụ sản xuất nông nghiệp, công nghiệp và sinh hoạt của hàng triệu người dân Đồng bằng sông Hồng. Đoạn sông Thái Bình chảy qua địa phận tỉnh Hải Dương có chiều dài 73 km với tổng lượng nước trung bình dao động từ 30 đến 40 tỷ m³/năm, tiếp nhận khoảng 22,9 triệu tấn phù sa mỗi năm từ sông Hồng qua sông Luộc và sông Đuống. Song song với đó, hệ thống thủy nông Bắc Hưng Hải với hơn 200 km kênh chính đảm nhiệm vai trò tưới tiêu cho diện tích canh tác khoảng 200.230 ha. Tuy nhiên, trước áp lực từ quá trình đô thị hóa nhanh chóng với quy mô dân số hơn 1,512 triệu người (mật độ 1.055 người/km²) cùng sự mở rộng của 18 khu công nghiệp tập trung và 38 cụm công nghiệp, chất lượng môi trường nước mặt tại Hải Dương đang đối mặt với nguy cơ suy thoái nghiêm trọng do ô nhiễm các hợp chất nitơ.

Vấn đề cốt lõi đặt ra là sự gia tăng mất kiểm soát của các dạng nitơ vô cơ như amoni, nitrit và nitrat trong nguồn nước, đe dọa trực tiếp đến hệ sinh thái thủy sinh và nguồn nước cấp sinh hoạt. Nghiên cứu của tác giả Lê Phú Đồng tập trung lượng hóa toàn diện các dòng vật chất nitơ, xác định chính xác nguồn phát thải then chốt từ các hoạt động sinh hoạt, chăn nuôi, trồng trọt và công nghiệp trên địa bàn tỉnh Hải Dương trong giai đoạn 2011 - 2013, đồng thời xây dựng kịch bản dự báo tải lượng ô nhiễm đến năm 2020. Bằng việc tích hợp đồng thời mô hình Phân tích Dòng chảy Vật chất (MFA) và Hệ thống Thông tin Địa lý (GIS), công trình mang lại giá trị thực tiễn cao khi giải quyết triệt để bài toán thiếu hụt dữ liệu quan trắc phân tán. Kết quả nghiên cứu đạt độ hội tụ thống kê cao với 1.000 vòng lặp mô phỏng xác suất, cung cấp luận cứ khoa học chuẩn xác để các cơ quan quản lý ban hành chính sách kiểm soát ô nhiễm nước mặt theo 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 được xây dựng dựa trên sự tích hợp giữa hai nền tảng lý thuyết và công cụ phân tích hiện đại trong quản lý môi trường lưu vực:

  • Lý thuyết Phân tích Dòng chảy Vật chất (Material Flow Analysis - MFA): Dựa trên nguyên lý bảo toàn khối lượng và phương pháp tiếp cận hệ thống do Baccini, Brunner và Steurer phát triển từ năm 1992, kết hợp các cải tiến ứng dụng của Montangero và Belevi. Mô hình MFA xem xét toàn bộ vòng tuần hoàn của nguyên tố nitơ thông qua phương trình cân bằng biến thiên tích lũy vật chất: hiệu số giữa tổng dòng vật chất đầu vào và tổng dòng vật chất đầu ra tại từng quá trình chuyển hóa xác định trong một không gian và thời gian cụ thể.
  • Lý thuyết Mô hình hóa Thủy lực và Phân phối Không gian GIS: Ứng dụng mô hình phân phối thông số dòng chảy của Lê Văn Chính và cộng sự trên cơ sở liên kết dữ liệu địa lý không gian với nồng độ quan trắc thực địa. Mô hình thiết lập hàm quan hệ phi tuyến giữa tỷ lưu tải lượng nitơ tiếp nhận với lưu lượng dòng chảy và tỷ lệ diện tích các loại hình sử dụng đất.
  • Các khái niệm chuyên ngành cốt lõi: Quá trình amoni hóa và nitrat hóa hợp chất hữu cơ; Tải lượng phát thải nguồn phân tán từ 68.542 ha đất canh tác nông nghiệp; Tải lượng phát thải nguồn điểm từ hệ thống trạm bơm thoát nước đô thị và khu công nghiệp; Độ không đảm bảo đo và phân bố xác suất tham số dòng vật chất.

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

  • Nguồn dữ liệu và mạng lưới thu thập: Dữ liệu nghiên cứu bao gồm số liệu quan trắc thực nghiệm tại 31 vị trí trọng điểm trên các hệ thống sông thuộc địa bàn Hải Dương (19 điểm trên hệ thống sông Thái Bình và 12 điểm trên hệ thống kênh Bắc Hưng Hải) với tần suất đo đạc 4 đợt/năm liên tục trong giai đoạn 2011 - 2013. Bên cạnh đó, các số liệu thứ cấp về hiện trạng sử dụng đất, dân số, công nghiệp và chăn nuôi được thu thập từ Niên giám Thống kê tỉnh Hải Dương năm 2013 và các báo cáo hiện trạng môi trường của Sở Tài nguyên và Môi trường.
  • Cỡ mẫu và phương pháp chọn mẫu: 31 điểm quan trắc được lựa chọn theo phương pháp phân tầng không gian dọc theo trục thủy văn chính, bao gồm các điểm hợp lưu nhạy cảm (ngã ba sông Lục Nam, sông Thương, sông Đuống), các mặt cắt tiếp nhận nước thải sinh hoạt đô thị (sông Sặt tại Cầu Cất, Cầu Ghẽ), khu vực xả thải công nghiệp (KCN Lai Vu, KCN Đại An) và vùng hạ lưu tưới tiêu nông nghiệp. Phương pháp lấy mẫu và bảo quản tuân thủ quy chuẩn quốc gia về chất lượng nước mặt.
  • Phương pháp phân tích hóa học và lý do lựa chọn: Các chỉ tiêu nitơ được phân tích bằng phương pháp trắc phổ hấp thụ phân tử UV-VIS trên máy Shimadzu (vùng bước sóng 190 - 800 nm). Cụ thể:
    • Chỉ tiêu nitrit (NO2--N) xác định theo TCVN 6178:1996 bằng phản ứng tạo màu điazô hóa với thuốc thử Griess A (axit sulfanilic) và Griess B (naphthylamin).
    • Chỉ tiêu nitrat (NO3--N) xác định theo TCVN 6180:1996 bằng phản ứng tạo màu với axit sunfosalixylic trong môi trường kiềm.
    • Chỉ tiêu amoni (NH4+-N) xác định theo TCVN 6179-1:1996 bằng phương pháp tạo phức màu xanh indophenol với natri salixylat và natri diclorosoxyanurat. Lý do lựa chọn nhóm phương pháp này là nhờ độ nhạy quang học cao, giới hạn phát hiện thấp, độ lặp lại ổn định và đã được kiểm định đạt chuẩn ISO/IEC 17025:2005 tại phòng thí nghiệm chuyên môn.
  • Kỹ thuật xử lý số liệu: Ứng dụng kỹ thuật Macro trên Microsoft Excel kết hợp thuật toán mô phỏng Monte Carlo với 1.000 vòng lặp ngẫu nhiên để đánh giá độ nhạy tham số, kiểm soát sai số và xác định các luồng phân phối nitơ có độ tin cậy vượt trên 68%.

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

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

  • Ô nhiễm Amoni (NH4+-N) ở mức báo động trên diện rộng: Nồng độ NH4+-N tại hầu hết các trạm quan trắc trên sông Thái Bình và sông Bắc Hưng Hải đều vượt ngưỡng quy định cho mục đích cấp nước sinh hoạt loại A1 và A2 theo QCVN 08:2008/BTNMT. Đặc biệt trên hệ thống kênh Bắc Hưng Hải, nồng độ amoni tại trạm Sđ25 (Cầu Sặt, huyện Bình Giang), Sđ26 (Cầu Cất) và Sđ27 (hạ lưu điểm xả TP. Hải Dương) thường xuyên vượt ngưỡng B1 gấp 2 đến 4 lần. Vào đợt quan trắc mùa khô năm 2013, các điểm hợp lưu phía thượng nguồn sông Thương và sông Lục Nam (S1, S2, S3) ghi nhận nồng độ amoni tăng đột biến, vượt qua cả giới hạn cho phép đối với giao thông thủy (mức B2).
  • Biến động phức tạp của Nitrit (NO2--N) và Nitrat (NO3--N): Hàm lượng nitrit tại nhiều vị trí lấy mẫu dao động từ 0,02 mg/l đến trên 0,05 mg/l, vượt ngưỡng tiêu chuẩn nước mặt loại A2 và chạm ngưỡng B1. Sự xuất hiện đồng thời của nồng độ amoni và nitrit cao cho thấy quá trình oxy hóa tự nhiên các hợp chất nitơ trong nước sông diễn ra chưa triệt để do điều kiện thiếu hụt oxy hòa tan cục bộ.
  • Lượng hóa nguồn thải qua mô hình MFA: Kết quả cân bằng vật chất chỉ ra rằng trên 60% tổng lượng nitơ tiêu thụ từ nguồn thực phẩm của các hộ gia đình bị đào thải trực tiếp ra môi trường nước mặt qua hệ thống vệ sinh và cống rãnh. Đáng chú ý, các bể tự hoại tại khu vực dân cư chỉ giữ lại được khoảng 5% đến 14% lượng nitơ dưới dạng bùn phân tích lũy, phần lớn lượng nitơ hòa tan còn lại bị xả thẳng ra kênh mương nội đồng.
  • Tác động từ hoạt động nông nghiệp và chăn nuôi: Hoạt động canh tác trên 68.542 ha đất trồng lúa kết hợp với nước thải từ hàng triệu gia súc, gia cầm chưa qua xử lý là nguồn phát sinh nitơ phân tán lớn nhất, chiếm trên 55% tổng tải lượng nitơ thất thoát vào các lưu vực sông nội tỉnh vào mùa mưa lũ.

Thảo luận kết quả

Nguyên nhân chính dẫn đến tình trạng ô nhiễm amoni nghiêm trọng tại hạ lưu sông Sặt là do tiếp nhận toàn bộ nước thải sinh hoạt chưa qua xử lý của khu vực nội thành Hải Dương thông qua trạm bơm Âu Thuyền với công suất xả thải lên tới 32.000 m³/giờ. Tại các huyện ngoại thành như Bình Giang, Cẩm Giàng và Gia Lộc, sự kết hợp giữa nước thải làng nghề sản xuất thực phẩm (chế biến nông sản, giết mổ gia súc) và dòng chảy tràn bề mặt mang theo phân bón hóa học dư thừa từ đồng ruộng đã làm trầm trọng thêm mức độ phú dưỡng hóa nguồn nước.

Dữ liệu quan trắc 12 đợt đo đạc từ năm 2011 đến năm 2013 có thể được minh họa trực quan thông qua các biểu đồ chuỗi thời gian phân bố nồng độ chất ô nhiễm và bảng tổng hợp tải lượng phát thải từng nguồn. Khi biểu diễn dòng chuyển hóa nitơ bằng sơ đồ dòng chảy vật chất Sankey, sự tương quan giữa dòng dinh dưỡng đầu vào từ phân bón, thức ăn chăn nuôi và dòng thất thoát ra nguồn nước mặt được thể hiện rõ ràng với độ chính xác cao. So sánh với các nghiên cứu của Montangero tại khu vực Việt Trì (Phú Thọ) và Hasan Belevi tại Kumasi (Ghana), cấu trúc thất thoát nitơ tại Hải Dương phản ánh quy luật điển hình của các vùng đô thị - nông nghiệp chuyển dịch tại các quốc gia đang phát triển: tỷ lệ xử lý nước thải tại chỗ thấp và cơ sở hạ tầng thoát nước chưa hoàn thiện là nguyên nhân cốt lõi gây ô nhiễm nước mặt.

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

  • Quy hoạch và xây dựng hệ thống xử lý nước thải đô thị tập trung: Xây dựng mới và nâng cấp đồng bộ mạng lưới thu gom nước thải sinh hoạt tại thành phố Hải Dương và các thị trấn trung tâm huyện, đưa vào vận hành các nhà máy xử lý nước thải sinh hoạt áp dụng công nghệ bùn hoạt tính hiếu khí kết hợp khử nitrat hóa sinh học. Mục tiêu giảm thiểu từ 40% đến 50% tải lượng nitơ hòa tan phát thải trực tiếp vào sông Sặt và sông Thái Bình trong giai đoạn 2015 - 2020. Chủ thể thực hiện: Ủy ban nhân dân tỉnh Hải Dương phối hợp cùng Sở Xây dựng và Sở Tài nguyên và Môi trường.
  • Kiểm soát và chuyển giao công nghệ phân bón trong sản xuất nông nghiệp: Triển khai quy trình canh tác lúa tiên tiến (nông nghiệp thông minh thích ứng biến đổi khí hậu), hướng dẫn nông dân cân đối tỷ lệ bón phân đạm urê trên 68.542 ha đất nông nghiệp, áp dụng các chế phẩm sinh học giữ đạm trong đất nhằm hạn chế thất thoát do dòng chảy bề mặt. Đặt mục tiêu cắt giảm từ 20% đến 25% lượng nitơ rửa trôi vào hệ thống kênh mương thủy nông đến năm 2018. Chủ thể thực hiện: Sở Nông nghiệp và Phát triển Nông thôn cùng hệ thống Hợp tác xã Nông nghiệp địa phương.
  • Giám sát nghiêm ngặt nước thải công nghiệp và cụm công nghiệp: Yêu cầu 100% trong số 18 khu công nghiệp và 38 cụm công nghiệp trên địa bàn tỉnh phải hoàn thiện trạm xử lý nước thải tập trung đạt tiêu chuẩn xả thải QCVN 40:2011/BTNMT trước khi xả ra lưu vực sông. Lắp đặt trạm quan trắc tự động liên tục đối với các thông số lưu lượng, pH, COD và amoni tại các cửa xả lớn. Hoàn thành kiểm tra và xử phạt vi phạm trước năm 2017. Chủ thể thực hiện: Ban Quản lý các Khu công nghiệp tỉnh Hải Dương và Cảnh sát Môi trường.
  • Ứng dụng hệ thống số hóa GIS trong quản lý cấp phép xả thải: Tích hợp mô hình tính toán tải lượng MFA-GIS vào phần mềm quản lý môi trường lưu vực cấp tỉnh, số hóa dữ liệu xả thải của 31 trạm quan trắc nhằm phục vụ công tác cấp phép xả nước thải vào nguồn nước dựa trên sức chịu tải thực tế của từng đoạn sông. Lộ trình triển khai định kỳ cập nhật dữ liệu hàng quý từ năm 2015 đến năm 2020. Chủ thể thực hiện: Trung tâm Quan trắc và Phân tích Môi trường tỉnh Hải Dương.

Đố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: Sở Tài nguyên và Môi trường, Chi cục Bảo vệ Môi trường và Phòng Tài nguyên Môi trường các huyện trên địa bàn tỉnh Hải Dương có thể ứng dụng kết quả mô hình MFA-GIS để thiết lập bản đồ nhạy cảm môi trường, xác định hạn ngạch xả thải nitơ cho từng tiểu lưu vực và xây dựng đề án bảo vệ chất lượng nước sông Thái Bình.
  • Kỹ sư và nhà quy hoạch hạ tầng đô thị - nông thôn: Các đơn vị tư vấn thiết kế hạ tầng cấp thoát nước có thể khai thác số liệu dòng thải của luận văn để tính toán dung tích hợp lý cho các trạm bơm, lựa chọn công nghệ xử lý nước thải vệ sinh phân tán phù hợp với đặc thù mật độ dân số 1.055 người/km² tại các thị tứ nông thôn.
  • Nhà nghiên cứu, giảng viên và sinh viên ngành Môi trường/Hóa học: Cung cấp tài liệu tham khảo giá trị về quy trình phân tích trắc phổ xác định các hợp chất của nitơ chuẩn ISO/IEC 17025:2005, đồng thời đóng vai trò là case-study mẫu mực trong việc áp dụng phương pháp mô phỏng Monte Carlo để lượng hóa sai số dữ liệu môi trường trong điều kiện cơ sở dữ liệu hạn chế.
  • Ban quản lý hạ tầng các Khu công nghiệp và doanh nghiệp xử lý nước thải: Sử dụng mô hình chuyển hóa vật chất để tính toán chính xác tải lượng amoni đầu vào, từ đó tối ưu hóa tỷ lệ tuần hoàn bùn và bổ sung nguồn cơ chất carbon trong các bể xử lý sinh học Anoxic-Oxic, giúp tiết kiệm từ 15% đến 20% chi phí vận hành hóa chất.

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

  • Sự kết hợp giữa mô hình MFA và GIS mang lại ưu thế vượt trội gì so với phương pháp quan trắc truyền thống? Phương pháp quan trắc truyền thống chỉ phản ánh nồng độ tức thời tại 31 điểm đo riêng lẻ mà không chỉ ra được nguồn gốc xuất phát của chất ô nhiễm. Việc tích hợp mô hình MFA giúp định lượng chính xác dòng vật chất nitơ từ khâu tiêu dùng đến các điểm xả, trong khi công cụ GIS trực quan hóa sự phân bố tải lượng ô nhiễm theo không gian lưu vực sông.
  • Tại sao nồng độ Amoni tại hệ thống kênh Bắc Hưng Hải thường xuyên cao hơn dòng chính sông Thái Bình? Hệ thống sông Bắc Hưng Hải chủ yếu là các kênh tiêu nông nghiệp có vận tốc dòng chảy chậm, độ dốc thủy lực nhỏ và tiếp nhận lượng lớn nước thải sinh hoạt đô thị chưa qua xử lý (như trạm bơm Âu Thuyền công suất 32.000 m³/giờ). Khả năng tự làm sạch và hòa loãng của hệ thống này thấp hơn rất nhiều so với lưu lượng 30 - 40 tỷ m³/năm của dòng chính sông Thái Bình.
  • Các phương pháp phân tích quang phổ UV-VIS trong luận văn có đảm bảo độ tin cậy khoa học không? Toàn bộ quy trình phân tích mẫu nước được thực hiện trên máy trắc phổ UV-VIS Shimadzu theo đúng các tiêu chuẩn Việt Nam TCVN 6178:1996, TCVN 6180:1996 và TCVN 6179-1:1996. Phòng thí nghiệm thực hiện đã được Văn phòng Công nhận Chất lượng chứng nhận đạt tiêu chuẩn quốc tế ISO/IEC 17025:2005, bảo đảm độ nhạy và tính lặp lại của số liệu.
  • Thuật toán mô phỏng Monte Carlo đóng vai trò gì trong việc xử lý số liệu của luận văn? Do dữ liệu nguồn thải tại các nước đang phát triển thường phân tán và thiếu đồng bộ, thuật toán Monte Carlo với 1.000 vòng lặp ngẫu nhiên giúp xác định phân bố xác suất của các tham số đầu vào. Kỹ thuật này giúp loại bỏ độ không đảm bảo đo và đảm bảo trên 68% giá trị dự báo tải lượng nitơ nằm trong khoảng tin cậy chấp nhận được.
  • Mô hình tính toán tải lượng nitơ trong nghiên cứu có thể nhân rộng cho các địa phương khác không? Phương pháp luận tích hợp MFA và GIS phân phối thông số hoàn toàn có thể chuyển giao áp dụng cho các lưu vực sông khác tại Đồng bằng Bắc Bộ như lưu vực sông Cầu, sông Đáy hoặc sông Nhuệ. Các địa phương chỉ cần hiệu chỉnh các tham số đầu vào về diện tích đất canh tác, mật độ dân số và hệ số phát thải đặc thù.

Kết luận

  • Luận văn đã phân tích toàn diện hiện trạng ô nhiễm các hợp chất nitơ tại 31 vị trí quan trắc trên sông Thái Bình và kênh Bắc Hưng Hải trong giai đoạn 2011 - 2013, xác định nồng độ amoni và nitrit tại nhiều khu vực đã vượt quy chuẩn QCVN 08:2008/BTNMT.
  • Ứng dụng thành công mô hình Phân tích Dòng chảy Vật chất (MFA) để lượng hóa các nguồn thải nitơ trên địa bàn tỉnh Hải Dương, chỉ ra rằng nước thải sinh hoạt đô thị và hoạt động nông nghiệp là hai nguồn gây ô nhiễm nitơ lớn nhất.
  • Tích hợp hiệu quả công nghệ GIS với mô hình cân bằng thủy lực nhằm mô phỏng quy luật phân bố tải lượng chất ô nhiễm theo không gian và thời gian với 1.000 lần lặp kiểm định Monte Carlo.
  • Đóng góp luận cứ khoa học thực tiễn giúp địa phương xây dựng kịch bản kiểm soát môi trường nước mặt đến năm 2020 và tối ưu hóa quy hoạch hệ thống trạm xử lý nước thải tập trung.
  • Đề xuất lộ trình kiểm soát ô nhiễm giai đoạn 2015 - 2020 với các mục tiêu định lượng cụ thể nhằm giảm thiểu tải lượng nitơ xả thải vào lưu vực sông.

Các nhà quản lý môi trường, nhà khoa học và đơn vị phát triển hạ tầng cần đẩy mạnh khai thác cơ sở dữ liệu và phương pháp luận của nghiên cứu này để xây dựng các giải pháp kiểm soát ô nhiễm nguồn nước mặt bền vững cho toàn bộ lưu vực sông Thái Bình.