Tổng quan nghiên cứu

Tài nguyên nước mặt giữ vai trò then chốt đối với sự phát triển kinh tế và an sinh xã hội của vùng đồng bằng sông Hồng, đặc biệt là tại tỉnh Hải Dương với tổng diện tích tự nhiên 1.656 km2 và quy mô dân số đạt trên 1,51 triệu người (mật độ bình quân 1.055 người/km2). Hải Dương sở hữu mạng lưới thủy văn dày đặc gồm 14 sông lớn có tổng chiều dài khoảng 500 km, trong đó dòng chính sông Thái Bình chảy qua địa bàn tỉnh dài 73 km với tổng lưu lượng hàng năm dao động từ 30 đến 40 tỷ m3 nước và tiếp nhận khoảng 22,9 triệu tấn phù sa từ sông Hồng. Cùng 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 phục vụ tưới tiêu cho hơn 200.000 ha đất canh tác nông nghiệp. Tuy nhiên, áp lực từ quá trình công nghiệp hóa với 18 khu công nghiệp tập trung và 38 cụm công nghiệp, kết hợp với 104.238 ha đất sản xuất nông nghiệp và nước thải sinh hoạt đô thị chưa qua xử lý triệt để đã làm suy giảm nghiêm trọng chất lượng nước mặt.

Vấn đề nghiên cứu trọng tâm là sự gia tăng nồng độ các hợp chất ô nhiễm Nitơ (gồm Amoni, Nitrit, Nitrat) tại nhiều điểm quan trắc, trong khi việc bóc tách và định lượng cụ thể tỷ trọng đóng góp từ từng nguồn thải vào lưu vực sông còn gặp nhiều thách thức. Mục tiêu cụ thể của luận văn là:

  • Đánh giá toàn diện hiện trạng ô nhiễm các dạng hợp chất Nitơ trên hệ thống sông Thái Bình và sông Bắc Hưng Hải giai đoạn 2011 - 2013.
  • Thiết lập và ứng dụng mô hình Phân tích Dòng chảy Vật chất (Material Flow Analysis - MFA) nhằm định lượng và truy xuất nguồn gốc dòng vật chất Nitơ từ các khâu phát sinh đến môi trường tiếp nhận.
  • Tích hợp Hệ thống Thông tin Địa lý (GIS) xây dựng mô hình phân phối thông số để tính toán tải lượng Nitơ và dự báo các kịch bản ô nhiễm đến năm 2020.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở định lượng có độ tin cậy đạt từ 68% đến 95% qua phân tích độ nhạy, tạo tiền đề kỹ thuật giúp các nhà quản lý địa phương xây dựng chiến lược cắt giảm từ 20% đến 30% tải lượng Nitơ phát thải vào nguồn nước mặt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn vận dụng đồng bộ hai nền tảng lý thuyết và công cụ không gian 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 (MFA): Được phát triển từ các nguyên lý nhiệt động học và phương pháp tiếp cận hệ thống, MFA vận hành trên nền tảng định luật bảo toàn khối lượng (Mass Balance). Mọi quá trình biến đổi chất được biểu diễn qua phương trình cân bằng: độ thay đổi lượng tích lũy trong hệ thống bằng tổng dòng vật chất đầu vào trừ tổng dòng vật chất đầu ra theo chu kỳ thời gian xác định. MFA đóng vai trò kết nối giữa nguồn thải, quá trình trung gian và môi trường tiếp nhận cuối cùng.
  • Lý thuyết mô hình hóa không gian tích hợp Hệ thống Thông tin Địa lý (GIS): GIS cho phép số hóa, phân tích đa lớp và trực quan hóa các dữ liệu không gian liên quan đến mạng lưới sông ngòi, hiện trạng sử dụng đất và tọa độ phân bố nguồn thải. Mô hình phân phối thông số GIS kết nối dữ liệu tải lượng xả thải với lưu lượng dòng chảy để mô phỏng sự lan truyền chất ô nhiễm.

