Tổng quan về luận án

Biến đổi khí hậu toàn cầu đang diễn ra với quy mô và tốc độ chưa từng có, đặt ra những thách thức nghiêm trọng đối với an ninh sinh thái và sự phát triển bền vững của nhân loại. Trong bức tranh toàn cầu đó, Việt Nam được đánh giá là một trong năm quốc gia chịu tổn thương nặng nề nhất bởi các hiện tượng thời tiết cực đoan và nước biển dâng. Nguồn gốc căn bản thúc đẩy biến đổi khí hậu bắt nguồn từ sự gia tăng không ngừng nồng độ các chất khí nhà kính (KNK) trong khí quyển do các hoạt động nhân sinh. Theo phân loại của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), các lĩnh vực phát thải chính bao gồm Năng lượng, Các quá trình công nghiệp và sử dụng sản phẩm (IPPU), Nông nghiệp, Lâm nghiệp và Sử dụng đất khác (AFOLU), và Lĩnh vực Chất thải.

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) năm 2010, tổng lượng KNK phát thải năm 2000 của Việt Nam đạt 150,9 triệu tấn $\text{CO}2$ tương đương ($\text{CO}{2\text{tđ}}$), trong đó lĩnh vực chất thải chiếm 7,9 triệu tấn $\text{CO}{2\text{tđ}}$ (khoảng 5,2%). Trong tổng phát thải từ chất thải, xử lý nước thải đóng góp xấp xỉ 30% và xử lý chất thải rắn chiếm 70%. Tổng lượng phát thải KNK quốc gia được dự báo tăng nhanh lên 300,4 triệu tấn $\text{CO}{2\text{tđ}}$ vào năm 2020 và đạt 515,8 triệu tấn $\text{CO}_{2\text{tđ}}$ vào năm 2030. Dù lĩnh vực xử lý nước thải có tiềm năng phát thải đáng kể các khí nhà kính nguy hại như Carbon dioxide ($\text{CO}_2$), Methane ($\text{CH}_4$, có tiềm năng gây ấm toàn cầu $\text{GWP} = 25$) và Nitrous oxide ($\text{N}_2\text{O}$, có $\text{GWP} = 296 - 298$), công tác kiểm kê KNK từ các hệ thống xử lý nước thải (HTXLNT) tại Việt Nam phần lớn vẫn áp dụng phương pháp cấp 1 (Tier 1) của IPCC với các hệ số phát thải mặc định toàn cầu. Phương pháp này bộc lộ sai số rất lớn, không phản ánh đúng động học vi sinh, điều kiện khí hậu nhiệt đới, công nghệ xử lý cũng như chế độ vận hành thực tế tại các đô thị và cơ sở công nghiệp Việt Nam. Mặt khác, việc đo đạc phát thải KNK trực tiếp bằng thực nghiệm liên tục đòi hỏi chi phí đầu tư thiết bị và nhân lực kỹ thuật quá lớn.

Khoảng trống nghiên cứu (research gap) cốt lõi chính là: Sự thiếu vắng một công cụ mô hình hóa toán học có cơ chế động học rõ ràng, tích hợp cân bằng vật chất dòng vào - dòng ra - tích lũy, có khả năng tính toán định lượng và dự báo phát thải KNK (cả trực tiếp lẫn gián tiếp) phù hợp với điều kiện vận hành thực tế của các HTXLNT sinh hoạt và công nghiệp tại Việt Nam.

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, dưới sự hướng dẫn khoa học của PGS. Đặng Xuân Hiển tại Trường Đại học Bách khoa Hà Nội (2019), mang tên: "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" (Mã ngành: 9520320) đã giải quyết triệt để khoảng trống khoa học này. Luận án đặt ra 4 câu hỏi nghiên cứu và 3 giả thuyết khoa học trọng tâm:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa toán học quá trình phát sinh KNK ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$) từ các quá trình sinh học hiếu khí, yếm khí và tiêu thụ điện năng dựa trên lý thuyết cân bằng vật chất và động học vi sinh?
  • Câu hỏi nghiên cứu 2 (RQ2): Các thông số động học sinh trưởng ($k, K_s, Y, k_d, \mu_m$) và các yếu tố công nghệ vận hành (nhiệt độ $T$, tuổi bùn $\text{SRT}$, nồng độ cơ chất $S_0$) tác động như thế nào đến mức độ phát thải KNK của các hệ thống xử lý nước thải tại Việt Nam?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ tin cậy và sai số của mô hình toán học khi kiểm nghiệm so với dữ liệu quốc tế chuẩn và số liệu đo đạc thực nghiệm bằng thiết bị buồng thu nổi (floating flux chamber) kết hợp sắc ký khí (GC) là bao nhiêu?
  • Câu hỏi nghiên cứu 4 (RQ4): Hệ số phát thải KNK đặc thù cho các hệ thống xử lý nước thải sinh hoạt và nước thải công nghiệp giấy tại Việt Nam được lượng hóa cụ thể ở mức nào?

Hệ thống giả thuyết nghiên cứu:

  • Giả thuyết 1 (H1): Phát thải KNK trực tiếp từ các bể sinh học tuân theo phương trình động học tăng trưởng tế bào vi sinh vật Herbert mở rộng kết hợp với bảo toàn khối lượng COD, BOD và Nitơ chuyển hóa.
  • Giả thuyết 2 (H2): Phát thải gián tiếp từ tiêu thụ điện năng phục vụ sục khí, bơm và khuấy trộn chiếm tỷ trọng chi phối trong tổng lượng phát thải $\text{CO}_{2\text{tđ}}$ của các hệ thống xử lý sinh học hiếu khí đô thị tại Việt Nam do hệ số phát thải lưới điện quốc gia tương đối cao.
  • Giả thuyết 3 (H3): Sự hình thành $\text{N}_2\text{O}$ trong quá trình nitrat hóa và khử nitrat chịu sự kiểm soát phi tuyến tính mạnh mẽ bởi nồng độ oxy hòa tan (DO), nồng độ cơ chất và thời gian lưu bùn ($\text{SRT}$).

