Tổng quan nghiên cứu

Theo báo cáo của Bộ Tài nguyên và Môi trường năm 2014, có trên 80% tổng lượng nước thải sinh hoạt tại Việt Nam chưa qua xử lý hoặc chỉ mới qua xử lý sơ bộ bằng bể tự hoại trước khi thải trực tiếp ra nguồn tiếp nhận như sông, hồ và kênh mương. Thực trạng này làm gia tăng nồng độ chất rắn lơ lửng, nhu cầu oxy sinh học (BOD5 đạt mức 178,67 ± 61,37 mg/l), tổng nitơ (TN đạt 74,30 ± 25,60 mg/l) và tổng photpho (TP đạt 15,98 ± 12,04 mg/l), dẫn đến hiện tượng suy giảm oxy hòa tan và gây phú dưỡng hóa nghiêm trọng các thủy vực tự nhiên. Trong khi các công nghệ hóa lý hiện đại như lọc màng siêu lọc hay thẩm thấu ngược có chi phí đầu tư và vận hành quá cao, công nghệ ao nuôi tảo hiệu suất cao (HRAPs) sử dụng vi tảo Chlorella vulgaris mở ra hướng tiếp cận sinh học bền vững với chi phí thấp và khả năng thu hồi tài nguyên sinh khối.

Mục tiêu cốt lõi của nghiên cứu là thiết kế, chế tạo và vận hành hệ thống xử lý nước thải HRAPs quy mô pilot; đồng thời xây dựng, tối ưu hóa và kiểm định mô hình toán học mô phỏng động học sinh trưởng của tảo Chlorella vulgaris cùng hiệu quả loại bỏ dinh dưỡng. Phạm vi nghiên cứu được thực hiện tại phòng thí nghiệm và khuôn viên Học viện Nông nghiệp Việt Nam trong mốc thời gian từ tháng 5/2015 đến tháng 8/2016. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn khi cung cấp bộ thông số công nghệ chuẩn xác, giúp giảm thiểu nồng độ các hợp chất ô nhiễm đạt chuẩn xả thải theo QCVN 14:2008/BTNMT Cột B, đồng thời tạo ra nguồn sinh khối giàu protein chiếm 40% đến 60% và lipid đạt 10% đến 15% trọng lượng khô phục vụ sản xuất nhiên liệu sinh họ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 hai nền tảng lý thuyết chính: lý thuyết sinh thái học vi tảo và lý thuyết mô hình hóa bể phản ứng quang sinh học. Cơ chế sinh học cốt lõi dựa trên hoạt động quang hợp tự dưỡng của vi tảo lục đơn bào Chlorella vulgaris (kích thước 2 đến 5 µm). Tảo sử dụng năng lượng ánh sáng thông qua hệ sắc tố diệp lục Chlorophyll-a (chiếm 2,2% trọng lượng khô) và Chlorophyll-b (0,58%) để đồng hóa nguồn carbon vô cơ, amoni (NH4+), nitrat (NO3-) và photphat (PO43-) thành sinh khối tế bào.

Bên cạnh đó, mô hình ao nuôi tảo hiệu suất cao (HRAPs) dạng đường đua (raceway) vận hành dựa trên nguyên lý kéo dài pha quang hợp, rút ngắn pha hô hấp và đảo trộn liên tục để tối ưu hóa khả năng tiếp xúc ánh sáng của tế bào. Về mặt mô hình hóa, tác giả áp dụng khung phân loại ba cấp độ của Bechet và cộng sự (2013) kết hợp mô hình động học quang hợp của Eilers và Peeters (1988) với ba trạng thái quang hợp: mở, đóng và ức chế quang học. Các khái niệm trung tâm bao gồm: bể phản ứng quang sinh học, đơn vị quang hợp (PSU), tỷ lệ dinh dưỡng N:P Redfield và chỉ số sai số căn bậc hai trung bình (RMSE).

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

Nghiên cứu kết hợp chặt chẽ giữa thực nghiệm trong phòng thí nghiệm và mô phỏng toán học trên máy tính. Dữ liệu thực nghiệm được thu thập từ ba chuỗi thí nghiệm độc lập lặp lại 3 lần cho mỗi công thức (tổng cộng 33 mẻ nuôi cấy) với mật độ tảo ban đầu 100.000 tế bào/ml trong thể tích 3,5 lít dưới điều kiện chiếu sáng 60 µmol/m²/s (tương đương 5000 lux) và chu kỳ sáng tối 16:8. Phương pháp chọn mẫu nước thải sau bể tự hoại được thực hiện ngẫu nhiên qua 20 đợt thu mẫu thực tế để sàng lọc ra 4 dải tỷ lệ N:P đặc trưng (5:1, 12:1, 16:1 và 30:1).