Các khái niệm chuyên ngành cốt lõi được chuẩn hóa bao gồm: Chu trình chuyển hóa Nitơ trong nước mặt, Tải lượng Nitơ phát thải (tính bằng tấn N/năm), Dạng tồn tại của Nitơ vô cơ (NH4+-N, NO2--N, NO3--N), Phân tích độ nhạy mô hình và Phân phối xác suất Monte Carlo.

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm và dữ liệu thứ cấp đồng bộ:

  • Nguồn dữ liệu quan trắc: Số liệu quan trắc chất lượng nước thực hiện liên tục trong 3 năm (2011 - 2013) với tần suất 4 đợt/năm tại 31 điểm quan trắc cố định trên 2 hệ thống sông chính (19 điểm trên sông Thái Bình và 12 điểm trên hệ thống Bắc Hưng Hải). Ngoài ra, nghiên cứu bổ sung số liệu tại 3 điểm trên sông Luộc và 21 điểm kênh mương nội đồng.
  • Phương pháp phân tích hóa nghiệm: Mẫu nước được thu thập, bảo quản và phân tích tại phòng thí nghiệm đạt chuẩn quốc tế ISO/IEC 17025:2005 thuộc Trung tâm Quan trắc và Phân tích Môi trường Hải Dương. Các chỉ tiêu được xác định bằng phương pháp trắc phổ hấp thụ phân tử trên máy quang phổ UV-VIS Shimadzu (dải bước sóng 190 - 800 nm), tuân thủ nghiêm ngặt các tiêu chuẩn quốc gia: TCVN 6179-1:1996 (xác định amoni bằng thuốc thử natri salixylat), TCVN 6178:1996 (xác định nitrit bằng thuốc thử Griess A và Griess B), TCVN 6180:1996 (xác định nitrat bằng axit sunfosalixylic).
  • Timeline nghiên cứu: Kế hoạch triển khai chia làm 4 giai đoạn logic từ năm 2011 đến năm 2014, gồm: thu thập dữ liệu cơ sở; khảo sát thực địa đợt 1 để xây dựng mô hình sơ bộ; khảo sát thực địa đợt 2 để thu thập thông số chi tiết và phân tích độ nhạy; hoàn thiện mô hình và mô phỏng kịch bản đến năm 2020.
  • Cỡ mẫu và phương pháp chọn mẫu: Cỡ mẫu thực nghiệm gồm 372 lượt mẫu nước mặt được lựa chọn theo phương pháp phân tầng không gian tại các vị trí nhạy cảm như hợp lưu sông, thượng lưu trạm cấp nước sinh hoạt, hạ lưu các khu công nghiệp và cống xả thải đô thị.
  • Lý do lựa chọn phương pháp phân tích: Sự kết hợp giữa thuật toán lặp Monte Carlo với 1.000 vòng lặp trên công cụ Macro Excel và mô hình GIS cho phép xử lý tối ưu các khoảng sai số dữ liệu thứ cấp, đảm bảo mô hình phản ánh sát thực tế mà không đòi hỏi chi phí đo đạc quá lớn.

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

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

Phân tích chuỗi số liệu quan trắc từ 31 điểm lấy mẫu và kết quả chạy mô hình MFA đã mang lại những phát hiện quan trọng về cấu trúc ô nhiễm Nitơ:

  • Ô nhiễm Amoni (NH4+-N) ở mức báo động trên cả hai hệ thống sông: Trên hệ thống sông Bắc Hưng Hải, nồng độ NH4+-N tại hầu hết các điểm đều vượt ngưỡng B1 của QCVN 08:2008/BTNMT. Điển hình tại điểm Sđ25 (Cầu Sặt), Sđ26 (Cầu Cất) và Sđ27 (hạ lưu cống xả nước thải thành phố Hải Dương), nồng độ amoni dao động từ 1,5 mg/L đến trên 4,5 mg/L, vượt từ 3 đến 9 lần quy chuẩn nguồn nước phục vụ sinh hoạt (mức A1 và A2). Trên sông Thái Bình, các điểm hợp lưu thượng nguồn như S1, S2 (sông Lục Nam) và S3 (sông Thương) ghi nhận nồng độ tăng đột biến vào đợt 2 năm 2013, vượt ngưỡng 1,2 mg/L do tiếp nhận tải lượng ô nhiễm từ các tỉnh phía thượng nguồn đổ về.
  • Nồng độ Nitrit (NO2--N) tăng cao bất thường: Hàm lượng NO2--N trên hệ thống sông Thái Bình và Bắc Hưng Hải có xu hướng dao động phức tạp và không tuân theo chu kỳ mùa cố định. Đa số các điểm quan trắc có hàm lượng vượt mức A2 (0,02 mg/L), nhiều thời điểm vượt ngưỡng 0,05 mg/L (vượt mức B1), cho thấy quá trình oxy hóa sinh hóa các hợp chất hữu cơ diễn ra mạnh mẽ nhưng không triệt để, tiềm ẩn nguy cơ cao đối với chất lượng nước đầu vào của các nhà máy nước sạch Cẩm Thượng và Phú Thái.
  • Định lượng nguồn phát thải qua mô hình MFA: Kết quả mô hình hóa chỉ ra rằng hoạt động sinh hoạt hộ gia đình và sản xuất nông nghiệp (trồng trọt trên 68.542 ha lúa và chăn nuôi gia súc) là hai nguồn đóng góp lớn nhất, chiếm hơn 60% tổng lượng Nitơ phát thải vào môi trường nước mặt. Đặc biệt, hệ thống công trình vệ sinh tại chỗ (bể tự hoại truyền thống) chỉ lưu giữ được từ 5% đến 14% lượng Nitơ trong bùn phân, trong khi đó khoảng 86% đến 95% lượng Nitơ còn lại tồn tại dưới dạng hòa tan trong nước thải và xả thẳng ra hệ thống thoát nước chung.
  • Mô hình phân phối tải lượng GIS: Đối với 4 địa bàn trọng điểm (huyện Bình Giang, Cẩm Giàng, Gia Lộc và thành phố Hải Dương), mô hình tính toán tải lượng Nitơ thiết lập hàm tương quan phi tuyến giữa tải lượng xả thải và lưu lượng dòng chảy với hệ số tương quan xác định R2 đạt trên 0,85, khẳng định tính tương thích cao của mô hình mô phỏng.

