Tổng quan nghiên cứu

Theo thống kê của ngành nông nghiệp, mỗi năm hoạt động chăn nuôi gia súc, gia cầm tại Việt Nam phát sinh khoảng 75 đến 85 triệu tấn chất thải rắn cùng hàng trăm triệu mét khối nước thải. Tuy nhiên, chỉ có khoảng 8,7% số hộ chăn nuôi xây dựng công trình khí sinh học và trên 23% cơ sở xả thải trực tiếp ra môi trường tự nhiên không qua xử lý. Tại tỉnh Bắc Giang, địa phương có tổng đàn lợn đứng thứ 3 toàn quốc với quy mô đạt trên 1,13 triệu con vào năm 2016, sức ép ô nhiễm môi trường nông thôn ngày càng trở nên gay gắt.

Huyện Yên Dũng có tổng đàn lợn xấp xỉ 80.000 con phân bố trên 19 xã và 2 thị trấn, chủ yếu theo phương thức nông hộ và gia trại nhỏ lẻ. Sự phân tán chất thải mang tính chất nguồn diện, dòng thải tự nhiên tích lũy cục bộ theo địa hình dòng chảy, trong khi công tác quản lý môi trường của địa phương vẫn còn mang tính dàn trải.

Luận văn thạc sĩ chuyên ngành Khoa học Môi trường của tác giả Ngô Phương Lan tại Học viện Nông nghiệp Việt Nam giải quyết trọn vẹn bài toán này thông qua việc đánh giá thực trạng phát sinh, xử lý chất thải và ứng dụng công nghệ GIS kết hợp mô hình hóa thủy văn để mô phỏng dòng vận chuyển chất thải chăn nuôi lợn.

Phạm vi nghiên cứu bao quát toàn bộ 19.093,04 ha diện tích tự nhiên của huyện Yên Dũng với chuỗi số liệu thực địa thu thập từ tháng 2/2017 đến tháng 4/2018. Kết quả nghiên cứu cung cấp bộ công cụ định lượng tải lượng ô nhiễm COD và BOD5 cho từng tiểu lưu vực, đạt hệ số biến thiên sai số kiểm chứng CV chỉ 8,7%, tạo luận cứ khoa học chuẩn xác phục vụ quy hoạch môi trường và quản lý rủi ro ô nhiễm tại đị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 tích hợp 3 khung lý thuyết và mô hình khoa học chuyên sâu:

  • Khung phân tích sinh kế bền vững của DFID: Đánh giá tác động của 5 nguồn vốn sinh kế hộ gia đình đến quyết định đầu tư công nghệ xử lý môi trường.
  • Lý thuyết mô hình hóa thủy văn và lan truyền chất ô nhiễm nguồn diện: Ứng dụng nguyên lý trọng lực và hướng dòng chảy mặt để xác định ranh giới lưu vực gom nước cùng điểm tích tụ ô nhiễm.
  • Mô hình cân bằng dòng vật chất và lưu chuyển chất dinh dưỡng của Vũ Chí Cương: Định lượng phát thải chất hữu cơ, nitơ và phốt pho qua các dòng thải lỏng, bùn cặn và khí sinh học.

Các khái niệm chuyên ngành then chốt bao gồm:

  • Chất ô nhiễm nguồn diện: Nguồn thải phân tán từ hàng nghìn chuồng nuôi hộ gia đình không có điểm xả cố định.
  • Tiểu lưu vực thoát nước: Đơn vị không gian thủy văn khép kín được phân định từ mô hình độ cao số, nơi toàn bộ dòng chảy mặt đổ về một cửa thoát chung.
  • Tải lượng ô nhiễm COD và BOD5: Khối lượng nhu cầu oxy hóa học và nhu cầu oxy sinh hóa tích lũy trên một hecta diện tích trong một năm (kg/ha/năm).
  • Hiệu suất xử lý khí sinh học: Tỷ lệ làm sạch các chất ô nhiễm sau hầm biogas, đạt 81% đối với COD và 86% đối với BOD5 theo thực nghiệm kiểm toán môi trường.

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