Phạm vi nghiên cứu bao quát các hệ thống XLNT sinh hoạt điển hình tại Hà Nội (Tòa nhà Goldmark City, 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 giấy quy mô lớn (Nhà máy Giấy Bãi Bằng, Nhà máy Giấy Corelex). Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp bộ công cụ phần mềm mô phỏng số chuyên dụng và xác lập bộ hệ số phát thải KNK cấp quốc gia (Tier 2/Tier 3), phục vụ trực tiếp cho Đóng góp do quốc gia tự quyết định (NDC) và Báo cáo kiểm kê KNK định kỳ của Việt Nam.

+-----------------------------------------------------------------------------------------------+
|                                 SƠ ĐỒ CẤU TRÚC TỔNG QUAN LUẬN ÁN                             |
+-----------------------------------------------------------------------------------------------+
|  BỐI CẢNH KHOA HỌC: BĐKH toàn cầu - Kiểm kê KNK lĩnh vực nước thải theo UNFCCC/IPCC          |
|  KHOẢNG TRỐNG: Thiếu mô hình động học vi sinh tính toán KNK đặc thù tại Việt Nam            |
+---------------------------------------------------------------+-------------------------------+
                                                                |
                                                                v
+-----------------------------------------------------------------------------------------------+
|  KHUNG LÝ THUYẾT & PHƯƠNG PHÁP:                                                              |
|  - Phương trình động học Herbert & Monod                                                      |
|  - Cân bằng vật chất C-N-P & Phương trình lượng hóa vi sinh                                  |
|  - Đo đạc thực địa: Buồng thu nổi Flux Chamber + Sắc ký khí GC                                |
+---------------------------------------------------------------+-------------------------------+
                                                                |
                                                                v
+-----------------------------------------------------------------------------------------------+
|  MÔ HÌNH HÓA & PHẦN MỀM:                                                                      |
|  - Hệ thống A (Hiếu khí), Hệ thống B (Yếm khí), Hệ thống C (Yếm khí - Hiếu khí kết hợp)       |
|  - Tính toán phát thải trực tiếp (CO2, CH4, N2O) & gián tiếp (Điện năng tiêu thụ)            |
+---------------------------------------------------------------+-------------------------------+
                                                                |
                                                                v
+-----------------------------------------------------------------------------------------------+
|  KIỂM NGHIỆM & ỨNG DỤNG THỰC TIỄN:                                                           |
|  - Kiểm nghiệm chuẩn: Hệ thống Windsor (Úc)                                                   |
|  - Thực nghiệm hiện trường: HTXLNT Goldmark City                                              |
|  - Ứng dụng quy mô lớn: Đô thị (Yên Sở, Trúc Bạch, Kim Liên) & Công nghiệp (Bãi Bằng, Corelex)|
+-----------------------------------------------------------------------------------------------+

Literature Review và Positioning

Các công trình nghiên cứu trên thế giới về phát thải KNK từ lĩnh vực xử lý nước thải đã trải qua nhiều giai đoạn phát triển với những bước tiến quan trọng về mặt lý thuyết và phương pháp luận:

                  TIẾN TRÌNH LÝ THUYẾT VỀ MÔ HÌNH HÓA KNK TỪ NƯỚC THẢI
========================================================================================
[1996 - 2006] IPCC Guidelines (Tier 1/2): Tiếp cận Top-down, sử dụng hệ số phát thải tĩnh
      │
      ▼
[1995 - 2000] Czepiel et al.; Henze et al. (ASM models): Cơ chế sinh học vi mô & cân bằng C-N
      │
      ▼
[2007] Bridle Consulting Model: Phân định 6 đường biên phát thải tĩnh trong HTXLNT
      │
      ▼
[2009 - 2018] Kampschreur et al.; Ni et al.: Cơ chế hình thành N2O trung gian trong nitrat hóa
      │
      ▼
[2019] Luận án Vân Anh (HUST): Mô hình hóa động học Herbert kết hợp điều kiện biên Việt Nam
========================================================================================

Dòng nghiên cứu thứ nhất dựa trên hướng dẫn của IPCC (1996, 2000, 2006, 2013). Phương pháp tiếp cận của IPCC mang tính vĩ mô (top-down), phân cấp từ Cấp 1 đến Cấp 3. Công thức cơ bản của IPCC xác định lượng phát thải methane từ nước thải sinh hoạt thông qua tổng lượng hữu cơ có khả năng phân hủy $\text{TOW} = P \times \text{BOD} \times 0{,}001 \times I \times 365$ và hệ số phát thải $\text{EF}_{i,j} = B_o \times \text{MCF}_j$, trong đó $B_o$ mặc định là $0{,}6\text{ kg CH}_4/\text{kg BOD}$ hoặc $0{,}25\text{ kg CH}_4/\text{kg COD}$. Đối với $\text{N}_2\text{O}$, IPCC (2006) áp dụng hệ số phát thải trực tiếp mặc định từ trạm xử lý tập trung là $3{,}2\text{ g N}_2\text{O}/\text{người/năm}$. Tuy nhiên, IPCC hoàn toàn bỏ qua phát thải $\text{CO}_2$ sinh học (biogenic $\text{CO}_2$) trong kiểm kê do giả định tính trung hòa carbon sinh khối ngắn hạn, đồng thời không tích hợp các thông số động học vận hành của từng loại bể phản ứng.

Dòng nghiên cứu thứ hai tập trung vào mô hình hóa biên xử lý công nghệ, tiêu biểu là mô hình của Bridle Consulting (2007). Bridle đã xây dựng mô hình tính toán phát thải KNK tổng thể cho trạm xử lý nước thải sinh hoạt dựa trên 6 ranh giới phát thải: (1) Quá trình xử lý sinh học (oxi hóa $\text{BOD}_5$, hô hấp nội bào, nitrat hóa, khử nitơ sinh $\text{CO}_2$ và $\text{N}_2\text{O}$); (2) Tiêu thụ năng lượng điện cho sục khí, khuấy và bơm; (3) Phân hủy bùn kỵ khí tạo $\text{CH}_4$ và $\text{CO}_2$; (4) Xử lý và vận chuyển bùn thải; (5) Thu hồi và sử dụng khí sinh học; (6) Sử dụng hóa chất trong keo tụ. Dù mô hình Bridle đã chi tiết hóa các nguồn thải nhưng vẫn mang tính chất tĩnh, các hệ số chuyển hóa chủ yếu dựa trên số liệu vận hành cố định tại các nước ôn đới.

Dòng nghiên cứu thứ ba đi sâu vào cơ chế vi sinh động học bậc cao, như các nghiên cứu của Czepiel et al. (1995), Kampschreur et al. (2009), và nhóm nghiên cứu mô hình bùn hoạt tính IWA (ASM1, ASM2d, ASM3) của Henze et al. (2000). Czepiel et al. (1995) đã chứng minh bể hiếu khí là nguồn phát thải $\text{N}_2\text{O}$ trực tiếp lớn nhất trong trạm xử lý nước thải đô thị, trung bình đạt $3{,}2\text{ g N}_2\text{O}/\text{người/năm}$. Kampschreur et al. (2009) chỉ ra rằng $\text{N}_2\text{O}$ phát sinh từ hai con đường chính: sản phẩm phụ của quá trình oxy hóa hydroxylamine ($\text{NH}_2\text{OH}$) thành $\text{NO}_2^-$ do vi khuẩn tự dưỡng Nitrosomonas khi thiếu DO, và quá trình khử nitrat không hoàn toàn của vi khuẩn dị dưỡng trong điều kiện ức chế oxy hoặc thiếu hụt nguồn carbon hữu cơ dễ phân hủy.

Trong bối cảnh học thuật quốc tế tồn tại sự tranh luận gay gắt giữa hai luồng quan điểm:

  • Quan điểm 1 (IPCC, Doorn et al., 2006): Cho rằng phát thải trực tiếp $\text{N}_2\text{O}$ và $\text{CO}_2$ từ trạm xử lý chỉ là nguồn thứ yếu, có thể bỏ qua hoặc áp dụng hệ số mặc định cố định theo đầu người nhằm đơn giản hóa công tác kiểm kê quốc gia.
  • Quan điểm 2 (Kampschreur et al., 2009; Ni et al., 2013; Bridle, 2007): Khẳng định phát thải KNK biến động rất mạnh theo chế độ thủy lực, cấu trúc dòng khí sục và tuổi bùn. Việc áp dụng hệ số tĩnh sẽ gây sai lệch từ $50%$ đến hơn $300%$ tổng lượng phát thải thực tế của trạm.

Công trình của Nguyễn Thị Vân Anh (2019) đã định vị chính xác vào giao điểm của hai luồng quan điểm trên: kế thừa tính bao quát theo ranh giới công nghệ của mô hình Bridle (2007), đồng thời nhúng lý thuyết động học vi sinh Herbert vào mô hình toán học để lượng hóa phát thải động. Tác giả so sánh trực tiếp kết quả với hệ thống chuẩn quốc tế Windsor (Australia) và các chỉ số mặc định của IPCC (2006), qua đó chứng minh rằng mô hình mới giải quyết triệt để sự thiếu hụt dữ liệu đầu vào đặc thù của Việt Nam mà vẫn đảm bảo độ chính xác cơ chế sinh học.

Đó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 nền tảng lý thuyết tăng trưởng và phân hủy sinh khối của Herbert (Herbert’s kinetic formulation) trong kỹ thuật môi trường, vốn trước đây chỉ dùng để tính toán nồng độ bùn và chất nền đầu ra, nay được chuyển hóa thành công cụ định lượng phát thải KNK đa thành phần ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$).

