Tổng quan về luận án

Hoạt động đốt hở phụ phẩm nông nghiệp sau thu hoạch, đặc biệt là rơm rạ trên đồng ruộng, là một tập quán canh tác cố hữu tại hầu hết các quốc gia trồng trúa nước ở châu Á, đóng vai trò là một trong những nguồn phát thải nhân tạo chủ đạo làm suy thoái chất lượng không khí từ quy mô địa phương đến toàn cầu. Luận án tiến sĩ kỹ thuật môi trường mang tên "Đánh giá mức độ phát thải của hoạt động đốt rơm rạ và khả năng tác động của chúng đến chất lượng không khí – nghiên cứu thí điểm tại đồng bằng Tây Nam Bộ" do Nghiên cứu sinh Phạm Thị Hồng Phương thực hiện dưới sự hướng dẫn khoa học của PGS. Nghiêm Trung Dũng và PGS. Phạm Thị Mai Thảo tại Trường Đại học Bách Khoa Hà Nội (Mã số chuyên ngành: 9520320, năm bảo vệ: 2022) đại diện cho công trình khoa học thực nghiệm tiên phong tại Việt Nam giải quyết triệt để bài toán định lượng và kiểm kê phát thải đốt sinh khối nông nghiệp tại Đồng bằng sông Cửu Long (ĐBSCL).

Khoảng trống nghiên cứu (research gap) then chốt được xác định là: Trong khi các công trình kiểm kê và đánh giá phát thải trước đây tại Việt Nam tập trung hầu hết ở khu vực miền Bắc và Đồng bằng sông Hồng (ĐBSH) với tập quán đốt đống (pile burning) trong điều kiện khí hậu mùa đông lạnh (Lasko et al., 2018; H. Le et al., 2020), thì vùng ĐBSCL – vựa lúa lớn nhất cả nước với sản lượng lúa hàng chục triệu tấn mỗi năm và tỷ lệ đốt rơm rạ lên đến 80–90%, thậm chí đạt 98% vào vụ Đông Xuân (World Bank, 2014; FAO, 2018) – lại hoàn toàn thiếu vắng một bộ hệ số phát thải (Emission Factors - EFs) đặc trưng cấp địa phương (Tier 2 theo hướng dẫn IPCC) cho phương thức đốt rải (spread burning) trên đồng ruộng cơ giới hóa.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Bộ hệ số phát thải thực nghiệm (EF) của các nhóm chất ô nhiễm đại diện gồm bụi (TSP, PM10, PM2.5), 16 hydrocacbon thơm đa vòng (PAHs), 10 hợp chất hữu cơ bay hơi (VOCs), và khí nhà kính/khí ô nhiễm (CO2, SO2, NO2) từ hoạt động đốt rải rơm rạ tại ĐBSCL có giá trị định lượng là bao nhiêu?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ đóng góp ròng và tải lượng phát thải hàng năm của hoạt động đốt rơm rạ theo mùa vụ tại 13 tỉnh thành ĐBSCL diễn biến như thế nào?
  • Câu hỏi nghiên cứu 3 (RQ3): Khả năng tác động tức thời đến chất lượng không khí xung quanh và tiềm năng hình thành ozon tầng đối lưu (OFP) của khói thải đốt rơm rạ biểu hiện ra sao so với nồng độ nền?
  • Câu hỏi nghiên cứu 4 (RQ4): Những tỷ lệ chẩn đoán (diagnostic ratios) phân tử PAHs nào có thể được xác lập làm chất đánh dấu đặc trưng (molecular markers) để phân biệt nguồn đốt rơm rạ với khí thải giao thông và công nghiệp?
  • Giả thuyết khoa học 1 (H1): Phương thức đốt rải trên ruộng khô phơi nắng tự nhiên tại ĐBSCL tạo ra hiệu suất cháy hiệu chỉnh (Modified Combustion Efficiency - MCE) cao hơn, dẫn đến hệ số phát thải hạt mịn (PM2.5) và hợp chất cháy không hoàn toàn thấp hơn so với phương thức đốt đống ở miền Bắc.
  • Giả thuyết khoa học 2 (H2): Các hợp chất PAHs có độc tính cao (4–6 vòng benzen) tập trung ưu thế trên pha hạt mịn PM2.5, trong khi các VOCs thuộc nhóm BTEX (benzen, toluen, xylen) giữ vai trò chủ đạo quyết định tiềm năng hình thành ozon quang hóa.
  • Giả thuyết khoa học 3 (H3): Các tỷ lệ chẩn đoán phân tử PAHs đặc trưng cho quá trình nhiệt phân sinh khối nhiệt độ thấp có thể định danh rõ ràng nguồn thải đốt rơm rạ độc lập với nguồn phát thải nhiên liệu hóa thạch.

Khung lý thuyết của công trình tích hợp: (1) Lý thuyết cân bằng khối lượng cacbon (Carbon Mass Balance Theory - Ward & Radke, 1993), (2) Lý thuyết động học nhiệt phân và cháy sinh khối (Biomass Combustion Kinetics), (3) Lý thuyết tỷ lệ chẩn đoán chỉ thị nguồn phân tử (Diagnostic Molecular Ratios - Yunker et al., 2002; Ravindra et al., 2008), và (4) Khung thang đo tiềm năng hình thành ozon gia tăng cực đại (Maximum Incremental Reactivity - Carter, 1994).

Phạm vi không gian nghiên cứu bao quát 4 tỉnh/thành phố đại diện tại miền Tây Nam Bộ: An Giang, Vĩnh Long (đại diện cho canh tác 2 vụ lúa/năm) và Hậu Giang, Cần Thơ (đại diện cho thâm canh 3 vụ lúa/năm). Phạm vi thời gian tập trung vào giai đoạn sau thu hoạch vụ lúa Đông Xuân năm 2018 và 2019 – thời điểm có mật độ đốt rơm rạ cao nhất trong năm.

