Việc tính toán hệ số phát thải của các chất [x] được dựa trên phương pháp cân bằng các bon. Công thức tính toán như sau [31, 70, 116]
- Hệ số phát thải của CO2
+ Phương trình cân bằng cacbon: Co = Cb − Ca (2.7)
Trong đó:
+ Co: Lượng C trong rơm đã bị cháy, kgC/kg;
+ Cb: Lượng C có trong rơm trước khi đốt, kgC/kg rơm; + Ca: Lượng C còn lại trong tro sau khi đốt, kgC/kg tro.
Khi bị đốt cháy, carbon (C) trong sinh khối đươc phát thải dưới dạng các hợp chất chứa C như CO2, CO, CH4, NMHC và C trên bụi. Tuy nhiên, các nghiên cứu đã chỉ ra rằng, phát thải dưới dạng CO2 và CO là chủ yếu, chiếm khoảng 90% tổng lượng cacbon đã phát thải [31, 116]
Khi đó, hệ số phát thải của CO2 được xác định theo công thức: EFCO2 = 0,9 x ∆CO2 ∆CO+∆C O2 x Co x44 12 (2.8)
Trong đó 0,9 tức là 90% lượng C thải ra là ở dạng CO2 và CO [31, 46]
- Hệ số phát thải của chất X (PM, CO, NO2, SO2, VOCs) được tính theo các phương trình sau:
+ Trường hợp CO là chất tham chiếu (quá trình cháy âm ỉ):
EFX = ERX/CO x ERCO/CO2 x EFCO2 (2.9) + Trường hợp CO2 là chất tham chiếu (quá trình cháy ngọn lửa):
EFX = ERX/CO2 x EFCO2 (2.10)
Trong đó:
EFCO2 : Hệ số phát thải của CO2 (g/kg)
ERCO/CO2 : Tỷ lệ phát thải của CO so với CO2
- Hệ số phát thải của PAHs được tính như sau:
EFPAHs = CPAHs x EFPM (2.11)
Trong đó :
CPAHs Hàm lượng trên PM (μg/g PM)
EFPM Hệ số phát thải của PM (g/ kg rơm rạ khô). 2.5. Kiểm kê phát thải
Lượng chất thải phát sinh từ hoạt động đốt rơm rạ được kiểm kê bằng lượng rơm rạ đem đốt và hệ số phát thải tương ứng của mỗi chất thải tương ứng.
2.5.1. Xác định lượng rơm rạ được đốt
(1) Công thức tính
Lượng rơm rạ đốt cháy được ước tính dựa trên sản lượng lúa, các tỉ lệ phụ phẩm và hiệu suất đốt theo công thức như sau [1]
M = P x N x Dx B x MCE (2.12)
Trong đó:
M: Lượng rơm rạ đốt cháy (kg/năm) P: Sản lượng lúa (ha/ năm)
N: Tỷ lệ rơm rạ phát sinh và sản lượng lúa D: Tỷ trọng khô của rơm rạ
B: Tỷ lệ rơm rạ được đốt trên đồng ruộng (0 -1); MCE: hiệu suất cháy (%).
(2) Các thông số phục vụ tính toán
• Sản lượng lúa gạo cả năm (P)
o Số liệu về sản lượng lúa hàng năm tại miền Tây Nam Bộ được kế thừa từ Báo cáo Niên giám thông kê của Tổng cục thống kê từ năm 2016 tới năm 2020.
o Số liệu về sản lượng lúa theo mùa vụ bao gồm: vụ Đông Xuân, Hè Thu và Thu Đông được sử dụng từ số liệu của Niên giám thống kê năm 2019 của Miền Tây Nam Bộ.
Tỷ lệ rơm rạ phát sinh theo sản lượng lúa là 1,3. Tỷ lệ này là giá trị trung bình của rơm rạ với sản lượng lúa được áp dụng theo kết quả được công bố mới nhất của nhóm tác giả Gurraj Singh và cộng sự (2021) áp dụng cho Việt Nam [9]
• Tỷ lệ đốt rơm rạ (B)
Số liệu về tỷ lệ đốt rơm rạ cho việc kiểm kê theo mùa tại miền Tây Nam Bộ là giá trị trung bình của kết quả phỏng vấn 120 hộ dân tại miền Tây Nam Bộ của năm 2019 (kế thừa số liệu của đề tài “Nghiên cứu xác định hệ số phát thải khí nhà kính từ hoạt động đốt hở các phụ phẩm nông nghiệp (trấu, rơm rạ) vùng Tây Nam Bộ” chủ nhiệm đề tài PGS.TS. Phạm Thị Mai Thảo. Kết quả kháo sát về phương thức sử dụng rơm rạ sau thu hoạch vào các mục đích khác nhau của các hộ dân tại miền Tây Nam Bộ (năm 2019) và kết quả kế thừa từ các công bố trước đây [10, 124] cho thấy tỷ lệ đốt rơm rạ đốt ngoài đồng ruộng tại miền Tây Nam Bộ chiếm 95%, 45% và 25% tương ứng với vụ đông xuân, hè thu và vụ thu đông [132]. Mặc dù nhiệt độ không khí trong khu vực nghiên cứu không có sự khác biệt lớn theo mùa, nhưng lượng mưa có xu hướng cao hơn và số giờ nắng có xu hướng ngắn hơn trong vụ hè thu so với các vụ lúa khác. Vì tỷ lệ đốt rơm rạ thu hoạch (trọng lượng khô của rơm đốt/trọng lượng khô của rơm rạ thu hoạch) có xu hướng thấp hơn vào vụ hè thu khi lượng mưa lúc thu hoạch cao. Điều đó chứng tỏ lượng mưa lúc thu hoạch làm gián đoạn việc đốt rơm rạ và nông dân trồng lúa không thích đốt rơm rạ vào vụ hè thu. Mặc khác, vụ đông xuân được biết đến là vụ khô hạn nhất trong năm, điều kiện thuận lợi nhất cho việc đốt rơm rạ ngay tại đồng ruộng. Do đó, vụ đông xuân được xác định là vụ có tỷ lệ đốt cao nhất trong ba vụ thu hoạch lúa tại miền Tây Nam Bộ.
Ngoài ra, các giống lúa truyền thống được gieo trong mỗi vụ lúa của hệ thống lúa ba vụ đã nghiên cứu trong đó vụ đông xuân (IR50404 chiếm 91% tổng số diện tích lúa), vụ xuân hè (OM4218 chiếm 80% tổng số diện tích lúa) và vụ
(OM4218 chiếm 77% tổng số diện tích lúa). Ngoài các giống lúa này, các giống OM551 và MTL480 cũng được sử dụng trong các diện tích canh tác nhỏ hơn. Theo thông tin có được từ nông dân, IR50404 cần nhiều bức xạ mặt trời hơn và thời gian sinh trưởng dài hơn (khoảng 105-115 ngày/vụ) so với các giống lúa không nhạy cảm với ánh sáng như OM4218 (khoảng 85-95 ngày/vụ). Vì thời vụ trồng IR50404 thường được coi là dài đối với hệ thống canh tác lúa ba vụ, nên trong hầu hết các trường hợp, giống cây này chỉ được trồng vào vụ đông xuân, mùa khô và hạn nhất. OM4218 và các giống lúa khác đã được sử dụng ở khu vực này cho hai vụ mùa vụ còn lại trong năm. Do đó, hầu hết rơm rạ thu hoạch trong vụ đông xuân đều bị đốt cháy để tiết kiệm thời gian chuẩn bị cho vụ gieo trồng tiếp theo.
Tuy nhiên, trong phạm vi nghiên cứu này, chỉ thực hiện tại một số tỉnh có hệ thống canh tác lúa ba vụ đại diện ở Đồng bằng sông Cửu Long, do đó, kết quả có thể không đủ rộng để ước tính việc quản lý rơm rạ ở toàn bộ Đồng bằng sông Cửu Long, bao gồm các hệ thống canh tác lúa một hoặc hai vụ hay bốn vụ/ năm. Tuy nhiên, kết quả của luận án xác định tỷ lệ đốt rơm rạ của ba vụ thu hoạch đối với hệ thống canh tác lúa ba vụ đại diện trong khu vực nghiên cứu.
Số liệu về tỷ lệ đốt rơm rạ dùng để kiểm kê phát thải hàng năm tại miền Tây Nam Bộ là giá trị trung bình của số liệu phỏng vấn theo mùa vụ của năm 2019 và số liệu của các công trình đã được công bố, cụ thể là 95%, 45% và 25% tương ứng với vụ đông xuân, hè thu và vụ thu đông [9, 10, 124].
• Hiệu suất cháy hiệu chỉnh MCE.
Căn cứ vào kết quả khảo sát thực tế và thực hiện thí nghiệm trên đồng ruộng, rơm rạ được rải đều trên mặt ruộng, phơi khô tự nhiên và đốt từ 3-5 ngày sau thu hoạch. Nhờ các lý do trên dẫn đến độ ẩm tương đối thấp, ngoài ra do đốt hở nên việc
cung cấp oxy khá lớn hơn so với quá trình đốt đống và quá trình cháy là cháy ngọn lửa. Vì vậy MCE = 0,9 được chọn áp dụng cho việc tính toán.
2.5.2. Xác định mức độ phát thải
(1) Công thức tính
Công thức tính phát thải từ đốt rơm rạ phục vụ kiểm kê phát thải theo hàng năm và kiểm kê phát thải theo mùa vụ được tỉnh như sau [15, 16]:
EAi = M x EFi (2.13)
Trong đó:
i: chất ô nhiễm i
EAi : lượng chất thải của chất ô nhiễm i từ hoạt động đốt rơm rạ Mj: sản lượng rơm rạ được đốt cháy (kg/ năm)
EFi: hệ số phát thải của chất ô nhiễm i (g/ kg)
Hệ số phát thải được dùng tính toán trong nghiên cứu này là giá trị trung bình của bộ hệ số phát thải từ kết quả nghiên cứu thực nghiệm tại 4 tỉnh miền Tây Nam Bộ. 2.6. Đánh giá tác động
2.6.1. Đánh giá tác động trực tiếp
Phương pháp so sánh được áp dụng trong nghiên cứu này nhằm đánh giá sự gia tăng nồng độ phát thải các chất ô nhiễm từ quá trình đốt rơm rạ so với môi trường nền tại cùng một vị trí, từ đó xác định được sự gia tăng nông độ của chất cần quan tâm. Trên cơ sở đó có thể đánh giá được tác động trực tiếp ở quy mô địa phương của chất thải phát sinh từ hoạt động đốt rơm rạ:
∆(i) = C(i) đốt – C(i) nền (2.14)
Trong đó:
∆(i) là sự gia tăng ( đóng góp ròng) của chất thải i
C(i) đốt : Nồng độ của chất i xác định được trong quá trình đốt C(i) nền : Nồng độ của chất i xác định được trong môi trường nền. 2.6.2. Đánh giá tiềm năng hình thành ozon
Các giá trị nồng độ trung bình của VOCs được sử dụng để tính toán tiềm năng hình thành ozon ở tầng đối lưu và rủi ro sức khỏe để đánh giá ảnh hưởng của VOCs trong khí quyển và sức khỏe con người. Trong nghiên cứu này sẽ sử dụng chỉ số MIR để đánh giá tiềm năng hình thành ozon trong khí quyển từ hoạt động đốt hở rơm rạ. Tiềm năng hình thành ozon (Ozon potential formation – OFP) được tính bằng công thức sau [133]:
OFP = MIR × [VOCs] (μg/m3 ) (2.15) Trong đó:
OFP: tiềm năng hình thành ozon (μg/m3 O3)
MIR: tỉ lệ gia tăng phản ứng tối đa, đặc trưng cho từng chất (g O3/gVOC) [VOCs]: nồng độ VOCs đo được (μg/m3).
Chương 3. KẾT QUẢ VÀ THẢO LUẬN
3.1. Hệ số phát thải các chất ô nhiễm không khí
3.1.1. Đặc tính của quá trình cháy và thông số tính toán
3.1.1.1. Đặc tính của quá trình cháy và đóng góp ròng của các chất ô nhiễm
Hàm lượng cacbon của rơm rạ liên quan đến từng giống lúa trước khi đốt và hàm lượng tro sau khi đốt trong mỗi thí nghiệm được xác định để tính toán hiệu suất oxy hóa, và xác định hệ số phát thải. Các giống lúa cho các thí nghiệm trong khu vực nghiên cứu này là các giống được sử dụng phổ biến hiện nay ở miền Tây Nam Bộ như OM4218, IR50404, OM551 và MTL480. Hàm lượng cacbon trong rơm rạ trung bình trong nghiên cứu này là 43,8%, trong đó giống lúa có hàm lượng C cao nhất lên tới 49,9% (Bảng 3.1). Từ hàm lượng C bị đốt cháy có thể tính được hệ số phát thải của CO2, là thông số quan trọng để tính tất cả các hệ số phát thải trong các thí nghiệm đốt hở tại đồng theo phương pháp cân bằng cacbon đã trình bày ở Chương 2. Hàm lượng cacbon trung bình trong nghiên cứu này nằm trong dải nghiên cứu của Kim Oanh và cộng sự, (2011) (49,0 ± 2,7%).
Nồng độ đóng góp ròng của các chất ô nhiễm từ quá trình đốt rơm rạ được xác định là nồng độ phát thải của các chất đó xác định trong quá trình đốt trừ đi nồng độ trong môi trường nền được xác định trong cùng một ngày tại cùng một vị trí (Bảng 3.1). Kết quả quan trắc đối với TSP và PM10 chỉ được xác định tại Cần Thơ và Vĩnh Long (năm 2019) và không được thực hiện tại An Giang và Hậu Giang (năm 2018) do thời điểm đó nghiên cứu hạn chế về thiết bị quan trắc. Ngoài ra, trong nghiên cứu này thông số CO cũng được lấy mẫu nhưng do trong quá trình vận chuyển mẫu gặp sự cố kĩ thuật nên mẫu này bị loại bỏ khỏi kết quả nghiên cứu.
bằng máy gặt đập liên hợp sẽ được rải đều trên mặt ruộng và phơi khô tự nhiên dưới ánh nắng mặt trời sau 2-5 ngày sẽ được người nông dân đốt ngay tại ruộng. Tuy nhiên, tùy thuộc vào khoảng thời gian lưu giữ hay thời gian từ khi gặt đến khi đốt mà độ ẩm giữa các loại rơm sẽ khác nhau. Hơn nữa tùy thuộc vào điều kiện thời tiết của từng địa điểm khác nhau sẽ cho độ ẩm của rơm là khác nhau.
Thời gian của 3 giai đoạn của quá trình đốt trong nghiên cứu này cũng tương tự như các nghiên cứu trước đây, bao gồm: bắt cháy, cháy bùng phát và cháy âm ỉ. Thời gian của quá trình cháy bùng phát từ 47-75 phút (trung bình là 57 phút - 1 tiếng). Giai đoạn bắt cháy từ 5s-5 phút là bắt đầu lấy mẫu. Kết thúc quá trình cháy âm ỉ là khoảng 2-3 tiếng. Tổng thời gian cho 1 đám cháy khoảng 2-5 tiếng. Do đó số lượng mẫu cho mỗi đợt nghiên cứu không lấy được nhiều (2 cỡ mẫu/ngày). Thí nghiệm được thực hiện thực tế theo đám cháy của người dân của địa phương, thời gian đốt thường tiến hành từ 2-5 ngày sau thu hoạch.
Bảng 3.1. Đặc tính của quá trình cháy và đóng góp ròng của các chất ô nhiễm
KHM % Đóng góp ròng (net contribution) Độ ẩm rơm rạ C trong rơm rạ C trong tro CO2, mg/m3 SO2 µg/m3 NO2 µg/m3 PM2,5 µg/m3 PM10 µg/m3 TSP µg/m3 F.HG1 22 34,4 12,9 207 340 225 2655 - - F.HG2 21,3 35,9 13,9 191 155 244 2740 - - F.HG3 20 40,2 13 245 224 279 2310 - - F.AG1 25,4 47,4 12,5 272 331 232 2383 - - F.AG2 32,3 48,0 13 289 266 251 2513 - - F.AG3 37,2 44,2 16,6 242 315 323 2656 - - F.VL1 34,1 34,1 14,2 273 339 287 2566 3175 3416 F.VL2 22,4 45,1 12 308 291 264 2482 3080 3506 F.VL3 31 46,9 10,1 250 297 253 2302 2534 3450 F.VL4 30 46,9 15,1 230 389 314 2294 3052 3171 F.CT1 40,1 48,1 18,2 278 190 334 1980 2255 2319 F.CT2 22,4 48,1 22 243 270 328 2938 3086 3595 % Đóng góp ròng (net contribution)
KHM Độ ẩm rơm rạ C trong rơm rạ C trong tro CO2, mg/m3 SO2 µg/m3 NO2 µg/m3 PM2,5 µg/m3 PM10 µg/m3 TSP µg/m3 F.CT3 21 49,9 13,1 344 304 182 2592 2853 3244 TB 27,6 43,8 14,3 259,3 285,5 270,5 2493,2 2862,1 3243
3.1.1.2. Tỷ lệ phát thải và hiệu suất cháy
Tỷ lệ phát thải trung bình của các chất ô nhiễm từ mỗi đám cháy và giá trị trung bình của các thí nghiệm đốt được trình bày trong Bảng 3.2. Vì NO2, SO2 và CO2 chủ yếu được thải ra trong quá trình cháy có ngọn lửa nên các giá trị đóng góp ròng (net contribution) của NO2 và SO2 được tính toán theo lượng đóng góp ròng của CO2 [134]. Hơn nữa, trong tất cả các thí nghiệm đốt rơm rạ, hiệu suất cháy dao động từ 89,9% tới 94,2% đều được đốt rải trên mặt ruộng sau 2-5 ngày phơi khô tự nhiên tại miền Tây Nam Bộ. Điều này chứng tỏ rằng quá trình cháy có ngọn lửa chi phối chính trong tất cả các đám cháy ngoài đồng ruộng. Điều đó có nghĩa là CO2 sẽ được sử dụng làm chất tham chiếu trong quá trình tính toán hệ số phát thải của các chất ô nhiễm khác. Kết luận này có ý nghĩa rất quan trọng trong việc xác định hệ số phát thải các chất ô nhiễm được trình bày trong mục tiếp theo.
Bảng 3.2. Tỷ lệ phát thải của các chất ô nhiễm và hiệu suất cháy
KH SO2/CO2 NO2/CO2 PM2,5/CO2 PM10/CO2 TSP/CO2 MCE (%)
F.HG1 0,002 0,001 0,013 - - 93,1 F.HG2 0,001 0,001 0,014 - - 92,6 F.HG3 0,001 0,001 0,009 - - 94,2 F.AG1 0,001 0,001 0,009 - - 91,7 F.AG2 0,001 0,001 0,009 - - 91,1 F.AG3 0,001 0,001 0,011 - - 90,0 F.VL1 0,001 0,001 0,009 0,012 0,012 90,2 F.VL2 0,001 0,001 0,008 0,010 0,011 93 F.VL3 0,001 0,001 0,009 0,010 0,014 90,7
F.VL4 0,002 0,001 0,010 0,013 0,014 89,9 F.CT1 0,001 0,001 0,007 0,008 0,008 90,7 F.CT2 0,001 0,001 0,012 0,013 0,015 92,3 F.CT3 0,001 0,001 0,008 0,008 0,009 91,4 TB 0,001 0,001 0,01 0,01 0,01 91,6 3.1.2. Hệ số phát thải của PM
Hệ số phát thải trung bình của PM (EFTSP, EFPM10 và EFPM2,5) từ đốt rơm rạ ngoài đồng ruộng ở Miền Tây Nam Bộ được trình bày trong Bảng 3.3. Có thể thấy từ Bảng 3.3, các giá trị EFPM (EFTSP, EFPM10, EFPM2,5) trong nghiên cứu này tương tự như những kết quả đã được công bố với cùng điều kiện đốt rải ở Thái Lan [116] và ở Mỹ [135]. Tuy nhiên, kết quả của nghiên cứu này thấp hơn so với kết quả ở miền Bắc Việt Nam [70]. Điều này có thể được giải thích do sự khác biệt về điều kiện khí