Nghiên cứu sử dụng hệ thống dữ liệu đa tầng bao gồm mô hình số độ cao DEM từ ảnh vệ tinh ALOS PALSAR-2, ảnh vệ tinh quang học Sentinel-2 chụp ngày 03/06/2017 để trích xuất lớp phủ khu dân cư, cùng số liệu thống kê chăn nuôi chính thức năm 2017 của 19 xã và 2 thị trấn do Phòng Nông nghiệp và Phát triển nông thôn huyện Yên Dũng cung cấp. Định mức phát thải cơ sở được tính toán với hệ số xả thải 6,01 m3 nước thải/năm/con nái, tương ứng 8,66 kg COD/năm và 5,28 kg BOD5/năm.

Cỡ mẫu điều tra thực địa gồm 90 hộ chăn nuôi được lựa chọn theo phương pháp chọn mẫu phân tầng ngẫu nhiên tại 3 xã đại diện cho 3 cấp độ mật độ chăn nuôi: xã Quỳnh Sơn (mật độ cao), xã Tiến Dũng (mật độ trung bình) và xã Tân Liễu (mật độ thấp). Cỡ mẫu này đảm bảo thu thập đầy đủ các biến số sinh kế, quy mô đàn và hiện trạng công trình khí sinh học.

Phương pháp phân tích hồi quy Logistic đa biến trên phần mềm SPSS 16.0 được lựa chọn nhằm dự báo chính xác xác suất áp dụng hầm biogas của hộ nuôi dựa trên quy mô đàn và điều kiện kinh tế. Trên phần mềm ArcGIS 10.3, công cụ Hydrology phân tích hướng dòng (Flow Direction) và tích lũy dòng (Flow Accumulation) từ dữ liệu DEM để vạch ranh giới các tiểu lưu vực. Mô hình được kiểm chứng độ tin cậy thông qua tiêu chuẩn thông tin Bayes (BIC) và đối chiếu với dữ liệu quan trắc nồng độ nước mặt tại 15 điểm tiếp nhận thực tế. Toàn bộ quá trình nghiên cứu diễn ra liên tục trong 14 tháng (từ tháng 2/2017 đến tháng 4/2018).

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

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

Mô hình hóa không gian đã mang lại 4 phát hiện quan trọng về quy luật phân tán chất thải chăn nuôi lợn trên địa bàn huyện Yên Dũng:

  • Tích tụ ô nhiễm cục bộ tại vùng trũng thấp: Chất thải từ các xã phía trên dồn về cuối các dòng chảy thuộc địa phận xã Đồng Phúc và xã Tân Liễu. Tải lượng ô nhiễm tại hai xã này ở mức cao nhất toàn huyện, với chỉ số COD dao động từ 10 đến 60 kg/ha/năm và BOD5 đạt từ 7 đến 34 kg/ha/năm.
  • Xác định các vùng an toàn sinh thái: Các tiểu lưu vực thuộc địa bàn xã Tân Phong, Nham Sơn, Đồng Việt, Trí Yên và thị trấn Neo có mức phát tán rất thấp, tải lượng COD và BOD5 đều duy trì dưới ngưỡng 2 kg/ha/năm.
  • Tác động nghịch đảo giữa mật độ nuôi và tỷ lệ xử lý chất thải: Tại các tiểu lưu vực có mật độ nuôi tập trung cao trên 50 con/hộ nhưng tỷ lệ hộ xây dựng bể biogas chỉ đạt dưới 35%, nguy cơ ô nhiễm môi trường nước mặt tăng gấp 4 lần. Ngược lại, những khu vực có tỷ lệ hộ lắp đặt biogas đạt trên 70% thì mức tải lượng ô nhiễm phát tán ra môi trường giảm tới 80 đến 85%, đưa môi trường nước về ngưỡng an toàn.
  • Độ chính xác vượt trội của mô hình: Kết quả mô phỏng trên phần mềm ArcGIS khi kiểm chứng với dữ liệu thực địa cho hệ số biến thiên sai số CV đạt 8,7%, nằm hoàn toàn trong khoảng tin cậy 95% của phân phối thống kê Bayes.