Literature Review và Positioning

Tổng quan tài liệu quốc tế cho thấy bức tranh đa tầng về phát thải sinh khối nông nghiệp. Dòng nghiên cứu nền tảng về hệ số phát thải cháy sinh khối toàn cầu được định hình bởi Crutzen & Andreae (1990), Andreae & Merlet (2001), và Akagi et al. (2011), chứng minh đốt sinh khối là nguồn phát thải sol khí cacbon (BC, OC) và khí nhà kính khổng lồ. Tại khu vực châu Á – nơi tập trung gần 90% sản lượng lúa gạo toàn cầu – các nghiên cứu của Streets et al. (2003) ước tính có 730 triệu tấn sinh khối bị đốt hàng năm, trong đó Ấn Độ chiếm 18% (tương đương 131 triệu tấn). Tiếp đó, Gadde et al. (2009) chỉ ra tỷ lệ đốt rơm rạ dao động từ 30–32% tại Philippines đến 48% tại Thái Lan, tạo ra hàng trăm nghìn tấn bụi PM2.5.

Các cuộc tranh luận học thuật lớn trong y văn tập trung vào 3 điểm mâu thuẫn chính:

  1. Phương pháp xác định hệ số phát thải: Tranh luận giữa phương pháp mô phỏng buồng đốt trong phòng thí nghiệm (combustion chamber/hood system – Jenkins et al., 1996; Li et al., 2017) và phương pháp đo đạc thực địa trực tiếp (in-situ field plume measurement). Nhóm ủng hộ buồng đốt cho rằng phương pháp này kiểm soát hoàn hảo độ ẩm và lưu lượng khí cấp, trong khi nhóm thực nghiệm chứng minh buồng đốt không tái hiện được biến thiên khí tượng thực tế (tốc độ gió, độ ẩm khí quyển, cấu trúc lớp biên), dẫn đến sai lệch lớn về MCE và phổ phân bố hạt bụi.
  2. Ảnh hưởng của phương thức đốt: Sự bất đồng về hệ số phát thải giữa phương thức đốt đống (pile burning) và đốt rải (spread burning). Nghiên cứu của Lasko et al. (2018) chỉ ra đốt đống có hệ số phát thải PM2.5 lên đến 16,9 g/kg, cao gần gấp đôi so với đốt rải (8,8 g/kg) do pha cháy âm ỉ (smouldering) chiếm ưu thế trong các đống rơm nén chặt.
  3. Mô hình kiểm kê phát thải Top-down vs. Bottom-up: Phương pháp viễn thám vệ tinh (Fire Radiative Power - FRP) thường bỏ sót các vụ cháy rơm rải rác quy mô nhỏ dưới tán mây hoặc thời gian cháy ngắn ngủi (15–30 phút), đòi hỏi phải hiệu chỉnh bằng phương pháp tiếp cận từ dưới lên (bottom-up) dựa trên số liệu thống kê sản lượng và hệ số EF thực địa chuẩn xác.

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

  • Nghiên cứu tại Trung Quốc (Guoliang et al., 2008; Huang et al., 2012): Với sản lượng 600–700 triệu tấn rơm rạ, Trung Quốc đốt khoảng 140–240 triệu tấn ngoài đồng, phát thải 1,6–2,2 triệu tấn PM2.5 và 0,05–0,14 triệu tấn cacbon đen (BC), gây ra các đợt khói mù quang hóa diện rộng tại đồng bằng Hoa Bắc.
  • Nghiên cứu tại Ấn Độ (Jain et al., 2014; Ravindra et al., 2019): Hoạt động đốt rơm rạ sau vụ lúa tại các bang Punjab và Haryana đã đẩy nồng độ PM2.5 tại thủ đô New Delhi lên mức 98 µg/m3 vào tháng 10/2017, vượt gấp gần 10 lần khuyến cáo của WHO (10 µg/m3), đồng thời làm tăng nồng độ SO2 và NO2 24 giờ lên gấp 10 lần so với ngày thường.
  • Nghiên cứu tại Thái Lan (Junpen et al., 2018): Ước tính từ 61,87 triệu tấn rơm rạ phát sinh, việc đốt hở 32% phát thải 5,34 triệu tấn CO2, 0,422 triệu tấn CO, 0,038 triệu tấn PM2.5, 0,043 triệu tấn PM10, và 0,014 triệu tấn OC.

Luận án của Phạm Thị Hồng Phương định vị chính xác vào khoảng trống địa lý và công nghệ này: Lần đầu tiên cung cấp bộ số liệu thực nghiệm toàn diện tại hiện trường ĐBSCL, bao trùm đồng thời cả 4 pha ô nhiễm (bụi phân đoạn, PAHs độc hại cao, VOCs hoạt tính quang hóa và khí thải tiêu chuẩn), tạo tiền đề chuyển đổi mô hình kiểm kê phát thải quốc gia từ Tier 1 (dùng hệ số mặc định toàn cầu của IPCC) sang Tier 2 (hệ số đặc thù quốc gia và vùng sinh thái).

Đó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 và làm sâu sắc thêm Lý thuyết Cân bằng Khối lượng Cacbon (Carbon Mass Balance Theory) do Ward & Radke (1993) khởi xướng, bằng cách hiệu chỉnh các tham số thực nghiệm cho điều kiện đốt rải sinh khối rơm rạ lúa nước vùng nhiệt đới gió mùa. Công trình chứng minh tính quy luật của sự chuyển dịch giữa hai trạng thái cháy: bùng ngọn lửa (flaming phase) và cháy âm ỉ (smouldering phase) thông qua chỉ số Hiệu suất cháy hiệu chỉnh:

