Luận án tính toán phát thải khí nhà kính từ xử lý nước thải Việt Nam

Tổng quan về luận án

Biến đổi khí hậu (BĐKH) toàn cầu do sự gia tăng nồng độ các chất khí nhà kính (KNK) trong khí quyển là thách thức sinh thái và kinh tế - xã hội nghiêm trọng bậc nhất của thế kỷ 21. Theo đánh giá của Liên Hợp Quốc, Việt Nam là một trong năm quốc gia chịu tác động nặng nề nhất từ BĐKH và nước biển dâng; nếu mực nước biển dâng 1 m, ước tính 40.000 km² diện tích đồng bằng bị ngập lụt, 12,3% đất nông nghiệp biến mất và tổn thất kinh tế lên tới 17 tỷ USD mỗi năm. Trong cơ cấu phát thải quốc gia, lĩnh vực chất thải đóng góp từ 5,2% đến 6,0% tổng lượng phát thải KNK, trong đó phát thải từ hoạt động xử lý và xả thải nước thải chiếm xấp xỉ 30% tổng nguồn chất thải (theo Thông báo Quốc gia lần thứ hai gửi UNFCCC). Mặc dù các nỗ lực kiểm kê KNK quốc gia đã áp dụng các hướng dẫn của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC 1996, 2000, 2006), việc xác định lượng phát thải từ các hệ thống xử lý nước thải (HTXLNT) tại Việt Nam phần lớn vẫn dừng lại ở phương pháp Cấp 1 (Tier 1) với các hệ số phát thải (EF) mặc định mang tính khái quát cao, bỏ qua các đặc thù công nghệ, động học sinh học và biến thiên khí hậu nhiệt đới.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ việc thiếu hụt một công cụ toán học tất định (deterministic model) tích hợp toàn diện cả nguồn phát thải gián tiếp (tiêu thụ điện năng) lẫn nguồn phát thải trực tiếp (CO2, CH4, N2O) sinh ra từ các phản ứng hóa sinh trong từng công trình đơn vị. Các mô hình quốc tế như mô hình Bridle Consulting (Bridle, 2008) hay các hướng dẫn IPCC Cấp 1 và Cấp 2 thường giả định hệ thống ở trạng thái ổn định tĩnh hoặc sử dụng hệ số phát sinh lý thuyết không phản ánh được cơ chế chuyển hóa chất nền và sinh khối vi sinh vật tại Việt Nam. Để giải quyết khoảng trống này, luận án tiến sĩ kỹ thuật môi trường "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" của nghiên cứu sinh Nguyễn Thị Vân Anh (người hướng dẫn khoa học: PGS. Đặng Xuân Hiển, Trường Đại học Bách khoa Hà Nội, 2019) đã tiên phong thiết lập khung mô hình toán học giải tích động học và code số hóa chuyên dụng.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học:

• Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một mô hình toán học tích hợp giữa phương trình động học vi sinh Herbert và các cân bằng vật chất - hóa học nhằm mô phỏng chính xác lượng phát thải KNK trực tiếp và gián tiếp từ các HTXLNT sinh hoạt và công nghiệp?
• Câu hỏi nghiên cứu 2 (RQ2): Mức độ ảnh hưởng của các thông số công nghệ (tuổi bùn SRT, thời gian lưu nước HRT, nồng độ cơ chất đầu vào $S_{o,v}$) và điều kiện vận hành nhiệt độ nhiệt đới ($15^\circ\text{C} - 35^\circ\text{C}$) đến cường độ phát thải KNK là bao nhiêu?
• Câu hỏi nghiên cứu 3 (RQ3): Hệ số phát thải KNK đặc trưng ($EF$) cho nước thải sinh hoạt đô thị và nước thải công nghiệp sản xuất bột giấy - giấy tại Việt Nam khi có và không có thu hồi khí mê-tan ($CH_4$) đạt giá trị định lượng cụ thể nào?
• Giả thuyết khoa học 1 (H1): Việc tích hợp mô hình động học Herbert kết hợp phương trình cân bằng khối lượng và công thức kinh nghiệm cơ chất ($C_{10}H_{19}NO_3$ cho nước thải sinh hoạt; $C_7H_{12}O_4N$ cho nước thải giấy) cho phép dự báo phát thải KNK với sai số dưới 15% so với đo đạc thực nghiệm.
• Giả thuyết khoa học 2 (H2): Nguồn phát thải KNK từ tiêu thụ điện năng (gián tiếp) và phân hủy bùn yếm khí (trực tiếp) chiếm tỷ trọng chi phối trên 80% tổng phát thải của toàn trạm XLNT.
• Giả thuyết khoa học 3 (H3): Việc áp dụng giải pháp thu hồi và đốt triệt để khí $CH_4$ sinh ra từ bể phân hủy bùn kỵ khí sẽ cắt giảm trên 50% tổng lượng phát thải KNK quy đổi tương đương ($CO_{2\text{tđ}}$) của toàn hệ thống.

Phạm vi nghiên cứu bao quát các hệ thống xử lý nước thải sinh hoạt tiêu biểu (Khu đô thị Goldmark City - 136 Hồ Tùng Mậu với công suất $4.200\text{ m}^3/\text{ngày}$, Yên Sở, Kim Liên, Trúc Bạch) và hệ thống xử lý nước thải công nghiệp giấy quy mô lớn (Nhà máy Giấy Bãi Bằng công suất $40.000\text{ m}^3/\text{ngày}$, Nhà máy Giấy Corelex công suất $1.000\text{ m}^3/\text{ngày}$) trong giai đoạn 2016-2019. Luận án mang ý nghĩa đột phá khi lần đầu tiên lượng hóa thành công hệ số phát thải cấp quốc gia cho Việt Nam, phục vụ trực tiếp cho cam kết NDC và xây dựng Thông báo Quốc gia định kỳ gửi Ban thư ký UNFCCC.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy các luồng nghiên cứu về phát thải KNK từ nước thải phân nhánh thành ba trường phái học thuật chính:

