Tổng quan về luận án
Biến đổi khí hậu toàn cầu đang là thách thức an ninh phi truyền thống nghiêm trọng nhất đối với nhân loại thế kỷ 21. Theo các báo cáo đánh giá của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và Thỏa thuận Paris (COP21), việc kiềm chế mức tăng nhiệt độ bề mặt Trái Đất dưới ngưỡng $1,5^\circ\text{C} - 2,0^\circ\text{C}$ so với thời kỳ tiền công nghiệp đòi hỏi mọi quốc gia phải lượng hóa chính xác và cắt giảm triệt để các nguồn phát thải khí nhà kính (KNK). Việt Nam là một trong năm quốc gia chịu tổn thương nặng nề nhất bởi biến đổi khí hậu và nước biển dâng; các kịch bản dự báo chỉ ra rằng nếu mực nước biển dâng 1 m, Việt Nam sẽ đối mặt với mức thiệt hại ước tính lên tới 17 tỷ USD mỗi năm, khiến 1/5 dân số mất nhà ở, 12,3% diện tích đất nông nghiệp bị ngập mặn và khoảng $40.000\text{ km}^2$ đồng bằng sông Cửu Long chịu rủi ro thiên tai khó lường.
Trong bức tranh phát thải quốc gia, lĩnh vực quản lý chất thải đóng vai trò then chốt nhưng lại thiếu hụt nghiêm trọng các công cụ tính toán chi tiết. Theo Thông báo quốc gia lần thứ hai của Việt Nam gửi Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC), tổng lượng phát thải KNK năm 2000 của Việt Nam là 150,9 triệu tấn $\text{CO}_2\text{tđ}$, trong đó nguồn chất thải đóng góp 7,9 triệu tấn $\text{CO}_2\text{tđ}$ (chiếm 5,2%), với tỷ trọng từ xử lý nước thải chiếm khoảng 30% tổng lượng phát thải của ngành chất thải và được dự báo sẽ tăng nhanh theo tiến trình đô thị hóa và công nghiệp hóa. Mặc dù vậy, công tác kiểm kê KNK từ các hệ thống xử lý nước thải (XLNT) tại Việt Nam trong các kỳ kiểm kê trước đây hầu như áp dụng hệ số mặc định cấp độ 1 (Tier 1) của IPCC, vốn có độ không chắc chắn rất lớn và bỏ qua đặc thù động học vi sinh, thành phần cơ chất nhiệt đới cũng như các nguồn phát thải gián tiếp từ tiêu thụ năng lượng.
Xuất phát từ khoảng trống khoa học này, luận án tiến sĩ kỹ thuật môi trường của Nghiên cứu sinh Nguyễn Thị Vân Anh (ngành Kỹ thuật Môi trường, mã số 9520320, Trường Đại học Bách khoa Hà Nội, dưới sự hướng dẫn khoa học của PGS. TS. Đặng Xuân Hiển) với đề tài: "Nghiên cứu ứng dụng mô hình tính toán phát thải khí nhà kính từ một số hệ thống xử lý nước thải ở Việt Nam" đã tạo ra một bước đột phá học thuật mang tính tiên phong.
Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu cốt lõi:
- RQ1: Làm thế nào để thiết lập một mô hình toán học tích hợp động học vi sinh Herbert, cân bằng khối lượng và cân bằng hóa học nhằm lượng hóa đồng thời phát thải KNK trực tiếp ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$) và gián tiếp (từ điện năng) từ các công nghệ XLNT sinh học?
- RQ2: Mức độ tương thích và độ tin cậy của mô hình toán khi kiểm nghiệm trên trạm thực nghiệm chuẩn quốc tế và các hệ thống XLNT quy mô công nghiệp thực tế tại Việt Nam đạt ngưỡng sai số bao nhiêu?
- RQ3: Các hệ số phát thải KNK đặc trưng cục bộ ($\text{kg CO}_2\text{tđ/kg BOD}_5$) cho nước thải sinh hoạt và nước thải công nghiệp sản xuất giấy ở điều kiện nhiệt đới Việt Nam là bao nhiêu khi so sánh giữa kịch bản có và không thu hồi khí sinh học?
Hai giả thuyết nghiên cứu chính được đặt ra và kiểm chứng:
- H1: Mô hình toán học động học bậc cao tích hợp phương trình Herbert và cơ chế chuyển hóa cơ chất - sinh khối có khả năng mô phỏng chính xác phát thải KNK từ hệ thống XLNT với sai số thấp hơn đáng kể so với phương pháp kiểm kê hệ số tĩnh của IPCC Tier 1 và Tier 2.
- H2: Việc thu hồi và đốt/tái sử dụng khí $\text{CH}_4$ từ bể phân hủy bùn kỵ khí kết hợp tối ưu hóa tuổi bùn ($SRT$) và nhiệt độ vận hành có thể giảm thiểu trên 50% tổng lượng phát thải KNK quy đổi tương đương trên toàn trạm xử lý.
Phạm vi nghiên cứu bao quát các hệ thống XLNT sinh hoạt đại diện cho các quy mô đô thị (Tòa nhà Goldmark City, Nhà máy XLNT Yên Sở, Trúc Bạch, Kim Liên) và các hệ thống XLNT công nghiệp sản xuất bột giấy và giấy (Nhà máy Giấy Bãi Bằng, Nhà máy Giấy Corelex) cùng trạm XLNT Windsor (Canada) làm hệ chuẩn đối sánh quốc tế.
Literature Review và Positioning
Khảo cứu lịch sử phát triển của các mô hình đánh giá phát thải KNK từ lĩnh vực nước thải trên thế giới cho thấy sự tiến hóa qua ba dòng lý thuyết chủ đạo:
Dòng lý thuyết 1 tập trung vào các phương pháp kiểm kê vĩ mô do Ủy ban Liên chính phủ về Biến đổi Khí hậu khởi xướng (IPCC 1996, 2000, 2006, 2013). Phương pháp IPCC Tier 1 và Tier 2 tiếp cận thông qua tải lượng hữu cơ phân hủy sinh học tổng số ($TOW$), tiềm năng sinh methane tối đa ($B_0 = 0,6\text{ kg CH}_4\text{/kg BOD}$ hoặc $0,25\text{ kg CH}_4\text{/kg COD}$) và hệ số điều chỉnh methane ($MCF$). Dòng tiếp cận này cung cấp bức tranh bao quát cấp quốc gia nhưng tồn tại điểm mù lý thuyết lớn: giả định hệ thống hoạt động ở trạng thái tĩnh, đồng nhất và bỏ qua bản chất các phản ứng sinh hóa nội tại trong từng bể phản ứng đơn vi phân.
Dòng lý thuyết 2 là các mô hình động học bùn hoạt tính (Activated Sludge Models - ASM1, ASM2, ASM3) do Hiệp hội Nước Quốc tế (IWA) phát triển, cùng các biến thể mở rộng phương trình Monod (1949) và Herbert (1958). Các nghiên cứu của Henze et al. (2000) và Grady et al. (1999) đặt nền tảng cho việc tính toán chuyển hóa carbon, nitơ và sinh khối vi sinh dị dưỡng ($X$) và tự dưỡng ($X_{nit}$), nhưng chủ yếu nhắm vào hiệu quả xử lý nước đầu ra mà chưa liên kết trực tiếp với thông lượng khí nhà kính phát tán vào khí quyển.
Dòng lý thuyết 3 là các mô hình dấu chân carbon tích hợp hiện đại, tiêu biểu là mô hình Bridle Consulting (Bridle, 2008; Monteith et al., 2005; Cicerone, 1989; Czepiel et al., 1995; Kampschreur et al., 2009). Mô hình Bridle đã bước đầu liên kết cân bằng khối lượng với phát thải $\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$ và lượng phát thải gián tiếp do tiêu thụ điện năng từ lưới.
[TIẾP CẬN KIỂM KÊ VĨ MÔ] [ĐỘNG HỌC XỬ LÝ NƯỚC THẢI]
IPCC Guidelines (1996, 2006) IWA Activated Sludge Models (ASM1-3)
- Tiếp cận tĩnh cấp quốc gia - Động học Monod / Herbert vi sinh
- Hệ số Bo, MCF mặc định - Tối ưu hóa chất lượng dòng ra
- Bỏ qua cơ chế sinh hóa nội tại - Chưa lượng hóa luồng phát thải KNK
\ /
\ /
v v
[MÔ HÌNH DẤU CHÂN CARBON TÍCH HỢP]
Bridle Consulting (2008), Monteith et al.
- Kết nối cân bằng vật chất & điện năng
- Giới hạn: Dữ liệu ôn đới, thiếu linh hoạt
|
v
+---------------------------------------+
| POSITIONING CỦA LUẬN ÁN TIẾN SĨ |
| GS. Nguyễn Thị Vân Anh & Đặng Xuân Hiển|
+---------------------------------------+
| • Động học Herbert tích hợp Stoichiometry |
| • Mô hình hóa 3 khí: CO2, CH4, N2O |
| • Khớp nối trực tiếp & gián tiếp điện |
| • Hiệu chuẩn đặc thù nhiệt đới VN |
+---------------------------------------+
Trong bức tranh học thuật này, các cuộc tranh luận khoa học diễn ra gay gắt xung quanh hai luồng quan điểm đối lập:
- Tranh luận về kiểm kê khí $\text{CO}_2$ sinh học (Biogenic $\text{CO}_2$): Hướng dẫn IPCC mặc định loại bỏ $\text{CO}_2$ sinh ra từ quá trình xử lý sinh học vì coi đây là nguồn carbon trung tính chu chuyển ngắn hạn trong tự nhiên. Ngược lại, trường phái mô hình hóa công nghệ môi trường (Bridle, 2008; Doorn et al., 2006) lập luận rằng việc tính toán đầy đủ $\text{CO}_2$ từ hô hấp nội bào vi sinh vật và oxy hóa chất hữu cơ là bắt buộc nhằm đánh giá chính xác cân bằng carbon tổng thể, cho phép tối ưu hóa công nghệ sục khí và lựa chọn giữa xử lý hiếu khí hay kỵ khí.
- Tranh luận về cơ chế phát tán $\text{N}_2\text{O}$: Quan điểm truyền thống của IPCC xem phát thải $\text{N}_2\text{O}$ trực tiếp tại trạm là thứ yếu (chỉ áp dụng hệ số phát thải tĩnh $3,2\text{ g N}_2\text{O/người/năm}$). Tuy nhiên, các nghiên cứu chuyên sâu của Kampschreur et al. (2009) và Ahn et al. (2010) chỉ ra rằng $\text{N}_2\text{O}$ có thể bùng phát rất mạnh trong điều kiện nồng độ oxy hòa tan (DO) thấp hoặc tải lượng nitrit tích lũy cao trong quá trình nitrat hóa do vi khuẩn Nitrosomonas và quá trình khử nitrat không hoàn toàn.
So sánh với hai nghiên cứu quốc tế điển hình là mô hình Bridle áp dụng tại nhà máy XLNT Windsor (Canada) và nghiên cứu của IPCC (2006), luận án đã định vị chính xác khoảng trống nghiên cứu: thiết lập một khung mô hình toán học giải tích số thích ứng hoàn hảo với đặc trưng nước thải có tỷ lệ chất rắn lơ lửng và tải lượng hữu cơ biến động mạnh ở các nước đang phát triển, bổ khuyết tính năng mô phỏng động học thời gian thực trong điều kiện khí hậu nhiệt đới tại Việt Nam.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã mở rộng khung lý thuyết động học tăng trưởng vi sinh vật thông qua việc cải biên phương trình Herbert (1958) kết hợp mô hình Monod, giải quyết triệt để mối quan hệ giữa tốc độ sử dụng chất nền ($r_s$), tốc độ sinh trưởng thực của sinh khối ($r_x$) và tốc độ phân hủy nội bào sinh khối dị dưỡng ($k_d$).
Hệ thống phương trình động học cốt lõi được xây dựng trên nền tảng:
$$r_s = \frac{k \cdot S \cdot X}{K_s + S}$$
$$r_x = Y \cdot r_s - k_d \cdot X$$
Trong đó:
- $k$: Tốc độ sử dụng chất nền riêng cực đại ($\text{mg/mg}\cdot\text{ngày}$).
- $S$: Nồng độ cơ chất trong bể sinh học ($\text{mg/L}$).
- $X$: Nồng độ sinh khối vi sinh vật dị dưỡng ($\text{mg/L}$).
- $K_s$: Hằng số bán bão hòa đối với cơ chất ($\text{mg/L}$).
- $Y$: Hệ số sản lượng sinh khối cực đại ($\text{mg sinh khối/mg cơ chất}$).
- $k_d$: Hệ số phân hủy nội bào ($\text{ngày}^{-1}$).
Luận án đã đóng góp đột phá khi thiết lập các phương trình hóa học và công thức kinh nghiệm đại diện chuẩn xác cho cơ chất và sinh khối trong điều kiện thực nghiệm Việt Nam:
- Nước thải sinh hoạt: Cơ chất được đại diện bởi công thức hóa học kinh nghiệm $\text{C}{10}\text{H}{19}\text{NO}_3$, sinh khối tế bào vi sinh vật đại diện bằng $\text{C}_5\text{H}_7\text{O}_2\text{N}$. Quá trình oxy hóa sinh học hiếu khí hoàn toàn diễn ra theo phương trình stoichiometry:
$$\text{C}{10}\text{H}{19}\text{NO}_3 + 4,5\text{O}_2 \rightarrow \text{C}_5\text{H}_7\text{O}_2\text{N} + 5\text{CO}_2 + 6\text{H}_2\text{O} + \text{NH}_3$$
- Nước thải công nghiệp sản xuất giấy: Cơ chất giàu hợp chất hữu cơ mạch vòng và cellulose được mô hình hóa dưới công thức $\text{C}7\text{H}{12}\text{O}_4\text{N}$, tạo nên bước chuyển hóa sinh hóa đặc thù trong các công trình xử lý yếm khí và hiếu khí kết hợp.
Mô hình lý thuyết đã chuyển dịch mô thức từ việc ước tính phát thải KNK dựa trên các chỉ số ô nhiễm đầu vào một cách tĩnh tại sang mô hình hóa cơ chế cân bằng vật chất phân đoạn: lượng carbon dòng vào được phân bổ chính xác thành 3 dòng chuyển hóa: (1) Carbon tích hợp vào sinh khối bùn mới hình thành ($P_X$), (2) Carbon chuyển hóa thành khí nhà kính trực tiếp ($\text{CO}_2, \text{CH}_4$), và (3) Carbon dư thừa thoát ra theo dòng thải sau xử lý ($S_r$).
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp toàn diện ba cấu phần độc lập nhưng tương hỗ mật thiết:
+---------------------------------------------------------------------------------------------------+
| TỔNG PHÁT THẢI KHÍ NHÀ KÍNH |
| (M_Total, kg CO2td/ngày) |
+-------------------------------------------------+-------------------------------------------------+
|
+--------------------------------+-------------------------------+
| |
v v
+---------------------------------+ +---------------------------------+
| PHÁT THẢI TRỰC TIẾP | | PHÁT THẢI GIÁN TIẾP |
| (M_KNK,qtxl) | | (M_CO2,điện) |
+---------------------------------+ +---------------------------------+
| • Bể hiếu khí (M_CO2,bexl, | | • Điện năng tiêu thụ: |
| M_N2O,bexl) | | M_CO2,điện = Q_E * |
| • Bể kỵ khí & Phân hủy bùn | | sum(PF_i * EF_i) |
| (M_CH4,bexlyk, M_CO2,bephanhuy)| | • Cơ cấu nguồn điện lưới VN: |
| • Dòng nước thải đầu ra | | Thủy điện, Than, Khí, Dầu |
| (M_CO2,dongra, M_CH4,dongra) | | |
+---------------------------------+ +---------------------------------+
- Cấu phần 1: Khối phát thải trực tiếp từ quá trình xử lý sinh học ($M_{\text{KNK,qtxl}}$): Tính toán chi tiết lượng $\text{CO}_2$ sinh ra từ bể hiếu khí dựa trên động học Monod-Herbert; lượng $\text{CH}_4$ và $\text{CO}_2$ từ bể phân hủy bùn yếm khí; lượng $\text{N}_2\text{O}$ sinh ra từ hai con đường nitrat hóa (oxy hóa amoni thành nitrit bởi vi khuẩn Nitrosomonas và oxy hóa nitrit thành nitrat bởi Nitrobacter) và khử nitrat dị hóa.
- Cấu phần 2: Khối phát thải gián tiếp từ tiêu thụ năng lượng điện ($M_{\text{CO}_2\text{,điện}}$):
$$M_{\text{CO}_2\text{,điện}} = Q_E \cdot \sum \left( PF_i \cdot EF_i \right)$$
Trong đó:
- $Q_E$: Tổng công suất tiêu thụ điện năng của trạm ($\text{kWh/ngày}$), phục vụ sục khí, bơm hồi lưu, khuấy trộn và xử lý bùn.
- $PF_i$: Tỷ lệ phần trăm đóng góp của nguồn phát điện thứ $i$ trong cơ cấu lưới điện quốc gia Việt Nam (nhiệt điện than, thủy điện, tuabin khí, dầu).
- $EF_i$: Hệ số phát thải KNK tương ứng của từng loại hình phát điện ($\text{kg CO}_2\text{tđ/kWh}$).
- Cấu phần 3: Khối phân rã bùn và BOD dòng ra: Lượng hóa lượng KNK phát sinh tiềm tàng khi nước thải sau xử lý xả ra các thủy vực tiếp nhận tự nhiên.
Điều kiện biên (boundary conditions) của mô hình được xác lập chặt chẽ: bao quát từ cửa nhận nước thải đầu vào trạm xử lý, cụm công trình xử lý cơ học sơ cấp, sinh học thứ cấp, lắng, xử lý và khử nước bùn hoạt tính, cho đến điểm xả thải cuối cùng ra môi trường tiếp nhận.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án theo đuổi lập trường nhận thức luận thực chứng (Positivism) kết hợp chặt chẽ giữa mô hình hóa toán giải tích và đo đạc thực nghiệm tại hiện trường. Thiết kế nghiên cứu đa tầng bậc bao gồm:
- Thiết lập hệ thống phương trình vi phân mô tả trạng thái biến đổi nồng độ cơ chất ($S$), sinh khối ($X$), bùn trơ không phân hủy sinh học ($X_{nb}$), sinh khối vi khuẩn nitrat hóa ($X_{nit}$), lượng bùn xả ($Q_x$), bùn tuần hoàn ($Q_T$) và thời gian lưu thủy lực ($HRT$), tuổi bùn ($SRT$).
- Thiết kế mô phỏng số trên nền tảng ngôn ngữ lập trình Scala và MATLAB, xây dựng phần mềm giao diện đồ họa trực quan (GUI) cho phép nhập dữ liệu công nghệ và xuất kết quả phát thải theo thời gian thực.
- Đo đạc thực địa và lấy mẫu thực nghiệm để kiểm chứng mô hình trên các hệ thống xử lý nước thải đại diện:
- Hệ thống A: Xử lý hiếu khí truyền thống kết hợp phân hủy bùn.
- Hệ thống B: Xử lý yếm khí.
- Hệ thống C: Xử lý kết hợp kỵ khí - thiếu khí - hiếu khí (A2O).
Quy trình nghiên cứu rigorous
Quy trình thu thập và phân tích dữ liệu thực nghiệm tuân thủ nghiêm ngặt các tiêu chuẩn quốc tế:
- Phương pháp lấy mẫu khí bề mặt: Sử dụng thiết bị buồng thu khí nổi chuyên dụng Chamber QUT (do Đại học Công nghệ Queensland - Brisbane, Australia thiết kế và chuyển giao cho Viện Môi trường Nông nghiệp Việt Nam). Buồng thu được thả nổi trên mặt thoáng các bể xử lý sinh học (bể hiếu khí, bể phân hủy bùn, bể lắng) với quy trình kiểm soát áp suất và nhiệt độ buồng kín nghiêm ngặt để thu mẫu khí theo chuỗi thời gian ($t_0, t_5, t_{10}, t_{15}, t_{20}\text{ phút}$).
- Phương pháp phân tích công cụ hiện đại: Các mẫu khí thu thập được phân tích trên hệ thống sắc ký khí cao cấp Shimadzu GC-2014 (Nhật Bản) trang bị đồng thời 3 đầu dò chuyên dụng:
- Khí $\text{CO}_2$: Xác định bằng Đầu dò dẫn nhiệt (TCD) vận hành ở nhiệt độ $200^\circ\text{C}$.
- Khí $\text{CH}_4$: Xác định bằng Đầu dò ion hóa ngọn lửa (FID) vận hành ở nhiệt độ $300^\circ\text{C}$.
- Khí $\text{N}_2\text{O}$: Xác định bằng Đầu dò bắt điện tử ($\text{ECD - } \mu\text{ECD}$) vận hành ở nhiệt độ cực cao $350^\circ\text{C}$.
| Thông số thiết bị / Quy trình |
Quy chuẩn kỹ thuật áp dụng trong luận án |
| Thiết bị sắc ký khí |
Shimadzu GC-2014 (Shimadzu Corporation, Japan) |
| Định lượng $\text{CO}_2$ |
Đầu dò dẫn nhiệt (TCD), nhiệt độ buồng dò: $200^\circ\text{C}$ |
| Định lượng $\text{CH}_4$ |
Đầu dò ion hóa ngọn lửa (FID), nhiệt độ buồng dò: $300^\circ\text{C}$ |
| Định lượng $\text{N}_2\text{O}$ |
Đầu dò bắt giữ điện tử (ECD), nhiệt độ buồng dò: $350^\circ\text{C}$ |
| Kỹ thuật thu mẫu mặt thoáng |
Floating Chamber QUT (Queensland University of Technology, Australia) |
| Ngôn ngữ giải thuật số |
Scala & MATLAB Numerical Computation Environment |
Quy trình bảo đảm chất lượng (QA/QC) được thực hiện qua các đường chuẩn hồi quy tuyến tính với $R^2 > 0,999$, lặp mẫu 3 lần cho mỗi điểm đo nhằm triệt tiêu sai số ngẫu nhiên.
Data và phân tích
Tập dữ liệu thực nghiệm và dữ liệu công nghệ được xử lý và phân tích thống kê sâu sắc:
- Dữ liệu vận hành hệ thống Goldmark City: Lưu lượng trung bình $2.700\text{ m}^3\text{/ngày}$, nồng độ $\text{BOD}_5$ đầu vào $220\text{ mg/L}$, $\text{COD} = 380\text{ mg/L}$, $\text{SS} = 200\text{ mg/L}$, $\text{TN} = 45\text{ mg/L}$.
- Dữ liệu nhà máy Giấy Bãi Bằng: Lưu lượng cực lớn $30.000\text{ m}^3\text{/ngày}$, $\text{COD}$ đầu vào $1.200\text{ mg/L}$, $\text{BOD}_5 = 650\text{ mg/L}$.
- Phân tích độ nhạy (Sensitivity Analysis): Luận án thực hiện khảo sát độ nhạy đa biến của mô hình đối với các thông số: nhiệt độ nước thải ($T = 15^\circ\text{C} - 35^\circ\text{C}$), tuổi bùn ($SRT = 5 - 30\text{ ngày}$), và nồng độ cơ chất đầu vào ($S_0$). Kết quả kiểm tra độ ổn định số học cho thấy mô hình đạt độ hội tụ vững chắc và phản ánh đúng bản chất lý sinh của các quá trình chuyển hóa vi sinh.
Phát hiện đột phá và implications
Những phát hiện then chốt
Luận án mang lại 5 phát hiện mang tính bước ngoặt cho ngành kỹ thuật môi trường và kiểm kê khí nhà kính tại Việt Nam:
Thứ nhất, lượng hóa và công bố chính thức bộ hệ số phát thải KNK đặc trưng quốc gia:
- Đối với hệ thống XLNT sinh hoạt:
- Kịch bản 1 (Không thu hồi khí $\text{CH}_4$, phóng không): Hệ số phát thải là $2,40\text{ kg CO}_2\text{tđ/kg BOD}_5$.
- Kịch bản 2 (Thu hồi và đốt/tận dụng năng lượng $\text{CH}_4$): Hệ số phát thải giảm mạnh xuống còn $1,16\text{ kg CO}_2\text{tđ/kg BOD}_5$.
- Đối với hệ thống XLNT công nghiệp sản xuất giấy:
- Kịch bản 1 (Không thu hồi $\text{CH}_4$): Hệ số phát thải đạt $3,58\text{ kg CO}_2\text{tđ/kg BOD}_5$.
- Kịch bản 2 (Thu hồi và đốt $\text{CH}_4$): Hệ số phát thải chỉ còn $0,93\text{ kg CO}_2\text{tđ/kg BOD}_5$.
HỆ SỐ PHÁT THẢI KNK ĐẶC TRƯNG QUỐC GIA (kg CO2tđ/kg BOD5)
===================================================================
Nước thải sinh hoạt:
• Không thu hồi CH4 (Phóng không) : [ 2.40 ] kg CO2tđ/kg BOD5
• Có thu hồi & đốt CH4 : [ 1.16 ] kg CO2tđ/kg BOD5 (-51.7%)
Nước thải công nghiệp giấy:
• Không thu hồi CH4 (Phóng không) : [ 3.58 ] kg CO2tđ/kg BOD5
• Có thu hồi & đốt CH4 : [ 0.93 ] kg CO2tđ/kg BOD5 (-74.0%)
===================================================================
Thứ hai, phát hiện mối tương quan phi tuyến tính mang tính nghịch lý giữa nhiệt độ và phát thải KNK:
Khi nhiệt độ nước thải tăng từ $20^\circ\text{C}$ lên $30^\circ\text{C}$, tốc độ phân hủy cơ chất tăng nhanh làm giảm $\text{BOD}_5$ dòng ra (hiệu quả xử lý nước tăng), nhưng đồng thời thúc đẩy mãnh liệt tốc độ hô hấp nội bào của vi sinh vật ($k_d$ tăng từ $0,05\text{ ngày}^{-1}$ lên $0,09\text{ ngày}^{-1}$), dẫn đến lượng phát thải $\text{CO}_2$ trực tiếp tại bể hiếu khí tăng hơn 28,5% và lượng $\text{CH}_4$ sinh ra tại bể phân hủy bùn tăng 34,2%.
Thứ ba, xác lập vai trò chi phối của phát thải gián tiếp từ điện năng:
Trong các hệ thống XLNT sinh hoạt đô thị hiếu khí như Goldmark City và Yên Sở, phát thải gián tiếp từ tiêu thụ điện năng cho máy thổi khí và máy bơm chiếm từ 35% đến 48% tổng lượng phát thải KNK toàn trạm. Điều này chứng minh rằng việc áp dụng công nghệ sục khí hiệu suất thấp sẽ triệt tiêu lợi ích giảm phát thải sinh học của hệ thống.
Thứ tư, xác thực độ tin cậy vượt trội của mô hình khi kiểm nghiệm thực chứng:
Khi chạy kiểm chứng với dữ liệu thực nghiệm tại trạm Goldmark City và trạm chuẩn quốc tế Windsor, sai số tương đối giữa giá trị mô hình tính toán và số liệu phân tích thực nghiệm trên sắc ký khí Shimadzu GC-2014 dao động dưới ngưỡng 8,5%, khẳng định mô hình có tính tương thích và độ vững rất cao.
Implications đa chiều
- Đóng góp lý thuyết: Luận án đã tích hợp thành công động học Herbert vào bài toán mô hình hóa phát thải KNK, lấp đầy khoảng trống giữa lý thuyết công nghệ môi trường cổ điển và khoa học biến đổi khí hậu hiện đại.
- Đổi mới phương pháp luận: Cung cấp một phương pháp luận nhất quán giải quyết bài toán phát thải KNK cho các quốc gia đang phát triển thuộc khu vực Đông Nam Á, nơi có đặc tính nước thải và điều kiện khí tượng tương đồng.
- Ứng dụng thực tiễn: Bộ phần mềm lập trình trên Scala và MATLAB là công cụ hỗ trợ đắc lực cho các kỹ sư tư vấn thiết kế trạm XLNT nhằm đánh giá vòng đời carbon (LCA), tối ưu hóa chế độ cấp khí, lựa chọn tuổi bùn ($SRT$) và thời gian lưu thủy lực ($HRT$) hợp lý nhằm đồng thời đạt hai mục tiêu: nước thải đạt chuẩn QCVN và phát thải KNK thấp nhất.
- Hàm ý chính sách quản lý: Kết quả nghiên cứu là cơ sở khoa học then chốt để Bộ Tài nguyên và Môi trường cập nhật Thông báo Quốc gia (National Communications), Báo cáo Cập nhật Hai năm một lần (BUR), và Đóng góp do Quốc gia Tự quyết định (NDC) gửi UNFCCC, tạo cơ sở cho việc phân bổ hạn ngạch phát thải và phát triển thị trường tín chỉ carbon trong ngành cấp thoát nước.
Limitations và Future Research
Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra 3 giới hạn nghiên cứu:
- Giả định thành phần cơ chất kinh nghiệm: Việc mô hình hóa cơ chất nước thải sinh hoạt dưới dạng đơn phân tử $\text{C}{10}\text{H}{19}\text{NO}_3$ và nước thải giấy dưới dạng $\text{C}7\text{H}{12}\text{O}_4\text{N}$ tuy bảo đảm tính khả thi cho việc giải hệ phương trình cân bằng nhưng chưa phản ánh trọn vẹn sự đa dạng của các hợp chất hữu cơ khó phân hủy, các chất hoạt động bề mặt và các chất vi ô nhiễm hữu cơ mới phát sinh (micropollutants).
- Ranh giới phát thải gián tiếp: Nghiên cứu mới tính toán phát thải gián tiếp từ điện năng tiêu thụ, chưa bao quát phát thải gián tiếp từ quá trình sản xuất và vận chuyển hóa chất keo tụ (như $\text{FeCl}_3\cdot 6\text{H}_2\text{O}$, PAC, Polymer) và công đoạn vận chuyển bùn thải ra ngoài bãi chôn lấp tập trung.
- Độ bao phủ dữ liệu thực địa: Dữ liệu đo đạc bằng sắc ký khí GC chủ yếu tập trung tại khu vực miền Bắc Việt Nam, cần mở rộng mạng lưới quan trắc tại các vùng khí hậu Tây Nam Bộ và Nam Trung Bộ để tăng cường tính đại diện không gian.
Các hướng nghiên cứu phát triển tiếp theo bao gồm:
- Mở rộng ứng dụng mô hình cho các loại hình nước thải công nghiệp có tải lượng hữu cơ và độ mặn cao như: chế biến thủy hải sản, thuộc da, dệt nhuộm và nước thải chăn nuôi lợn quy mô công nghiệp.
- Tích hợp mô hình động học phát thải KNK với mô hình động lực học chất lưu tính toán (Computational Fluid Dynamics - CFD) để mô phỏng sự phân bố không gian của các túi khí kỵ khí cục bộ trong bể hiếu khí.
- Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) kết hợp mạng lưới cảm biến IoT đo trực tuyến nồng độ $\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$ để điều khiển tự động hệ thống cấp khí thông minh, hướng tới các nhà máy XLNT Net-Zero Energy và Net-Zero Carbon.
Tác động và ảnh hưởng
Công trình luận án mang lại những tác động sâu rộng và đa tầng:
- Ảnh hưởng học thuật: Mở ra hướng tiếp cận học thuật mới trong các trường đại học khối kỹ thuật tại Việt Nam, thiết lập tiêu chuẩn mới cho các nghiên cứu đánh giá dấu chân carbon trong công nghệ môi trường; tạo lập nguồn trích dẫn học thuật quan trọng cho các công bố quốc tế về lĩnh vực kiểm kê khí nhà kính đô thị nhiệt đới.
- Chuyển đổi ngành công nghiệp: Hướng dẫn các doanh nghiệp sản xuất giấy lớn (như Bãi Bằng, Corelex) và các tập đoàn bất động sản quản lý đô thị chuyển đổi công nghệ xử lý bùn sang thu hồi khí sinh học $\text{CH}_4$ phát điện, giúp giảm hàng ngàn tấn $\text{CO}_2\text{tđ}$ phát thải mỗi năm và tiết kiệm chi phí năng lượng đáng kể.
- Đóng góp xã hội và quốc tế: Hỗ trợ Chính phủ Việt Nam hiện thực hóa cam kết đạt mức phát thải ròng bằng "0" (Net Zero) vào năm 2050 được tuyên bố tại COP26, nâng cao tính minh bạch và độ chính xác của cơ sở dữ liệu phát thải quốc gia trước cộng đồng quốc tế.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận mô hình toán học tích hợp hoàn chỉnh, kế thừa mã nguồn phần mềm trên Scala/MATLAB và phương pháp luận đo đạc buồng Chamber kết hợp sắc ký khí Shimadzu GC-2014.
- Kỹ sư thiết kế và Các đơn vị vận hành trạm XLNT: Sở hữu công cụ số trực quan để dự báo chính xác tải lượng KNK, từ đó lựa chọn phương án công nghệ tối ưu hóa năng lượng và kinh tế.
- Cơ quan quản lý Nhà nước (Bộ TN&MT, Cục Biến đổi Khí hậu, Sở TN&MT): Có được bộ hệ số phát thải chuẩn xác cấp quốc gia phục vụ công tác xây dựng thể chế, ban hành quy chuẩn kỹ thuật và định giá hạn ngạch phát thải carbon trong ngành hạ tầng kỹ thuật đô thị.
Câu hỏi chuyên sâu
-
Đóng góp lý thuyết độc đáo nhất của luận án là gì?
Đó là việc mở rộng thành công lý thuyết động học Herbert kết hợp stoichiometry thực nghiệm để xây dựng hệ phương trình liên tục mô tả sự biến thiên đồng thời của cơ chất, sinh khối và các luồng phát thải KNK ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$), khắc phục hoàn toàn tính tĩnh tại và độ không chắc chắn của mô hình IPCC truyền thống.
-
Đổi mới phương pháp luận nổi bật khi so sánh với các nghiên cứu trước đây?
Khác với mô hình Bridle (2008) chỉ áp dụng cho khí hậu ôn đới hoặc phương pháp IPCC chỉ tính toán vĩ mô, luận án đã tích hợp đo đạc thực nghiệm bằng kỹ thuật buồng Chamber QUT tiên tiến và phân tích đa đầu dò sắc ký khí GC (TCD, FID, ECD) với mô hình số hóa trên Scala/MATLAB, đạt mức tương thích cao khi hiệu chuẩn trên hệ thống thực tế tại Việt Nam.
-
Phát hiện gây bất ngờ nhất được dữ liệu thực nghiệm chứng minh là gì?
Phát hiện về tác động vượt trội của việc thu hồi khí $\text{CH}_4$: thu hồi và đốt khí sinh học tại bể bùn kỵ khí giúp giảm tới 51,7% tổng phát thải KNK ở nước thải sinh hoạt và giảm tới 74,0% ở nước thải nhà máy sản xuất giấy. Đồng thời, lượng phát thải gián tiếp do điện năng tiêu thụ chiếm tỷ trọng cực lớn (lên tới 48% tổng phát thải ở các trạm xử lý hiếu khí đô thị).
-
Luận án có cung cấp quy trình tái lặp thực nghiệm (Replication Protocol) rõ ràng không?
Có. Luận án mô tả chi tiết từng bước quy trình chế tạo buồng thu khí Chamber QUT, các thông số cài đặt nhiệt độ đầu dò sắc ký khí Shimadzu GC-2014 (TCD: $200^\circ\text{C}$, FID: $300^\circ\text{C}$, ECD: $350^\circ\text{C}$), cùng hệ thống phương trình toán và thuật giải số, cho phép các nhà khoa học khác tái lặp hoàn toàn trên các trạm XLNT khác.
-
Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Hướng tới việc tích hợp mô hình động học Herbert với công nghệ mô phỏng thủy lực 3D (CFD), mở rộng hệ số phát thải cho toàn bộ 10 ngành công nghiệp trọng điểm phát sinh nước thải lớn tại Việt Nam, và xây dựng nền tảng giám sát KNK trực tuyến thời gian thực dựa trên cảm biến quang học và trí tuệ nhân tạo (AI-enabled Digital Twin).
Kết luận
Luận án tiến sĩ của Nghiên cứu sinh Nguyễn Thị Vân Anh đã xuất sắc hoàn thành trọn vẹn mục tiêu đề ra với 5 đóng góp học thuật và thực tiễn cốt lõi:
- Thiết lập thành công mô hình toán học giải tích số đầu tiên tại Việt Nam dựa trên phương trình động học Herbert và các cân bằng vật chất - hóa học để lượng hóa phát thải KNK trực tiếp và gián tiếp từ hệ thống XLNT.
- Xây dựng bộ mã nguồn và phần mềm tính toán có giao diện trực quan, giải quyết hệ phương trình vi phân phi tuyến nhanh chóng và chuẩn xác.
- Xác lập và công bố bộ hệ số phát thải KNK đặc trưng quốc gia cho nước thải sinh hoạt ($2,40\text{ kg CO}_2\text{tđ/kg BOD}_5$ khi không thu hồi $\text{CH}_4$; $1,16\text{ kg CO}_2\text{tđ/kg BOD}_5$ khi thu hồi $\text{CH}_4$) và nước thải sản xuất giấy ($3,58\text{ kg CO}_2\text{tđ/kg BOD}_5$ khi không thu hồi $\text{CH}_4$; $0,93\text{ kg CO}_2\text{tđ/kg BOD}_5$ khi thu hồi $\text{CH}_4$).
- Kiểm nghiệm thực chứng thành công trên hệ thống chuẩn quốc tế Windsor và hệ thống thực tế Goldmark City, đạt độ chính xác cao với sai số dưới 8,5%.
- Cung cấp luận cứ khoa học vững chắc phục vụ đắc lực cho công tác kiểm kê khí nhà kính quốc gia, xây dựng Thông báo Quốc gia cho UNFCCC và định hình các chính sách phát triển kinh tế tuần hoàn, carbon thấp trong tương lai.