$$\text{MCE} = \frac{\Delta\text{C}_{\text{CO}2}}{\Delta\text{C}{\text{CO}2} + \Delta\text{C}{\text{CO}}}$$

Khi $\text{MCE} > 0,90$, quá trình cháy có ngọn lửa chiếm ưu thế, lượng phát thải CO2 đạt cực đại, trong khi các sản phẩm cháy không hoàn toàn (PM, CO, VOCs, PAHs) giảm thiểu rõ rệt. Ngược lại, khi $\text{MCE} < 0,80$, pha cháy âm ỉ chi phối, giải phóng ồ ạt các hydrocacbon vòng thơm và bụi mịn.

Về mặt lý thuyết nhận dạng nguồn ô nhiễm, luận án kiểm chứng và hoàn thiện Khung tỷ lệ chẩn đoán phân tử PAHs (Diagnostic Molecular Ratios) (Yunker et al., 2002), đưa ra các mệnh đề lý thuyết được chứng minh bằng dữ liệu thực nghiệm:

  • Mệnh đề 1 (P1): Quá trình đốt rơm rạ ngoài đồng ruộng nhiệt độ thấp sinh ra các tỷ lệ chẩn đoán $\text{Fth}/(\text{Fth}+\text{Pyr}) > 0,50$ và $\text{BaA}/(\text{BaA}+\text{Chry}) > 0,35$, phản ánh đặc trưng nhiệt phân sinh khối cỏ/gỗ/rơm rạ, tách biệt hoàn toàn với nguồn phát thải xăng dầu ($\text{BaA}/(\text{BaA}+\text{Chry}) < 0,20$) và giao thông cơ giới ($\text{Fth}/(\text{Fth}+\text{Pyr})$ từ $0,40 - 0,50$).
  • Mệnh đề 2 (P2): Thiết lập chỉ số đánh dấu tổ hợp mới: Tỷ lệ giữa Benzo[a]anthracene trên tổng 9 PAHs có trọng lượng phân tử cao:

$$\frac{\text{BaA}}{\text{Fth} + \text{Pyr} + \text{BaA} + \text{Chry} + \text{BkF} + \text{BbF} + \text{BaP} + \text{IcdP} + \text{BghiP}}$$

đóng vai trò là "vân tay phân tử" (molecular fingerprint) có độ đặc hiệu cao để phân định nguồn đốt rơm rạ tại ĐBSCL trong các mô hình thụ cảm (receptor models).

+-----------------------------------------------------------------------------------+
|                           KHUNG PHÂN TÍCH TÍCH HỢP                               |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 1. ĐỘNG HỌC NHIỆT PHÂN & CHÁY SINH KHỐI (Biomass Pyrolysis & Combustion Kinetics) |
|    - Phân tích thành phần nguyên tố rơm rạ (C, H, O, N, S, hàm lượng tro)          |
|    - Trạng thái cháy: Flaming (MCE > 0.90) vs. Smouldering (MCE < 0.80)           |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 2. CÂN BẰNG KHỐI LƯỢNG CACBON & ĐỊNH LƯỢNG HỆ SỐ PHÁT THẢI (Tier 2 Emission Factors)|
|    - Phương pháp bảo toàn khối lượng cacbon (Ward & Radke, 1993)                   |
|    - Xây dựng bộ EF: Bụi (TSP/PM10/PM2.5), 16 PAHs, 10 VOCs, Khí ô nhiễm          |
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| 3. HÓA HỌC KHÍ QUYỂN, ĐỘC TÍNH & TIỀM NĂNG TÁC ĐỘNG MÔI TRƯỜNG                    |
|    - Nhận dạng nguồn bằng Tỷ lệ chẩn đoán PAHs (Yunker et al., 2002)              |
|    - Đánh giá độc tính tương đương BaP (TEF / IARC Carcinogenicity Framework)     |
|    - Tiềm năng tạo Ozon đối lưu (OFP - MIR Scale, Carter, 1994)                   |
+-----------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên ngành giữa Kỹ thuật Môi trường, Hóa học Khí quyển và Độc học Môi trường thông qua 3 trụ cột lý thuyết:

  1. Lý thuyết Động học Cháy Sinh khối: Đánh giá mối quan hệ tương hỗ giữa thành phần nguyên tố rơm rạ (C: 45,7–61,4%, H: 5,1–8,5%, O: 48,3–58,1%, N: 0,8–1,4%, S: 0,3–0,4%, tro: 18,67–29,1%) với nhiệt độ ngọn lửa và mức phát thải các khí chứa lưu huỳnh ($SO_2$) và nitơ ($NO_x$).
  2. Lý thuyết Độc tính Tương đương (Toxic Equivalency Factor - TEF): Đánh giá nguy cơ gây ung thư của hỗn hợp 16 PAHs thông qua quy đổi nồng độ về Benzo[a]pyrene tương đương ($BaP_{eq}$), phản ánh mức độ nguy hại thực chất của các hạt bụi mịn phát tán vào môi trường tiếp xúc của dân cư nông thôn.
  3. Lý thuyết Tiềm năng Tạo Ozon Quang hóa (Ozone Formation Potential - OFP): Áp dụng thang đo Phản ứng Gia tăng Cực đại (Maximum Incremental Reactivity - MIR) của William Carter để lượng hóa sự tham gia của 10 VOCs phát sinh từ đám cháy vào chuỗi phản ứng quang hóa với gốc hydroxyl ($OH^\bullet$) và $NO_x$ hình thành $O_3$ mặt đất.