Các phương pháp phân tích chất lượng nước tuân thủ nghiêm ngặt theo Tiêu chuẩn Quốc gia (TCVN):

  • Đo pH bằng điện cực thủy tinh theo TCVN 6492:2011;
  • Đo oxy hòa tan (DO) bằng đầu đo điện hóa theo TCVN 7325:2004;
  • Xác định COD bằng phương pháp dicromat theo TCVN 6491:1999;
  • Xác định BOD5 theo TCVN 6001-2:2008;
  • Nồng độ N-NH4+ và N-NO3- được định lượng bằng quang phổ hấp thụ phân tử UV-VIS ở bước sóng tương ứng 667 nm và 430 nm;
  • Mật độ tế bào tảo được đếm trực tiếp trên kính hiển vi quang học Olympus CX21 bằng buồng đếm Burker-Turk chuyên dụng nhằm đảm bảo độ chuẩn xác đối với vi sinh vật kích thước hiển vi.

Các thuật toán giải hệ phương trình vi phân mô tả tốc độ sinh trưởng và suy giảm cơ chất được lập trình hoàn toàn trên nền tảng phần mềm MATLAB. Sai số giữa mô hình và số liệu đo thực tế được kiểm chứng qua chỉ số sai số RMSE kết hợp phân tích phương sai ANOVA một nhân tố và hậu kiểm Tukey trên phần mềm R ở mức ý nghĩa thống kê p = 0,05. Toàn bộ timeline nghiên cứu kéo dài liên tục 15 tháng đảm bảo thu thập trọn vẹn diễn biến theo mùa.

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

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

Thứ nhất, vi tảo Chlorella vulgaris thể hiện khả năng thích nghi và tăng trưởng vượt trội trong môi trường nước thải sau bể tự hoại so với môi trường dinh dưỡng chuẩn BBM. Giai đoạn thích nghi diễn ra rất nhanh (dưới 24 giờ). Giai đoạn tăng trưởng trong nước thải kéo dài từ ngày thứ 5 đến ngày thứ 10 và đạt đỉnh mật độ tế bào lên tới 7,0 × 10^6 tế bào/ml, cao gấp 2,5 lần so với môi trường nhân tạo BBM (đạt đỉnh 2,8 × 10^6 tế bào/ml) với mức ý nghĩa thống kê p < 0,05. Quá trình quang hợp mạnh mẽ đã thúc đẩy nồng độ oxy hòa tan trong nước thải tăng vọt từ mức ban đầu 0,2 mg/l lên 5,4 mg/l và đẩy giá trị pH từ 6,2 lên 8,7.

Thứ hai, tỷ lệ dinh dưỡng N:P đóng vai trò quyết định đến hiệu suất đồng hóa dinh dưỡng. Trong 4 dải tỷ lệ khảo sát, công thức có tỷ lệ N:P = 16:1 mang lại hiệu quả hấp thụ nitơ cao nhất, đạt hiệu suất loại bỏ tổng nitơ lên đến 96%, đồng thời hiệu suất loại bỏ tổng photpho đạt mức 80% trong dải N:P từ 1:1 đến 20:1.

Thứ ba, điều kiện chiếu sáng ảnh hưởng trực tiếp đến động thái tích lũy sinh khối. Thí nghiệm nguồn sáng đơn sắc cho thấy ánh sáng trắng và dải bước sóng đỏ (650–750 nm) kích thích tổng hợp Chlorophyll-a và gia tăng mật độ tảo nhanh hơn 35% so với ánh sáng xanh lục. Khi thử nghiệm hệ thống HRAPs pilot ngoài trời dưới bức xạ tự nhiên yếu, mật độ tảo chỉ đạt khoảng 3,0 × 10^6 tế bào/ml, trong khi hệ thống HRAPs lắp đặt dàn đèn LED nhân tạo 60W duy trì mật độ ổn định trên 6,5 × 10^6 tế bào/ml.

Thứ tư, mô hình toán học trên MATLAB đã mô phỏng thành công đường cong sinh trưởng bốn pha của tảo và diễn biến suy giảm nồng độ amoni, nitrat. Giá trị RMSE tính toán giữa số liệu thực nghiệm và kết quả mô phỏng dao động ở mức rất thấp (dưới 0,15 đối với mật độ tế bào), khẳng định độ tin cậy vượt trội của mô hình.

Thảo luận kết quả

Cơ chế tăng vọt của pH từ 6,2 lên 8,7 và DO từ 0,2 lên 5,4 mg/l trong bể phản ứng được giải thích do vi tảo tiêu thụ mạnh mẽ nguồn CO2 hòa tan và ion HCO3- để phục vụ phản ứng quang hợp, từ đó làm dịch chuyển cân bằng carbonat và giải phóng một lượng lớn oxy phân tử vào nước. Hiện tượng sụt giảm mật độ sau ngày thứ 10 là do hiện tượng tự che tối khi mật độ quang học tăng cao (cường độ sáng xuyên thấu giảm xuống dưới 1200 lux) kết hợp tình trạng cạn kiệt photphat hòa tan.