Thảo luận kết quả

Quy luật tích tụ ô nhiễm tại Yên Dũng bắt nguồn trực tiếp từ đặc điểm địa hình và kinh tế - xã hội. Huyện có 82,90% diện tích đất có độ dốc dưới 3 độ, nằm kẹp giữa dãy núi Nham Biền cao 254 m và hệ thống 3 con sông lớn (sông Cầu, sông Thương, sông Lục Nam). Địa hình lòng chảo khiến chất thải từ vùng gò đồi nhanh chóng bị rửa trôi và lắng đọng tại các vùng trũng ven sông. Đồng thời, khu vực nông thôn chiếm tới 91,24% dân số toàn huyện (120.730 người trên tổng số 132.395 người năm 2016), tỷ lệ lao động nông nghiệp chiếm 44,39%, dẫn đến việc đầu tư hệ thống xử lý nước thải hoàn chỉnh sau biogas gặp khó khăn do hạn chế về nguồn vốn tự có.

Khi so sánh với các nghiên cứu mô hình hóa lưu vực tại Thái Lan của Thanapongtharm (2016) hay nghiên cứu tải lượng nitrate sông Raccoon tại Hoa Kỳ của Jha (2010), kết quả luận văn khẳng định tiếp cận theo ranh giới tiểu lưu vực tự nhiên mang lại độ chuẩn xác cao hơn nhiều so với việc thống kê theo ranh giới hành chính.

Dữ liệu nghiên cứu được trực quan hóa thành công qua hệ thống bản đồ số đa lớp (Raster Map) với thang cảnh báo 5 cấp màu sắc từ xanh lá cây (dưới 2 kg/ha/năm) đến đỏ đậm (50 - 60 kg/ha/năm), kết hợp bảng ma trận thống kê tương quan không gian giúp nhà quản lý nhận diện tức thì các điểm nóng ô nhiễm cần ưu tiên can thiệp.

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

Dựa trên kết quả mô phỏng không gian dòng thải, luận văn đưa ra 4 nhóm giải pháp chiến lược:

  1. Chuẩn hóa và nâng cấp công nghệ khí sinh học nông hộ:
  • Vận động và hỗ trợ kỹ thuật để các hộ chăn nuôi chuyển đổi sang hầm biogas composite hoặc phủ bạt HDPE quy mô lớn, kết hợp bể lắng kỵ khí UASB và hồ sinh thái lọc thực vật.
  • Mục tiêu: Nâng tỷ lệ hộ chăn nuôi có hệ thống xử lý đạt chuẩn từ dưới 40% lên 85% vào năm 2025, cắt giảm tải lượng COD phát tán ra các kênh mương xuống dưới 5 kg/ha/năm.
  • Đơn vị thực hiện: Phòng Nông nghiệp & PTNT phối hợp cùng Phòng Tài nguyên & Môi trường huyện Yên Dũng, triển khai giai đoạn 2018 - 2023.
  1. Quy hoạch phân vùng chăn nuôi theo sức tải môi trường tiểu lưu vực:
  • Di dời triệt để các cơ sở chăn nuôi quy mô từ 50 con trở lên nằm xen kẽ trong khu dân cư tại các xã vùng trũng Đồng Phúc và Tân Liễu vào 91 khu chăn nuôi tập trung xa khu dân cư theo quy hoạch của tỉnh Bắc Giang.
  • Mục tiêu: Giảm mật độ đàn tại các điểm trũng nhạy cảm xuống dưới 10 con/ha đất tự nhiên, chấm dứt 100% tình trạng xả tràn phân tươi chưa xử lý ra môi trường.
  • Đơn vị thực hiện: Ủy ban nhân dân huyện Yên Dũng và UBND các xã, thực hiện trong giai đoạn 2019 - 2025.
  1. Ứng dụng công nghệ đệm lót sinh học và ủ phân hữu cơ vi sinh:
  • Triển khai sử dụng chế phẩm vi sinh hữu hiệu (Balasa No1, EM) làm đệm lót chuồng kỵ nước và kỹ thuật ủ compost phân rắn kết hợp phụ phẩm rơm rạ, trấu.
  • Mục tiêu: Tái chế 75% lượng phân thải rắn thành phân bón hữu cơ chất lượng cao cung cấp cho 14 cánh đồng mẫu lớn của huyện, giảm trên 65% nồng độ phát thải các khí độc hại NH3 và H2S.
  • Đơn vị thực hiện: Trạm Khuyến nông và Trạm Chăn nuôi Thú y huyện hướng dẫn trực tiếp cho các hộ dân từ năm 2018.
  1. Xây dựng cơ sở dữ liệu WebGIS quản lý nguồn thải theo thời gian thực:
  • Thiết lập hệ thống bản đồ số trực tuyến định danh 100% hộ chăn nuôi trên nền tảng WebGIS và duy trì 10 trạm quan trắc nước mặt định kỳ tại các cửa thoát tiểu lưu vực trọng điểm.
  • Mục tiêu: Kiểm soát toàn diện chất lượng nước mặt, đảm bảo các chỉ tiêu BOD5 và COD đạt Quy chuẩn kỹ thuật quốc gia QCVN 08-MT:2015/BTNMT (cột B1).
  • Đơn vị thực hiện: Sở Tài nguyên và Môi trường tỉnh Bắc Giang chủ trì, cập nhật dữ liệu 6 tháng một lần trong giai đoạn 2020 - 2025.

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