Cơ chế toán học được xây dựng dựa trên phương trình động học tăng trưởng thực của sinh khối: $$r_x = \mu X - k_d X = \left(\frac{\mu_m S}{K_s + S}\right) X - k_d X$$ và tốc độ sử dụng cơ chất: $$r_s = -\frac{k S X}{K_s + S} = -\frac{\mu_m S}{Y(K_s + S)} X$$

Trong đó:

  • $\mu_m$: Tốc độ tăng trưởng riêng cực đại của vi sinh vật ($\text{ngày}^{-1}$).
  • $k$: Tốc độ sử dụng cơ chất riêng cực đại ($\text{mg/mg}\cdot\text{ngày}$).
  • $K_s$: Hằng số bán bão hòa đối với cơ chất ($\text{mg/L}$ hoặc $\text{g/m}^3$).
  • $k_d$: Hệ số phân hủy nội bào ($\text{ngày}^{-1}$).
  • $Y$: Hệ số sản lượng sinh khối cực đại ($\text{mg VSS/mg BOD}$ hoặc $\text{mg VSS/mg COD}$).
  • $X$: Nồng độ sinh khối hoạt tính dị dưỡng trong bể ($\text{mg/L}$).
  • $S$: Nồng độ cơ chất hòa tan trong bể phản ứng ($\text{mg/L}$).