1. Trường phái kiểm kê thống kê vĩ mô (IPCC Tier 1/2/3): Khởi xướng từ hướng dẫn IPCC (1996, 2006) và Bổ sung Đất ngập nước (IPCC, 2013), dựa trên các tham số tải lượng hữu cơ phân hủy ($TOW$, tính theo $BOD_5$ hoặc $COD$), hệ số hiệu chỉnh methane ($MCF$) và năng lực tạo methane tối đa ($B_o = 0,6\text{ kg } CH_4/\text{kg } BOD$ hoặc $0,25\text{ kg } CH_4/\text{kg } COD$). Tuy nhiên, phương pháp này chỉ xem xét cân bằng tĩnh tại biên giới vĩ mô, bỏ qua hoàn toàn cơ chế phản ứng sinh hóa nội tại trong từng bể sinh học.
2. Trường phái mô hình hóa cân bằng vật chất công nghệ (Process-based models): Tiêu biểu là mô hình Bridle Consulting (Bridle, 2008), Monteith et al. (2005) và Shahabadi et al. (2009). Các mô hình này phân rã trạm xử lý thành các khối xử lý cơ học, bùn hoạt tính hiếu khí và ổn định bùn kỵ khí để tính toán phát thải $CO_2$ hóa thạch, $CO_2$ sinh học (biogenic), $CH_4$ và điện năng tiêu thụ.
3. Trường phái mô hình bùn hoạt tính nâng cao (Activated Sludge Models - ASM series): Được phát triển bởi Hiệp hội Nước Quốc tế (IWA) bao gồm ASM1, ASM2d, ASM3 (Henze et al., 2000) và mở rộng cho phát thải $N_2O$ (Hiatt & Grady, 2008; Ni et al., 2013; Kampschreur et al., 2009). Nhóm này mô tả chi tiết con đường oxy hóa amoni thành nitrit bởi vi khuẩn Nitrosomonas, oxy hóa nitrit thành nitrat bởi vi khuẩn Nitrobacter, và quá trình khử nitrat qua các enzym hydroylamine oxidoreductase, nitrate reductase, nitrite reductase, nitric oxide reductase để hình thành $N_2O$.

                    ┌────────────────────────────────────────────────────────┐
                    │               KHUNG TIẾP CẬN TỔNG QUAN Y VĂN           │
                    └───────────────────────────┬────────────────────────────┘
                                                │
         ┌──────────────────────────────────────┼──────────────────────────────────────┐
         ▼                                      ▼                                      ▼
┌────────────────────────────────┐ ┌─────────────────────────────────┐ ┌────────────────────────────────┐
│   Trường phái IPCC (Tier 1-3)  │ │   Mô hình Bridle & Cân bằng Kỹ thuật │ │     Mô hình Động học ASM/IWA   │
│ - Tải lượng hữu cơ vĩ mô (TOW) │ │ - Tiếp cận cân bằng cơ học/bùn  │ │ - Động học phản ứng phức tạp   │
│ - Hệ số tĩnh Bo, MCF mặc định  │ │ - Tách dòng CO2, CH4, Điện năng │ │ - Tập trung chuyên biệt chu    │
│ - Thiếu chi tiết vi sinh nội tại│ │ - Giả định trạng thái ổn định   │ │   trình Nito và phát thải N2O  │
└────────────────┬───────────────┘ └────────────────┬────────────────┘ └────────────────┬───────────────┘
                 │                                  │                                   │
                 └──────────────────────────────────┼───────────────────────────────────┘
                                                    │
                                                    ▼
                    ┌────────────────────────────────────────────────────────┐
                    │               ĐỊNH VỊ ĐỘT PHÁ CỦA LUẬN ÁN              │
                    │ - Động học Herbert cải biên tích hợp Cân bằng chất      │
                    │ - Công thức cơ chất thực nghiệm: C10H19NO3 & C7H12O4N │
                    │ - Bao quát toàn diện: Điện năng + Trực tiếp + N2O      │
                    │ - Code số Scala/Matlab & Thực nghiệm Chamber-GC        │
                    └────────────────────────────────────────────────────────┘

Trong y văn tồn tại hai luồng quan điểm mâu thuẫn sâu sắc về cơ chế phát sinh và hạch toán KNK trong xử lý nước thải:

• Tranh luận 1 (Về tính chất phát thải của $CO_2$ sinh học): Hướng dẫn của IPCC (2006) loại bỏ hoàn toàn $CO_2$ phát thải từ quá trình hô hấp vi sinh hiếu khí và phân hủy sinh học chất hữu cơ vì quy ước đây là nguồn gốc sinh học ngắn hạn (biogenic carbon neutral). Ngược lại, các công trình của Bridle (2008) và Doorn et al. (2006) lập luận rằng đối với các hệ thống công nghiệp (như sản xuất giấy, dệt nhuộm) và các trạm đô thị có bổ sung nguồn cơ chất ngoài gốc hóa thạch (methanol, acetate), việc bỏ qua $CO_2$ sinh học dẫn đến đánh giá sai lệch tổng dấu chân carbon (carbon footprint) và cân bằng năng lượng thực tế của trạm. Luận án định vị quan điểm tính toán đầy đủ phát thải $CO_2$ sinh hóa để phục vụ tối ưu hóa công nghệ, nhưng tách biệt rõ ràng khi lập báo cáo kiểm kê theo chuẩn IPCC.
• Tranh luận 2 (Về nguồn gốc phát thải $N_2O$): Czepiel et al. (1995) và Kampschreur et al. (2009) khẳng định phát thải $N_2O$ xảy ra chủ yếu trong điều kiện thiếu khí (anoxic) do khử nitrat không hoàn toàn. Trong khi đó, các nghiên cứu của Wunderlin et al. (2012) và Law et al. (2012) chứng minh vi khuẩn oxy hóa amoni tự dưỡng (Ammonia-Oxidizing Bacteria - AOB) trong điều kiện hiếu khí có nồng độ oxy hòa tan thấp ($DO < 1,0\text{ mg/L}$) lại là nguồn phát sinh $N_2O$ chính thông qua con đường nitrifier denitrification.

So sánh với các nghiên cứu quốc tế điển hình:

• So với mô hình tính toán tại Trạm xử lý nước thải Windsor (Ontario, Canada) do Bridle Consulting thực hiện (2008), luận án kế thừa khung cân bằng khối lượng nhưng cải tiến căn bản bằng việc nhúng mô hình động học Herbert vào hệ thống phương trình vi phân phi tuyến, cho phép mô phỏng động học cả trong điều kiện không ổn định (unsteady-state).
• So với nghiên cứu của Foley et al. (2010) tại 15 nhà máy xử lý nước thải ở Úc (chỉ tập trung vào phát thải $N_2O$ bằng cảm biến trực tuyến), mô hình của luận án có tính toàn diện hơn khi kết hợp đồng thời mô phỏng và đo đạc thực nghiệm cả 3 loại khí ($CO_2, CH_4, N_2O$) và phát thải điện năng gián tiếp trong cùng một khung phần mềm thống nhất.

Đó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 lý thuyết động học tăng trưởng vi sinh vật cổ điển Monod thông qua việc ứng dụng và hoàn thiện phương trình động học Herbert trong kỹ thuật môi trường. Phương trình Monod truyền thống chỉ xét đến tốc độ tăng trưởng riêng phụ thuộc vào nồng độ chất nền giới hạn mà bỏ qua tốc độ phân hủy nội bào (endogenous decay) và sự hình thành các mảnh vụn tế bào không phân hủy sinh học (non-biodegradable cell debris).

          ┌────────────────────────────────────────────────────────────────────────┐
          │                    KHUNG LÝ THUYẾT ĐỘNG HỌC MỞ RỘNG                    │
          └───────────────────────────────────┬────────────────────────────────────┘
                                              │
                    ┌─────────────────────────┴─────────────────────────┐
                    ▼                                                   ▼
  ┌───────────────────────────────────┐               ┌───────────────────────────────────┐
  │      MÔ HÌNH MONOD TRUYỀN THỐNG   │               │       MÔ HÌNH HERBERT CẢI BIÊN    │
  │  - Tốc độ tăng trưởng riêng:      │               │  - Tốc độ tăng trưởng thực tế:    │
  │    µ = µmax * [S / (Ks + S)]      │               │    rx = (µ - kd) * X              │
  │  - Bỏ qua khoáng hóa nội bào      │               │  - Tích hợp hệ số phân hủy kd     │
  │  - Bỏ qua mảnh vụn tế bào trơ     │               │  - Tích hợp mảnh vụn trơ fd      │
  └─────────────────┬─────────────────┘               └─────────────────┬─────────────────┘
                    │                                                   │
                    └─────────────────────────┬─────────────────────────┘
                                              │
                                              ▼
                    ┌────────────────────────────────────────────────────────┐
                    │            CÂN BẰNG KHỐI LƯỢNG VÀ HÓA HỌC MỞ RỘNG      │
                    │                                                        │
                    │ Cơ chất Sinh hoạt:      C10H19NO3 + O2  ──► CO2 + X    │
                    │ Cơ chất Giấy:           C7H12O4N + O2   ──► CO2 + X    │
                    │ Sinh khối Vi sinh vật:  C5H7O2N         ──► CO2 + CH4  │
                    │                                                        │
                    │               Lượng bùn sinh học tạo thành:             │
                    │       Px = Qv * Y * (So,v - Sv) / [1 + kd * SRT]       │
                    └────────────────────────────────────────────────────────┘

Luận án thiết lập hệ phương trình mô tả mối quan hệ giữa tốc độ sử dụng cơ chất riêng ($r_s$) và tốc độ tăng trưởng sinh khối thực tế ($r_x$):

$$r_s = \frac{k \cdot S}{K_s + S} \cdot X$$

$$r_x = Y \cdot r_s - k_d \cdot X$$

Trong đó:

• $k$: Tốc độ sử dụng chất nền cực đại ($\text{mg}/\text{mg}\cdot\text{ngày}$).
• $K_s$: Hằng số bán bão hòa ($\text{mg}/L$).
• $Y$: Hệ số năng suất sinh khối cực đại ($\text{mg sinh khối}/\text{mg cơ chất}$).
• $k_d$: Hệ số phân hủy nội bào ($\text{ngày}^{-1}$).
• $X$: Nồng độ sinh khối vi sinh vật hoạt tính ($\text{mg}/L$).
• $S$: Nồng độ cơ chất hòa tan trong bể phản ứng ($\text{mg}/L$).

Sự chuyển đổi mô hình này mang tính đột phá khi liên kết trực tiếp tốc độ phân hủy nội bào $k_d$ với lượng khí $CO_2$ sinh ra do hô hấp nội bào và lượng sinh khối dư thừa $P_X$ cần xả bỏ sang bể phân hủy bùn kỵ khí:

$$P_X = \frac{Q_v \cdot Y \cdot (S_{o,v} - S_v)}{1 + k_d \cdot \text{SRT}} + f_d \cdot k_d \cdot X \cdot V$$

Trong đó $\text{SRT}$ là tuổi bùn (ngày), $f_d$ là phần sinh khối thu được từ mảnh vụn tế bào (thường lấy giá trị mặc định 0,15), $Q_v$ là lưu lượng dòng vào ($\text{m}^3/\text{ngày}$), $S_{o,v}$ và $S_v$ lần lượt là nồng độ cơ chất ban đầu và nồng độ cơ chất trong bể sinh học ($\text{mg}/L$).

Khung phân tích độc đáo

Khung phân tích của luận án là sự tích hợp đa ngành giữa:

1. Lý thuyết Động học Vi sinh vật (Herbert-Monod Kinetics): Mô tả tốc độ phản ứng tiêu thụ cơ chất hữu cơ và oxy hóa các hợp chất nitơ.
2. Lý thuyết Cân bằng Khối lượng và Cân bằng Hóa học Thực nghiệm (Stoichiometric Mass Balance): Sử dụng các công thức thực nghiệm đại diện chuẩn mực:
   • Nước thải sinh hoạt: $C_{10}H_{19}NO_3$
   • Nước thải công nghiệp sản xuất giấy: $C_7H_{12}O_4N$
   • Tế bào vi sinh vật (sinh khối): $C_5H_7O_2N$
3. Lý thuyết Hệ thống Cân bằng Năng lượng (Energy-Emission Coupling): Quy đổi điện năng tiêu thụ từ lưới điện quốc gia sang phát thải KNK dựa trên tỷ trọng cơ cấu nguồn phát điện tại Việt Nam (thủy điện, nhiệt điện than, nhiệt điện khí).

         ┌────────────────────────────────────────────────────────────────────────┐
         │              KHUNG PHÂN TÍCH TỔNG THỂ CỦA MÔ HÌNH TOÁN HỌC             │
         └───────────────────────────────────┬────────────────────────────────────┘
                                             │
      ┌──────────────────────────────────────┼──────────────────────────────────────┐
      ▼                                      ▼                                      ▼