Luận văn là tài liệu tham khảo học thuật và thực tiễn giá trị cho 4 nhóm đối tượng:

  • Cán bộ quản lý môi trường và nhà hoạch định chính sách cấp huyện, tỉnh: Cung cấp bản đồ phân vùng ô nhiễm chi tiết cho 19 xã và 2 thị trấn, giúp Sở TN&MT và UBND cấp huyện phân bổ kinh phí sự nghiệp môi trường đúng trọng tâm cho các tiểu lưu vực chịu tải COD trên 10 kg/ha/năm thay vì dàn trải ngân sách.
  • Nhà nghiên cứu, giảng viên và học viên cao học ngành Khoa học Môi trường, Địa tin học: Cung cấp quy trình tích hợp mô hình số độ cao DEM ALOS PALSAR-2, ảnh vệ tinh Sentinel-2 và phân tích hồi quy Logistic trên SPSS và ArcGIS với hệ số kiểm định sai số CV chỉ 8,7%, đóng vai trò như một bài giảng mẫu về mô phỏng ô nhiễm nguồn diện.
  • Kỹ sư môi trường và đơn vị tư vấn xử lý chất thải nông nghiệp: Khai thác bộ thông số thực nghiệm chuẩn xác về định mức xả thải (6,01 m3/năm/nái) cùng hiệu suất xử lý (81% COD, 86% BOD5) để thiết kế các modul xử lý nước thải phân tán chi phí thấp phù hợp địa hình nông thôn miền núi.
  • Chủ trang trại, gia trại và hợp tác xã nông nghiệp tuần hoàn: Nắm bắt quy trình kỹ thuật chuyển hóa chất thải rắn thành phân bón vi sinh phục vụ canh tác trên 15.784,37 ha đất nông nghiệp toàn huyện, giúp tiết kiệm 20 đến 30% chi phí phân bón hóa học và triệt tiêu rủi ro vi phạm quy định xả thải.

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

Phương pháp mô phỏng theo tiểu lưu vực có ưu thế gì vượt trội so với quản lý theo ranh giới hành chính truyền thống? Quản lý theo xã phường thường bỏ qua quy luật dòng chảy tự nhiên, dẫn đến việc không phát hiện được hiện tượng ô nhiễm tích tụ xuyên địa giới. Tiếp cận theo tiểu lưu vực thủy văn xác định chính xác các điểm trũng như xã Đồng Phúc và Tân Liễu đang chịu tải lượng COD lên tới 10 đến 60 kg/ha/năm, giúp tối ưu hóa nguồn lực xử lý đúng vị trí tích lũy.