+-----------------------------------------------------------------------------------------------+
|                       CƠ CHẾ ĐỘNG HỌC PHÁT THẢI KHÍ NHÀ KÍNH TRONG BỂ HIẾU KHÍ                |
+-----------------------------------------------------------------------------------------------+
|                                                                                               |
|   Cơ chất dòng vào (S0, Q) ──> [ BỂ HIẾU KHÍ ] ──> Cơ chất dòng ra (S, Q)                    |
|                                     │                                                         |
|         ┌───────────────────────────┼───────────────────────────┐                             |
|         ▼                           ▼                           ▼                             |
|   1. OXY HÓA CƠ CHẤT         2. PHÂN HỦY NỘI BÀO         3. NITRAT HÓA & PHÁT THẢI N2O        |
|   S + O2 ──> CO2 + H2O       Tế bào (X) ──> CO2 + NH3    NH4+ ──> NH2OH ──> NO2- ──> NO3-     |
|   (M_CO2,oxi)                (M_CO2,phan_huy)                       │ (Sản phẩm phụ: N2O)     |
|         │                           │                               ▼                         |
|         └───────────────────────────┴──────────────────────> [ M_CO2,be_xl + M_N2O,be_xl ]   |
|                                                                                               |
+-----------------------------------------------------------------------------------------------+

Mô hình lý thuyết đã thiết lập các mệnh đề toán học định lượng:

Mệnh đề 1: Phát sinh $\text{CO}_2$ từ oxy hóa chất hữu cơ và hô hấp nội bào trong bể hiếu khí: $$M_{\text{CO}2,\text{bexl}} = M{\text{CO}2,\text{oxi}} + M{\text{CO}2,\text{phanhuy}}$$ $$M{\text{CO}2,\text{oxi}} = Q_v (S_v - S) \cdot Y{\text{CO}2/\text{BOD}}$$ $$M{\text{CO}2,\text{phanhuy}} = V \cdot X \cdot k_d \cdot Y{\text{CO}_2/\text{VSS}}$$

Mệnh đề 2: Phát sinh $\text{CH}_4$ từ quá trình phân hủy bùn kỵ khí: $$M_{\text{CH}4,\text{bephanhuy}} = P_X \cdot (1 - f_d) \cdot Y{\text{CH}4/\text{VSS}} \cdot \eta{\text{yemkhi}}$$ Trong đó $P_X$ là lượng bùn dư sinh ra hàng ngày ($\text{kg/ngày}$), $f_d$ là phần sinh khối trơ từ mảnh vỡ tế bào ($0{,}10 - 0{,}15$), và $\eta_{\text{yemkhi}}$ là hiệu suất chuyển hóa kỵ khí.

Mệnh đề 3: Phát thải $\text{CO}_2$ gián tiếp từ tiêu thụ năng lượng điện: $$M_{\text{CO}2,\text{dien}} = Q_E \times \sum{i} (PF_i \times EF_i)$$ Trong đó $Q_E$ là tổng công suất tiêu thụ điện năng của trạm ($\text{kWh/ngày}$), $PF_i$ là tỷ trọng nguồn phát điện loại $i$ trong cơ cấu lưới điện quốc gia ($%$), và $EF_i$ là hệ số phát thải của nguồn điện tương ứng ($\text{kg CO}_{2\text{tđ}}/\text{kWh}$).

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

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết cốt lõi:

  1. Lý thuyết nhiệt động học và phương trình lượng hóa sinh học: Thiết lập tỷ lượng mol phản ứng hóa sinh chuẩn cho từng loại nước thải (nước thải sinh hoạt và nước thải sản xuất bột giấy/giấy tái chế).
  2. Lý thuyết động học sinh trưởng vi sinh vật: Tích hợp các hệ số $k, K_s, Y, k_d$ biến thiên theo nhiệt độ môi trường nhiệt đới dựa trên phương trình Arrhenius: $\theta^{(T - 20)}$.
  3. Lý thuyết cân bằng vật chất bảo toàn ranh giới hệ thống: Phân chia ranh giới thành ba cấu hình công nghệ điển hình:
    • Hệ thống A: Cấu hình hiếu khí hoàn toàn (Bể lắng sơ cấp $\rightarrow$ Bể Aerotank $\rightarrow$ Bể lắng thứ cấp $\rightarrow$ Bể phân hủy bùn kỵ khí).
    • Hệ thống B: Cấu hình yếm khí triệt để (Bể UASB/Bể phản ứng kỵ khí tiếp xúc $\rightarrow$ Xử lý bùn).
    • Hệ thống C: Cấu hình kết hợp Yếm khí - Thiếu khí - Hiếu khí ($\text{A}_2\text{O}$ / Anoxic-Aerobic).
Thông số động học / Hệ số chuyển đổi Ý nghĩa sinh học Giá trị nước thải sinh hoạt Giá trị nước thải công nghiệp giấy
$k$ ($\text{mg BOD/mg VSS}\cdot\text{ngày}$) Tốc độ sử dụng chất nền cực đại $2{,}0 - 5{,}0$ $1{,}2 - 3{,}5$
$K_s$ ($\text{mg/L}$) Hằng số bán bão hòa cơ chất $40 - 100$ $80 - 180$
$Y$ ($\text{mg VSS/mg BOD}$) Hệ số sản lượng sinh khối $0{,}4 - 0{,}6$ $0{,}3 - 0{,}5$
$k_d$ ($\text{ngày}^{-1}$) Hệ số phân hủy nội bào $0{,}04 - 0{,}08$ $0{,}02 - 0{,}05$
$Y_{\text{CO}_2/\text{BOD}}$ ($\text{kg CO}_2/\text{kg BOD}$) Hệ số phát sinh $\text{CO}_2$ do oxy hóa $0{,}95 - 1{,}10$ $1{,}15 - 1{,}35$
$Y_{\text{CH}_4/\text{VSS}}$ ($\text{kg CH}_4/\text{kg VSS}$) Hệ số sinh $\text{CH}_4$ kỵ khí $0{,}25 - 0{,}35$ $0{,}20 - 0{,}30$

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan khoa học thực chứng (positivism) kết hợp với chủ nghĩa duy thực phê phán (critical realism), khẳng định các quy luật phát sinh KNK là các quá trình hóa sinh khách quan có thể lượng hóa chính xác bằng mô hình toán và thực nghiệm kiểm chứng.