Kết quả loại bỏ 96% TN tại tỷ lệ N:P = 16:1 trong nghiên cứu này tương đồng với các công bố quốc tế của Choi (2014) và cao hơn mức 76% TN của Liang Wang và cộng sự (2009). Điều này chứng minh tỷ lệ 16:1 hoàn toàn trùng khớp với công thức phân tử tế bào tảo và là điều kiện tối ưu để tránh hiện tượng nghẽn dinh dưỡng cục bộ. Toàn bộ diễn biến động học sinh trưởng và phân hủy cơ chất có thể được trực quan hóa rõ nét thông qua biểu đồ đường thời gian đa trục biểu diễn đồng thời mật độ tế bào và hàm lượng diệp lục, kết hợp bảng ma trận so sánh hiệu suất suy giảm N-NH4+ và P-PO43- theo từng dải bước sóng ánh sáng.

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

Dựa trên kết quả thực nghiệm và mô phỏng, tác giả đề xuất 4 giải pháp công nghệ trọng tâm nhằm ứng dụng hiệu quả hệ thống HRAPs vào thực tế:

  1. Chuẩn hóa tỷ lệ dinh dưỡng đầu vào: Tiến hành phối trộn hoặc tiền xử lý nước thải sinh hoạt để đưa tỷ lệ N:P về ngưỡng tối ưu 16:1, hướng tới mục tiêu duy trì hiệu suất loại bỏ tổng nitơ ổn định trên 90% và tổng photpho trên 75%. Thời gian áp dụng trong vòng 1 đến 3 tháng đầu vận hành, do kỹ sư quản lý công nghệ xử lý nước thải trực tiếp thực hiện.

  2. Kiểm soát cơ học và chiếu sáng bổ sung: Thiết kế cánh khuấy dạng guồng quay duy trì vận tốc dòng chảy ổn định ở mức 20 cm/s để chống lắng cặn bùn sinh học và đảm bảo xáo trộn tuần hoàn; lắp đặt dàn đèn LED phụ trợ công suất 60W cung cấp cường độ sáng 5000 lux (60 µmol/m²/s) với chu kỳ 16 giờ sáng : 8 giờ tối nhằm khắc phục triệt để tình trạng thiếu sáng vào mùa đông. Hạng mục này cần hoàn thành trong vòng 6 tháng bởi đơn vị thi công kỹ thuật.

  3. Ứng dụng mô hình toán học MATLAB trong điều khiển tự động: Tích hợp module mô phỏng thuật toán vào hệ thống cảm biến đo online (pH, DO, nhiệt độ) để dự báo chính xác thời gian lưu thủy lực (HRT) dao động từ 4 đến 9 ngày theo mùa, kiểm soát sai số RMSE dưới 10%. Timeline triển khai trong vòng 12 tháng do các nhóm nghiên cứu công nghệ môi trường phối hợp cùng doanh nghiệp phần mềm.

  4. Tận thu và chế biến sinh khối tảo: Thiết lập quy trình thu hoạch sinh khối bằng máy ly tâm liên tục ở tốc độ 10.000 vòng/phút, thu hồi nguồn phụ phẩm hữu cơ chứa 40% đến 60% protein làm thức ăn thủy sản và chiết xuất 10% đến 15% chất béo phục vụ tổng hợp biodiesel. Hoạt động này cần thực hiện định kỳ hàng tuần bởi các hợp tác xã nông nghiệp công nghệ cao.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng cụ thể:

  • Học viên cao học và nghiên cứu sinh ngành Khoa học Môi trường, Kỹ thuật Môi trường: Tiếp cận phương pháp luận mẫu mực về ứng dụng ngôn ngữ MATLAB để giải hệ phương trình vi phân mô phỏng sinh thái học vi tảo, kế thừa bộ số liệu thực nghiệm về hằng số động học sinh trưởng của Chlorella vulgaris.
  • Kỹ sư thiết kế và vận hành hệ thống xử lý nước thải: Sử dụng trực tiếp các thông số kỹ thuật đã được tối ưu hóa (độ sâu mương 30 cm, tốc độ khuấy 20 cm/s, chu kỳ sáng tối 16:8) để thiết kế các trạm xử lý nước thải sinh hoạt phân tán quy mô vừa và nhỏ.
  • Doanh nghiệp năng lượng sinh học và sản xuất thức ăn chăn nuôi: Nắm bắt quy trình nuôi cấy vi tảo trong nước thải giàu dinh dưỡng để tận thu sinh khối giàu lipid (10–15%) sản xuất nhiên liệu sinh học biodiesel gốc methyl ester (C16:0 và C18:2) và bột tảo giàu protein (40–60%).
  • Cơ quan quản lý nhà nước về tài nguyên, môi trường và nông nghiệp: Tham khảo mô hình xử lý nước thải sinh học phát thải carbon thấp để xây dựng chính sách kinh tế tuần hoàn, quy chuẩn kỹ thuật cho các khu dân cư nông thôn mới.

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