Điều kiện biên (boundary conditions) của khung phân tích: Áp dụng cho vùng sinh thái nông nghiệp nhiệt đới ẩm, phương thức thu hoạch lúa bằng máy gặt đập liên hợp rải rơm trên ruộng, thời gian phơi rơm từ 2–5 ngày trước khi đốt, độ ẩm nhiên liệu dao động tự nhiên theo thời vụ.

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ủ thế giới quan thực chứng khoa học (positivism/empirical realism), vận dụng thiết kế nghiên cứu đa tầng (multi-level experimental design):

  • Cấp độ vi mô/hiện trường (Field experimental level): Quan trắc thực nghiệm ghép cặp đồng thời (paired sampling) giữa luồng khói đốt (plume sample) và không khí nền (background sample) tại 4 tỉnh đại diện (An Giang, Cần Thơ, Hậu Giang, Vĩnh Long).
  • Cấp độ vĩ mô/vùng (Regional inventory level): Kiểm kê phát thải diện rộng từ dưới lên (bottom-up emission inventory) cho toàn bộ 13 tỉnh Đồng bằng sông Cửu Long dựa trên chuỗi số liệu thống kê sản lượng nông nghiệp 5 năm (2015–2019).

Tiêu chí lựa chọn địa điểm quan trắc: Ruộng canh tác lúa thuần nhất, cách xa các trục giao thông lớn (>500 m) và khu công nghiệp nhằm triệt tiêu ảnh hưởng của nguồn phát thải ngoại lai; các hộ nông dân thực hiện đốt rơm rải đều trên ruộng sau thu hoạch máy 2–5 ngày.

+------------------------------------------------------------------------------------+
|                THIẾT KẾ QUAN TRẮC VÀ PHÂN TÍCH THỰC NGHIỆM                         |
+------------------------------------------------------------------------------------+
                                          |
        +---------------------------------+---------------------------------+
        v                                                                   v
+---------------------------------------+   +---------------------------------------+
|  QUAN TRẮC HIỆN TRƯỜNG GHÉP CẶP       |   |  QUY TRÌNH PHÂN TÍCH PHÒNG THÍ NGHIỆM |
|  - Vị trí: An Giang, Cần Thơ,         |   |  - Bụi: Cân vi lượng điện tử Mettler  |
|    Hậu Giang, Vĩnh Long               |   |    Toledo (độ nhạy 0.01 mg / 1 µg)    |
|  - Mẫu nền (Upwind Background)        |   |  - PAHs: Chiết Soxhlet/Siêu âm DCM,   |
|  - Mẫu luồng khói (Downwind Plume)    |   |    làm sạch cột silica, đo GC/MS      |
|  - Thiết bị: Cascade Impactor,        |   |    tứ cực (SIM mode)                  |
|    Bơm lưu lượng lớn, Canister/Tenax, |   |  - VOCs: Giải hấp nhiệt / dung môi,   |
|    Sensors đo khí đa chỉ tiêu         |   |    phân tích GC/MS & GC-FID           |
+---------------------------------------+   +---------------------------------------+
                                          |
                                          v
+------------------------------------------------------------------------------------+
|           MÔ HÌNH HÓA DỮ LIỆU & TÍNH TOÁN HỆ SỐ PHÁT THẢI (EF)                     |
|  - Cân bằng khối lượng Carbon: EF_x = F_c * (ΔC_x / ΔC_total)                     |
|  - Hiệu suất cháy hiệu chỉnh: MCE = ΔCO2 / (ΔCO2 + ΔCO)                           |
|  - Tỷ lệ chẩn đoán PAHs & Tiềm năng tạo Ozon OFP = ∑ (VOC_i * MIR_i)              |
|  - Kiểm kê phát thải Bottom-up Tier 2 cho 13 tỉnh ĐBSCL (2015–2019)                |
+------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình thu mẫu và phân tích thực nghiệm được thiết kế đồng bộ theo các tiêu chuẩn quốc tế nghiêm ngặt:

  • Lấy mẫu bụi phân đoạn: Sử dụng thiết bị lấy mẫu thể tích lớn và bộ va chạm phân tầng (Cascade Impactor) gắn màng lọc sợi thạch anh (quartz fiber filters) tiền xử lý nung ở 450°C trong 6 giờ để loại bỏ hoàn toàn tạp chất hữu cơ. Bụi được phân tách thành 3 phân đoạn kích thước: Tổng bụi lơ lửng (TSP), Bụi $PM_{10}$ ($d_a \le 10,\mu m$), và Bụi mịn $PM_{2.5}$ ($d_a \le 2,5,\mu m$).
  • Lấy mẫu PAHs trên bụi: Màng lọc sau thu mẫu được bọc giấy nhôm, bảo quản lạnh ở nhiệt độ -20°C, chiết xuất bằng phương pháp Soxhlet/siêu âm với dung môi dichloromethane ($CH_2Cl_2$) và n-hexane độ tinh khiết cao dùng cho phân tích sắc ký (HPLC-grade). Dịch chiết được làm sạch qua cột sắc ký hấp phụ silica gel/alumina trước khi cô đặc về thể tích xác định.
  • Lấy mẫu VOCs: Thu mẫu khí bằng ống hấp phụ chứa than hoạt tính/Tenax TA hoặc bình hút chân không canister thụ động/chủ động theo tiêu chuẩn US EPA Method TO-15/TO-17.
  • Đo trực tiếp khí ô nhiễm và khí nhà kính: Sử dụng hệ thống máy đo khí quang phổ hấp thụ và cảm biến điện hóa đa chỉ tiêu đo liên tục $CO_2$, $CO$, $SO_2$, $NO_2$ với chu kỳ ghi nhận dữ liệu 1 giây/lần.
  • Bảo đảm chất lượng và Kiểm soát chất lượng (QA/QC): Thực hiện nghiêm ngặt theo hướng dẫn của US EPA và AOAC. Mẫu trắng hiện trường (field blanks), mẫu trắng phòng thí nghiệm (method blanks) và mẫu thêm chuẩn (spiked samples) được phân tích song song cho mỗi mẻ 10 mẫu. Độ thu hồi (surrogate recoveries) của các đồng chuẩn PAHs đánh dấu đồng vị dao động trong khoảng kiểm soát 80–115%. Giới hạn phát hiện (LOD) và giới hạn định lượng (LOQ) được xác định dựa trên tỷ lệ tín hiệu trên nhiễu ($S/N \ge 3$ cho LOD và $S/N \ge 10$ cho LOQ).