Thiết kế nghiên cứu đa tầng bậc (multi-level design) được triển khai:

  • Tầng vi mô (Micro-level): Động học chuyển hóa cơ chất và enzym vi sinh vật (Nitrosomonas, Nitrobacter, vi khuẩn sinh methane kỵ khí).
  • Tầng trung mô (Meso-level): Cân bằng dòng vật chất và thủy lực trong các đơn vị bể phản ứng cụ thể (CSTR, bể lắng, bể phân hủy bùn).
  • Tầng vĩ mô (Macro-level): Lượng hóa phát thải KNK toàn diện của toàn bộ nhà máy, bao gồm cả phát thải gián tiếp từ cơ cấu nguồn điện lưới quốc gia Việt Nam (thủy điện, nhiệt điện than, nhiệt điện khí, năng lượng tái tạo).
+-----------------------------------------------------------------------------------------------+
|                             THIẾT KẾ NGHIÊN CỨU ĐA TẦNG CỦA LUẬN ÁN                           |
+-----------------------------------------------------------------------------------------------+
|                                                                                               |
|  [ TẦNG VĨ MÔ ]   ──> Tổng kiểm kê toàn nhà máy (Trực tiếp + Gián tiếp qua lưới điện quốc gia)|
|         │                                                                                     |
|         ▼                                                                                     |
|  [ TẦNG TRUNG MÔ ]──> Cân bằng vật chất C-N-P từng cụm bể (Hiếu khí, Yếm khí, Bể lắng, Bùn)   |
|         │                                                                                     |
|         ▼                                                                                     |
|  [ TẦNG VI MÔ ]   ──> Động học enzym & vi sinh vật (Herbert-Monod: k, Ks, Y, kd, Nitrosomonas)|
|                                                                                               |
+-----------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu kết hợp giữa đo đạc thực nghiệm thực địa và mô phỏng số học:

                            QUY TRÌNH THỰC NGHIỆM VÀ MÔ HÌNH HÓA
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  1. THU MẪU KHÍ THỰC ĐỊA:                                                                   |
|     - Thiết bị: Buồng thu nổi Flux Chamber (vật liệu composite, kín khí, diện tích chuẩn)   |
|     - Vị trí: Bề mặt bể Aerotank, bể lắng, bể chứa bùn tại HTXLNT Goldmark City             |
|     - Thời gian: Lấy mẫu định kỳ theo chu kỳ sục khí và các khung giờ vận hành khác nhau     |
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  2. PHÂN TÍCH KHÍ NHÀ KÍNH BẰNG SẮC KÝ KHÍ (GC):                                            |
|     - Thiết bị: Sắc ký khí tích hợp đầu dò GC-FID (CH4, CO2) & GC-ECD (N2O)                 |
|     - Khí mang: He/N2 siêu tinh khiết; kiểm định đường chuẩn bằng mẫu khí chuẩn quốc tế     |
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  3. XÂY DỰNG MÃ NGUỒN SỐ & GIAO DIỆN PHẦN MỀM:                                              |
|     - Lập trình thuật toán mô phỏng cân bằng khối lượng và phương trình vi phân động học    |
|     - Thiết kế giao diện GUI thân thiện cho phép nhập thông số công nghệ & đặc tính nước    |
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  4. HIỆU CHỈNH & KIỂM NGHIỆM MÔ HÌNH:                                                       |
|     - Bước 1: Kiểm nghiệm với số liệu chuẩn của hệ thống Windsor (Bridle, 2007)             |
|     - Bước 2: Kiểm nghiệm đối sánh với kết quả đo đạc thực tế tại HTXLNT Goldmark City      |
+─────────────────────────────────────────────────────────────────────────────────────────────+

Độ tin cậy và giá trị khoa học được bảo đảm qua nguyên tắc tam giác giác hóa (triangulation):

  • Tam giác hóa phương pháp: Đối chiếu chéo kết quả giữa tính toán theo phương pháp IPCC, mô hình Bridle, mô hình động học thiết lập của luận án và số liệu đo trực tiếp bằng sắc ký khí.
  • Tam giác hóa dữ liệu: Kiểm nghiệm trên nhiều đối tượng đa dạng từ nước thải sinh hoạt khu đô thị cao tầng (Goldmark City) đến hệ thống xử lý nước thải đô thị tập trung (Yên Sở công suất $200.000\text{ m}^3/\text{ngày}$, Trúc Bạch, Kim Liên) và các tổ hợp công nghiệp giấy quy mô lớn (Nhà máy Giấy Bãi Bằng $25.000\text{ m}^3/\text{ngày}$, Giấy Corelex).

Data và phân tích

Số liệu đầu vào bao gồm:

  1. Thông số chất lượng nước: $\text{BOD}_5$, $\text{COD}$, $\text{SS}$, Tổng Nitơ ($\text{TN}$), $\text{NH}_4^+$, $\text{NO}_3^-$, Tổng Phốt pho ($\text{TP}$) tại các điểm: Nước thải thô vào, sau lắng sơ cấp, trong bể sinh học, và dòng ra sau lắng thứ cấp.
  2. Thông số công nghệ và vận hành: Lưu lượng dòng vào ($Q$), tỷ lệ tuần hoàn bùn ($Q_T/Q$), lưu lượng xả bùn dư ($Q_x$), thể tích bể ($V$), thời gian lưu nước ($\text{HRT}$), tuổi bùn ($\text{SRT}$), lượng điện năng tiêu thụ ($Q_E$).
  3. Phân tích độ nhạy (Sensitivity Analysis): Đánh giá mức độ ảnh hưởng của sự biến thiên các thông số động học ($k, K_s, Y, k_d$), nhiệt độ môi trường ($15^\circ\text{C} - 35^\circ\text{C}$), và nồng độ chất nền đầu vào đối với lượng phát thải từng loại khí.

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

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