1. Vi tảo Chlorella vulgaris có bị ức chế sinh trưởng bởi các độc chất trong nước thải bể tự hoại không?
Hoàn toàn không. Kết quả thực nghiệm cho thấy thời gian thích nghi của vi tảo diễn ra dưới 24 giờ và mật độ tế bào đạt đỉnh 7,0 × 10^6 tế bào/ml, cao gấp 2,5 lần so với môi trường dinh dưỡng tổng hợp BBM, chứng minh nước thải bể tự hoại là môi trường dinh dưỡng lý tưởng cho loài tảo này phát triển.

2. Tại sao tỷ lệ N:P = 16:1 lại mang lại hiệu quả xử lý nước thải cao nhất?
Tỷ lệ 16:1 tương ứng chính xác với tỷ lệ cân bằng hóa học nội bào của vi tảo lục (Redfield ratio). Khi đạt tỷ lệ này, tế bào đồng hóa tối đa cả hai nguyên tố mà không gặp hiện tượng giới hạn chất dinh dưỡng, giúp hiệu suất loại bỏ tổng nitơ đạt tới 96% và tổng photpho đạt trên 80%.

3. Vai trò chính của bánh khuấy điều chỉnh tốc độ 20 cm/s trong hệ thống HRAPs là gì?
Bánh khuấy tạo vận tốc dòng chảy 20 cm/s giúp đảo trộn theo phương thẳng đứng, đưa các tế bào tảo từ đáy lên bề mặt để tiếp nhận ánh sáng liên tục, ngăn chặn hiện tượng tự che tối, chống lắng đọng bùn kỵ khí và tăng cường hòa tan oxy khắp tầng nước.

4. Mô hình toán xây dựng trên phần mềm MATLAB có khả năng ứng dụng thực tế ra sao?
Mô hình toán mô phỏng chính xác đường cong sinh trưởng và sự suy giảm dinh dưỡng với sai số RMSE rất nhỏ. Mô hình cho phép dự báo trước nồng độ đầu ra khi các yếu tố thời tiết thay đổi, giúp người vận hành điều chỉnh thời gian lưu nước linh hoạt từ 4 đến 9 ngày.

5. Sinh khối tảo thu gom sau quá trình xử lý có thể sử dụng làm nhiên liệu sinh học không?
Có. Sinh khối Chlorella vulgaris chứa 10% đến 15% hàm lượng lipid tổng số, phân bố chủ yếu ở các mạch acid béo C15 đến C20 (đặc biệt là C16:0 và C18:2). Đây là nguồn nguyên liệu chất lượng cao để thực hiện phản ứng transester hóa tổng hợp nhiên liệu sinh học biodiesel.

Kết luận

  • Luận văn đã chứng minh vi tảo Chlorella vulgaris có khả năng sinh trưởng xuất sắc trong nước thải sau bể tự hoại, đạt mật độ cực đại 7,0 × 10^6 tế bào/ml và nâng chỉ số DO từ 0,2 lên 5,4 mg/l.
  • Xác lập dải tỷ lệ tối ưu N:P = 16:1 giúp nâng hiệu suất xử lý tổng nitơ lên 96% và tổng photpho lên 80%, giải quyết triệt để nguy cơ phú dưỡng nguồn tiếp nhận.
  • Thiết kế hoàn chỉnh hệ thống HRAPs pilot dạng đường đua vận hành ổn định với vận tốc khuấy 20 cm/s và cường độ sáng 5000 lux (chu kỳ 16:8).
  • Ứng dụng thành công phần mềm MATLAB trong việc xây dựng mô hình toán học mô phỏng động thái sinh thái với sai số RMSE tin cậy dưới 0,15.
  • Đóng góp giải pháp công nghệ sinh học tuần hoàn, vừa làm sạch môi trường nước theo chuẩn QCVN 14:2008/BTNMT Cột B, vừa thu hồi sinh khối giàu đạm (40–60%) và dầu béo (10–15%).

Kế hoạch phát triển tiếp theo cần tập trung nâng cấp hệ thống điều khiển tự động hóa trong vòng 12 đến 24 tháng tới và thử nghiệm mở rộng tại các trang trại chăn nuôi quy mô lớn. Các đơn vị nghiên cứu và doanh nghiệp môi trường hãy cùng kết nối để chuyển giao và thương mại hóa công nghệ HRAPs bền vững này ngay hôm nay.