┌───────────────────────────┐  ┌───────────────────────────┐  ┌───────────────────────────┐
│     HỆ THỐNG A (HIẾU KHÍ) │  │    HỆ THỐNG B (YẾM KHÍ)   │  │   HỆ THỐNG C (KẾT HỢP)    │
│ - Bể Bùn hoạt tính (AER)  │  │ - Bể UASB / Anaerobic     │  │ - Thiếu khí - Hiếu khí    │
│ - Bể Lắng sơ cấp/thứ cấp  │  │ - Quá trình lên men CH4   │  │ - Phân hủy bùn kỵ khí     │
│ - Phân hủy bùn kỵ khí     │  │ - Thu hồi/phóng không CH4 │  │ - Xử lý Nito & BOD sâu    │
└─────────────┬─────────────┘  └─────────────┬─────────────┘  └─────────────┬─────────────┘
              │                              │                              │
              └──────────────────────────────┼──────────────────────────────┘
                                             │
                                             ▼
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│                           HAI DÒNG PHÁT THẢI CẤU THÀNH                                  │
├────────────────────────────────────────────┬────────────────────────────────────────────┤
│     1. NGUỒN PHÁT THẢI TRỰC TIẾP           │     2. NGUỒN PHÁT THẢI GIÁN TIẾP           │
│ - Khí CO2: Hô hấp vi sinh + Nội bào        │ - Điện năng tiêu thụ cho sục khí (Aeration)│
│ - Khí CH4: Phân hủy bùn kỵ khí / Bể yếm khí│ - Bơm, khuấy trộn, ép bùn, chiếu sáng      │
│ - Khí N2O: Quá trình Nitrat / Khử Nitrat   │ - Cơ cấu nguồn điện VN: Than, Khí, Thủy điện│
└────────────────────────────────────────────┴────────────────────────────────────────────┘

Phương trình hóa học cho quá trình hiếu khí đối với nước thải sinh hoạt được xác lập như sau:

$$C_{10}H_{19}NO_3 + a O_2 + b NH_4^+ + b HCO_3^- \longrightarrow c C_5H_7O_2N + d CO_2 + e H_2O$$

Đối với nước thải công nghiệp giấy:

$$C_7H_{12}O_4N + a' O_2 \longrightarrow c' C_5H_7O_2N + d' CO_2 + e' H_2O$$

Đối với quá trình phân hủy bùn kỵ khí:

$$C_5H_7O_2N + 2 H_2O \longrightarrow 2,5 CH_4 + 2,5 CO_2 + NH_3$$

Điều kiện biên và giới hạn áp dụng (boundary conditions): Mô hình áp dụng cho các HTXLNT xử lý bằng phương pháp sinh học với dòng chảy liên tục (continuous flow stirred-tank reactor - CSTR), nồng độ chất hữu cơ nằm trong khoảng $100 - 2.500\text{ mg } BOD_5/L$, nhiệt độ vận hành từ $15^\circ\text{C}$ đến $35^\circ\text{C}$, và nồng độ oxy hòa tan trong bể hiếu khí được duy trì $DO \ge 1,5\text{ mg } O_2/L$.

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

Thiết kế nghiên cứu

Nghiên cứu áp dụng quan điểm thực chứng (positivism) kết hợp mô hình hóa giải tích tất định (deterministic mathematical modeling) và kiểm chứng thực nghiệm tại hiện trường. Thiết kế nghiên cứu đa tầng bậc (multi-level design) gồm ba cấp độ:

• Cấp độ 1 (Lý thuyết và Giải thuật số): Xây dựng hệ phương trình cân bằng chất và động học vi sinh; mã hóa phần mềm bằng ngôn ngữ Scala và MATLAB R2017b.
• Cấp độ 2 (Kiểm chuẩn quốc tế): Chạy kiểm nghiệm mô hình bằng cách tái lập dữ liệu vận hành chuẩn của Trạm xử lý nước thải Windsor (Canada) từ công bố của Bridle (2008).
• Cấp độ 3 (Thực nghiệm hiện trường): Quan trắc trực tiếp, lấy mẫu và phân tích sắc ký khí tại Trạm xử lý nước thải sinh hoạt Goldmark City và Nhà máy Giấy Bãi Bằng.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn quốc tế:

• Giao thức lấy mẫu khí: Sử dụng kỹ thuật hộp kín (Chamber QUT) được chuyển giao từ Đại học Công nghệ Queensland (Úc) cho Viện Môi trường Nông nghiệp Việt Nam. Hộp chamber hình trụ kín được đặt nổi trên mặt thoáng của các bể xử lý hiếu khí, bể phân hủy bùn và mương dẫn nước thải. Khí tích tụ trong buồng kín được rút định kỳ tại các thời điểm $0, 10, 20, 30\text{ phút}$ bằng bơm tiêm khí kín khí chuyên dụng và bơm vào các ống nghiệm chân không thủy tinh (vacutainer) $20\text{ mL}$.
• Quy trình phân tích sắc ký khí (Gas Chromatography - GC): Phân tích đồng thời 3 thành phần KNK trên hệ thống sắc ký khí cao cấp GC-2014 của Shimadzu (Nhật Bản):
  • Khí $CO_2$: Xác định bằng Detector dẫn nhiệt ($TCD$), nhiệt độ buồng dò $200^\circ\text{C}$, cột phân tích Porapak Q.
  • Khí $CH_4$: Xác định bằng Detector ion hóa ngọn lửa ($FID$), nhiệt độ buồng dò $300^\circ\text{C}$, khí mang $N_2$ siêu sạch.
  • Khí $N_2O$: Xác định bằng Detector bắt giữ điện tử ($ECD$), nhiệt độ buồng dò $350^\circ\text{C}$, sử dụng khí mang $Ar/CH_4$ (95/5%).
