Tổng quan luận án

Luận án tiến sĩ kỹ thuật môi trường với đề tài "Nghiên cứu mức độ và phần đóng góp của các dạng nguồn thải chính tới nồng độ bụi nano trong không khí" được thực hiện bởi nghiên cứu sinh Nguyễn Thị Thu Thủy, dưới sự hướng dẫn khoa học của PGS.TS Nghiêm Trung Dũng tại Trường Đại học Bách khoa Hà Nội (bảo vệ năm 2019, mã số chuyên ngành: 9520320). Công trình sử dụng một phần kết quả của đề tài nghiên cứu khoa học cấp Bộ Giáo dục và Đào tạo mã số B2016-BKA-28.

Tính cấp thiết và khoảng trống nghiên cứu

Bụi nano, hay bụi siêu mịn (nanoparticles / ultrafine particles), được định nghĩa là các hạt bụi có đường kính khí động học (aerodynamic diameter - AED) nhỏ hơn hoặc bằng 100 nm (PM0,1). Khi hít vào đường dẫn khí, bụi nano có khả năng lắng đọng sâu tại tiểu phế quản và phế nang theo cơ chế khuếch tán, thay vì cơ chế va chạm quán tính hay lắng đọng trọng lực như các dải bụi kích thước lớn (PM10, PM2,5). Với diện tích bề mặt riêng lớn, các hạt này hấp phụ nhiều hợp chất độc hại như hydrocacbon thơm đa vòng (PAHs), kim loại nặng và vi sinh vật; có khả năng xâm nhập qua biểu mô, nội mô vào hệ tuần hoàn, hệ bạch huyết, tủy xương, lách, tim và hệ thần kinh trung ương. Khi tích điện, mức độ tích tụ trong phổi của bụi nano tăng từ 5 đến 6 lần so với hạt không tích điện cùng kích thước. Về môi trường, bụi nano làm thay đổi tính chất quang học khí quyển, ảnh hưởng tới sự hình thành mây và gây cưỡng bức bức xạ.

Tại Việt Nam đến thời điểm năm 2014, các nghiên cứu kiểm soát ô nhiễm bụi ngoài trời chủ yếu tập trung vào tổng bụi lơ lửng (TSP), PM10 và PM2,5; hoàn toàn chưa có công bố nào về nồng độ khối lượng, đặc trưng hóa lý và nguồn phát thải của bụi nano. Đến năm 2017, chỉ có hai công bố khảo sát nồng độ số lượng hạt PM0,1 trong thời gian ngắn (từ 6 đến 94 giờ) tại trường học và tòa nhà cao tầng. Trên phạm vi quốc tế, các cơ quan quản lý như US EPA và WHO vẫn chưa ban hành quy chuẩn nồng độ bụi nano do thiếu các dữ liệu quan trắc dài hạn và nghiên cứu dịch tễ tương ứng. Đề tài được triển khai nhằm giải quyết khoảng trống dữ liệu thực nghiệm và phương pháp luận nhận dạng nguồn thải PM0,1 tại đô thị Việt Nam.