Thảo luận kết quả

Các phát hiện thực nghiệm được minh chứng rõ rệt khi giải thích dưới góc độ phát triển kinh tế - xã hội và cơ sở hạ tầng. Tốc độ đô thị hóa nhanh tại Hải Dương nâng tỷ lệ dân số thành thị lên 22,1% (với 385.513 người), nhưng hạ tầng xử lý nước thải sinh hoạt chưa phát triển tương xứng. Trạm bơm Âu Thuyền tại thành phố Hải Dương với công suất vận hành 32.000 m3/giờ thực hiện bơm xả trực tiếp nước thải chưa qua xử lý sinh học ra sông Thái Bình và một phần vào sông Sặt, tạo nên các dải ô nhiễm cục bộ nghiêm trọng. Bên cạnh đó, hoạt động thâm canh lúa nước và chăn nuôi nông hộ với tập quán sử dụng phân bón hóa học vượt định mức dẫn đến lượng Nitơ dư thừa bị rửa trôi qua các đợt mưa lớn (lượng mưa trung bình năm đạt 1.650 mm, tập trung 80% - 85% vào mùa mưa từ tháng 5 đến tháng 10).

Khi so sánh với các nghiên cứu dòng vật chất tại Việt Trì và khu vực Hà Nội, kết quả nghiên cứu này hoàn toàn đồng nhất về quy luật phân bổ: nguồn thải sinh hoạt nông thôn và đô thị luôn chi phối phần lớn lượng dưỡng chất Nitơ và Phốt pho thoát vào lưu vực sông.

Trong báo cáo nghiên cứu, các tập dữ liệu được trình bày tối ưu qua hệ thống biểu đồ và bản đồ không gian:

  • Biểu đồ tần suất phân phối xác suất Monte Carlo thể hiện rõ độ lệch chuẩn và khoảng tin cậy 95% của các dòng tải lượng N từ hộ gia đình sang hệ thống thoát nước.
  • Bản đồ chuyên đề GIS thể hiện mật độ phát thải N (tính theo tấn N/km2/năm) phân tầng theo từng tiểu vùng địa hình từ cao trình 1,0 m đến 4,0 m.
  • Đồ thị đường hồi quy phi tuyến mô tả mối tương quan chặt chẽ giữa lưu lượng xả từ trạm bơm và nồng độ chất ô nhiễm tại các mặt cắt ngang sông.

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

