Trong đó, ERx và ER(tracer gas) lần lượt là tỷ lệ phát thải của chất ô nhiễm và của chất đánh dấu, trong khi [X] và [tracer gas] là nồng độ chất khí được đo.
Từ ERX thu được, sau đó có thể được chuyển đổi thành EFX bằng cách nhân nó với khoảng thời gian thí nghiệm hoặc chia nó cho tiết diện bề mặt của vùng thực hiện thí nghiệm.
Như vậy, giống phương pháp cân bằng cacbon, phương pháp này có ưu điểm là cho kết quả xác định EF tương đối chính xác, phản ảnh được điều kiện cháy ngoài hiện trường. Tuy nhiên, hạn chế chính của phương pháp này là giả định về động lực vận chuyển của các chất ô nhiễm là bằng nhau (thông qua các dòng đối lưu). Điều này có thể tạo ra những sai số nhất định so với thực tế [90].
Một số tác giả đã thực hiện nghiên cứu so sánh giữa việc xác định hệ số phát thải bằng kỹ thuật sử dụng chất đánh dấu khí quyển và sử dụng mô hình AERMOD (trên cùng một thí nghiệm). Kết quả cho thấy, EFs thu được bằng phương pháp sử dụng chất đánh dấu khí quyển có độ biến thiên thấp hơn không đáng kể so với phương pháp sử dụng mô hình [108].
Vì vậy, phương pháp sử dụng chất đánh dấu khí quyển có thể thay thế cho các phương pháp khác đã được đề cập, có tiềm năng mang lại kết quả tốt với chi phí và
( ) [ ] [ ] tracergas x X ER ER tracergas
39
kỹ thuật thực hiện không quá khó khăn. Với cách tiếp cận này, phương pháp này hoàn toàn khả thi trong điều kiện nghiên cứu của Việt Nam.
1.1.3.5. Phương pháp biên dạng thẳng đứng
Phương pháp biên dạng thẳng đứng dựa trên các phép đo hiện trường về tốc độ gió và nồng độ các chất ô nhiễm trong khí quyển từ hoạt động đốt rơm rạ ngoài đồng ruộng để suy ra biên dạng thẳng đứng của tốc độ gió và biên dạng thẳng đứng của nồng độ chất ô nhiễm (phân bố nồng độ chất ô nhiễm theo chiều cao) [109-112]. Có thể thu được biên dạng tốc độ gió bằng phương trình biên dạng của gió dạng logarit, với máy đo tốc độ gió dùng âm thanh 3D hoặc bằng cách đo tốc độ gió ở hai độ cao khác nhau [90]. Biên dạng nồng độ chất ô nhiễm theo chiều cao thu được bằng cách đo nồng độ chất ô nhiễm ở bốn độ cao khác nhau với thiết bị đo thẳng đứng. Độ cao được chọn phụ thuộc vào khoảng cách của thiết bị đo tới vùng phát thải. Hệ số phát thải của chất ô nhiễm sau đó được tính toán theo phương trình sau [113]
(1.2)
Trong đó, EFs là hệ số phát thải của chất ô nhiễm (mg/m2), z là chiều cao so với mặt đất (m), z0 là chiều dài độ nhám bề mặt, u(z) là tốc độ gió trung bình ở độ cao z (m/s) trong suốt quá trình đo, c(z) là nồng độ trung bình của chất cần quan tâm ở độ cao z (mg/m3), t là khoảng thời gian thí nghiệm, θ là góc giữa sức căng ma sát ở bề mặt và gió trong khí quyển, w là khoảng rộng gió dịch chuyển so với mặt đất trong thời gian thử nghiệm và zmax là chiều cao mà tại đó nồng độ chất ô nhiễm cần quan tâm được tính là 0.
40
Ưu điểm của phương pháp này cho phép tính toán EFs dựa hoàn toàn vào các phép đo tại hiện trường nên kết quả tương đối chính xác. Tuy nhiên, phương pháp biên dạng thẳng đứng cũng tồn tại một số hạn chế như sau:
- Cần rất nhiều dụng cụ thiết yếu để thực hiện các phép đo nồng độ và tốc độ gió ở các độ cao khác nhau;
- Ước tính phân bố nồng độ theo phương thẳng đứng, chiều cao luồng khói và phân bố nồng độ tốc độ gió là không cố định vì nó dựa trên các phép đo không trễ;
- Đối với khoảng cách để đo PM theo chiều gió, tỷ lệ PM2,5/PM10 giữa phép đo phát thải gần nguồn (PM2,5/PM10 khoảng 50%) khác hẳn phép đo phát thải xa nguồn (PM2,5/PM10 khoảng 10%). Sự khác biệt này có thể là do phần bụi mịn hơn (PM2,5) có xu hướng phân tán theo chiều dọc, khiến cho việc phát hiện trong các phép đo nồng độ xa nguồn trở nên khó khăn hơn [109].
Trong điều kiện nghiên cứu ở Việt Nam việc áp dụng phương pháp này là khó khăn do hạn chế về thiết bị đo và sai số trong quá trình đo không trễ.
1.2.3.6. Phương pháp sử dụng công nghệ LIDAR
Trong những năm gần đây, công nghệ khảo sát từ xa sử dụng nguồn kích thích bằng tia laser, công nghệ LIDAR (Light Detection And Ranging) đã và đang được các nhà khoa học nghiên cứu thử nghiệm để nghiên cứu quá trình phát thải PM từ các hoạt động nông nghiệp, đặc biệt để xác định được các thông số phát tán từ một đám cháy [113, 114]. Nguyên lý hoạt động của một hệ LIDAR tương tự nguyên lý hoạt động của một hệ radar, bao gồm một khối phát bức xạ điện từ kích thích và một khối thu tín hiệu tán xạ ngược. Về cơ bản hệ LIDAR có cấu trúc gồm một khối phát tia laser hướng về đối tượng cần quan trắc và một khối thu tín hiệu tán xạ ngược trở lại. Bức xạ laser hướng về phía đối tượng nghiên cứu, tương tác với đối tượng cần khảo
41
sát, bức xạ điện từ sẽ biến đổi tính chất trước khi trở về đầu thu. Bức xạ điện từ tán xạ trở về đầu thu sẽ mang các thông tin về đối tượng khảo sát, tuân theo lý thuyết tán xạ tùy thuộc vào bản chất của đối tượng tán xạ. Sự thay đổi tính chất của bức xạ trở về cho phép xác định các thông số đặc trưng của môi trường nghiên cứu như: đặc trưng tán xạ ngược, mật độ, sự phân bố, hình dạng và kích thước PM. Trên cơ sở đó, người ta có thể ước tính được ER của chất ô nhiễm cần quan tâm thông qua việc hiệu chỉnh các mẫu lọc thu được từ công nghệ LIDRA. Những ứng dụng này cũng cho phép đánh giá mức độ không chắc chắn của các ước tính chiều cao của khói bằng phương pháp biên dạng thẳng đứng [113]. Tương tự, công nghệ LIDAR cũng đã và đang được sử dụng để đánh giá mức độ không chắc chắn của các ước tính độ cao của thông số khói được thực hiện bằng các mô hình. Hệ số phát thải được ước tính bằng công nghệ LIDAR có khoảng không chắc chắn nhỏ hơn không đáng kể so với phương pháp sử dụng mô hình AERMOD [115].
Một nghiên cứu khác liên quan đến việc sử dụng công nghệ LIDAR trong việc quan sát và xác định phát thải từ một hoạt động cháy ngoài đồng ruộng cho thấy, dưới những điều kiện đối lưu, khói có xu hướng di chuyển theo phương thẳng đứng hơn là phương nằm ngang. Những phát hiện này rất có ý nghĩa trong việc lựa chọn vị trí đặt thiết bị xác định PM trên mặt đất nên đặt gần hay xa nguồn phát thải thì cho kết quả chính xác hơn. Một ưu điểm nữa của ứng dụng công nghệ LIDAR gần đây là có thể phân biệt các sol khí từ những nguồn khác nhau như phát thải từ động cơ, bụi đất hay từ hoạt động đốt [109]. Trong những năm gần đây, công nghệ LIDAR đã và đang trở thành một công cụ quan trọng để ước tính EFs của PM, đặc biệt là trong những thử nghiệm thí điểm như là một phương pháp tham chiếu để đánh giá tính chính xác của các mô hình. Nhược điểm chính của công nghệ này liên quan đến chi phí và độ phức
42
tạp của nó trong quá trình sử dụng thiết bị và hiệu chỉnh. Mặt khác, kỹ thuật này là một kỹ thuật cung cấp nhiều thông tin nhất về hình dạng và động học của khói bụi.
Trong khuôn khổ đề tài này, phương pháp cân bằng cacbon được áp dụng trong việc xác định hệ số phát thải các chất ô nhiễm không khí từ hoạt động đốt rơm rạ ngoài đồng ruộng thuộc bốn tỉnh miền Tây Nam Bộ là: An Giang, Hậu Giang (năm 2018) và Vĩnh Long, Cần Thơ (2019).
1.2.4. Tổng hợp nghiên cứu xác định hệ số phát thải của hoạt động đốt rơm rạ
Sử dụng hai cơ sở dữ liệu là Web of Science cho các ấn phẩm tiếng Anh và Cơ sở dữ liệu quốc gia tại Việt Nam để tổng hợp các nghiên cứu xác định hệ số phát thải từ hoạt động đốt rơm rạ. Kết quả tổng hợp có khoảng 284 bài nghiên cứu đã được công bố và trích dẫn để phân tích dữ liệu [91]. EFs của các chất ô nhiễm không khí điển hình, bao gồm CO2, CO, PM2.5, SO2, NOx và PAHs, được tổng hợp trong phần Phụ lục 1. EFs được tính bằng đơn vị đo khối lượng chất ô nhiễm trên khối lượng nhiên liệu đốt.
(A) Số lượng các công trình công bố EF hàng năm (B) Số lượng các nghiên cứu xác định EF của các quốc gia
Nguồn: [91]
43
Từ Hình 1.3 có thể thấy trong một vài năm gần đây, số lượng các nghiên cứu xác định hệ số phát thải tại đồng ruộng, tuy có tăng về số lượng nhưng vẫn còn hạn chế tại một vài quốc gia điển hình như: Trung Quốc, Ấn Độ, Nepal, Malawi, Benin, Rwanda, Zambia, Bangladesh, Kenya, Thái Lan và Việt Nam [91]. Trong các nghiên cứu tại đồng ruộng kể trên, hầu hết được thực hiện ở một làng hoặc một số làng trong một vùng và có quy mô mẫu nhỏ, ngoại trừ một vài nghiên cứu xuyên quốc gia. Các nghiên cứu kết hợp giữa phòng thí nghiệm và hiện trường hiện nay còn rất hạn chế. Tại Thái Lan, nghiên cứu xác định hệ số phát thải các chất ô nhiễm không khí từ hoạt động đốt rơm rạ được thực hiện bằng phương pháp cân bằng cacbon trên cả thí nghiệm tại đồng ruộng và thí nghiệm trong phòng thí nghiệm [32, 116]. Đối với nghiên cứu này, quan trắc hiện trường áp dụng đối với phương thức đốt rải được triển khai.
Tuy nhiên, như phần trên đã trình bày, hệ số phát thải của chất ô nhiễm không khí từ hoạt động đốt rơm rạ phụ thuộc vào nhiều yếu tố, như giống lúa, đặc điểm thành phần tính chất của rơm rạ, phương thức thu hoạch và phương thức đốt, độ ẩm của nhiên liệu, đặc điểm khí hậu khu vực đốt. Do đó, việc xác định một bộ hệ số phát thải đặc trưng cho điều kiện khí hậu riêng của mỗi quốc gia là điều cần thiết. Tại Việt Nam, gần đây một nghiên cứu được thực hiện tại miền Bắc bởi nhóm tác giả Pham.C.T và cộng sự (2021), trong đó các hệ số phát thải tại đồng ruộng và phòng thí nghiệm cũng được thực hiện. Kết quả nghiên cứu này đã xác định được hệ số phát thải của TSP, PM2,5, CO, CO2, SO2, NO2, 10 PAHs và không nghiên cứu hệ số phát thải của VOCs, một nhóm hợp chất quan trọng phát sinh từ hoạt động đốt rơm rạ trong việc ảnh hưởng đến khí quyển. Trong nghiên cứu này, rơm rạ sau thu hoạch được chất thành đống và tiến hành đốt.
Trong khi đó, có sự khác biệt đáng kể về khí hậu giữa miền Bắc và miền Tây Nam Bộ của Việt Nam. Ở miền Bắc có mùa đông lạnh với nhiệt độ trung bình dưới
44
20°C [117] trong khi ở Miền Tây Nam Bộ nắng nóng quanh năm với nhiệt độ trung bình hàng năm khoảng 20-35°C . Sự khác biệt này dẫn đến một số yếu tố khác nhau, bao gồm cả lịch thời vụ khác nhau giữa hai vùng. Miền Bắc chỉ có hai vụ lúa/năm, trong khi đó ở Miền Tây Nam Bộ một năm có ba vụ lúa/năm, thậm chí bảy vụ lúa chao hai năm. Do đó, thời gian thu hoạch, phương thức thu hoạch, và điều kiện không khí xung quanh (nhiệt độ và độ ẩm) cũng khác nhau. Hơn nữa, diện tích trồng lúa và cách thức thực hiện đốt rơm rạ giữa hai vùng cũng không giống nhau. Ở miền Bắc, rơm rạ chất thành từng đống nhỏ rồi đốt trong khi đó ở miền Tây Nam Bộ thì rơm rạ được phơi khô tự nhiên, rải trên ruộng rồi đốt. Miền Tây Nam Bộ được coi là vùng trồng lúa lớn nhất cả nước và cũng có tỷ lệ đốt rơm rạ cao nhất [9, 10], các nghiên cứu về đốt rơm rạ ngoài đồng ruộng ở vùng này vẫn còn rất khan hiếm. Theo hiểu biết tốt nhất của chúng tôi, một số công trình được tìm thấy trong nguồn cơ sở dữ liệu mở. Một nghiên cứu như vậy đã sử dụng phương pháp đốt trong chụp hút để xác định hệ số phát thải của các khí nhà kính như CH4 và N2O từ quá trình đốt rơm rạ ở Miền Tây Nam Bộ [118]. Do đó, nghiên cứu này được thiết kế để tiến hành các thí nghiệm đốt rơm rạ tại hiện trường (đốt hở) khu vực miền Tây Nam Bộ nhằm đánh giá phát thải các chất ô nhiễm không khí từ hoạt động đốt rơm rạ.
1.3. Kiểm kê phát thải từ hoạt động đốt hở sinh khối
1.3.1. Khái niệm chung
Kiểm kê phát thải (KKPT) là một công cụ hữu ích để định lượng tải lượng ô nhiễm tại một khu vực cụ thể tại một thời điểm cụ thể [119]. KKPT có thể được thiết kế và tiến hành ở quy mô địa phương, khu vực và quốc gia bằng cách sử dụng các phương pháp khác nhau tùy thuộc vào mục đích, tốc độ phát thải và sự sẵn có của dữ liệu đầu vào [120].
Kiểm kê phát thải của quá trình đốt sinh khối quy mô khu vực đến toàn cầu được tính bằng cách nhân tổng lượng sinh khối bị đốt cháy với các hệ số phát thải
45
tương ứng. Nhìn chung, việc ước tính phát thải chất ô nhiễm không khí nói chung và từ quá trình đốt sinh khối nói riêng có thể được tiến hành bằng hai cách tiếp cận: (i) từ trên xuống (top-down) và (ii) từ dưới lên (bottom-up) [3, 121]. Sự khác biệt giữa hai cách tiếp cận là tính toán lượng sinh khối bị đốt cháy.
1.3.2. Phương pháp tiếp cận từ trên xuống (top-down)
Cách tiếp cận từ trên xuống thường đề cập đến việc chia nhỏ một hệ thống tổng quan để có được cái nhìn sâu sắc về các hệ thống thành phần, nó giống như thiết kế ngược. Trong trường hợp kiểm kê, đó là việc sử dụng các quan sát vệ tinh trong không khí, thường kết hợp với mô hình khí quyển, để ước tính lượng phát thải. Cách tiếp cận từ trên xuống cung cấp một phương pháp đáng tin cậy để định lượng đốt sinh khối một cách nhất quán trên các khu vực rộng lớn, chẳng hạn như đốt rừng và đốt xavan hay đồng cỏ. Dữ liệu về diện tích của đám cháy (m2), chủ yếu được lấy từ dữ liệu MODIS [122, 123]. Dữ liệu ảnh MODIS có các kênh phổ có bước sóng trong khoảng cho phép tính toán được nhiệt độ của bề mặt đất, do đó các đám cháy ngoài thực tế lúc sắp hoặc mới xảy ra đều phát sinh nhiệt rất lớn, sensor của vệ tinh thu các ảnh trong khu vực có cháy và sử dụng các thuật toán để đoán đọc ảnh cho phép xác định khu vực có xảy ra hoạt động đốt cháy. Đối với hoạt động đốt hở rơm rạ ngoài đồng ruộng thì diện tích các đám cháy thường ở quy mô nhỏ (khoảng 1ha ở miền Tây Nam Bộ [124]) và ở các vị trí phân tán, nên khó đánh giá tốt qua vệ tinh. Do đó, mặc dù cách tiếp cận từ trên xuống cho kết quả kiểm kê tổng thể, nhanh gọn nhưng nó cũng có thể có những sai số nhất định và không phản ánh chính xác điều kiện địa phương. Do vậy, khi áp dụng cách tiếp cận này cho việc kiểm kê phát thải từ hoạt động đốt hở rơm rạ ngoài đồng ruộng người ta thường phải xác định độ không chắc chắn của phương pháp, để từ đó có những giá trị tính toán có độ tin cậy cao.
46
1.3.3. Phương pháp tiếp cận từ dưới lên (bottom-up)
Nếu như phương pháp tiếp cận từ trên xuống tiếp cận tính toán phát thải từ dữ liệu tổng hợp từ trên xuống thông qua quan sát vệ tinh thì phương pháp tiếp cận từ dưới lên thực hiện kiểm kê phát thải dựa trên tính toán chi tiết về lượng phát thải từ tất cả các nguồn riêng lẻ của một khu vực, sau đó được tổng hợp lại để thu được lượng phát thải từ tổng khu vực. Trong cách tiếp cận từ dưới lên, tổng lượng sinh khối bị đốt cháy được xác định dựa trên dữ liệu thống kê (ví dụ: sản lượng cây trồng, tỷ lệ chất thải trên sản lượng sinh khối, tỷ lệ sinh khối khô, tỷ lệ đốt, v.v.). Việc thu thập