• Hiệu chuẩn và độ tin cậy: Đường chuẩn phân tích khí được dựng bằng các mẫu khí chuẩn đa thành phần được chứng nhận quốc tế (Standard Gas Mixtures - Scott Gas / Air Liquide) với hệ số tương quan tuyến tính $R^2 \ge 0,999$. Mỗi điểm quan trắc được lấy lặp lại tối thiểu 3 lần ($n=3$) để tính độ lệch chuẩn và khoảng tin cậy.

       ┌────────────────────────────────────────────────────────────────────────┐
       │                   SƠ ĐỒ QUY TRÌNH THỰC NGHIỆM VÀ PHÂN TÍCH             │
       └───────────────────────────────────┬────────────────────────────────────┘
                                           │
         ┌─────────────────────────────────┴─────────────────────────────────┐
         ▼                                                                   ▼
┌─────────────────────────────────┐                         ┌─────────────────────────────────┐
│     LẤY MẪU NGOÀI HIỆN TRƯỜNG   │                         │     PHÂN TÍCH SẮC KÝ KHÍ (GC)   │
│ - Thiết bị: Chamber QUT (Úc)    │                         │ - Thiết bị: GC-2014 SHIMADZU    │
│ - Tần suất lấy: 0, 10, 20, 30 ph│                         │ - Khí CO2: Detector TCD (200°C) │
│ - Mẫu lưu trong Vacutainer chân │ ──────────────────────► │ - Khí CH4: Detector FID (300°C) │
│   không chuyên dụng             │                         │ - Khí N2O: Detector ECD (350°C) │
└─────────────────────────────────┘                         └────────────────┬────────────────┘
                                                                             │
                                                                             ▼
┌─────────────────────────────────────────────────────────────────────────────────────────────┐
│                            KIỂM CHỨNG VÀ ĐÁNH GIÁ MÔ HÌNH                                   │
│ - Đối chiếu mô hình giải tích lập trình (Scala/Matlab) với dữ liệu thực nghiệm GC           │
│ - Kiểm nghiệm độ nhạy: Nhiệt độ T (15-35°C), Tuổi bùn SRT (5-25 ngày), Tải lượng vào So,v │
└─────────────────────────────────────────────────────────────────────────────────────────────┘

Data và phân tích

Số liệu vận hành và đặc tính lý hóa của các trạm nghiên cứu:

• Trạm Goldmark City: Lưu lượng trung bình $Q = 4.200\text{ m}^3/\text{ngày}$; $BOD_5$ vào $= 250\text{ mg}/L$; $BOD_5$ ra $= 18\text{ mg}/L$; $SS$ vào $= 220\text{ mg}/L$; $TN$ vào $= 45\text{ mg}/L$; $TN$ ra $= 8,5\text{ mg}/L$. Công suất điện tiêu thụ trung bình $Q_E = 1.450\text{ kWh}/\text{ngày}$.
• Trạm Giấy Bãi Bằng: Lưu lượng $Q = 40.000\text{ m}^3/\text{ngày}$; $COD$ vào $= 1.200\text{ mg}/L$; $COD$ ra $= 110\text{ mg}/L$; $BOD_5$ vào $= 480\text{ mg}/L$; $BOD_5$ ra $= 35\text{ mg}/L$. Công suất điện tiêu thụ $Q_E = 28.500\text{ kWh}/\text{ngày}$.
• Hệ số phát thải lưới điện quốc gia ($EF_i$): Dựa trên cơ cấu nguồn điện năm 2016 (Thủy điện 42,3% - $EF = 0,02\text{ kg } CO_{2\text{tđ}}/\text{kWh}$; Than đá 34,4% - $EF = 0,98\text{ kg } CO_{2\text{tđ}}/\text{kWh}$; Khí đốt 21,8% - $EF = 0,55\text{ kg } CO_{2\text{tđ}}/\text{kWh}$), hệ số phát thải trung bình của điện lưới Việt Nam được tính toán chính xác là $0,562\text{ kg } CO_{2\text{tđ}}/\text{kWh}$.

Hệ số động học thực nghiệm xác định cho hệ thống hiếu khí ở $20^\circ\text{C}$:

• Tốc độ sử dụng chất nền cực đại $k = 4,5\text{ mg } BOD_5/\text{mg sinh khối}\cdot\text{ngày}$.
• Hằng số bán bão hòa $K_s = 60\text{ mg } BOD_5/L$.
• Hệ số năng suất sinh khối cực đại $Y = 0,6\text{ mg sinh khối}/\text{mg } BOD_5$.
• Hệ số phân hủy nội bào $k_d = 0,06\text{ ngày}^{-1}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả mô phỏng và đo đạc thực địa đã xác lập 5 phát hiện mang tính đột phá khoa học:

1. Lượng hóa chính xác hệ số phát thải KNK cấp quốc gia cho xử lý nước thải:

   • Nước thải sinh hoạt: Trong trường hợp không thu hồi khí $CH_4$ (phóng không từ bể tự hoại và bể phân hủy bùn), hệ số phát thải đạt $2,40\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$. Trong trường hợp thu hồi và đốt/tái sử dụng năng lượng khí $CH_4$, hệ số phát thải giảm sâu xuống còn $1,16\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ (mức cắt giảm phát thải đạt $51,67%$).
   • Nước thải công nghiệp sản xuất giấy: Trong trường hợp không thu hồi khí $CH_4$, hệ số phát thải đạt $3,58\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$. Khi thu hồi và đốt khí $CH_4$, hệ số phát thải chỉ còn $0,93\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ (mức cắt giảm phát thải đạt $74,02%$).

2. Tỷ trọng chi phối tuyệt đối của phát thải gián tiếp và phân hủy bùn: Tại hệ thống Goldmark City, khi phóng không $CH_4$, nguồn phát thải từ bể phân hủy bùn kỵ khí chiếm tới $53,4%$ tổng phát thải KNK toàn trạm, phát thải từ tiêu thụ điện năng sục khí chiếm $38,2%$, phát thải $CO_2$ sinh học tại bể hiếu khí chiếm $7,1%$, và phát thải $N_2O$ chiếm $1,3%$. Khi thu hồi và đốt $CH_4$, tỷ trọng phát thải do tiêu thụ điện năng vọt lên chiếm vị trí số 1 với $78,9%$ tổng phát thải của toàn trạm.

       CƠ CẤU PHÁT THẢI TẠI TRẠM XỬ LÝ NƯỚC THẢI SINH HOẠT (GOLDMARK CITY)
       
       Trường hợp 1: Phóng không CH4              Trường hợp 2: Thu hồi & Đốt CH4
       ┌─────────────────────────────────┐        ┌─────────────────────────────────┐
       │ Phân hủy bùn (CH4):       53.4% │        │ Tiêu thụ điện năng:       78.9% │
       │ Tiêu thụ điện năng:       38.2% │        │ CO2 sinh học bể hiếu khí: 14.6% │
       │ CO2 sinh học bể hiếu khí:  7.1% │        │ Phát thải N2O:             2.7% │
       │ Phát thải N2O:             1.3% │        │ Phân hủy bùn (dư):         3.8% │
       └─────────────────────────────────┘        └─────────────────────────────────┘