Dựa trên kết quả định lượng dòng vật chất và mô hình phân phối GIS, luận văn đề xuất 4 nhóm giải pháp kỹ thuật và quản trị chiến lược:

  • 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 tỉnh Hải Dương cùng Sở Xây dựng cần đẩy nhanh tiến độ triển khai các nhà máy xử lý nước thải sinh hoạt tập trung tại thành phố Hải Dương và các thị trấn trung tâm; đặt mục tiêu thu gom và xử lý sinh học đạt 70% tổng lượng nước thải đô thị trước năm 2020, nhằm cắt giảm ít nhất 25% tải lượng amoni xả trực tiếp vào lưu vực sông Sặt và sông Thái Bình.
  • Cải tạo công trình vệ sinh và áp dụng công nghệ tách dòng nước thải nông thôn: Sở Nông nghiệp và Phát triển nông thôn chủ trì phối hợp cùng Sở Khoa học và Công nghệ chuyển giao mô hình nhà vệ sinh tách nước tiểu và bể tự hoại cải tiến có ngăn lọc kỵ khí tại các khu dân cư nông thôn; đặt mục tiêu tăng tỷ lệ thu hồi và tái sử dụng dưỡng chất Nitơ làm phân bón hữu cơ lên mức 35% đến 40% trong giai đoạn 2015 - 2020.
  • Kiểm soát dư lượng phân bón và tối ưu hóa canh tác nông nghiệp: Chi cục Trồng trọt và Bảo vệ thực vật Hải Dương chỉ đạo các hợp tác xã áp dụng quy trình quản lý dinh dưỡng cây trồng tổng hợp trên toàn bộ 68.542 ha đất canh tác lúa; đặt mục tiêu cắt giảm từ 15% đến 20% lượng phân đạm hóa học bón dư thừa trước năm 2018, hạn chế tối đa lượng Nitơ rửa trôi vào hệ thống kênh mương thủy lợi Bắc Hưng Hải.
  • Hiện đại hóa mạng lưới giám sát và tích hợp cảnh báo sớm trên nền tảng GIS: Sở Tài nguyên và Môi trường tỉnh Hải Dương cần đầu tư lắp đặt bổ sung 10 trạm quan trắc tự động liên tục tại các điểm nóng ô nhiễm dọc sông Thái Bình và kênh chính Bắc Hưng Hải trước năm 2017; tích hợp toàn bộ dữ liệu thủy văn và tải lượng xả thải lên phần mềm quản lý GIS phục vụ công tác thanh kiểm tra xả thải của các khu, cụm công nghiệp.

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

Nội dung và phương pháp luận của công trình nghiên cứu mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước về môi trường và quy hoạch đô thị: Cán bộ thuộc Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn, Ban Quản lý các Khu công nghiệp tỉnh Hải Dương và các tỉnh thuộc lưu vực sông Thái Bình có thể ứng dụng trực tiếp mô hình MFA-GIS để xây dựng báo cáo sức chịu tải nguồn nước, thiết lập hạn ngạch cấp phép xả thải và quy hoạch phân vùng bảo vệ nguồn nước mặt.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Các nhà nghiên cứu thuộc chuyên ngành Hóa môi trường, Kỹ thuật môi trường, Quản lý tài nguyên nước và Khoa học Trái đất có thể sử dụng luận văn như tài liệu tham khảo chuẩn mực về phương pháp tích hợp mô hình dòng chảy vật chất với công nghệ thông tin địa lý và kỹ thuật mô phỏng Monte Carlo trong điều kiện dữ liệu hạn chế.
  • Các công ty cấp thoát nước và doanh nghiệp hạ tầng môi trường: Đội ngũ kỹ sư tại các công ty khai thác công trình thủy lợi Bắc Hưng Hải và các nhà máy nước sạch sinh hoạt có thể khai thác số liệu nồng độ các dạng Nitơ để điều chỉnh công nghệ tiền xử lý nước thô, đồng thời tối ưu hóa vị trí đặt cửa lấy nước an toàn.
  • Các tổ chức phi chính phủ và chuyên gia phát triển nông nghiệp tuần hoàn: Các chuyên gia môi trường có thể tham khảo mô hình phân tách dòng dưỡng chất để xây dựng các dự án thí điểm về vệ sinh sinh thái, xử lý phân bùn bể phốt và phục hồi dinh dưỡng Nitơ cho cây trồng tại khu vực nông thôn đồng bằng Bắc Bộ.

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

  • Tại sao mô hình dòng chảy vật chất (MFA) lại phù hợp để đánh giá ô nhiễm nitơ tại Hải Dương? Mô hình MFA dựa trên nguyên tắc cân bằng khối lượng giúp kết nối định lượng toàn bộ dòng vật chất từ khâu phát sinh đến khâu xả thải cuối cùng. Với đặc thù địa bàn Hải Dương có 68.542 ha đất lúa và hơn 1,51 triệu dân sinh sống đan xen, MFA cho phép phân lập rõ nguồn thải sinh hoạt và nông nghiệp chiếm ưu thế áp đảo mà không phụ thuộc vào chuỗi số liệu quan trắc dày đặc vốn rất tốn kém.

  • Nồng độ Amoni (NH4+-N) tại hệ thống sông Bắc Hưng Hải thực tế ghi nhận mức ô nhiễm ra sao? Chuỗi quan trắc liên tục giai đoạn 2011 - 2013 tại 12 điểm trên hệ thống Bắc Hưng Hải cho thấy nồng độ amoni thường xuyên vượt mức B1 theo chuẩn QCVN 08:2008/BTNMT. Điển hình tại khu vực Cầu Sặt (Sđ25) và sau điểm xả thải thành phố Hải Dương (Sđ27), hàm lượng amoni ghi nhận từ 3,0 mg/L đến trên 4,5 mg/L, vượt từ 3 đến 9 lần tiêu chuẩn cấp nước sinh hoạt.

  • Việc tích hợp công nghệ GIS mang lại ưu thế gì vượt trội so với cách tiếp cận truyền thống? Hệ thống GIS số hóa toàn diện bản đồ hiện trạng sử dụng 104.238 ha đất nông nghiệp và vị trí 18 khu công nghiệp, liên kết với dữ liệu chất lượng nước tại 31 điểm quan trắc. Nhờ đó, mô hình biểu diễn trực quan sự phân bố tải lượng ô nhiễm theo không gian và mô phỏng lan truyền chất với độ tương quan R2 đạt trên 0,85.

  • Đâu là nguồn phát sinh tải lượng Nitơ lớn nhất làm suy thoái nguồn nước mặt địa phương? Nghiên cứu khẳng định sinh hoạt hộ gia đình và canh tác nông nghiệp là hai nguồn thải lớn nhất, chiếm hơn 60% tổng lượng Nitơ phát tán vào nguồn nước mặt. Bể tự hoại truyền thống chỉ giữ lại 5% đến 14% lượng Nitơ trong bùn lắng, trong khi 86% đến 95% lượng Nitơ dạng lỏng hòa tan tiếp tục thoát ra hệ thống kênh mương.

  • Kịch bản mô phỏng tải lượng ô nhiễm Nitơ đến năm 2020 đưa ra những cảnh báo quan trọng nào? Nếu không có biện pháp can thiệp xử lý nước thải sinh hoạt và kiểm soát phân bón, tải lượng Nitơ xả vào nước mặt đến năm 2020 ước tính sẽ tăng thêm khoảng 15% đến 25% theo đà tăng dân số đô thị và sản xuất công nghiệp. Điều này gia tăng áp lực ô nhiễm lên các nhà máy nước sạch Cẩm Thượng và Phú Thái.