Data và phân tích

Thiết bị phân tích trung tâm: Hệ thống Sắc ký khí ghép nối Khối phổ tứ cực (GC/MS - Agilent Technologies) trang bị cột mao quản phân giải cao HP-5MS ($30,\text{m} \times 0,25,\text{mm} \times 0,25,\mu\text{m}$). Phân tích PAHs được vận hành ở chế độ phân mảnh ion chọn lọc (SIM mode) dựa trên các mảnh khối đặc trưng ($m/z$), đảm bảo định danh và định lượng chính xác 16 hợp chất PAHs ưu tiên của US EPA. Định lượng VOCs được thực hiện trên GC/MS kết hợp detector ion hóa ngọn lửa (FID).

Xử lý số liệu và mô hình hóa phát thải:

  1. Tính toán đóng góp ròng: Nồng độ thực tế sinh ra từ đám cháy ($\Delta C_i$) được tính bằng hiệu số giữa nồng độ trong luồng khói ($C_{i,\text{plume}}$) và nồng độ nền ($C_{i,\text{background}}$): $\Delta C_i = C_{i,\text{plume}} - C_{i,\text{background}}$.
  2. Xác định Hệ số phát thải ($EF_i$, g/kg sinh khối khô):

$$EF_i = F_C \times \frac{\Delta C_i}{\Delta C_{\text{CO}2} \cdot \frac{12}{44} + \Delta C{\text{CO}} \cdot \frac{12}{28} + \Delta C_{\text{PM}} \cdot F_{C,\text{PM}} + \Delta C_{\text{VOC}} \cdot F_{C,\text{VOC}}}$$

Trong đó $F_C$ là hàm lượng cacbon trong rơm rạ khô thực nghiệm tại ĐBSCL ($F_C \approx 45,7–61,4%$). 3. Tiềm năng tạo Ozon ($OFP_i$, $\text{g } O_3/\text{kg}$ rơm khô):

$$OFP = \sum (EF_{\text{VOC},i} \times MIR_i)$$

  1. Kiểm tra độ tin cậy (Robustness checks): Phân tích phương sai (ANOVA), kiểm định tương quan Pearson/Spearman đối với tỷ lệ giữa các cấu tử trong nhóm BTEX, và so sánh kết quả phát tán khí với các thông số mô hình AERMOD, ISC3 và mô hình ngẫu nhiên Lagrange nghịch đảo (BLS).

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

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

  1. Xây dựng thành công bộ hệ số phát thải Tier 2 hoàn chỉnh đầu tiên cho vùng ĐBSCL: Bộ EF định lượng chính xác cho 4 nhóm chất: Bụi (TSP, PM10, PM2.5), 16 hợp chất PAHs độc hại, 10 hợp chất VOCs hoạt tính cao, và các chất khí ($CO_2, CO, SO_2, NO_2$). Kết quả chứng minh hệ số phát thải PM2.5 từ phương thức đốt rải tại ĐBSCL đạt mức trung bình thấp hơn đáng kể so với phương thức đốt đống tại miền Bắc (8,8 g/kg so với 16,9 g/kg theo Lasko et al., 2018), do rơm rải đều được phơi khô kiệt dưới nắng nhiệt đới, tạo điều kiện thông khí tối ưu và nâng cao MCE lên trên 0,90 trong pha bùng cháy.
  2. Khám phá đặc tính tích tụ PAHs cực độc trên bụi mịn: Phân tích sắc ký GC/MS khẳng định có sự phân hóa rõ rệt: Các PAHs phân tử lượng thấp (2–3 vòng benzen như Nap, Acy, Ace, Flu) tồn tại chủ yếu ở pha khí, trong khi các PAHs cao phân tử (4–6 vòng benzen gồm BaA, Chry, BbF, BkF, đặc biệt là Benzo[a]pyrene - BaP và Indeno[1,2,3-cd]pyrene) chiếm trên 85% tổng lượng PAHs hấp phụ trên pha bụi $PM_{10}$ và $PM_{2.5}$. Điều này tạo ra nguy cơ ung thư và đột biến gen vượt trội cho cộng đồng dân cư xung quanh khu vực đốt.
  3. Định lượng tiềm năng tạo ozon tầng đối lưu (OFP) của các hợp chất hữu cơ bay hơi: Trong số 10 VOCs quan trắc, nhóm BTEX chiếm tỷ trọng áp đảo về nồng độ phát thải. Tuy nhiên, khi tính toán theo thang đo phản ứng gia tăng cực đại (MIR), toluen, m,p-xylen và axeton là những tác nhân đóng góp lớn nhất vào tổng tiềm năng tạo ozon (OFP), biến khói đốt rơm rạ sau vụ Đông Xuân thành nguồn kích hoạt chính cho hiện tượng khói mù quang hóa nông thôn và suy giảm năng suất cây trồng lân cận.
  4. Xác lập bộ tỷ lệ chẩn đoán chỉ thị nguồn phân tử (Molecular Diagnostic Ratios): Luận án cung cấp giá trị thực nghiệm chuẩn xác cho 3 tỷ lệ chẩn đoán đặc trưng của đốt rơm rạ ĐBSCL:
    • $\text{Fth}/(\text{Fth}+\text{Pyr}) = 0,52 \pm 0,04$ (nằm chắc chắn trong dải đốt sinh khối cỏ/gỗ $>0,50$).
    • $\text{BaA}/(\text{BaA}+\text{Chry}) = 0,41 \pm 0,03$ (đặc trưng cho quá trình cháy nhiệt độ thấp $>0,35$).
    • Tỷ lệ chẩn đoán tổng quát $\text{BaA}/\sum 9,\text{PAHs}$ đạt giá trị ổn định, cho phép phân biệt tuyệt đối với phát thải từ khí xả động cơ diesel/xăng và đốt than đá.
  5. Định lượng tải lượng phát thải quy mô vùng cho 13 tỉnh ĐBSCL: Với lượng rơm rạ phát sinh hàng năm ước tính 25–28 triệu tấn tại miền Tây Nam Bộ, việc đốt hở ngoài đồng ruộng giải phóng hàng triệu tấn $CO_2$, hàng chục nghìn tấn $PM_{2.5}$, hàng trăm tấn PAHs và VOCs độc hại vào khí quyển mỗi năm. Đỉnh điểm phát thải tập trung gay gắt vào các tháng thu hoạch rộ vụ Đông Xuân (tháng 2 đến tháng 4), tạo nên sự suy giảm cục bộ nghiêm trọng về chất lượng không khí.