Phát hiện 1: Cơ cấu phát thải KNK bị chi phối áp đảo bởi điện năng tiêu thụ gián tiếp trong các hệ thống hiếu khí đô thị. Tại hệ thống xử lý nước thải sinh hoạt tòa nhà Goldmark City, phát thải KNK từ tiêu thụ điện năng chiếm tới $68{,}4% - 72{,}1%$ tổng lượng phát thải $\text{CO}_{2\text{tđ}}$ toàn hệ thống. Điều này chứng minh rằng hiệu quả năng lượng của máy thổi khí và bơm là chìa khóa then chốt để giảm thiểu dấu chân carbon của các trạm XLNT sinh hoạt tại Việt Nam.

Phát hiện 2: Đóng góp tiềm tàng khổng lồ của $\text{N}_2\text{O}$ dù nồng độ phát thải thể tích rất nhỏ. Tại bể Aerotank của hệ thống Goldmark City, nồng độ $\text{N}_2\text{O}$ đo đạc qua sắc ký khí dao động từ $0{,}015$ đến $0{,}042\text{ mg/L}$ trong pha khí thoát ra bề mặt. Tuy nhiên, do có hệ số gây ấm toàn cầu $\text{GWP} = 296$, phát thải $\text{N}_2\text{O}$ đóng góp tới $18{,}2%$ tổng phát thải KNK trực tiếp từ các quá trình sinh học, khẳng định việc kiểm soát nitrat hóa/khử nitrat không hoàn toàn có ý nghĩa sống còn trong việc cắt giảm KNK.

Phát hiện 3: Ảnh hưởng phi tuyến tính của nhiệt độ và tuổi bùn ($\text{SRT}$) đến phát thải $\text{CO}_2$ sinh học. Khi nhiệt độ tăng từ $20^\circ\text{C}$ lên $30^\circ\text{C}$ (điều kiện mùa hè đặc trưng tại miền Bắc Việt Nam), hệ số phân hủy nội bào $k_d$ tăng từ $0{,}05\text{ ngày}^{-1}$ lên $0{,}082\text{ ngày}^{-1}$, dẫn đến lượng $\text{CO}_2$ phát sinh do hô hấp nội bào tăng $38{,}5%$, trong khi lượng bùn hoạt tính dư sinh ra giảm $19{,}4%$. Khi tăng tuổi bùn $\text{SRT}$ từ $5$ ngày lên $15$ ngày, lượng $\text{CO}_2$ sinh học tăng dần nhưng phát thải $\text{N}_2\text{O}$ giảm mạnh do quá trình nitrat hóa diễn ra hoàn toàn hơn.

Phát hiện 4: Sự khác biệt bản chất giữa nước thải sinh hoạt và nước thải công nghiệp giấy. Đối với hệ thống xử lý nước thải Nhà máy Giấy Bãi Bằng, tỷ lệ phát thải $\text{CH}_4$ từ các bể xử lý kỵ khí và hồ lắng bùn chiếm tỷ trọng cao hơn đáng kể ($42{,}6%$ tổng phát thải trực tiếp) so với nước thải sinh hoạt do hàm lượng lignocellulose và hợp chất hữu cơ khó phân hủy cao, tạo điều kiện thuận lợi cho các vi khuẩn sinh methane kỵ khí hoạt động kéo dài trong lớp bùn lắng.

+-----------------------------------------------------------------------------------------------+
|                       CƠ CẤU PHÁT THẢI KHÍ NHÀ KÍNH ĐẶC TRƯNG THEO CÔNG NGHỆ                  |
+-----------------------------------------------------------------------------------------------+
|                                                                                               |
|  [ NƯỚC THẢI SINH HOẠT - GOLDMARK CITY ]                                                      |
|  ┌────────────────────────────────────────────────────────┬───────────────────┬───────────┐   |
|  │ Phát thải gián tiếp từ điện năng tiêu thụ: ~70%        │ CO2 sinh học: ~18%│ N2O: ~12% │   |
|  └────────────────────────────────────────────────────────┴───────────────────┴───────────┘   |
|                                                                                               |
|  [ NƯỚC THẢI CÔNG NGHIỆP GIẤY - BÃI BẰNG ]                                                   |
|  ┌──────────────────────────┬───────────────────────────────┬─────────────────┬───────────┐   |
|  │ Điện năng: ~45%          │ CH4 từ xử lý yếm khí/bùn: ~38%│ CO2 sinh học:14%│ N2O: ~3%  │   |
|  └──────────────────────────┴───────────────────────────────┴─────────────────┴───────────┘   |
|                                                                                               |
+-----------------------------------------------------------------------------------------------+

Phát hiện 5: Độ chính xác vượt trội của mô hình thiết lập so với phương pháp IPCC Tier 1. Khi kiểm nghiệm với hệ thống thực tế Goldmark City, kết quả tính toán phát thải từ mô hình của luận án có độ lệch so với số liệu đo đạc thực nghiệm bằng sắc ký khí dao động trong khoảng $4{,}8% - 9{,}2%$. Trong khi đó, phương pháp tính toán theo hệ số mặc định của IPCC (2006) cho sai số vượt quá $+65%$ đối với $\text{CH}_4$ và $-48%$ đối với $\text{N}_2\text{O}$.

Implications đa chiều

                                  HỆ THỐNG IMPLICATIONS ĐA CHIỀU
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  1. ĐÓNG GÓP HỌC THUẬT (THEORETICAL ADVANCES):                                              |
|     - Bổ sung bộ thông số động học KNK phù hợp với điều kiện khí hậu nhiệt đới              |
|     - Cung cấp mô hình toán học giải thích sự tương tác giữa hô hấp nội bào và tạo khí     |
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  2. ĐỔI MỚI PHƯƠNG PHÁP LUẬN (METHODOLOGICAL INNOVATIONS):                                  |
|     - Chuyển đổi kiểm kê từ "hộp đen mặc định" (IPCC Tier 1) sang mô phỏng động (Tier 2/3)  |
|     - Quy trình kết hợp Buồng thu nổi + Sắc ký khí chuẩn hóa cho kiểm chứng mô hình        |
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  3. ỨNG DỤNG THỰC TIỄN & VẬN HÀNH (PRACTICAL APPLICATIONS):                                 |
|     - Tối ưu hóa chế độ sục khí DO và tuổi bùn SRT để giảm đồng thời điện năng và N2O       |
|     - Tận dụng triệt để khí sinh học (Biogas/CH4) trong công nghiệp giấy nhằm bù trừ năng lượng|
+───────────────────────────────────────────────┬─────────────────────────────────────────────+
                                                │
                                                ▼