3. Quy luật tác động phi tuyến của nhiệt độ môi trường: Khi nhiệt độ nước thải tăng từ $15^\circ\text{C}$ lên $35^\circ\text{C}$, tốc độ phân hủy vi sinh $k_d$ tăng từ $0,04\text{ ngày}^{-1}$ lên $0,095\text{ ngày}^{-1}$, dẫn đến lượng $CO_2$ sinh ra do hô hấp nội bào tăng $38%$. Tuy nhiên, lượng sinh khối dư thừa giảm $22%$, làm giảm tổng lượng bùn cần phân hủy kỵ khí, qua đó giảm phát thải $CH_4$ tại bể phân hủy bùn khoảng $15%$.

4. Phát hiện bất ngờ về động học tuổi bùn (SRT): Kéo dài tuổi bùn SRT từ 5 ngày lên 20 ngày giúp tăng cường hiệu quả xử lý nước và khoáng hóa bùn, làm tăng lượng $CO_2$ sinh học tại bể hiếu khí thêm $28%$, nhưng giúp giảm lượng bùn hoạt tính dư thừa sinh ra $41%$. Vì $CH_4$ có tiềm năng ấm lên toàn cầu $GWP = 25$, cao hơn nhiều so với $CO_2$ ($GWP = 1$), tổng lượng phát thải KNK quy đổi toàn hệ thống giảm $18,4%$ khi tăng SRT trong trường hợp không thu hồi $CH_4$.

5. Độ chính xác và tương thích cao của mô hình số thiết lập: So sánh kết quả tính toán của mô hình với số liệu đo đạc thực địa bằng Chamber-GC tại Goldmark City cho thấy sai số phát thải $CO_2$ tại bể hiếu khí chỉ là $6,8%$, sai số $CH_4$ tại bể bùn là $9,2%$, và sai số $N_2O$ là $13,5%$. Điều này chứng minh độ tin cậy vượt trội của mô hình giải tích so với sai số lên tới $45 - 80%$ khi dùng hệ số mặc định IPCC Tier 1.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp phương pháp tiếp cận mới tích hợp giữa sinh học phân hủy chất ô nhiễm và khí động học môi trường, mở rộng lý thuyết cân bằng vật chất trong các hệ sinh thái nhân tạo.
• Về mặt phương pháp luận: Cung cấp bộ công cụ mã nguồn số hóa và giao diện phần mềm trực quan, cho phép tính toán phát thải KNK nhanh chóng cho cả kỹ sư vận hành lẫn các nhà tư vấn thiết kế môi trường.
• Về mặt ứng dụng thực tiễn: Chỉ rõ con đường tối ưu hóa trạm XLNT theo hướng "Carbon trung tính":
  1. Thay thế máy thổi khí truyền thống bằng thiết bị biến tần hiệu suất cao và máy khuấy chìm tiết kiệm năng lượng để giảm thiểu tiêu thụ điện (vốn chiếm gần $80%$ phát thải khi đã thu hồi khí).
  2. Bắt buộc thu hồi khí biogas từ bể phân hủy bùn kỵ khí để đốt phát điện hoặc tận dụng nhiệt, triệt tiêu nguồn phát thải $CH_4$.
  3. Kiểm soát nồng độ $DO$ tối ưu ($1,8 - 2,2\text{ mg}/L$) để ngăn chặn triệt để sự tạo thành khí $N_2O$ do thiếu khí cục bộ.
• Về mặt chính sách quản lý: Là căn cứ kỹ thuật chuẩn xác để Bộ Tài nguyên và Môi trường cập nhật cơ sở dữ liệu hệ số phát thải quốc gia cho báo cáo BUR (Biennial Update Report) và Báo cáo Quốc gia Định kỳ gửi UNFCCC, hỗ trợ xây dựng thị trường hạn ngạch phát thải carbon tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật và thực tiễn nổi bật, luận án thừa nhận 4 hạn chế nghiên cứu mang tính điều kiện biên:

1. Giới hạn phạm vi phát thải gián tiếp: Mô hình chỉ tính toán phát thải gián tiếp từ điện năng tiêu thụ, chưa bao gồm dấu chân carbon từ quá trình sản xuất và vận chuyển hóa chất keo tụ (như $Al_2(SO_4)_3$, $FeCl_3$, polymer), chất kiềm hóa chất điều chỉnh pH (NaOH, vôi), hoặc nguồn carbon bổ sung (methanol) do hạn chế về chuỗi dữ liệu chu trình sống (LCA).
2. Động học phát sinh $N_2O$ chưa được mô hình hóa chi tiết theo cấp enzym: Việc tính toán $N_2O$ dựa trên hệ số chuyển đổi nitơ dòng vào và các biểu thức kinh nghiệm, chưa tích hợp đầy đủ hệ phương trình vi phân mô tả hoạt tính enzym hydroxylamine oxidoreductase và nitric oxide reductase trong điều kiện $DO$ biến thiên liên tục.
3. Quy mô loại hình nước thải: Luận án mới tập trung vào nước thải sinh hoạt và nước thải công nghiệp sản xuất giấy; các ngành có tải lượng ô nhiễm hữu cơ và độc chất cao khác như dệt nhuộm, chế biến thủy hải sản, thuộc da, tinh bột sắn chưa được xây dựng hệ số động học riêng.
4. Dữ liệu kiểm chứng thời gian thực: Các đợt đo đạc thực nghiệm bằng hộp chamber mới thực hiện theo các chiến dịch quan trắc định kỳ, chưa triển khai được hệ thống cảm biến quang phổ hồng ngoại (NDIR) hoặc laser đo liên tục 24/7 để theo dõi đỉnh phát thải theo chu kỳ ngày - đêm.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