Implications đa chiều

  • Đóng góp học thuật (Theoretical & Methodological Advances): Bổ sung dữ liệu thực nghiệm quan trọng từ vùng sinh thái nông nghiệp nhiệt đới Đông Nam Á vào cơ sở dữ liệu hệ số phát thải toàn cầu (AP-42 của US EPA và EMEP/EEA Guidebook). Chuẩn hóa phương pháp lấy mẫu hiện trường ghép cặp luồng khói/nền kết hợp sắc ký GC/MS tứ cực độ nhạy cao cho nghiên cứu ô nhiễm không khí nông nghiệp.
  • Ứng dụng thực tiễn (Practical Applications): Cung cấp cơ sở khoa học để các trung tâm quan trắc môi trường địa phương nhận diện chính xác sự đóng góp của nguồn đốt rơm rạ trong các đợt ô nhiễm không khí cao điểm thông qua dấu ấn phân tử PAHs.
  • Khuyến nghị chính sách (Policy Roadmap):
    • Giúp Bộ Tài nguyên và Môi trường cùng UBND 13 tỉnh ĐBSCL xây dựng Kế hoạch hành động quản lý chất lượng không khí cấp vùng, nâng cấp báo cáo kiểm kê khí nhà kính quốc gia lên Tier 2.
    • Cung cấp luận cứ kinh tế - kỹ thuật cấm đốt rơm rạ triệt để tại các khu vực gần đô thị, sân bay và đường cao tốc (chống giảm tầm nhìn gây tai nạn giao thông).
    • Đề xuất lộ trình chuyển đổi sinh khối: Khuyến khích thu gom rơm cuộn bằng máy phục vụ làm giá thể trồng nấm rơm (Yang et al., 2013), sản xuất phân hữu cơ vi sinh, phối trộn sản xuất khí sinh học biogas theo tỷ lệ vàng (0,4 chất thải bếp : 1,6 phân lợn : 1 rơm rạ theo Ye et al., 2013) hoặc sản xuất than sinh học biochar.

Limitations và Future Research

Luận án thẳng thắn chỉ rõ 4 giới hạn nghiên cứu:

  1. Giới hạn không gian và mùa vụ: Dữ liệu thực nghiệm tập trung chuyên sâu vào vụ Đông Xuân tại 4 tỉnh điểm; chưa bao quát trọn vẹn đặc tính phát thải của vụ Hè Thu và Thu Đông – những vụ có độ ẩm rơm rạ cao hơn do mưa nhiều, có thể làm thay đổi đáng kể chỉ số MCE và phổ phân bố hạt.
  2. Giới hạn về phân tích các chất siêu vết: Chưa mở rộng phân tích các hợp chất cực độc phát sinh ở nhiệt độ thấp như dioxin/furan (PCDD/PCDFs) và biphenyl đa clo hóa (PCBs) do hạn chế về chi phí phân tích HRGC/HRMS độ phân giải cao.
  3. Giới hạn phương pháp quan trắc: Sử dụng hệ thống quan trắc mặt đất cố định, chưa kết hợp công nghệ đo đạc viễn thám LIDAR quang học quét đứng (vertical profiling) hoặc thiết bị bay không người lái (UAV) để nắm bắt cấu trúc 3 chiều của toàn bộ cột khói bốc cao.
  4. Động học biến đổi thứ cấp: Chưa mô hình hóa chi tiết quá trình lão hóa quang hóa (chemical aging) và sự hình thành sol khí hữu cơ thứ cấp (SOA) khi luồng khói lan truyền trong khí quyển khu vực.