+─────────────────────────────────────────────────────────────────────────────────────────────+
|  4. KHUYẾN NGHỊ CHÍNH SÁCH QUỐC GIA (POLICY RECOMMENDATIONS):                               |
|     - Xây dựng hệ số phát thải đặc thù quốc gia phục vụ Báo cáo Đóng góp NDC gửi UNFCCC     |
|     - Ban hành tiêu chuẩn kỹ thuật về kiểm toán dấu chân carbon trong thiết kế HTXLNT mới   |
+─────────────────────────────────────────────────────────────────────────────────────────────+

Limitations và Future Research

Nhằm bảo đảm tính khách quan và liêm chính học thuật, các giới hạn nội tại của luận án cần được làm rõ:

  1. Giới hạn về mô hình hóa cơ chế vi sinh $\text{N}_2\text{O}$: Mô hình luận án tiếp cận $\text{N}_2\text{O}$ thông qua hệ số chuyển đổi tỷ lệ dựa trên nồng độ nitơ dòng vào và hiệu suất nitrat hóa, chưa tích hợp chi tiết động học enzym nội bào của hai nhánh vi khuẩn oxy hóa amoniac (AOB) và vi khuẩn khử nitrat dị dưỡng (heterotrophic denitrifiers) như trong các mô hình sinh học phân tử chuyên sâu (ví dụ mô hình hai bước AOB của Hiệp hội Nước Quốc tế IWA).
  2. Giới hạn ranh giới đánh giá vòng đời (LCA Boundary): Mô hình mới chỉ tập trung vào phát thải trực tiếp tại các bể sinh học và phát thải gián tiếp từ tiêu thụ điện năng vận hành. Phát thải gián tiếp thượng nguồn từ quá trình sản xuất hóa chất keo tụ (PAC, Phèn, Polymer), phát thải trong quá trình vận chuyển bùn đến bãi chôn lấp tập trung chưa được tích hợp đầy đủ.
  3. Giả định đồng nhất thủy lực: Mô hình giả định các bể phản ứng sinh học hoạt động theo chế độ khuấy trộn hoàn chỉnh (CSTR) lý tưởng, chưa tính đến sự phân tầng nồng độ oxy hòa tan và hiện tượng dòng chảy tắt (dead zones) thường xảy ra trong các bể Aerotank kích thước lớn tại thực tế Việt Nam.

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Hướng 1: Mở rộng mô hình toán học tích hợp động lực học chất lưu tính toán (CFD - Computational Fluid Dynamics) với mô hình động học sinh hóa để mô phỏng sự phân bố không gian 3D của quá trình phát thải KNK trong các bể phản ứng quy mô lớn.
  • Hướng 2: Tích hợp mô hình vào hệ thống điều khiển tự động SCADA/IoT theo thời gian thực (Real-time Carbon Footprint Monitoring), cho phép tự động điều chỉnh lưu lượng khí cấp dựa trên cảm biến online để tối thiểu hóa phát thải KNK.
  • Hướng 3: Mở rộng nghiên cứu sang các loại hình nước thải công nghiệp phức tạp khác như nước thải chế biến cao su, rỉ rác, dệt nhuộm và chế biến thủy hải sản.
  • Hướng 4: Đánh giá toàn diện vòng đời (LCA) kết hợp phân tích chi phí - lợi ích kinh tế kỹ thuật (Techno-Economic Analysis) cho các giải pháp thu hồi biogas và đồng phát nhiệt - điện (CHP) từ bùn thải.

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

                                  MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG
=================================================================================================
 LĨNH VỰC                      TÁC ĐỘNG ĐỊNH TÍNH                          CHỈ SỐ ĐỊNH LƯỢNG
─────────────────────────────────────────────────────────────────────────────────────────────────
 Học thuật           - Xác lập nền tảng mô phỏng KNK nước thải     - Đóng góp bài báo khoa học
 (Academic)            nhiệt đới đầu tiên tại Việt Nam               chuyên ngành, trích dẫn cao
                     - Đào tạo thế hệ nghiên cứu sinh chuyên sâu

 Công nghiệp         - Giúp doanh nghiệp tối ưu chi phí điện năng  - Cắt giảm 15-25% điện tiêu thụ
 (Industry)          - Kiểm soát rò rỉ khí Biogas kỵ khí           - Giảm 20-30% tổng lượng CO2tđ

 Chính sách          - Cung cấp hệ số Tier 2/3 chuẩn xác           - Phục vụ báo cáo UNFCCC,
 (Policy)            - Căn cứ xây dựng thị trường tín chỉ carbon    hướng tới Net Zero 2050

 Xã hội & Môi trường - Cải thiện chất lượng môi trường không khí    - Giảm thiểu khí mùi độc hại
 (Societal/Global)   - Hạn chế phát sinh khí nhà kính gây ấm       - Tối ưu hóa chu trình C-N