Mục tiêu nghiên cứu

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Xác định mức độ bụi nano trong không khí và đặc trưng lý hóa của chúng.
  2. Nhận dạng và xác định phần đóng góp của các dạng nguồn chính tới nồng độ bụi nano trong không khí.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Bụi nano trong môi trường không khí ngoài trời và mô hình nơi tiếp nhận nhân tố hóa ma trận dương (Positive Matrix Factorization - PMF).
  • Phạm vi không gian: Thực hiện tại hai điểm quan trắc ở thành phố Hà Nội:
    1. Khuôn viên Đại học Bách khoa Hà Nội (HUST): Đặt trên sân thượng tầng 3 Trung tâm Ngoại ngữ (tọa độ 21°00.37'B, chiều cao đầu lấy mẫu 8 m, cách trục đường chính 200–300 m), đại diện cho điểm chịu tác động hỗn hợp của các nguồn thải đô thị (giao thông, dân sinh).
    2. Tòa nhà Vinacomin: Đặt trên sân thượng tầng 5 tại đường Nguyễn Văn Cừ, quận Long Biên (tọa độ 21°03.04'B, chiều cao đầu lấy mẫu 18 m, cách mép đường 15 m), đại diện cho điểm chịu tác động trực tiếp của giao thông trục chính, bến xe, ga tàu và hoạt động đốt rơm rạ ngoại thành.
  • Phạm vi thời gian: Quan trắc trong năm 2015, bao gồm mùa mưa (tháng 8/2015) và mùa khô (tháng 10 đến tháng 12/2015).
  • Phạm vi nhận dạng nguồn thải: Áp dụng mô hình PMF cho tập dữ liệu đo đạc liên tục 56 ngày trong mùa khô (tháng 11 và tháng 12/2015) tại vị trí HUST.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Thuật ngữ "bụi siêu mịn" được sử dụng trong các nghiên cứu khí quyển từ thập niên 1950 và phát triển mạnh từ thập niên 1980; thuật ngữ "bụi nano" phổ biến từ thập niên 1990 gắn với khoa học vật liệu và y học. Trong quản lý chất lượng không khí, US EPA từng bước quy chuẩn hóa từ TSP (1971), PM10 (1987) đến PM2,5 (1997); Liên minh châu Âu đưa ra quy chuẩn PM10 năm 1999 và PM2,5 năm 2005.

Các công trình điểm luận được tác giả hệ thống hóa theo các nhóm nội dung:

  • Nồng độ số lượng hạt (number concentration): Morawska và cộng sự tổng hợp 71 nghiên cứu, ghi nhận nồng độ trung bình tăng dần theo tính chất khu vực: nền sạch (2,6 × 10³ hạt/cm³), nông thôn (4,8 × 10³ hạt/cm³), nền đô thị (7,3 × 10³ hạt/cm³), đô thị (10,8 × 10³ hạt/cm³), đường hẻm (42,1 × 10³ hạt/cm³), lề đường (48,2 × 10³ hạt/cm³), trên đường (71,5 × 10³ hạt/cm³) và đường hầm (167,7 × 10³ hạt/cm³). Kumar và cộng sự (2014) khảo sát tại 42 thành phố chỉ ra rằng nồng độ số lượng ven đường tại châu Á cao nhất thế giới (trung bình 1,17 ± 1,04 × 10⁵ hạt/cm³, trong đó Delhi đạt 3,08 × 10⁵ hạt/cm³, Thượng Hải đạt 1,20 × 10⁵ hạt/cm³), trong khi châu Âu thấp nhất (3,15 ± 1,60 × 10⁴ hạt/cm³). Tại Atlanta (Hoa Kỳ), các hạt dưới 100 nm chiếm 89% về nồng độ số lượng dù chỉ chiếm hơn 20% khối lượng.
  • Nồng độ khối lượng (mass concentration): Thống kê từ các nghiên cứu dao động từ 0,49 đến 4,6 µg/m³; tại các thành phố lớn của Hoa Kỳ (Los Angeles, San Francisco, Houston, Miami, New York), nồng độ mùa hè thường vượt 2 µg/m³.
  • Thành phần hóa học: Cass và cộng sự nghiên cứu tại 7 thành phố phía bắc California (Hoa Kỳ) xác định cacbon hữu cơ (OC) chiếm 50% khối lượng PM0,1 (dao động 32–67%), oxit kim loại chiếm 14% (chủ yếu Fe, Ti, Cr, Zn, Ce), cacbon nguyên tố (EC) 8,7%, sunfat 8,2%, nitrat 6,8%, amoni 3,7%, Na⁺ 0,6%, Cl⁻ 0,5%. Hughes và Cass (1998) tại Los Angeles ghi nhận tổng OC và EC chiếm 46–62%. Tỷ lệ áp đảo của OC được củng cố bởi Geller và cộng sự (2002), Kim và cộng sự (2002), Sardar (2005). Trái lại, Miguel và cộng sự (2004) thấy EC chiếm 43% và OC chiếm 24%, nhận diện các hợp chất phenanthrene/fluoranthene và benzo[a]anthracene/indeno[1,2,3-cd]pyrene từ khí thải động cơ. Herner và cộng sự (2005, 2006) ghi nhận 98% khối lượng hạt nano là cacbon đen (BC) và OC. Fine và cộng sự (2004) dùng levoglucosan và cholesterol làm chất chỉ thị cho thấy đun nấu sinh hoạt đóng góp không đáng kể tại bắc California. Tại Houston, silic chiếm tới 70% ở dải hạt 35–50 nm. Chen và cộng sự (2009) tại Đài Loan so sánh tỷ lệ các thành phần tại ba môi trường:
Thành phần (%) Ven đường Đường hầm Khu vực rừng
Cacbon hữu cơ (OC) 38,7 31,4 30,4
Cacbon nguyên tố (EC) 16,1 27,8 13,0
Ion hòa tan 23,3 3,0 11,5
Nguyên tố kim loại 11,3 11,6 7,9
Thành phần khác 10,6 26,2 37,2
  • Ảnh hưởng sức khỏe và khí hậu: Báo cáo của Viện Nghiên cứu Sức khỏe Hoa Kỳ (HEI, 2013) tổng hợp 86 công trình dịch tễ học ngắn hạn (11 nghiên cứu tử vong, 15 tim mạch cấp, 17 hô hấp, 43 tim mạch nói chung), phân bố chủ yếu ở Tây Âu và Bắc Mỹ, khu vực châu Á chỉ có 5 công trình. Về khí hậu, bụi trong khí quyển tạo cưỡng bức bức xạ chung là -1,2 W/m² (trực tiếp -0,5 W/m², gián tiếp -0,7 W/m²), song hạt muội than khi phủ thêm sulfate hoặc chất hữu cơ có thể chuyển thành cưỡng bức dương lên tới +0,6 W/m².
  • Dự án quan trắc tiêu biểu: Dự án UFIPOLNET (châu Âu, 2004–2008) tại Stockholm, Dresden, Prague, Augsburg; dự án kiểm kê phát thải mode2sec của BAAQMD (Hoa Kỳ, 2012); mạng lưới quan trắc bụi nano Đông Á (EA Nanonet, thành lập 2013) gồm Nhật Bản, Hàn Quốc, Trung Quốc, Thái Lan, Việt Nam, Campuchia và Singapore.

Vị trí và khoảng trống luận án lựa chọn giải quyết

Luận án được định vị trong khuôn khổ mạng lưới EA Nanonet. Luận án khắc phục tình trạng thiếu hụt dữ liệu thực nghiệm về nồng độ khối lượng và thành phần hóa học chi tiết của bụi nano tại Việt Nam, đồng thời tiên phong thiết lập bộ hồ sơ nguồn thải thực nghiệm bằng mô hình toán thống kê đa biến PMF kết hợp phân tích hướng gió CPF và mô hình quỹ đạo khối khí HYSPLIT.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết mô hình nơi tiếp nhận

Mô hình nơi tiếp nhận hoạt động dựa trên định luật bảo toàn khối lượng. Tổng khối lượng mẫu $x$ tại nơi tiếp nhận là tổng đóng góp của $p$ nguồn riêng rẽ $g_k$: $$x = \sum_{k=1}^{p} g_k$$ Nồng độ khối lượng của cấu tử $j$ trong mẫu thứ $i$ ($x_{ij}$) biểu diễn qua ma trận: $$x_{ij} = \sum_{k=1}^{p} g_{ik} f_{kj}$$ trong đó $g_{ik}$ là tỷ trọng đóng góp của nguồn $k$ vào mẫu $i$; $f_{kj}$ là phần đóng góp của cấu tử $j$ trong hồ sơ nguồn $k$.

Luận án phân tích hạn chế của hai phương pháp truyền thống:

  • Mô hình cân bằng khối lượng hóa học (CMB): Đòi hỏi phải biết trước số lượng nguồn và cấu trúc hồ sơ nguồn $f_{kj}$. Do Việt Nam chưa có bộ số liệu nguồn cho hạt nano nên không thể áp dụng CMB.
  • Phân tích thành phần chính (PCA/FA): Dựa trên phân tích giá trị riêng SVD ($X = GF + E = U \bar{S} \bar{V}' + E$), có thể tạo ra giá trị âm không phù hợp thực tế môi trường, nhạy cảm với việc chuẩn hóa dữ liệu, không xử lý được giá trị khuyết (MV) hoặc giá trị dưới ngưỡng phát hiện (BDL).

Mô hình nhân tố hóa ma trận dương (PMF)

Mô hình PMF (Paatero và Tapper, 1993) phân tách ma trận quan trắc $X(n \times m)$ thành hai ma trận không âm $G(n \times p)$, $F(p \times m)$ và ma trận phần dư $E(n \times m)$ thông qua tối thiểu hóa hàm mục tiêu $Q(E)$: $$Q(E) = \sum_{i=1}^{n} \sum_{j=1}^{m} \left[ \frac{e_{ij}}{s_{ij}} \right]^2 = \sum_{i=1}^{n} \sum_{j=1}^{m} \left[ \frac{x_{ij} - \sum_{k=1}^{p} g_{ik} f_{kj}}{s_{ij}} \right]^2$$ với điều kiện ràng buộc $g_{ik} \ge 0$ và $f_{kj} \ge 0$, trong đó $s_{ij}$ là độ không đảm bảo đo của cấu tử $j$ trong mẫu $i$. Thuật toán PMF2 kiểm soát điều kiện không âm thông qua việc cực tiểu hóa hàm mở rộng $\bar{Q}(E, G, F)$ chứa các hàm phạt $P(G), P(F)$ và hàm chuẩn hóa $R(G), R(F)$. Nghiên cứu sử dụng phần mềm EPA-PMF phiên bản 5.0 của Cục Bảo vệ Môi trường Hoa Kỳ.

Các công cụ mô hình khí tượng bổ trợ

  • Hàm xác suất có điều kiện (CPF): Xác định hướng của các nguồn thải cục bộ dựa trên dữ liệu khí tượng bề mặt: $$CPF = \frac{m_\theta}{n_\theta}$$ trong đó $m_\theta$ là số mẫu có nồng độ nguồn vượt ngưỡng phân vị quy định tại hướng gió $\theta$; $n_\theta$ là tổng số mẫu tại hướng gió $\theta$.
  • Mô hình quỹ đạo lùi HYSPLIT: Tính toán quỹ đạo lùi 3 ngày (backward trajectory) của khối khí lan truyền tới Hà Nội ở độ cao 300–500 m nhằm đánh giá khả năng vận chuyển tầm xa.

Phương pháp thực nghiệm và thiết bị quan trắc

  • Thu mẫu bụi nano: Sử dụng bộ va chạm quán tính nhiều tầng do Đại học Kanazawa (Nhật Bản) chế tạo (Model 26A57C), lưu lượng hút 40 L/phút. Thiết bị gồm 4 tầng va chạm thu các dải hạt kích thước > 10 µm; 2,5–10 µm; 1,0–2,5 µm và 0,5–1,0 µm; một cartridge chứa sợi thép không gỉ (SUS) thu hạt 0,1–0,5 µm; và tầng cuối cùng thu bụi nano (< 0,1 µm) trên giấy lọc quartz đường kính 55 mm (Pallflex 2500 QAT-UP, Mỹ). Thiết bị Model 25A137F được vận hành song song để kiểm chuẩn chất lượng (QA/QC).
  • Thu mẫu bụi PM2,5 và PM10: Thiết bị Sibata LV40B (Nhật Bản, đầu xyclon, lưu lượng 16,7 L/phút) lấy mẫu PM2,5; thiết bị Minivol (Airmetrics, Hoa Kỳ, lưu lượng 5,0 L/phút) lấy mẫu PM10. Cả hai dùng giấy lọc quartz đường kính 47 mm (Pallflex).
  • Đo trực tiếp phân bố số lượng: Thiết bị phân tích kích thước hạt quét linh động điện Nanoscan (SMPS TSI Model 3910), dải đo 13 kênh kích thước từ 10 đến 420 nm, lưu lượng 0,6 L/phút.
  • Xử lý giấy lọc: Giấy lọc quartz được nung tại 900°C trong 4 giờ để loại bỏ cacbon hữu cơ nền trước khi cân và lấy mẫu.

Kế hoạch lấy mẫu năm 2015

Số lượng mẫu được tổng hợp trong bảng số liệu thực nghiệm:

Đợt lấy mẫu Vị trí Thời gian (năm 2015) Mẫu PM0,1 (Thực / Trắng) Mẫu PM2,5 (Thực / Trắng) Mẫu PM10 (Thực / Trắng) Đặc trưng thời tiết
Đợt 1 HUST 04/08 – 20/08 17 / 3 16 / 3 16 / 3 Mùa mưa, 8 ngày mưa, hướng E, SE
Đợt 2 Vinacomin 22/08 – 28/08 7 / 3 7 / 3 7 / 3 Mùa mưa, 2 ngày mưa, hướng E, SE
Đợt 3 Vinacomin 13/10 – 19/10 7 / 3 7 / 3 7 / 3 Mùa khô, vụ đốt rơm rạ, hướng NE, NW
Đợt 4 HUST 04/11 – 30/12 56 / 7 38 / 7 39 / 7 Mùa khô, 23 ngày mưa nhỏ, hướng NE, E
Tổng cộng 87 / 16 68 / 16 69 / 16

Phân tích trong phòng thí nghiệm

  • Khối lượng: Cân vi lượng tại Phòng Kiểm chuẩn Môi trường, Trung tâm Quan trắc Môi trường Miền Bắc (Tổng cục Môi trường).
  • Thành phần cacbon (OC, EC): Phân tích bằng phương pháp nhiệt - quang học tại Đại học Saitama (Nhật Bản), phân tách thành phần char-EC và soot-EC.
  • Thành phần ion vô cơ: Phân tích 9 ion hòa tan trong nước (SO₄²⁻, NO₃⁻, Cl⁻, NH₄⁺, Na⁺, K⁺, Ca²⁺, Mg²⁺...) bằng sắc ký ion (IC) tại Đại học Saitama.
  • Thành phần nguyên tố: Định lượng 22 nguyên tố kim loại và á kim bằng phương pháp khối phổ plasma cảm ứng cao tần (ICP-MS) tại Đại học Georgia Gwinnett và Viện Công nghệ Georgia (Hoa Kỳ).

Nội dung chính theo từng chương

Chương 1. Tổng quan về bụi nano

Chương này tổng kết cơ sở lý thuyết về động học sol khí, các tiêu chuẩn môi trường và hiện trạng nghiên cứu trong nước và quốc tế.

  • Đặc tính vật lý và cơ chế lắng đọng: Bụi khí quyển được phân đoạn thành bụi siêu mịn (AED ≤ 0,1 µm), bụi mịn (0,1–2,5 µm), bụi thô (2,5–10 µm) và bụi siêu thô (> 10 µm). Bụi nano có thời gian lưu khí quyển từ vài phút đến vài ngày; cơ chế loại bỏ chủ yếu là khuếch tán và va chạm kết tụ. Về mặt độc học, xác suất lắng đọng của hạt bụi trong phế nang đạt cực đại tại khoảng kích thước 20–30 nm. Tốc độ đào thải bụi nano khỏi phế nang đạt khoảng 20%, thấp hơn nhiều so với mức 80% của các hạt thô lớn hơn 500 µm.
  • Tác động khí hậu: Bụi nano hoạt động như hạt nhân ngưng tụ mây (CCN). Trong môi trường đô thị ô nhiễm, sự dư thừa hạt nano làm tăng mật độ hạt mây nhưng làm giảm kích thước giọt, khiến mây đặc hơn, tăng độ phản xạ bức xạ sóng ngắn mặt trời, ức chế giáng thủy và giảm khả năng tự làm sạch của khí quyển.
  • Tổng quan phương pháp mô hình hóa nguồn thải: So sánh chi tiết nguyên lý, cấu trúc toán học, điều kiện biên và giới hạn ứng dụng của ba mô hình nơi tiếp nhận: CMB, PCA và PMF; khẳng định tính ưu việt của thuật toán bình phương tối thiểu có trọng số và ràng buộc không âm trong PMF khi xử lý số liệu môi trường thực tế.

Chương 2. Phương pháp nghiên cứu

Chương này mô tả chi tiết quy trình thu thập mẫu hiện trường, quy chuẩn kỹ thuật lắp đặt trạm đo, vận hành thiết bị và quy trình phân tích trong phòng thí nghiệm.

  • Quy chuẩn bố trí trạm lấy mẫu: Tuân thủ tiêu chuẩn vị trí quan trắc theo chỉ dẫn 2008/50/EC: chiều cao đầu lấy mẫu từ 2 đến 15 m; cách xa tường chắn ít nhất 2 m; góc mở không bị che chắn tối thiểu 180°; khoảng cách tới cây xanh tối thiểu 10–20 m.
  • Vận hành hệ thống va chạm tầng Kanazawa: Lưu lượng dòng khí qua hệ thống được kiểm soát ở mức 40 L/phút để duy trì đúng vận tốc cắt tách các phân đoạn hạt tại từng tầng va chạm quán tính.
  • Quy trình kiểm soát chất lượng (QA/QC): Đánh giá độ chệch giữa hai máy lấy mẫu nano Model 26A57C và Model 25A137F; kiểm tra tương quan nồng độ tổng khối lượng các chất phân tích được (hóa học) so với khối lượng cân vi lượng trực tiếp; thiết lập đường chuẩn phân tích sắc ký ion và phổ khối ICP-MS; phân tích mẫu trắng hiện trường song song để hiệu chỉnh nền.

Chương 3. Kết quả và bàn luận

Chương này trình bày toàn bộ số liệu đo đạc thực nghiệm, quy luật biến thiên theo không gian - thời gian, tương quan giữa các dải hạt và kết quả mô hình hóa nguồn thải.

  • Độ tin cậy của chuỗi dữ liệu thực nghiệm: Kết quả so sánh song song giữa thiết bị Model 26A57C và Model 25A137F cho thấy hệ số tương quan tuyến tính cao ($R^2$ đạt từ 0,85 đến 0,95 đối với nồng độ khối lượng và các chỉ tiêu cacbon, ion). Tổng khối lượng các hợp phần hóa học phân tích được (OC, EC, ion, nguyên tố) tương thích với khối lượng bụi xác định bằng cân vi lượng. Cân bằng điện tích giữa tổng cation và tổng anion thể hiện tính trung hòa điện tích của mẫu bụi.
  • Đặc trưng nồng độ khối lượng và phân bố số lượng:
    • Nồng độ khối lượng bụi nano có sự phân hóa rõ rệt giữa mùa khô và mùa mưa: nồng độ mùa khô tại HUST và Vinacomin cao hơn đáng kể so với mùa mưa do điều kiện nghịch nhiệt bức xạ và tần suất mưa rửa trôi thấp.
    • Tỷ lệ khối lượng giữa bụi nano với các phân đoạn bụi lớn hơn được lượng hóa cụ thể: tỷ lệ PM0,1/PM2,5 và PM0,1/PM10 thay đổi phụ thuộc vào điều kiện khí tượng và nguồn phát thải cục bộ.
    • Nồng độ số lượng đo bằng thiết bị Nanoscan SMPS 3910 phản ánh sự ưu thế vượt trội của các hạt có đường kính nhỏ hơn 50 nm về mặt số lượng hạt, đối lập hoàn toàn với phân bố nồng độ khối lượng (nơi phân đoạn hạt thô và hạt tích tụ chiếm ưu thế).
    • Phân tích quỹ đạo lùi HYSPLIT chỉ ra rằng khi các khối khí di chuyển từ lục địa phía Bắc và Đông Bắc qua các vùng công nghiệp phát triển tràn về Hà Nội trong mùa khô, nồng độ các phân đoạn bụi tăng cao đột biến.
  • Đặc trưng thành phần hóa học của PM0,1:
    • Cacbon hữu cơ (OC) và cacbon nguyên tố (EC): OC luôn là thành phần chiếm tỷ trọng lớn nhất trong khối lượng bụi nano. Tỷ lệ OC/EC được xác định chi tiết cho từng mùa và từng vị trí; nồng độ cacbon hữu cơ thứ cấp (SOC) được ước tính thành công dựa trên phương pháp tỷ lệ (OC/EC) nhỏ nhất. Sự tương quan giữa char-EC và soot-EC khẳng định sự hiện diện đồng thời của khí thải phương tiện giao thông và đốt sinh khối/than.
    • Thành phần ion hòa tan: Các ion chiếm ưu thế tuyệt đối là SO₄²⁻, NH₄⁺ và NO₃⁻, chứng minh sự đóng góp của các phản ứng hóa học thứ cấp pha khí - hạt. Mối tương quan chặt chẽ giữa SO₄²⁻ và NH₄⁺ cho thấy sự tồn tại chủ yếu dưới dạng muối (NH₄)₂SO₄ và NH₄HSO₄. K⁺ xuất hiện rõ nét trong các đợt quan trắc trùng với vụ đốt rơm rạ ngoại thành (Đợt 3 tại Vinacomin).
    • Thành phần nguyên tố: 22 nguyên tố được định lượng chính xác ở đơn vị ng/m³. Các nguyên tố có nguồn gốc vỏ Trái Đất (Al, Fe, Ca, Mg) có tỷ lệ thấp hơn hẳn so với trong bụi PM10, trong khi các kim loại vết có nguồn gốc đốt nhiên liệu và công nghiệp (Pb, Zn, Cu, Ni, Cr, Cd) tập trung đáng kể trong phân đoạn hạt nano.
  • Kết quả nhận dạng và xác định phần đóng góp nguồn thải bằng mô hình PMF:
    • Áp dụng trên tập số liệu mùa khô tại trạm HUST (56 mẫu 24 giờ, Đợt 4). Dựa trên phân tích các chỉ số độ không đảm bảo (IM, IS), các tham số xoay trục (Fpeak), sự suy giảm hàm mục tiêu $Q$ và tương quan giữa nồng độ mô phỏng với quan trắc, mô hình xác lập tối ưu 6 nhân tố nguồn thải chính:
      1. Nguồn thứ cấp (Secondary aerosol): Đặc trưng bởi hàm lượng cao của các ion SO₄²⁻, NH₄⁺, NO₃⁻ và một phần OC thứ cấp.
      2. Nguồn động cơ chạy xăng (Gasoline vehicles): Đặc trưng bởi tỷ lệ OC cao, các kim loại vết chỉ thị phụ gia nhiên liệu và tỷ lệ soot-EC thấp.
      3. Nguồn động cơ diesel (Diesel vehicles): Đặc trưng bởi hàm lượng soot-EC chiếm ưu thế cùng với EC, Pb, Zn.
      4. Nguồn đốt sinh hoạt và dịch vụ (Residential/service combustion): Đặc trưng bởi ion K⁺, char-EC và OC phát sinh từ đun nấu dân sinh, bếp than tổ ong và đốt chất thải sinh hoạt.
      5. Nguồn công nghiệp (Industrial emissions): Chứa hàm lượng đặc trưng của các kim loại nặng và nguyên tố vết công nghiệp (Fe, Zn, Cu, Pb, Mn, Cr).
      6. Nguồn bụi đất và bụi đường tái cuốn (Soil/Road dust): Nhận diện qua các nguyên tố khoáng vật vỏ Trái Đất (Al, Ca, Fe, Mg).
    • Kết hợp hàm xác suất có điều kiện (CPF) với hoa gió thực nghiệm cho phép định vị hướng không gian phát thải chủ yếu của từng nguồn tới vị trí tiếp nhận HUST.

Kết quả và những đóng góp mới

Đóng góp mới về phương pháp luận

  • Luận án là công trình đầu tiên tại Việt Nam nghiên cứu và thiết lập hoàn chỉnh quy trình quan trắc thực nghiệm bụi nano ngoài trời: từ tiêu chuẩn lựa chọn vị trí trạm đo, vận hành thiết bị phân tách hạt va chạm nhiều tầng (Cascade impactor Kanazawa), quy chuẩn hóa xử lý nhiệt giấy lọc quartz, đến quy trình phân tích tích hợp đồng bộ các chỉ tiêu hóa lý siêu vết (cân vi lượng, IC, nhiệt - quang học OC/EC, ICP-MS).
  • Là công trình đầu tiên ứng dụng thành công mô hình nơi tiếp nhận đa biến nhân tố hóa ma trận dương (PMF 5.0) kết hợp với hàm xác suất có điều kiện (CPF) để định lượng tỷ lệ đóng góp của từng nguồn thải tới nồng độ bụi nano PM0,1 tại môi trường đô thị Việt Nam.

Đóng góp mới về mặt dữ liệu thực nghiệm và thực tiễn

  • Thiết lập bộ dữ liệu thực nghiệm đồng bộ đầu tiên tại Hà Nội với 6 dải kích thước hạt (> 10 µm; 2,5–10 µm; 1,0–2,5 µm; 0,5–1,0 µm; 0,1–0,5 µm và < 0,1 µm), bao gồm nồng độ khối lượng, nồng độ số lượng, thành phần 9 ion vô cơ, 22 nguyên tố kim loại và các phân đoạn cacbon (OC, EC, char-EC, soot-EC, SOC).
  • Lượng hóa các tỷ lệ khối lượng thực nghiệm giữa PM0,1 với PM2,5 và PM10 theo mùa và theo tính chất vị trí tiếp nhận.
  • Phân lập thành công 6 dạng nguồn thải chính đóng góp vào nồng độ bụi nano tại khu vực nội thành Hà Nội trong mùa khô (thứ cấp, phương tiện chạy xăng, phương tiện chạy diesel, đốt sinh hoạt - dịch vụ, công nghiệp, bụi đất - đường).
  • Bước đầu xây dựng bộ hồ sơ nguồn thải (source profiles) đầu tiên cho bụi nano tại Hà Nội, cung cấp dữ liệu cơ sở cho các mô hình tính toán ô nhiễm không khí tiếp theo khi thiếu điều kiện lấy mẫu dài ngày.

Đề xuất và kiến nghị

  • Sử dụng các tương quan tỷ lệ khối lượng giữa PM0,1/PM2,5 và chỉ số chỉ thị OC/EC, char-EC/soot-EC để ước tính sơ bộ mức độ phơi nhiễm bụi nano từ mạng lưới trạm đo tự động PM2,5 sẵn có.
  • Cung cấp luận cứ khoa học định hướng chính sách kiểm soát chất lượng không khí đô thị: ưu tiên kiểm soát khí thải phương tiện giao thông (đặc biệt là động cơ diesel), cấm triệt để bếp than tổ ong, kiểm soát hoạt động đốt rơm rạ ngoại thành và xây dựng lộ trình quy chuẩn kỹ thuật quốc gia đối với bụi siêu mịn tại Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế của luận án:
    • Do hạn chế về kinh phí, chi phí phân tích phòng thí nghiệm tại nước ngoài và thời gian vận hành thiết bị, việc chạy mô hình PMF để nhận dạng nguồn thải mới chỉ thực hiện được cho một mùa (mùa khô) tại một vị trí trạm đo (khuôn viên HUST, 56 mẫu 24 giờ).
    • Chưa thiết lập được chuỗi quan trắc liên tục đa năm để đánh giá xu thế biến đổi dài hạn của bụi nano tại nhiều dạng sinh cảnh đô thị khác nhau (khu công nghiệp, khu đô thị mới, vùng nông thôn lân cận).
    • Các đánh giá tác động sức khỏe trong luận án dừng lại ở mức tổng thuật tài liệu dịch tễ quốc tế, chưa triển khai nghiên cứu độc học hoặc dịch tễ học thực địa trực tiếp trên quần thể dân cư Hà Nội.
  • Hướng nghiên cứu tiếp:
    • Mở rộng quan trắc bụi nano liên tục cho cả mùa mưa và các thời điểm chuyển mùa; tăng cường mật độ điểm quan trắc trên các vùng địa lý khác nhau của Việt Nam.
    • Tiếp tục cập nhật, hiệu chỉnh và hoàn thiện bộ hồ sơ nguồn thải bụi nano phục vụ chạy mô hình CMB và PMF ổn định hơn.
    • Phối hợp liên ngành y tế - môi trường triển khai các nghiên cứu phơi nhiễm ngắn hạn và dài hạn nhằm đánh giá rủi ro sức khỏe cộng đồng do bụi nano gây ra.

Giá trị tham khảo

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên:
    • Cung cấp tài liệu tham khảo về phương pháp luận quan trắc hạt kích thước nano, nguyên lý kỹ thuật của thiết bị va chạm kiểu tầng Kanazawa và quy trình tiền xử lý giấy lọc quartz nung nhiệt độ cao.
    • Cung cấp quy trình chi tiết thiết lập ma trận nồng độ và ma trận độ không đảm bảo (uncertainty matrix) để xử lý các số liệu khuyết (MV) và dưới ngưỡng phát hiện (BDL) trong mô hình EPA-PMF v5.0.
    • Tài liệu tham khảo về phương pháp phân tách nhiệt - quang học xác định char-EC, soot-EC và ước tính cacbon hữu cơ thứ cấp (SOC).
  • Đối với cơ quan quản lý và hoạch định chính sách:
    • Giúp Tổng cục Môi trường, Sở Tài nguyên và Môi trường Hà Nội nhận diện mức độ đóng góp của từng nhóm phát thải (giao thông chạy xăng/diesel, sinh hoạt, công nghiệp, thứ cấp) để ưu tiên nguồn lực quản lý.
    • Cung cấp cơ sở khoa học để các nhà làm chính sách nghiên cứu, chuẩn bị luận cứ phục vụ việc xem xét xây dựng quy chuẩn kỹ thuật quốc gia về giới hạn nồng độ bụi nano trong không khí xung quanh trong tương lai.

Câu hỏi thường gặp

1. Bụi nano được định nghĩa như thế nào và cơ chế lắng đọng sinh học có gì khác biệt so với bụi PM2,5?

Bụi nano (PM0,1) là những hạt bụi khí quyển có đường kính khí động học nhỏ hơn hoặc bằng 100 nm (0,1 µm). Trong khi bụi PM2,5 chủ yếu lắng đọng tại khí quản, phế quản nhờ cơ chế va chạm quán tính và lắng trọng lực, bụi nano di chuyển theo cơ chế khuếch tán Brown, có khả năng đi sâu và lắng đọng với xác suất cao nhất tại vùng tiểu phế quản và phế nang (đạt đỉnh tại dải kích thước 20–30 nm). Bụi nano có tốc độ làm sạch sinh học chậm hơn (chỉ khoảng 20% bị đào thải khỏi phế nang), có thể hấp thu vào biểu mô, nội mô để đi vào hệ tuần hoàn máu và hệ bạch huyết tới các cơ quan nội tạng.

2. Vì sao luận án lựa chọn mô hình PMF thay vì CMB hay PCA để xác định nguồn thải bụi nano?

Mô hình CMB đòi hỏi phải biết trước số lượng nguồn và cấu trúc hồ sơ nguồn thải ($f_{kj}$), trong khi tại Việt Nam chưa có bất kỳ hồ sơ nguồn thực nghiệm nào cho hạt nano. Mô hình PCA có nhược điểm tạo ra kết quả âm (phi thực tế môi trường), nhạy cảm với việc chuẩn hóa dữ liệu và không xử lý được các mẫu khuyết (MV) hoặc nồng độ dưới giới hạn phát hiện (BDL). Mô hình PMF khắc phục toàn bộ các nhược điểm trên nhờ thuật toán bình phương tối thiểu có trọng số từng điểm (dựa trên ma trận độ không đảm bảo $s_{ij}$), tích hợp ràng buộc toán học không âm ($G \ge 0, F \ge 0$) và không đòi hỏi hồ sơ nguồn ban đầu.

3. Thiết bị lấy mẫu bụi nano sử dụng trong nghiên cứu có nguyên lý và cấu tạo ra sao?

Thiết bị lấy mẫu bụi nano (Model 26A57C) do Đại học Kanazawa (Nhật Bản) phát triển, hoạt động với lưu lượng dòng khí 40 L/phút. Thiết bị gồm 4 tầng va chạm kiểu tầng phía trên sử dụng các vòi phun giảm dần diện tích tiết diện để tăng vận tốc dòng khí, qua đó tách và thu các hạt bụi có kích thước > 10 µm; 2,5–10 µm; 1,0–2,5 µm và 0,5–1,0 µm theo nguyên lý va chạm quán tính. Tiếp theo là cartridge chứa màng lọc sợi kim loại không gỉ (SUS) thu dải hạt 0,1–0,5 µm. Tầng đáy cuối cùng thu gom các hạt bụi nano (< 0,1 µm) trên giấy lọc quartz đường kính 55 mm.

4. Kế hoạch lấy mẫu thực nghiệm tại Hà Nội được triển khai tại những địa điểm nào và số lượng mẫu ra sao?

Nghiên cứu bố trí lấy mẫu tại 2 địa điểm: điểm hỗn hợp đô thị HUST (tầng 3 Trung tâm Ngoại ngữ ĐH Bách khoa Hà Nội, cao 8 m) và điểm giao thông/đốt sinh khối Vinacomin (tầng 5 tòa nhà Vinacomin đường Nguyễn Văn Cừ, Long Biên, cao 18 m). Tổng số mẫu thu thập trong năm 2015 gồm 87 mẫu thực PM0,1 (16 mẫu trắng), 68 mẫu thực PM2,5 (16 mẫu trắng) và 69 mẫu thực PM10 (16 mẫu trắng). Việc lấy mẫu chia làm 4 đợt, trong đó đợt quan trắc phục vụ mô hình PMF được tiến hành liên tục 56 ngày (từ 04/11 đến 30/12/2015) trong mùa khô tại vị trí HUST.

5. Mô hình PMF đã tách được những dạng nguồn thải chính nào đóng góp vào bụi nano tại Hà Nội trong mùa khô?

Từ bộ số liệu 56 mẫu mùa khô tại HUST, mô hình EPA-PMF 5.0 đã nhận diện và định lượng được 6 nguồn thải chính đóng góp vào nồng độ khối lượng bụi PM0,1 gồm:

  1. Nguồn aerosol thứ cấp (đặc trưng bởi SO₄²⁻, NH₄⁺, NO₃⁻).
  2. Nguồn khí thải động cơ chạy xăng (đặc trưng bởi tỷ lệ OC cao và các kim loại vết).
  3. Nguồn khí thải động cơ diesel (đặc trưng bởi hàm lượng soot-EC cao).
  4. Nguồn đốt sinh hoạt và dịch vụ dân sinh (đặc trưng bởi K⁺, char-EC, OC).
  5. Nguồn khí thải sản xuất công nghiệp (đặc trưng bởi Fe, Zn, Cu, Pb, Mn, Cr).
  6. Nguồn bụi đất và bụi đường tái cuốn (đặc trưng bởi các nguyên tố khoáng Al, Ca, Fe, Mg).

Kết luận

Luận án của nghiên cứu sinh Nguyễn Thị Thu Thủy đã giải quyết hoàn chỉnh bài toán thực nghiệm về xác định mức độ ô nhiễm, thành phần hóa học và nhận diện các nguồn đóng góp chính của bụi nano (PM0,1) tại Hà Nội. Bằng việc kết hợp hệ thống lấy mẫu phân tầng Kanazawa, các kỹ thuật phân tích hóa lý hiện đại và mô hình receptor PMF/CPF, công trình đã xác lập bộ số liệu cơ sở đồng bộ đầu tiên tại Việt Nam về bụi PM0,1 với 6 nguồn thải chính được lượng hóa trong mùa khô. Kết quả nghiên cứu cung cấp luận cứ khoa học phục vụ quản lý chất lượng môi trường không khí đô thị và định hướng xây dựng quy chuẩn kỹ thuật quốc gia về bụi nano trong tương lai.