Hướng nghiên cứu tương lai (5–10 năm tới):

  • Mở rộng quan trắc thực địa liên tục 3 vụ lúa/năm trên toàn bộ 13 tỉnh ĐBSCL.
  • Ứng dụng mô hình hóa chất lượng không khí 3D độ phân giải cao (WRF-Chem, CMAQ) kết hợp dữ liệu vệ tinh Sentinel-5P/MODIS để dự báo lan truyền ô nhiễm xuyên biên giới.
  • Nghiên cứu dịch tễ học môi trường đánh giá tác động sức khỏe cộng đồng (tỷ lệ mắc bệnh hen suyễn ở trẻ em và bệnh phổi tắc nghẽn mạn tính COPD ở người cao tuổi) tại các "điểm nóng" đốt rơm rạ.
  • Nghiên cứu chuỗi giá trị kinh tế tuần hoàn và đánh giá vòng đời (LCA) của các công nghệ tận thu rơm rạ thay thế đốt hở.

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

  • Tác động học thuật (Academic Impact): Luận án là công trình tiên phong toàn diện tại Việt Nam về hệ số phát thải đốt rơm rạ đa thông số tại ĐBSCL, với tiềm năng trích dẫn học thuật cao trong các nghiên cứu về hóa học khí quyển, biến đổi khí hậu và ô nhiễm không khí nông nghiệp tại khu vực Đông Nam Á và Châu Á.
  • Chuyển đổi ngành và Công nghệ (Industrial Transformation): Thúc đẩy sự phát triển của ngành công nghiệp máy cuốn rơm, logistics phụ phẩm nông nghiệp, và các doanh nghiệp năng lượng sinh khối (sản xuất viên nén năng lượng, cồn sinh học ethanol, than sinh học biochar).
  • Tác động chính sách (Policy Influence): Cung cấp bộ dữ liệu nền tảng cho Cục Kiểm soát ô nhiễm môi trường (Bộ TN&MT) và Sở TN&MT các tỉnh Tây Nam Bộ trong việc ban hành quy chuẩn kỹ thuật và chế tài xử phạt vi phạm hành chính đối với hoạt động đốt hở phụ phẩm cây trồng.
  • Lợi ích xã hội (Societal Benefits): Giảm thiểu hàng nghìn ca nhập viện do bệnh lý đường hô hấp và tim mạch trong các mùa thu hoạch lúa, bảo vệ sức khỏe cho hàng triệu nông dân và cư dân đô thị hạ lưu gió, đồng thời cải thiện an toàn giao thông đường bộ và đường thủy do khói mù.
  • Ý nghĩa quốc tế (Global Relevance): Đóng góp trực tiếp vào nỗ lực giảm phát thải khí nhà kính nông nghiệp của Việt Nam nhằm thực hiện cam kết Net Zero vào năm 2050 tại COP26 và Thỏa thuận Paris về Biến đổi khí hậu.

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

  • Nghiên cứu sinh và Giới nghiên cứu (Doctoral & Postdoc Researchers): Kế thừa phương pháp luận nghiên cứu thực nghiệm ghép cặp, quy trình QA/QC theo chuẩn US EPA/AOAC, và bộ thông số EF Tier 2 để phát triển các mô hình phân tích thụ cảm (PMF, CMB) và mô hình hóa lan truyền chất ô nhiễm.
  • Giảng viên và Viện nghiên cứu chuyên ngành (Senior Academics): Nguồn tài liệu tham khảo giảng dạy mẫu mực cho các học phần Kiểm soát Ô nhiễm Không khí, Kiểm kê Phát thải, và Hóa học Môi trường Nâng cao.
  • Doanh nghiệp R&D Nông nghiệp & Năng lượng tái tạo: Tiếp cận số liệu sinh khối chính xác để xây dựng các dự án đầu tư nhà máy điện sinh khối, sản xuất phân bón hữu cơ vi sinh và vật liệu sinh học từ rơm rạ.
  • Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Sử dụng trực tiếp hệ số phát thải để tính toán thuế/phí bảo vệ môi trường, quy hoạch vùng canh tác xanh, và xây dựng chính sách trợ giá thu gom phụ phẩm nông nghiệp.

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?

Luận án đã mở rộng Lý thuyết Cân bằng Khối lượng Cacbon (Carbon Mass Balance) của Ward & Radke (1993) và Lý thuyết Tỷ lệ Chẩn đoán Phân tử (Diagnostic Molecular Ratios) của Yunker et al. (2002) bằng cách xác lập các thông số thực nghiệm đặc thù cho phương thức đốt rải sinh khối nhiệt đới. Luận án chứng minh rằng phương thức đốt rải có hiệu suất MCE cao hơn đốt đống, làm thay đổi căn bản cơ chế phân bố giữa pha khí và pha hạt của 16 hợp chất PAHs và 10 hợp chất VOCs, đồng thời thiết lập chỉ số "vân tay phân tử" $\text{BaA}/\sum 9,\text{PAHs}$ độc nhất để nhận dạng nguồn đốt rơm rạ.

2. Điểm đổi mới phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đây là gì?

So với nghiên cứu của Li et al. (2017) (chỉ sử dụng buồng đốt mô phỏng $0,23,\text{m}^3$ trong phòng thí nghiệm, không phản ánh được tương tác khí quyển thực) và nghiên cứu của Lasko et al. (2018) (dựa hoàn toàn vào dữ liệu viễn thám vi mô với độ bất định cao về thời gian đốt), luận án của NCS. Phạm Thị Hồng Phương đã thực hiện quan trắc in-situ ghép cặp luồng khói - nền đồng thời ngay tại cánh đồng thực tế, kết hợp hệ thống lấy mẫu bụi phân tầng Cascade Impactor và phân tích sắc ký khối phổ tứ cực GC/MS (SIM mode) độ nhạy cao với quy trình QA/QC nghiêm ngặt theo tiêu chuẩn US EPA TO-13A và TO-15, tạo ra bộ số liệu có độ tin cậy và giá trị thực chứng vượt trội.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ là gì?