Mô hình phân tán dòng thải trên phần mềm ArcGIS đạt độ tin cậy khoa học ở mức độ nào? Mô hình đạt độ chuẩn xác rất cao khi được kiểm định bằng tiêu chuẩn thông tin Bayes (BIC) và đối chiếu với dữ liệu quan trắc nồng độ thực địa tại 15 điểm tiếp nhận. Kết quả cho thấy hệ số biến thiên sai số CV chỉ là 8,7% và toàn bộ số liệu đều nằm trọn vẹn trong khoảng tin cậy 95%.

Tại sao nhiều hộ chăn nuôi đã lắp đặt hầm biogas nhưng nguồn nước xung quanh vẫn bị ô nhiễm? Dù hầm biogas đạt hiệu suất xử lý COD 81% và BOD5 86%, nhưng nhiều hộ nuôi quy mô trên 30 đến 50 con có lượng nước thải xả ra vượt quá dung tích hầm ủ. Nước thải sau biogas chưa qua bể lắng UASB hoặc hồ sinh học vẫn mang hàm lượng hữu cơ rất cao, khi xả thẳng ra kênh mương thủy lợi sẽ gây hiện tượng phú dưỡng tức thì.

Giải pháp đệm lót sinh học có thực sự phát huy hiệu quả tại các vùng nông thôn miền núi? Công nghệ đệm lót sinh học sử dụng men vi sinh tuyển chọn (như Balasa No1, EM) rất phù hợp cho quy mô nông hộ và gia trại tại các vùng đất cao ráo. Phương pháp này phân giải phân ngay tại chuồng, giảm trên 65% phát thải khí mùi hôi NH3 và H2S, đồng thời tạo ra nguồn phân hữu cơ giàu dinh dưỡng cho 14 cánh đồng mẫu lớn.

Ảnh vệ tinh Sentinel-2 giải quyết bài toán thiếu hụt tọa độ hộ chăn nuôi như thế nào? Do số liệu thống kê không lưu tọa độ GPS của từng hộ, ảnh vệ tinh Sentinel-2 chụp ngày 03/06/2017 được giải đoán để bóc tách chính xác ranh giới lớp phủ khu dân cư. Kết hợp thuật toán tạo điểm ngẫu nhiên có trọng số sinh kế từ 90 mẫu khảo sát, mô hình đã định vị chính xác vị trí xả thải của xấp xỉ 80.000 con lợn trên toàn huyện.

Kết luận

Công trình nghiên cứu của tác giả Ngô Phương Lan đã mang lại những đóng góp khoa học và giá trị thực tiễn nổi bật:

  • Định lượng chính xác bức tranh phân tán ô nhiễm với tải lượng COD từ 10 đến 60 kg/ha/năm và BOD5 từ 7 đến 34 kg/ha/năm tại các tiểu lưu vực điểm nóng thuộc xã Đồng Phúc và Tân Liễu.
  • Xây dựng thành công quy trình công nghệ tích hợp viễn thám Sentinel-2, mô hình số độ cao ALOS PALSAR-2 và hồi quy Logistic trên ArcGIS với độ tin cậy cao, hệ số biến thiên CV chỉ 8,7%.
  • Làm sáng tỏ quy luật tương tác giữa mật độ chăn nuôi, tỷ lệ che phủ công trình biogas (hiệu suất xử lý COD 81%, BOD5 86%) và tải lượng tích lũy ô nhiễm nguồn diện.
  • Đề xuất đồng bộ 4 nhóm giải pháp kỹ thuật và quy hoạch khả thi, hướng tới mục tiêu nâng tỷ lệ xử lý chất thải đạt chuẩn lên 85% và tái chế 75% phụ phẩm rắn đến năm 2025.
  • Mở ra định hướng ứng dụng mở rộng giai đoạn 2023 - 2030 trong việc tích hợp mô hình thủy văn SWAT và giám sát phát thải trực tuyến cho toàn bộ lưu vực sông Cầu và sông Thương.

Hãy tham khảo và ứng dụng ngay mô hình phân tích không gian từ luận văn thạc sĩ này để nâng cao hiệu quả quản trị môi trường và phát triển ngành chăn nuôi bền vững!