=================================================================================================

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

  • Nghiên cứu sinh và Giới học thuật (Academic Researchers): Tiếp cận được một khung lý thuyết hoàn chỉnh kết hợp giữa động học Herbert, cân bằng vật chất và mô hình hóa KNK; thừa hưởng bộ mã nguồn và phương pháp luận để mở rộng nghiên cứu sang các đối tượng nước thải đặc thù khác.
  • Kỹ sư vận hành và Doanh nghiệp quản lý trạm XLNT (Industry Operators): Sử dụng phần mềm được xây dựng để mô phỏng dự báo trước tải lượng phát thải KNK, tối ưu hóa quá trình sục khí, giảm hóa đơn tiền điện vận hành từ $15% - 25%$ và hạn chế phát thải khí mùi độc hại.
  • Các Cơ quan Quản lý Nhà nước và Hoạch định Chính sách (Policy Makers - Bộ TN&MT): Sở hữu bộ cơ sở dữ liệu và hệ số phát thải KNK cấp quốc gia có độ tin cậy cao, thay thế cho các hệ số mặc định của IPCC, phục vụ trực tiếp công tác đàm phán biến đổi khí hậu quốc tế và thẩm định hạn ngạch phát thải KNK trong tương lai.

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 độc đáo nhất là việc tích hợp thành công phương trình động học tăng trưởng và phân hủy sinh khối nội bào của Herbert với định luật bảo toàn khối lượng các nguyên tố $\text{C}, \text{N}$ để thiết lập hệ phương trình vi phân mô tả tốc độ phát sinh đồng thời ba loại khí nhà kính ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$) trong điều kiện vận hành động học không ổn định (unsteady-state). Luận án đã mở rộng lý thuyết Herbert truyền thống từ một công cụ tính toán nồng độ bùn lắng trở thành một khung lý thuyết định lượng dấu chân carbon vi sinh vật học.

2. Sự đổi mới về mặt 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 phương pháp kiểm kê của IPCC (2006) vốn chỉ là mô hình tĩnh (black-box) dựa trên hệ số mặc định vĩ mô, luận án đã xây dựng mô hình cơ chế vi sinh động học chi tiết. So với mô hình 6 ranh giới của Bridle Consulting (2007) vốn áp dụng các hệ số chuyển hóa tĩnh tại các nước ôn đới, luận án đã bản địa hóa và tích hợp phương trình phụ thuộc nhiệt độ Arrhenius và cơ cấu phát thải lưới điện đặc thù của Việt Nam, đồng thời kiểm chứng nghiêm ngặt bằng sắc ký khí thực địa.

3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là tại hệ thống xử lý nước thải sinh hoạt đô thị hiện đại (Goldmark City), lượng $\text{CO}2$ gián tiếp sinh ra do tiêu thụ điện năng cho máy thổi khí sục khí chiếm tới hơn $70%$ tổng lượng $\text{CO}{2\text{tđ}}$, vượt trội hoàn toàn so với tổng phát thải trực tiếp từ các phản ứng hóa sinh tại các bể phản ứng. Điều này chỉ ra rằng chiến lược giảm phát thải KNK hiệu quả nhất trong xử lý nước thải đô thị tại Việt Nam không phải là thay đổi công nghệ sinh học mà là nâng cao hiệu suất sử dụng năng lượng của hệ thống cấp khí.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không?

Có. Luận án trình bày chi tiết toàn bộ thuật toán tính toán, các phương trình lượng hóa cân bằng vật chất cho từng loại bể, danh mục đầy đủ các hệ số động học ($k, K_s, Y, k_d$), bản vẽ cấu tạo và quy trình vận hành thiết bị buồng thu nổi (Flux chamber), quy trình phân tích sắc ký khí đạt chuẩn ISO, cũng như cấu trúc sơ đồ khối phần mềm, cho phép các nhà nghiên cứu khác dễ dàng tái lập và hiệu chỉnh mô hình cho các trạm xử lý khác.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào: (1) Số hóa và nhúng mô hình vào các nền tảng điều khiển thông minh IoT/AI để quản trị phát thải KNK thời gian thực; (2) Tích hợp mô hình động học với mô phỏng thủy lực CFD 3D; (3) Mở rộng kiểm kê phát thải KNK sang toàn bộ các ngành công nghiệp phát thải lớn tại Việt Nam (hóa chất, dệt nhuộm, chế biến thực phẩm); (4) Xây dựng sàn giao dịch tín chỉ carbon từ các dự án tối ưu hóa trạm xử lý nước thải.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thị Vân Anh đã đạt được các đóng góp khoa học và thực tiễn mang tính đột phá:

  1. Thiết lập thành công mô hình toán học động học đầu tiên tại Việt Nam có khả năng tính toán định lượng chính xác phát thải KNK ($\text{CO}_2, \text{CH}_4, \text{N}_2\text{O}$) cho cả nguồn trực tiếp và gián tiếp trong các hệ thống xử lý nước thải sinh hoạt và công nghiệp giấy.
  2. Xây dựng bộ công cụ phần mềm chuyên dụng với giao diện đồ họa thân thiện, hỗ trợ đắc lực cho các kỹ sư và nhà quản lý trong việc dự báo và tối ưu hóa phát thải KNK.
  3. Xác lập bộ thông số động học vi sinh và hệ số phát thải KNK đặc thù cho điều kiện nhiệt đới tại Việt Nam, thu hẹp khoảng cách sai số từ hơn $65%$ của phương pháp IPCC Tier 1 xuống dưới $10%$.
  4. Chứng minh vai trò chi phối của điện năng tiêu thụ trong cơ cấu phát thải của trạm xử lý hiếu khí và đề xuất các giải pháp kỹ thuật tối ưu hóa $\text{DO} - \text{SRT}$ để giảm đồng thời năng lượng và phát sinh $\text{N}_2\text{O}$.
  5. Đóng góp cơ sở dữ liệu khoa học chuẩn mực phục vụ trực tiếp cho công tác kiểm kê KNK quốc gia và thực thi các cam kết giảm phát thải của Việt Nam theo Hiệp định Paris về Biến đổi Khí hậu.

Công trình không chỉ mở ra hướng nghiên cứu liên ngành kết hợp giữa kỹ thuật môi trường, toán tin học và khoa học khí hậu, mà còn để lại một di sản học thuật quan trọng, góp phần thúc đẩy lộ trình chuyển dịch xanh và phát triển bền vững của quốc gia.