Kết luận

  • Luận văn đã xây dựng thành công mô hình Phân tích Dòng chảy Vật chất (MFA) đối với chu trình Nitơ trên quy mô toàn tỉnh Hải Dương, tích hợp 1.000 vòng lặp mô phỏng ngẫu nhiên Monte Carlo để kiểm soát độ không đảm bảo đo.
  • Đánh giá toàn diện hiện trạng chất lượng nước qua 372 lượt mẫu tại 31 trạm quan trắc trên hệ thống sông Thái Bình và Bắc Hưng Hải giai đoạn 2011 - 2013, chỉ rõ sự vượt chuẩn nghiêm trọng của amoni (từ 1,5 đến 4,5 mg/L) và nitrit (vượt 0,02 mg/L).
  • Xác định rõ cơ cấu phát thải: Hoạt động sinh hoạt dân cư và trồng trọt nông nghiệp đóng góp trên 60% tổng lượng Nitơ vào nguồn nước mặt, trong đó bể tự hoại để thất thoát tới 86% - 95% lượng Nitơ ra môi trường.
  • Thiết lập mô hình phân phối thông số GIS tính toán tải lượng Nitơ cho 4 địa bàn trọng điểm (Bình Giang, Cẩm Giàng, Gia Lộc, TP Hải Dương) đạt hệ số tương quan R2 trên 0,85, mở ra hướng ứng dụng hiệu quả trong quy hoạch bảo vệ môi trường nước.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật, hạ tầng và chính sách với lộ trình thực hiện cụ thể từ năm 2015 đến năm 2020 nhằm cắt giảm từ 20% đến 30% tải lượng ô nhiễm.

Lộ trình triển khai các giải pháp kiểm soát nguồn thải Nitơ cần được ưu tiên thực hiện ngay trong giai đoạn 2015 - 2020 nhằm bảo đảm an ninh nguồn nước cho toàn vùng hạ lưu châu thổ sông Hồng. Các cơ quan quản lý, viện nghiên cứu và doanh nghiệp hãy chủ động ứng dụng công cụ MFA-GIS vào công tác quản trị lưu vực sông để kiến tạo môi trường nước bền vững cho tương lai!