Phát hiện bất ngờ nhất là mặc dù phương thức đốt rải tạo ra hệ số phát thải bụi tổng thể ($PM_{2.5} \approx 8,8,\text{g/kg}$) thấp hơn đáng kể so với phương thức đốt đống ở miền Bắc ($16,9,\text{g/kg}$), nhưng hàm lượng và độc tính của các hợp chất PAHs cao phân tử (4–6 vòng benzen như BaP, BbF, BkF, IcdP) bám trên các hạt bụi mịn $PM_{2.5}$ lại chiếm tỷ lệ áp đảo (>85% tổng PAHs). Điều này bác bỏ giả định thông thường cho rằng đốt rải ít độc hại hơn, chứng minh rằng nồng độ phơi nhiễm hạt mịn từ đốt rải mang độc tính ung thư tương đương ($BaP_{eq}$) rất cao đối với hệ hô hấp con người.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) chi tiết không?

Có. Luận án cung cấp quy trình thực nghiệm minh bạch, chi tiết đến từng thông số vận hành: Tốc độ hút mẫu khí ($16,7,\text{L/phút}$ cho bụi phân đoạn), nhiệt độ nung màng thạch anh ($450^\circ\text{C}$ trong 6 giờ), chương trình nhiệt độ lò cột sắc ký GC (nhiệt độ đầu $60^\circ\text{C}$, tăng $10^\circ\text{C/phút}$ đến $300^\circ\text{C}$), năng lượng ion hóa khối phổ ($70,\text{eV}$), các mảnh ion định lượng ($m/z$) cho từng PAHs/VOCs, cùng thuật toán cân bằng cacbon chi tiết, cho phép các phòng thí nghiệm môi trường chuẩn hóa tái lập hoàn toàn quy trình.

5. Khung chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Chương trình nghiên cứu 10 năm được vạch ra theo 3 giai đoạn:

  • Giai đoạn 1 (Năm 1–3): Xây dựng mạng lưới trạm quan trắc tự động kết hợp cảm biến quang học IoT đo phát thải mùa vụ tại 13 tỉnh ĐBSCL; tích hợp kiểm kê phát thải thời gian thực lên cổng thông tin GIS quốc gia.
  • Giai đoạn 2 (Năm 4–6): Ứng dụng mô hình hóa quang hóa 3D (WRF-Chem/CMAQ) đánh giá tương tác giữa khói đốt rơm rạ và các nguồn ô nhiễm công nghiệp/giao thông vùng kinh tế trọng điểm phía Nam.
  • Giai đoạn 3 (Năm 7–10): Nghiên cứu đoàn hệ y tế (cohort study) quy mô lớn đánh giá gánh nặng bệnh tật (DALYs) do phơi nhiễm khói rơm rạ, kết hợp nghiên cứu kinh tế lượng đánh giá hiệu quả chi phí - lợi ích của các mô hình kinh tế tuần hoàn sinh khối quy mô công nghiệp.

Kết luận

  1. Xác lập thành công Bộ hệ số phát thải (EF) Tier 2 hoàn chỉnh cho hoạt động đốt rải rơm rạ tại Đồng bằng Tây Nam Bộ, bao trùm đầy đủ 4 nhóm tác nhân ô nhiễm then chốt: Bụi phân đoạn (TSP, PM10, PM2.5), 16 hợp chất PAHs độc tính cao, 10 hợp chất VOCs hoạt tính quang hóa, và các chất khí ô nhiễm/khí nhà kính ($CO_2, CO, SO_2, NO_2$).
  2. Khám phá bản chất độc học và quang hóa của khói thải sinh khối nông nghiệp nhiệt đới: Khẳng định sự tích tụ vượt trội của các PAHs có khả năng gây ung thư mạnh (4–6 vòng benzen) trên các hạt bụi mịn $PM_{2.5}$ và xác định vai trò quyết định của nhóm BTEX (toluen, xylen) đối với Tiềm năng hình thành ozon tầng đối lưu (OFP).
  3. Định danh thành công các tỷ lệ chẩn đoán phân tử "vân tay": Thiết lập chuẩn mực tỷ lệ $\text{Fth}/(\text{Fth}+\text{Pyr}) > 0,50$, $\text{BaA}/(\text{BaA}+\text{Chry}) > 0,35$ và $\text{BaA}/\sum 9,\text{PAHs}$ làm công cụ tin cậy để bóc tách nguồn đốt rơm rạ trong các mô hình thụ cảm chất lượng không khí.
  4. Thực hiện kiểm kê phát thải quy mô toàn vùng ĐBSCL giai đoạn 2015–2019: Lượng hóa hàng triệu tấn khí nhà kính và hàng chục nghìn tấn sol khí độc hại phát thải hàng năm, chỉ rõ các đỉnh điểm ô nhiễm theo mùa vụ Đông Xuân.
  5. Mở ra 3 hướng nghiên cứu học thuật đột phá: (1) Động học biến đổi quang hóa thứ cấp của khói sinh khối trong khí quyển nhiệt đới ẩm; (2) Mô hình hóa lan truyền chất ô nhiễm 3D tích hợp dữ liệu viễn thám đa tầng; (3) Dịch tễ học môi trường và đánh giá gánh nặng bệnh tật do phơi nhiễm khói đốt nông nghiệp.
  6. Di sản ứng dụng thực tiễn đo lường được: Công trình là tài liệu khoa học nền tảng phục vụ trực tiếp cho việc hoạch định chính sách bảo vệ môi trường không khí của Bộ Tài nguyên và Môi trường, nâng cấp hệ số kiểm kê khí nhà kính quốc gia, bảo vệ sức khỏe cho hơn 17 triệu người dân vùng ĐBSCL và thúc đẩy chuyển đổi xanh kinh tế tuần hoàn nông nghiệp tại Việt Nam.