• Hướng nghiên cứu 1: Mở rộng mô hình toán học tích hợp đánh giá vòng đời toàn diện (Life Cycle Assessment - LCA), bao gồm toàn bộ phát thải từ thượng nguồn sản xuất hóa chất và hạ nguồn tái sử dụng bùn thải làm vật liệu xây dựng.
• Hướng nghiên cứu 2: Nhúng phân hệ mô phỏng vi sinh vật học phân tử (Metagenomics and Enzyme Kinetics) để mô hình hóa chính xác sự biểu hiện của các gen vi sinh (amoA, nirS, nirK, nosZ) trong phát thải $N_2O$.
• Hướng nghiên cứu 3: Mở rộng thư viện hệ số phát thải và cơ sở dữ liệu động học vi sinh cho 10 ngành công nghiệp trọng điểm tại Việt Nam.
• Hướng nghiên cứu 4: Phát triển module tối ưu hóa kinh tế - môi trường đa mục tiêu (Multi-objective Genetic Algorithm - MOGA) nhằm tìm kiếm điểm cân bằng tối ưu giữa chi phí vận hành trạm XLNT và tổng lượng phát thải KNK.

Tác động và ảnh hưởng

• Tác động học thuật: Luận án đã công bố các bài báo khoa học chất lượng trên các tạp chí chuyên ngành uy tín trong nước và quốc tế; mở ra một nhánh nghiên cứu chuyên sâu về "Xử lý nước thải phát thải carbon thấp" tại Viện Khoa học và Công nghệ Môi trường - Đại học Bách khoa Hà Nội.
• Chuyển đổi ngành công nghiệp xử lý môi trường: Cung cấp giải pháp cho các tập đoàn bất động sản (như Vingroup, TNR Holdings - chủ đầu tư Goldmark City) và các tổng công ty giấy (Tổng công ty Giấy Việt Nam - Bãi Bằng) phương pháp luận để kiểm toán năng lượng - phát thải, từ đó đầu tư cải tạo hệ thống sục khí và thu hồi biogas, giúp tiết kiệm hàng tỷ đồng chi phí năng lượng mỗi năm.
• Đóng góp cho quản trị khí hậu quốc gia: Hỗ trợ Cục Biến đổi Khí hậu và Tổng cục Môi trường (Bộ Tài nguyên và Môi trường) nâng cấp phương pháp kiểm kê khí nhà kính lĩnh vực chất thải từ Cấp 1 (Tier 1) lên Cấp 2/3 (Tier 2/3), nâng cao uy tín quốc tế của Việt Nam trong việc thực thi Thỏa thuận Paris về Biến đổi Khí hậu (COP21) và cam kết Net Zero tại COP26.
• Lợi ích xã hội và toàn cầu: Cắt giảm hàng trăm nghìn tấn $CO_{2\text{tđ}}$ phát thải vào khí quyển mỗi năm thông qua các giải pháp thu hồi mê-tan và tối ưu hóa vận hành, góp phần cải thiện chất lượng môi trường không khí đô thị và sức khỏe cộng đồng.

Đối tượng hưởng lợi

                              CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI TRỌNG TÂM
                                                │
         ┌──────────────────────────────────────┼──────────────────────────────────────┐
         ▼                                      ▼                                      ▼
┌────────────────────────────────┐ ┌─────────────────────────────────┐ ┌────────────────────────────────┐
│      GIỚI HỌC THUẬT & NCS      │ │    DOANH NGHIỆP VÀ KỸ SƯ R&D    │ │      CƠ QUAN QUẢN LÝ NHÀ NƯỚC  │
│ - Kế thừa mô hình toán Herbert │ │ - Ứng dụng code số mô phỏng trạm│ │ - Cơ sở dữ liệu kiểm kê Tier 2 │
│ - Tham chiếu hệ số động học k, │ │ - Tối ưu hóa năng lượng sục khí │ │ - Xây dựng khung hạn ngạch KNK │
│   Ks, Y, kd chuẩn thực nghiệm  │ │ - Giảm chi phí vận hành 20-30%  │ │ - Cập nhật báo cáo BUR/UNFCCC  │
└────────────────────────────────┘ └─────────────────────────────────┘ └────────────────────────────────┘

1. Nghiên cứu sinh và Giới nghiên cứu học thuật: Kế thừa khung lý thuyết động học Herbert cải biên, phương pháp cân bằng hóa học cơ chất đặc thù ($C_{10}H_{19}NO_3$, $C_7H_{12}O_4N$) và bộ số liệu động học thực nghiệm chuẩn hóa tại Việt Nam làm tiền đề mở rộng nghiên cứu.
2. Các chuyên gia và Kỹ sư thiết kế R&D Môi trường: Sử dụng phần mềm mô hình hóa để mô phỏng, dự báo chính xác tải lượng KNK ngay từ giai đoạn thiết kế cơ sở (FEED), lựa chọn phương án công nghệ tối ưu giữa bùn hoạt tính truyền thống, mương oxy hóa (OD), hay công nghệ màng sinh học (MBBR/MBR).
3. Doanh nghiệp và Ban quản lý trạm xử lý nước thải: Nhận diện rõ điểm nghẽn tiêu thụ năng lượng và nguồn rò rỉ khí $CH_4$, triển khai giải pháp kỹ thuật cụ thể giúp giảm $20 - 30%$ hóa đơn tiền điện và đăng ký nhận tín chỉ carbon (carbon credits) từ các cơ chế thị trường như CDM hoặc Gold Standard.
4. Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Cung cấp bộ hệ số phát thải thực chứng để ban hành các quy chuẩn kỹ thuật quốc gia về ngưỡng phát thải KNK đối với các công trình xử lý chất thải, hỗ trợ thực thi Luật Bảo vệ Môi trường 2020.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc hoàn thiện và tích hợp mô hình động học Herbert với phương trình cân bằng chất cơ chất - sinh khối đặc thù cho điều kiện nhiệt đới Việt Nam. Luận án đã mở rộng lý thuyết tăng trưởng vi sinh vật cổ điển Monod bằng cách đưa hệ số phân hủy nội bào $k_d$ và mảnh vụn tế bào trơ $f_d$ vào phương trình cân bằng hóa học stoichiometry, liên kết trực tiếp động học khoáng hóa sinh khối với tốc độ phát thải khí $CO_2$ sinh học và sinh khối bùn xả thải $P_X$.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?

• So với mô hình Bridle Consulting (Bridle, 2008): Mô hình của Bridle chủ yếu áp dụng phương pháp cân bằng tĩnh ở trạng thái ổn định (steady-state) và sử dụng hệ số phát thải cố định cho vùng ôn đới. Luận án đã nâng cấp thành hệ phương trình vi phân giải tích động học, cho phép mô phỏng các biến thiên nhiệt độ, tải lượng cơ chất đầu vào không ổn định (unsteady-state).
• So với phương pháp kiểm kê IPCC Tier 1/2 (IPCC, 2006): Phương pháp IPCC dựa vào tải lượng hữu cơ tổng thể ($TOW$) và hệ số $MCF$ giả định mà không mô tả cấu trúc bể phản ứng. Luận án phân rã chi tiết từng công trình đơn vị (bể hiếu khí, bể kỵ khí, bể phân hủy bùn) và tích hợp phát thải gián tiếp từ lưới điện Việt Nam ($0,562\text{ kg } CO_{2\text{tđ}}/\text{kWh}$), giảm sai số ước tính từ $\pm 50%$ xuống dưới $10%$.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì và giải thích cơ chế khoa học?

Phát hiện bất ngờ nhất là việc tăng tuổi bùn $\text{SRT}$ (từ 5 ngày lên 20 ngày) tuy làm tăng lượng $CO_2$ sinh học phát thải tại bể hiếu khí thêm $28%$, nhưng lại làm giảm $18,4%$ tổng phát thải KNK quy đổi toàn hệ thống trong điều kiện không thu hồi khí $CH_4$. Cơ chế khoa học: Khi $\text{SRT}$ dài, vi sinh vật trải qua quá trình tự phân hủy nội bào sâu sắc, biến đổi sinh khối thành $CO_2$ (có $GWP = 1$). Điều này làm giảm mạnh lượng bùn hữu cơ xả sang bể phân hủy kỵ khí, từ đó triệt tiêu đáng kể lượng khí $CH_4$ sinh ra (có $GWP = 25$). Hiệu ứng giảm phát thải $CH_4$ vượt trội hơn hẳn sự gia tăng phát thải $CO_2$.

4. Giao thức nhân bản nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Có. Luận án cung cấp tường minh toàn bộ giao thức nhân bản:

• Hệ thống phương trình giải tích toán học đầy đủ và thuật toán giải phương trình vi phân.
• Cấu trúc mã code số hóa bằng ngôn ngữ Scala và script MATLAB.
• Quy trình lấy mẫu chi tiết bằng hộp Chamber QUT (kích thước, thời gian lấy mẫu, cách bảo quản ống vacutainer).
• Điều kiện vận hành thiết bị sắc ký khí GC-2014 Shimadzu ($TCD$ $200^\circ\text{C}$, $FID$ $300^\circ\text{C}$, $ECD$ $350^\circ\text{C}$).

5. Lộ trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?

Lộ trình 10 năm gồm 3 giai đoạn kế tiếp:

• Giai đoạn 1 (Năm 1 - 3): Mở rộng kiểm chứng mô hình trên 20 trạm XLNT thuộc các ngành dệt nhuộm, chế biến thủy sản, hóa chất; nâng cấp phần mềm thành nền tảng điện toán đám mây SaaS.
• Giai đoạn 2 (Năm 4 - 6): Tích hợp cảm biến IoT quan trắc KNK trực tuyến thời gian thực và mô hình hóa mạng nơ-ron nhân tạo (AI/Machine Learning) để dự báo tự động phát thải theo thời gian thực.
• Giai đoạn 3 (Năm 7 - 10): Phát triển hệ thống điều khiển tối ưu hóa tự động SCADA nhằm điều chỉnh lưu lượng sục khí và tuần hoàn bùn, hiện thực hóa các nhà máy xử lý nước thải phát thải ròng bằng Không (Net Zero WWTPs) trên quy mô toàn quốc.

Kết luận

Công trình nghiên cứu luận án tiến sĩ của NCS. Nguyễn Thị Vân Anh đã xác lập các đóng góp khoa học và thực tiễn cốt lõi:

1. Thiết lập thành công mô hình toán học giải tích tất định đầu tiên tại Việt Nam tích hợp phương trình động học Herbert, cân bằng vật chất cơ chất - sinh khối và cân bằng năng lượng để mô phỏng toàn diện phát thải KNK ($CO_2, CH_4, N_2O$, Điện năng).
2. Xây dựng thành công công cụ code số trên nền tảng Scala và MATLAB với giao diện trực quan, cho phép tính toán phát thải KNK với độ chính xác cao (sai số dưới $10%$ so với thực nghiệm đo đạc bằng sắc ký khí GC-2014).
3. Lượng hóa và công bố bộ hệ số phát thải KNK cấp quốc gia chuẩn xác cho nước thải sinh hoạt ($2,40\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ khi phóng không $CH_4$; $1,16\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ khi thu hồi $CH_4$) và nước thải công nghiệp giấy ($3,58\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ khi phóng không $CH_4$; $0,93\text{ kg } CO_{2\text{tđ}}/\text{kg } BOD_5$ khi thu hồi $CH_4$).
4. Làm sáng tỏ định lượng cơ chế tác động của nhiệt độ nhiệt đới ($15 - 35^\circ\text{C}$), tuổi bùn $\text{SRT}$ và nồng độ cơ chất đầu vào tới sự biến thiên phát thải của từng công trình đơn vị.
5. Chứng minh vai trò quyết định của giải pháp thu hồi khí biogas và tối ưu hóa năng lượng sục khí, mở ra tiềm năng cắt giảm từ $51%$ đến $74%$ dấu chân carbon cho các trạm xử lý nước thải tại Việt Nam.

Luận án đánh dấu bước tiến tiên phong trong kỹ thuật môi trường, tạo dựng cơ sở khoa học vững chắc để Việt Nam chủ động kiểm kê KNK chính xác và xây dựng lộ trình phát triển ngành nước bền vững hướng tới mục tiêu phát thải ròng bằng Không.