PM2.5 làm gia tăng tử vong do ung thư hệ hô hấp tại thành phố Hồ Chí Minh năm 2018

Nghiên cứu nhằm xây dựng bản đồ phân bố PM2.5 và đánh giá tác động của PM2.5 đến tử vong do ung thư phổi và ung thư hệ hô hấp tại TPHCM năm 2018. Bản đồ PM2.5 được xây dựng dựa trên quan trắc 96 điểm vào mùa mưa và khô bằng thiết bị AirBeam2 và thuật toán nội suy IDW.

TẠP CHÍ NGHIÊN CỨU Y HỌC PM2.5 LÀM GIA TĂNG TỬ VONG DO UNG THƯ HỆ HÔ HẤP TẠI THÀNH PHỐ HỒ CHÍ MINH NĂM 2018 Nguyễn Trường Viên1, Nguyễn Ngọc Nhật Thanh1, Phan Hoàng Thùy Dung1, Nguyễn Đào Thiên Ân1, Trương Thị Thùy Dung1, Đinh Thị Giang1, Trần Ngọc Đăng1, Đại học Y Dược Thành phố Hồ Chí Minh Nghiên cứu nhằm xây dựng đồ phân bố PM2.5 đánh giá tác động PM2.5 đến tử vong ung thư phổi ung thư hệ hô hấp TPHCM năm 2018 Bản đồ PM2.5 xây dựng dựa quan trắc 96 điểm vào mùa mưa khơ thiết bị AirBeam2 thuật tốn nội suy IDW Phân tích tác động PM2.5 đến tử vong ung thư phổi hệ hô hấp dựa liệu tử vong A6/YTCS Bộ Y tế mơ hình BenMAP Ở tất quận/huyện nồng độ PM2.5 trung bình năm cao tiêu chuẩn an tồn sức khỏe WHO (PM2.5 < 10 µg/m3) PM2.5 đóng góp 6,3% (35/557) tử vong ung thư phổi 6,5% (41/629) tử vong ung thư hệ hơ hấp Mỗi 10 µg/m3 PM2.5 tăng thêm số ca tử vong ung thư phổi gia tăng thêm 56 ca ung thư hệ hô hấp 64 ca Cần có biện pháp ngắn hạn phòng tránh tác hại PM2.5 biện pháp dài hạn giảm thiểu nhiễm PM2.5 Từ khóa: PM2.5, ung thư hệ hô hấp, tử vong I ĐẶT VẤN ĐỀ Ơ nhiễm khơng khí yếu tố nguy lớn sức khỏe, đặc biệt ô nhiễm khơng khí ngồi trời Trong đó, phơi nhiễm chất dạng hạt kích thước nhỏ < 2,5μm (PM2.5) xếp thứ nguy gây tử vong toàn cầu vào năm 2018.1 PM2.5 hạt bụi có đường kính 2.5 micromet chứa chất rắn giọt chất lỏng siêu nhỏ mang chất.2, Với kích thước nhỏ lơ lửng khơng khí, PM2.5 dễ dàng xâm nhập vào sâu phổi gây kích thích, ăn mòn thành phế nang dẫn đến suy giảm chức phổi làm tăng nguy bệnh đường hô hấp, sau thời gian dài tiếp xúc dẫn đến ung thư đường hô hấp, đặc biệt ung thư phổi.2, Theo Song cộng sự, giai đoạn 2014-2016 PM2.5 đóng góp 23,9% số ca tử vong ung thư phổi 15,5% cho tất trường hợp tử vong.5 Theo Báo cáo Đo lường gánh nặng bệnh tật toàn cầu, phơi nhiễm Tác giả liên hệ: Trần Ngọc Đăng Đại học Y Dược Thành phố Hồ Chí Minh Email: ngocdangytcc@gmail.com Ngày nhận: 01/03/2021 dài hạn với PM2.5 đóng góp 16% số ca tử vong ung thư phổi năm 2016.6 Chất lượng khơng khí Việt Nam từ năm 2010 đến mức báo động có diễn tiến theo chiều hướng xấu cách so sánh kết đo quan trắc bụi PM2.5 Tại TPHCM, liệu từ WHO cho thấy mức độ ô nhiễm hạt mịn cao hàng đầu Việt Nam sau Hà Nội Nồng độ hạt mịn PM2.5 trung bình năm 20167 42 µg/m3 cao nhiều lần so với tiêu chuẩn WHO với giới hạn 10 µg/m3.8 Trước tình hình nhiễm khơng khí, đặc biệt nhiễm PM2.5 hữu, người dân TPHCM đối diện với nguy cao sức khỏe liên quan đến ung thư phổi hệ hô hấp Tuy nhiên, chưa có đánh giá tác động nhiễm PM2.5 đến sức khỏe, cụ thể tác động PM2.5 lên gánh nặng tử vong ung thư phổi hệ hơ hấp Do đó, chúng tơi thực nghiên cứu với mục tiêu đánh giá tác động PM2.5 đến tử vong ung thư phổi ung thư hệ hô hấp người dân sống TPHCM năm 2018 Ngày chấp nhận: 22/04/2021 108 TCNCYH 142 (6) - 2021 TẠP CHÍ NGHIÊN CỨU Y HỌC II ĐỐI TƯỢNG VÀ PHƯƠNG PHÁP Đối tượng Dữ liệu quan trắc Nồng độ PM2.5, nhiệt độ, độ ẩm quan trắc trực tiếp nhóm nghiên cứu sau tập huấn lựa chọn địa điểm kỹ thuật quan trắc Tại quận/huyện, điểm quan trắc lựa chọn theo hướng dẫn chọn mẫu xếp hạng Cơ quan bảo vệ mô trường Hoa Kỳ.9 Giao thông xác định nguồn phát thải ô nhiễm quan trọng TPHCM, đó, nghiên cứu chọn điểm đo cách lề 5m đoạn đường ứng với mức mật độ giao thông vào cao điểm (7h00 18h00) theo liệu Google Map Ứng với 24 quận/huyện, có tổng cộng 96 điểm quan trắc tồn TPHCM Q trình quan trắc chia thành hai đợt ứng với mùa mưa (tháng 12 - năm sau) mùa khô (tháng - 11) Trong đợt, điểm quan trắc ngày gồm ngày tuần ngày cuối tuần ô nhiễm thường khác biệt ngày tuần với cuối tuần, ngày lần (từ 6h00 8h00 17h00 - 19h00), lần kéo dài 2h Thiết bị AirBeam2 sử dụng để quan trắc nghiên cứu Đây thiết bị giá rẻ với mức giá tham khảo $249 USD (khoảng triệu đồng Việt Nam), có sẵn thị trường, có độ tin cậy cao so sánh với thiết bị quan trắc tự động cố định tiêu chuẩn với hệ số tương quan đạt 0,84.10 Dữ liệu tử vong Hồi cứu liệu thống kê nguyên nhân tử vong năm 2018 theo sổ A6/YTCS Bộ Y tế Nguyên nhân tử vong mã hóa lại dựa phân loại bệnh tật quốc tế ICD-10 Dựa mã ICD-10, trường hợp ung thư phổi (C33 - C34.9, C39.8, C45.7), ung thư hệ hô hấp xác định (C11, C32, C33 - C34.9, C39.8, C45.7) Những ca báo cáo tử vong trước phải có năm liên tục sinh sống địa TCNCYH 142 (6) - 2021 phương Các ca tử vong thiếu số kiện tuổi, giới, ngày tử vong, nguyên nhân tử vong, địa (quận) bị loại khỏi nghiên cứu Dữ liệu dân số, địa lý Dữ liệu tổng dân số cấp độ quận/huyện năm 2018 thu thập dựa theo báo cáo Cục Thống kê TPHCM Dữ liệu địa lý trích xuất lớp địa giới hành từ Trung tâm Ứng dụng Hệ thống Thơng tin Địa lý TPHCM Phương pháp Thiết kế nghiên cứu Nghiên cứu sinh thái đánh giá tác động sức khỏe dựa liệu hồi cứu tử vong, dân số quan trắc PM2.5 sử dụng thiết bị cảm biến giá rẻ AirBeam2 Xử lý số liệu Xây dựng đồ phân bố nồng độ PM2.5: Thuật toán nội suy IDW sử dụng để xây dựng đồ phân bố PM2.5 với liệu đầu vào nồng độ nhiễm trung bình năm 96 điểm đo, tọa độ điểm đo, lớp đồ địa giới Thuật toán nội suy IDW sử d quận huyện TPHCM IDW thực độ ô nhiễm trung bình năm 96 điểm phần mềm ứng dụng cơng nghệ hệ thống thông IDW phần mềm ứn tin địa lý ArcGIS10.3.1 củathực Viện nghiên cứu hệ nghiên hệ công thống môi trường Hoa K thống môi trường Hoa cứu Kỳ 11 với thức: ∑+ %& (&)* 𝑍𝑍" = ∑&,+ )* &,-.(& Trong đó: Trong đó: (cơng thức 1) 𝑍𝑍0: giá trị điểm cần ước tính Z0: giá trị điểm cần ước tính z: giá trị z 𝑧𝑧: mẫu điểm I giátạitrị z mẫu điểm I d: khoảng cách từ điểm mẫu đến điểm ước 𝑑𝑑: khoảng cách từ điểm mẫu tính đến điểm N: hệ số xác định trọng số dựa khoảng N: hệ số xác định trọng số dựa cách n: tổng số n: điểm dự đoán tổng số điểm kho dự đoán Kết đầu IDW phần mềm 109 xác định nồng độ PM2.5 dự đố Ước tính tác động PM2.5 đến TẠP CHÍ NGHIÊN CỨU Y HỌC Kết đầu IDW phần mềm ArcGIS10.3.1 đồ phân bố khơng gian nồng độ PM2.5, từ xác định nồng độ PM2.5 dự đoán tọa độ TPHCM.11, 12 Ước tính tác động PM2.5 đến tử vong ung thư phổi hệ hơ hấp: Mơ hình BenMAP phát triển Cơ quan Bảo vệ Môi trường Hoa Kỳ sử dụng để đánh giá tác động ô nhiễm PM2.5 đến tử vong ung thư phổi ung thư hệ hô hấp với liệu đầu vào lớp đồ nồng độ PM2.5, địa giới quận/huyện, số ca tử vong ung thư phổi ung thư hệ hô hấp quận huyện, tổng dân số cấp quận/huyện Ước tính tác động sức khỏe từ mơ hình BenMAP dựa công thức: 13 Tử vong liên quan đến PM2.5 = A x B x C x D (công thức 2) Trong đó: (A) Nồng độ PM2.5: Chênh lệch nồng độ PM2.5 trung bình quận/huyện dựa kết nội suy IDW so với ngưỡng tham chiếu 10µg/m3 (B) Ước tính ảnh hưởng sức khỏe: Dựa nghiên cứu dịch tễ, tử vong tăng 6,2% nồng độ PM2.5 tăng 10μg/m3 (C) Tỷ suất tử vong nền: Sử dụng tỉ suất tử vong ung thư phổi ung thư hệ hô hấp năm 2018 (D) Dân số phơi nhiễm: Sử dụng dân số trung bình quận huyện năm 2018 từ Cục thống kê TPHCM.13 Kết đầu BenMAP ước tính số ca tử vong liên quan đến ô nhiễm PM2.5 theo quận huyện Đạo đức nghiên cứu Đề tài số 66/2019/HĐ- ĐHYD ngày 10/01/2020 Hội đồng Đạo đức Trường Đại học Y Dược Tp Hồ Chí Minh chấp thuận theo chứng nhận chấp thuận số 01/HĐĐĐ, ngày 06/01/2020 III KẾT QUẢ Dựa liệu quan trắc được, đồ ô nhiễm PM2.5 xây dựng giúp hình dung tình trạng nhiễm khơng khí mùa mưa mùa khô TPHCM Các kết liệu cung cấp đồ phân bố nồng độ PM2.5 kết hợp với liệu tử vong, dân số, địa lý để ước tính tác động PM2.5 lên tử vong ung thư phổi ung thư hệ hô hấp PHỤ LỤC Phân bố nồng độ PM2.5 theo TPHCM theo IDW Hình Bảntrung đồ phân bốbình nồng độ PM2.5 trungmùa bình theo mùa TPHCM theo IDW Biểu đồ Phân bố nồng độ PM2.5 mùa mưa, mùa khơ trung bình năm Hình Bản đồ phân bố nồng độ PM2.5 trung bình theo mùa TPHCM theo IDW 110 TCNCYH 142 (6) - 2021 TẠP CHÍ NGHIÊN CỨU Y HỌC Phân bố nhiễm PM2.5 tồn thành phố có khác biệt lớn quận thay đổi theo mùa mùa mưa thường nhiễm mùa khơ Mặc dù huyện có mức độ thị hóa thấp so với quận, ô nhiễm PM2.5 huyện Củ Chi, Hóc Mơn cao TPHCM mùa khô mùa mưa Phân bố ô nhiễm PM2.5 xu hướng liên quan đến thị hóa (hình 1) Nồng độ PM2.5 trung bình ước tính giao động quận khoảng 10 µg/m3 đến 33,9 µg/m3 Vào mùa mưa nồng độ PM2.5 cao mùa khô ghi nhận 22 số 24 quận huyện Đặc biệt quận Tân Phú, Nhà Bè, Bình Thạnh, Bình Chánh, Bình Tân, Tân Bình, Quận có chênh lệch lớn mùa mưa mùa khô theo xu hướng mùa mưa ô nhiễm cao (biển đồ 1) Quận Huyện Nhà Bè 35 Huyện Hóc Môn 30 Huyện Củ Chi 25 20 Huyện Cần Giờ 15 10 Huyện Bình Chánh Gị Vấp Quận Quận Quận Quận Quận Quận Thủ Đức Quận Tân Phú Quận Bình Tân Tân Bình Phú Nhuận Quận 10 Bình Thạnh Cả năm Quận 11 Quận 12 Mùa khô Mùa mưa Biểu đồ Phân bố nồng độ PM2.5 mùa mưa, mùa khơ trung bình năm Tác động PM2.5 lên tử vong ung thư phổ ung thư hệ hô hấp 2.1 Đặc điểm tử vong TPHCM năm 2018 Bảng Đặc điểm tử vong ung thư phổi ung thư hệ hô hấp TPHCM năm 2018 Đặc điểm Tử vong ung thư phổi Tử vong ung thư hệ hô hấp Tần số % Tần số % Nam 372 66,79 424 67,41 Nữ 185 33,21 205 32,59 Dưới 15 0,00 0,00 Từ 15 - 60 281 50,45 319 50,72 Trên 60 276 49,55 310 49,28 Giới Độ tuổi Thông qua liệu 29.696 trường hợp tử vong năm 2018, ghi nhận 598 trường hợp tử vong ung thư phổi 690 trường hợp tử vong ung thư hệ hơ hấp; đó, số trường hợp đảm bảo tiêu chí đưa vào phân tích 557 (93%) 629 (91%) tương ứng Đặc điểm dân số tử vong ung thư phổi, ung thư hệ hô hấp năm 2018 ghi nhận tỉ lệ nam cao gấp đôi nữ Không ghi nhận trường hợp tử vong 15 tuổi Khoảng 49% ca tử vong có độ tuổi 60 Khơng có khác biệt nhiều tử vong ung thư phổi với tử vong ung thư hệ hô hấp (bảng 1) TCNCYH 142 (6) - 2021 111 TẠP CHÍ NGHIÊN CỨU Y HỌC 2.2 Tỉ suất tử vong Bảng Tỉ suất tử vong ung thư phổi ung thư hệ hô hấp TPHCM năm 2018 Đặc điểm Tử vong ung thư phổi Quy trách PM2.5 Tử vong Tử vong ung thư hệ hô hấp Quy trách PM2.5 Tử vong Tần suất Tỉ suất (/10.000) Tần suất Tỉ suất (/10.000) Tần suất Tỉ suất (/10.000) Tần suất Tỉ suất (/10.000) Quận 22 1,561 0,071 22 1,561 0,071 Quận 18 1,057 0,059 20 1,174 0,059 Quận 21 1,098 0,105 22 1,150 0,105 Quận 11 0,626 0,057 13 0,740 0,057 Quận 13 0,812 0,062 17 1,062 0,062 Quận 31 1,319 0,043 34 1,447 0,085 Quận 0,258 0,000 10 0,286 0,000 Quận 31 0,729 0,024 41 0,964 0,047 Quận 12 0,312 0,026 13 0,338 0,026 Quận 10 15 0,639 0,043 16 0,682 0,043 Quận 11 20 0,946 0,047 27 1,277 0,095 Quận 12 20 0,332 0,017 21 0,349 0,017 Bình Thạnh 36 0,731 0,041 40 0,812 0,041 Phú Nhuận 16 0,975 0,061 18 1,097 0,061 Tân Bình 21 0,448 0,021 21 0,448 0,021 Bình Tân 18 0,233 0,013 30 0,389 0,013 Tân Phú 23 0,481 0,021 24 0,502 0,042 Thủ Đức 38 0,655 0,017 41 0,707 0,017 Gò Vấp 39 0,587 0,030 39 0,587 0,030 Huyện Bình Chánh 34 0,492 0,029 37 0,535 0,029 Huyện Cần Giờ 11 1,569 0,000 13 1,854 0,000 Huyện Củ Chi 68 1,528 0,202 73 1,640 10 0,225 Huyện Hóc Mơn 24 0,452 0,056 29 0,546 0,075 Huyện Nhà Bè 0,310 0,000 0,413 0,000 Toàn TPHCM Dân số (8.831.865) 557 0,631 35 0,040 629 0,712 41 0,046 Quận/Huyện Nghiên cứu ghi nhận tỉ suất tử vong ung thư phổi ung thư hệ hô hấp quy trách PM2.5 0,631/10.000 dân 0,712/10.000 dân Ơ nhiễm PM2.5 đóng góp 6,3% (35/557) số ca tử vong ung thư phổi 6,5% (41/629) số ca tử vong ung thư hệ hô hấp (bảng 2) 112 TCNCYH 142 (6) - 2021 TẠP CHÍ NGHIÊN CỨU Y HỌC Bảng Mối liên quan nồng độ PM2.5 tỉ suất tử vong ung thư phổi ung thư hệ hô hấp Liên quan nồng độ PM2.5 với Tử vong ung thư phổi Liên quan nồng độ PM2.5 với Tử vong ung thư hệ hô hấp p R b (KTC 95%) (/10000 dân) p 0,733 0,007 (0,004 – 0,010) < 0,001 0,752 0,008 (0,005 – 0,011) < 0,001 25 15 05 0 05 05 1 15 15 2 25 15 05 15 25 b (KTC 95%) (/10000 dân) 25 R 15 20 20 25 pm25 25 pm25 30 30 35 35 Tỉ suất tử vong ung phổi trách cho PM2.5 Đường hồi quy Tỉ suất tử vong ung thư phổi quy thư trách choquy PM2.5 Đường hồi quy 15 15 20 20 25 pm25 25 pm25 30 30 35 35 Tỉ suất tử vong dohấp ung quy thư trách hô hấp quy trách cho Đường PM2.5 hồi quy Đường hồi quy Tỉ suất tử vong ung thư hô cho PM2.5 Hệ số quan; tươngbquan; hệ số góc phương trình quy tính; tuyếnKTC tính; KTCKhoảng 95%, Khoảng cậy 95% R, Hệ sốR,tương hệ sốbgóc phương trình hồi quyhồi tuyến 95%, tin cậytin 95% R, Hệ số tương quan; b hệ số góc phương trình hồi quy tuyến tính; KTC 95%, Khoảng tin cậy 95% Có mối tương quan thuận mạnh nồng độ PM2.5 tỉ suất tử vong quy thuộc ung thư mối quan quan mạnh mạnhnồng giữađộ nồng độ PM2.5 tỉtửsuất tửquy vong quy thuộc ung thưvàphổi Cóung mốiCó tương thuận PM2.5 tỉ suất vong ungdo thư phổi phổi thư hệtương hô hấp vớithuận hệ số tương quan 0,733 thuộc 0,752 (bảng 3) Ứng với ung thư hệ hô hấp với hệ số tương quan 0,733 0,752 (bảng 3) Ứng với 10µg/m3 ung thư3 hệ hơ hấp với hệ số tương quan 0,733 0,752 (bảng 3) Ứng với 10µg/m3 10μg/m PM2.5 gia tăng, trực tiếp làm gia tăng tỉ suất tử vong ung thư phổi lên 0,07/10.000 giasẽtăng, trực tiếp vong ung thưlên phổi lên 0,07/10.000 tỉtửsuất tử PM2.5 PM2.5 gia tăng, trựcsẽ tiếp làm gialàm tănggiatỉ tăng suất tỉ tửsuất vongtửdo ungdo thư phổi 0,07/10.000 dân vàdân tỉ suất dân tỉ suất tử vong ung thư hệ hô hấp lên 0,08/10.000 dân Như vậy, với dân số 8.831.865 vong ung thư hệ hô hấp lên 0,08/10.000 dân Như vậy, với dân số 8.831.865 TPHCM năm 2018, vong ung thư hệ hô hấp lên 0,08/10.000 dân Như vậy, với dân số 8.831.865 TPHCM năm 2018, TPHCM năm 2018, độ tăng thêm 10 số catăng tửgia vong ung phổi gia nồng độ PM2.5 tăngnồng thêm 10PM2.5 µg/m3 sốvong ca tửdo vong ung thưsẽ phổi tăngdo thêm 56thư ca ung ung thư nồng độ PM2.5 tăng thêm 10 µg/m3 số ca tử ungμg/m thư phổi gia thêm 56 ca thư hệ hô hấp hệ hô hấp 64 tăng thêm 56 caca và64doca.ung thư hệ hô hấp 64 ca IV BÀN LUẬN BÀN LUẬN BÀN LUẬN Tình trạng nhiễm khơng khí, nhiễm phù hợp điều kiện nguồn lực nghiên Tìnhơtrạng nhiễm khơng khí,tâm nhiễm PM2.5 TPHCM đángCác quan tâmphản phản ánh qua cácđã cứu chế chứng trước PM2.5Tình tạitrạng TPHCM đángkhí, quan nhiễm khơng nhiễm PM2.5 TPHCM rấthạn đáng quan tâm ánh qua kết đo lường trước vượt QCVN 05: 2013/BTNMT ngưỡng an toàn sức khỏe Tổ chức Y tế kết lường trước QCVN 05: 2013/BTNMT ngưỡng an toàn sức khỏe Tổ chức chovàthấy độ tin cậy AirBeam2 Khi Ysotếsánh phản ánhđoqua kết quảvượt đo lường trước 7, 8, 14 7, 8, 14 PM2.5 có phân bố theo khơng gian địa lý, đó, việc lập đồ PM2.5 với độ phân Thế giới (WHO) Thế giới có phân bốvà theo khơng gian địa lý, việc lập PM2.5 với với độ phân đầu dịđó, PMS7003 củađồAirBeam2 thiết bị TSI vượt quá(WHO) QCVN 05:PM2.5 2013/BTNMT ngưỡng giải gian không gian cao tảng nềncơtảng đểgiá đánh nguy cơ.tại, Hiện tại, lưới mạng lướitrắc quan trắc môi giải không cao bảncơ đểbản đánh yếugiá tố yếu nguytốcơ Hiện mạng quan mơi DustTrak DRX Aerosol Monitor 8533 cho thấy an tồntrường sức khỏe Tổ chức Y tế Thế giới vềgiải phân giải khơng gian chỉcộng tổng30 cộng vị trítrắc quanTrong trắc nghiên Trong nghiên trường TPHCMTPHCM hạncòn chếhạn chế phân khơng gian với chỉvới tổng vị trí30quan cứu cứu 7, 8, 14 nồng độ PM2.5 đo tương quan mạnh (WHO) PM2.5 cógiá theođể khơng tơi, chúng tôi,bịthiết bịphân đođược giábố rẻdùng dùng đểtrắc quan 96đo điểm đonăm năm 2020 vàomưa mùavà mưa mùa chúng thiết đo rẻ quan 96trắc điểm 2020 vào mùa mùavà 15 Sửthiết dụngviệc bị rẻ đo giáđồ rẻ PM2.5 AirBeam2 đãđộ cảiđược thiệnmẽ độgiải phân giải không chosựvới thấy phùbiến hợp với R=0,94 Khivàgian so gian lý, đó, với khơ.địa Sửkhơ dụng bịthiết đo lập giá AirBeam2 cải thiện độ phân không gian chosánh thấy phùcảm hợp điều kiện nguồn lực nghiên cứu hạn chế Các chứng trước cho thấy độ tin cậy điều kiện nguồn lực nghiên cứu hạn chế Các chứng trước cho thấy độ tin cậy tiêu chuẩn tham chiếu, AirBeam2 cho hệ phân giải không gian cao tảng để AirBeam2 Khi so sánh đầu dò PMS7003 AirBeam2 vớibịthiết bị TSI DustTrak DRX Aerosol Monitor 8533 AirBeam2 Khi tố so sánh PMS7003 AirBeam2 thiết TSI DustTrak DRX10Aerosol Monitor 8533 tương quan Để tiếp tục cải thiện độ đánh giá yếu nguyđầu cơ.dòHiện tại, mạng lưới với số 15 >0,8 15 cho thấy nồng độ PM2.5 đo tương quan mạnh mẽ với R=0,94 Khi so sánh với cảm biến tiêu chuẩn cho thấy nồng độ PM2.5 đo tương quan mạnh mẽ với R=0,94 Khi so sánh với cảm biến tiêu chuẩn 10 phân giải khơng gian,độthuật tốn nội suy quan trắc mơi trường TPHCM cịnhệhạn chế tham chiếu, AirBeam2 số tương quan10>0,8 cảiđộ thiện phân giải không gian,IDW thuật tham chiếu, AirBeam2 cho hệcho số tương quan >0,8 Để tiếpĐể tụctiếp cảitục thiện phân giải không gian, thuật dùng liệu quan trắc để tính tốn nồng phân giải không gian với tổng cộng 30 vị toán nội suy IDW dùng liệu quan trắc để tính tốn nồng độ khu vực khơng quan trắc toán nội suy IDW dùng liệu quan trắc để tính tốn nồng độ khu vực không quan trắc theo giácàng trị đo gầnchúng điểm dự đốn thìđộ ảnh hưởng Phương pháp có ưutheo có khu vực khơng quan trắc trítheo quan trắc.ngun Trong cứu tơi, ngun tắc giá tắc trịnghiên đo gần điểm dự đốn cócàng ảnh hưởng Phương pháp cónày ưu điểm đơn giản, nhanh chóng, tiết kiệm sử dụng tốt điểm đo có phân bố rộng khắp diện tích tính điểmbịđơn nhanh chóng, tiếtđể kiệm sửtrắc dụng96 tốt điểm đo phân khắp diệnđiểm tích tính ngun tắccógiá trịbố đorộng gần dự đốn thiết đo giản, giá rẻ dùng quan tốn Độ xác IDW đánh giá nghiên cứu năm 2015 Trung Quốc cho thấy có tốn.đo Độtrong chínhnăm xác IDWvào đánh giá nghiênthì cứucàng năm có 2015 hưởng Trung Quốc cho thấy IDWpháp cóIDWnày ảnh Phương điểm 2020 mùa mưa mùa 16 độ xác cao với hệ số tương quan 16 0,99 Như vậy, việc kết hợp quan trắc phạm vi rộng với độ xác cao với hệ số tương quan 0,99 Như vậy, việc kết hợp quan trắc phạm vi rộng với ưu điểm đơn giải, giản, nhanh chóng, tiếtcao kiệm khơ SửIDW, dụng thiết cứu bị đo giá rẻ AirBeam2 cải xâyđược dựng đồbố phân bố có PM2.5 độgiải, phân xác IDW, nghiênnghiên cứu xây dựng đồbản phân PM2.5 với độ với phân độ tin độ cậytin cậy xác cao sử dụng tốt điểm đo có phân bố thiện độ phân giải không gian cho thấy nghiên tơi chúng ghinồng nhậnđộ nồng độ PM2.5 trungước bình ước tínhgiao đượcđộng giaogiữa động Kết quảKết nghiên cứu củacứu chúng ghi tơi nhận PM2.5 trung bình tính 3 quận khoảng 10 µg/m đến 33,93 µg/m ; mùa mưa có xu hướng cao mùa khô Phần lớn các quận khoảng 10 µg/m đến 33,9 µg/m ; mùa mưa có xu hướng cao mùa khô Phần lớn 14 TCNCYH 142 (6) - đảm 2021 113 quận/huyện bảo 05: QCVN 05: 2013/BTNMT < 253)µg/m , nhiên sốvượt quậnngưỡng vượt ngưỡng 14 quận/huyện đảm bảo QCVN 2013/BTNMT (PM2.5(PM2.5 < 25 µg/m , tuy) nhiên sốmột quận Hóc mơn, Củ Chi Đáng ý, tất quận cho thấy nồng độ PM2.5 trung bình năm cao Hóc mơn, Củ Chi Đáng ý, tất quận cho thấy nồng độ PM2.5 trung bình năm cao tiêu an chuẩn tồn sứccủa khỏe của(PM2.5 WHO (PM2.5 < 103).µg/m Cácchứng chứng trước đâynêu lên tình tiêu chuẩn tồnansức khỏe WHO < 10 µg/m Các ).bằng trước lênnêu tình trạng nhiễm PM2.5 tương tự Nồng độ PM2.5 ghi nhận Sở Tài nguyên Mơi trường TPHCM TẠP CHÍ NGHIÊN CỨU Y HỌC rộng khắp diện tích tính tốn Độ xác IDW đánh giá nghiên cứu năm 2015 Trung Quốc cho thấy IDW có độ xác cao với hệ số tương quan 0,99.16 Như vậy, việc kết hợp quan trắc phạm vi rộng với IDW, nghiên cứu xây dựng đồ phân bố PM2.5 với độ phân giải, độ tin cậy xác cao Kết nghiên cứu chúng tơi ghi nhận nồng độ PM2.5 trung bình ước tính giao động quận khoảng 10 µg/m3 đến 33,9 µg/m3; mùa mưa có xu hướng cao mùa khô Phần lớn quận/huyện đảm bảo QCVN 05: 2013/BTNMT (PM2.5 < 25 µg/ m3),14 nhiên số quận vượt ngưỡng Hóc Mơn, Củ Chi Đáng ý, tất quận cho thấy nồng độ PM2.5 trung bình năm cao tiêu chuẩn an tồn sức khỏe WHO (PM2.5 < 10 µg/m3).8 Các chứng trước nêu lên tình trạng nhiễm PM2.5 tương tự Nồng độ PM2.5 ghi nhận Sở Tài nguyên Môi trường TPHCM giai đoạn 2013 đến 2017 thường xuyên vượt tiêu chuẩn WHO nồng độ ô nhiễm 24 (< 25 µg/m3) Kết đo lường năm 2016 thành phố Hồ Chí Minh ghi nhận nồng độ PM2.5 trung bình năm 42 µg/m.37 Theo báo cáo Trung tâm Phát triển Sáng tạo Xanh (GreenID), năm 2017 nồng độ PM2.5 trung bình TPHCM 29,6 μg/m3; số ngày thành phố có nồng độ PM2.5 vượt tiêu chuẩn trung bình 24 WHO (25 μg/m3) 222 ngày vượt tiêu chuẩn trung bình 24 QCVN 14 ngày (50 μg/m3).17 Các chứng trước kết hợp với kết nghiên cứu cho thấy tình trạng nhiễm PM2.5 TPHCM gây ảnh hưởng đến sức khỏe người dân Xây dựng đồ phân bố PM2.5 bước tảng để đánh giá tác động sức khỏe PM2.5; nhờ nhà hoạch định có cảnh báo biện pháp phòng ngừa giảm thiểu tác động sức khỏe 114 Nghiên cứu kết hợp đồ phân bố PM2.5 với liệu dân số tử vong sẵn có TPHCM để ước tính tác động sức khỏe Mơ hình BenMAP sử dụng để đánh giá tác động lên tử vong PM2.5.13 Mơ hình sử dụng rộng rãi Hoa Kỳ, Nigeria, Chi Lê, Trung Quốc, Thái Lan với ưu điểm ước tính số người tử vong ô nhiễm PM2.5 khu vực trường hợp có thơng tin dân số phơi nhiễm tỷ suất tử vong năm khu vực nghiên cứu Tính xác ước lượng tử vong PM2.5 phụ thuộc vào tham số đầu vào Nghiên cứu chọn tham số đầu vào dựa chứng nhằm đảm bảo tính xác kết Nồng độ PM2.5 (tham số A, cơng thức 2) trích xuất từ đồ phân bố toàn TPHCM giá trị tham chiếu chọn ngưỡng 10mg/m3 ngưỡng an toàn sức khỏe WHO đề xuất.18 Ước tính tác động tăng 6,2% tử vong PM2.5 tăng 10 μg/m3 tham khảo từ nghiên cứu Tallin nguồn liệu xác hợp (tham số B, công thức 2).19 Tỉ suất tử vong ung thư phổi, hệ hô hấp thu thập từ sổ A6/YTCS Bộ Y tế (tham số C, công thức 2) Các liệu dân số (tham số D, công thức 2) thu thập dựa vào thống kê năm 2018 Cục thống kê TPHCM, nguồn liệu phù hợp phản ánh phân bố dân cư toàn thành phố.13, 20 Dữ liệu địa lý trích xuất lớp địa giới hành từ Trung tâm Ứng dụng Hệ thống Thơng tin Địa lý TPHCM Như vậy, liệu sử dụng ước tính tác động sức khỏe có tính xác hợp thống cao Tuy nhiên, hạn chế đặc thù liệu thứ cấp tồn nghiên cứu đánh giá đầy đủ chất lượng liệu Nghiên cứu ghi nhận tỉ suất tử vong ung thư phổi ung thư hệ hô hấp 0,631/10.000 dân 0,712/10.000 dân Theo thống kê từ liệu năm 2018, tỉ TCNCYH 142 (6) - 2021 TẠP CHÍ NGHIÊN CỨU Y HỌC suất tử vong chung toàn dân số năm 2018 33,031/10000 dân Như vậy, tử vong ung thư phổi, ung thư hệ hô hấp chiếm khoảng 2% số ca tử vong nguyên nhân TPHCM năm 2018 Kết nghiên cứu chúng tơi thể rõ đóng góp PM2.5 vào 6,3% (35/557) số ca tử vong ung thư phổi 6,5% (41/629) số ca tử vong ung thư hệ hô hấp năm 2018 Tỉ suất tử vong ung thư phổi ung thư hệ hô hấp cho thấy mối liên quan thuận mạnh với nồng độ PM2.5 với hệ số tương quan 0,733 0,752 với p < 0,001 Với dân số 8.831.865 TPHCM năm 2018, nồng độ PM2.5 tăng thêm 10 μg/m3 số ca tử vong ung thư phổi gia tăng thêm 56 ca ung thư hệ hô hấp 64 ca Các kết nghiên cứu chúng tơi đóng góp thêm chứng tác động PM2.5 đến sức khỏe Nhiều nghiên cứu trước mối liên quan PM2.5 với giảm tuổi thọ,19, 21 tử vong chung,22-24 nhập viện,23, 24 thiệt hại kinh tế.20, 23 Các nghiên cứu tử vong không xác định đóng góp PM2.5 vào nguyên nhân cụ thể, đó, kết bị ảnh hưởng trường hợp tử vong nguyên nhân độc lập với PM2.5 Mối liên quan PM2.5 ung thư phổi, hệ hô hấp củng cố dựa nhiều chứng.25-28 Do đó, việc tập trung phân tích tử vong ung thư phổi, hệ hô hấp loại trừ nguyên nhân tử vong độc lập với PM2.5 từ cung cấp ước tính tin cậy Bên cạnh đó, cảnh báo dựa bệnh lý cụ thể giúp tăng mức độ quan tâm đem lại hiệu cao so với cảnh báo tử vong chung, đặc biệt bệnh lý ung thư có xu hướng gia tăng người dân quan tâm Nghiên cứu ghi nhận số hạn chế ảnh hưởng đến phát Thiết bị quan trắc giá rẻ có độ xác thấp so với thiết bị quan trắc cố định tiêu chuẩn Dù vậy, nghiên cứu cho thấy thiết bị TCNCYH 142 (6) - 2021 giá rẻ AirBeam có độ tin cậy cao so sánh với thiết bị quan trắc tự động cố định tiêu chuẩn (được xem tiêu chuẩn vàng); hệ số tương quan ghi nhận đạt 0,84 Thêm vào đó, thiết bị giá rẻ AirBeam cho phép thu thập liệu có độ phân giải lớn với mức chi phí thấp Tiêu chuẩn năm sinh sống địa phương chưa đủ để làm tăng tỉ lệ mắc bệnh Tuy nhiên, việc chọn giới hạn năm phụ thuộc vào hệ thống thu thập liệu tử vong A6/YTCS Bộ Y tế Nghiên cứu khơng tìm liệu có giới hạn dài hơn, hạn chế mà nhóm nghiên cứu khơng khắc phục Sử dụng thiết kế sinh thái, số liệu phơi nhiễm với yếu tố gây nhiễu không tiếp cận cấp độ quận/huyện dẫn đến khơng thể kiểm sốt yếu tố gây nhiễu tiềm ẩn Đây điểm hạn chế điển hình nghiên cứu sinh thái Trong tương lai, nghiên cứu ước lượng gánh nặng tử vong PM2.5 cần tiếp tục mở rộng kỹ thuật đo lường dựa thiết bị giá rẻ để đạt độ phân giải không gian cao phục vụ phân tích cấp độ cá nhân cấp hành thấp xã/phường Cùng với liệu yếu tố gây nhiễu cần tiếp cận cấp độ phù hợp để loại trừ tác động gây nhiễu cung cấp ước lượng xác V KẾT LUẬN Ơ nhiễm PM2.5 TPHCM thường xuyên ghi nhận vượt giới hạn an tồn sức khỏe WHO đóng góp vào tử vong ung thư phổi ung thư hệ hơ hấp năm 2018 Cần có biện pháp ngắn hạn phòng tránh tác hại PM2.5 biện pháp dài hạn giảm thiểu ô nhiễm PM2.5 Những kết nguồn thông tin quan trọng cho chương trình truyền thơng sức khỏe, bảo vệ mơi trường, nguồn thông tin tham khảo cho nhà hoạch định sách đề mục tiêu giảm thiểu ô nhiễm môi trường để bảo vệ sức khỏe 115 TẠP CHÍ NGHIÊN CỨU Y HỌC LỜI CẢM ƠN Nghiên cứu nhận tài trợ từ Đề tài tiềm Trường Đại học Y Dược Tp.HCM hợp đồng số 66/2019/HĐ- ĐHYD TS Trần Ngọc Đăng làm chủ nhiệm Ngồi ra, nhóm nghiên cứu xin bày tỏ lòng biết ơn chân thành đến: Sở Tài nguyên Môi trường Tp.HCM, Cục thống kê TP.HCM, Cục quản lý môi trường y tế - Bộ Y tế, Trung tâm Ứng dụng Hệ thống Thông tin Địa lý TPHCM tạo điều kiện giúp đỡ chúng tơi q trình thu thập số liệu Global Burden of Disease Project, The Health Effects Institute The State of Global Air https:// www.stateofglobalair.org/sites/default/files/ soga-2018-report.pdf, Accessed Feb 12, 2019 Chúng chân thành cảm ơn PGS.TS Tô Thị Hiền, trường Đại học Khoa Học Tự Nhiên TpHCM PGS.TS Hồ Quốc Bằng, Viện Môi Trường Tài Nguyên, Đại học Quốc Gia TpHCM cố vấn chuyên môn cho nghiên cứu World Health Organization WHO outdoor air quality guidelines http://www.euro.who int/en/health-topics/environment-and-health/ airquality/policy/who-outdoor-air-qualityguidelines, Accessed 01/05/2020 TÀI LIỆU THAM KHẢO Environmental Protection Agency Guidance on Choosing a Sampling Design for Environmental Data Collection: for Use in Developing a Quality Assurance Project Plan 2002 Environmental Pollution Centers Air Pollution Causes https://www environmentalpollutioncenters.org/air/ causes/, Accessed Dec19, 2018 Xing Y-F, Xu Y-H, Shi M-H, Lian Y-X The impact of PM2.5 on the human respiratory system Journal of Thoracic Disease 2016;8(1):69-74 Brook RD, Franklin B, Cascio W, et al Air Pollution and Cardiovascular Disease: A Statement for Healthcare Professionals From the Expert Panel on Population and Prevention Science of the American Heart Association Circulation 2004;109:2655-2671 Environmental Pollution Centers What Is Air Pollution? https://www environmentalpollutioncenters.org/air/, Accessed Feb 2, 2020 Song C, He J, Wu L, et al Health burden attributable to ambient PM2.5 in China Environmental Pollution 2017;223:575-586 Institute for Health Metrics, Evaluation’s 116 World Health Organization Western Pacific Region Frequently Asked Questions about Ambient and Household Air Pollution and Health http://www.wpro.who.int/vietnam/ mediacentre/features/air_pollution_QandA/en/, Accessed Jun 30, 2019 10 Feinberg S, Williams R, Hagler GSW, et al Long-term evaluation of air sensor technology under ambient conditions in Denver, Colorado Atmospheric Measurement Techniques 2018;11(8):4605-4615 doi:10.5194/amt-114605-2018 11 Setianto A, Triandini T Comparison of Kriging and Inverse Distance Weighted (IDW) interpolation methods in lineament extraction and analysis Journal of Applied Geology 2013;5(1):21-29 12 Li L, Losser T, Yorke C, Piltner R Fast Inverse Distance Weighting-Based Spatiotemporal Interpolation: A Web-Based Application of Interpolating Daily Fine Particulate Matter PM2.5 in the Contiguous U.S Using Parallel Programming and k-d Tree Environmental Research and Public Health 2014;11:9101-9141 13 United States Environmental TCNCYH 142 (6) - 2021 TẠP CHÍ NGHIÊN CỨU Y HỌC Protection Agency BENMAP Environmental Benefits Mapping and Analysis Program: User’s Manual Appendices 2008:1-14 term exposure to PM (2.5) in 23 European cities Eur J Epidemiol 2006;21(6):449-58 doi:10.1007/s10654-006-9014-0 14 Bộ Tài Nguyên Môi Trường QCVN 05:2013/BTNMT Quy chuẩn kỹ thuật quốc gia chất lượng khơng khí xung quanh 2013 22 Chen T, Deng S, Li M Spatial Patterns of Satellite-Retrieved PM2.5 and Long-Term Exposure Assessment of China from 1998 to 2016 Int J Environ Res Public Health Dec 2018;15(12) doi:10.3390/ijerph15122785 15 Michael H, Chris Chaeha Lim AirBeam2 Technical Specifications, Operation & Performance http://www.takingspace.org/ airbeam2-technical-specifications-operationperformance/, Accessed access on June 1st 2019 16 Zhang P, Shen T Comparison of different spatial interpolation methods for atmospheric pollutant PM2.5 by using GIS and Spearman correlation Journal of Chemical and Pharmaceutical Research 2015;7(12):452-469 17 Nguyen Thi Anh Thu, Blume L Báo cáo chất lượng khơng khí năm 2017 2017:1415 18 World Health Organization WHO Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide: global update 2005 : summary of risk assessment 2005 19 Orru H, Erik Teinemaa TL, Tamm T, Kaasik M, Kimmel V, al e Health impact assessment of particulate pollution in Tallinn using fine spatial resolution and modeling techniques Environmental Health Perspectives 2009;8:7 20 Ho Quoc Bang Modeling PM10 in Ho Chi Minh City, Vietnam and evaluation of its impact on human health Sustainable Environment Research 2017;27(2):95-102 21 Boldo E, Medina S, LeTertre A, et al Apheis: Health impact assessment of long- TCNCYH 142 (6) - 2021 23 Abe KC, Miraglia SG Health Impact Assessment of Air Pollution in Sao Paulo, Brazil Int J Environ Res Public Health Jul 11 2016;13(7):694 doi:10.3390/ijerph13070694 24 Le LTP, Leung A Associations between urban road-traffic emissions, health risks, and socioeconomic status in Ho Chi Minh City, Vietnam: a cross-sectional study The Lancet 2018;2:5 25 Huang F, Pan B, Wu J, Chen E, Chen L Relationship between exposure to PM2.5 and lung cancer incidence and mortality: A metaanalysis Oncotarget 2017;8(26):4332243331 doi:10.18632/oncotarget.17313 26 Yang S, Liucun Z, Fei Y, Xiangyin K, Tao H, Yu-Dong C Analysis of the relationship between PM2.5 and lung cancer based on protein-protein interactions Combinatorial Chemistry & High Throughput Screening 2016;19(2):100-108 doi:http://dx.doi.org/10 2174/1386207319666151110123345 27 Li R, Zhou R, Zhang J Function of PM2.5 in the pathogenesis of lung cancer and chronic airway inflammatory diseases Oncol Lett 2018;15(5):7506-7514 doi:10.3892/ ol.2018.8355 28 Xing Y-F, Xu Y-H, Shi M-H, Lian Y-X The impact of PM2.5 on the human respiratory system Journal of thoracic disease 2016;8(1): E69-E74 doi:10.3978/j issn.2072-1439.2016.01.19 117 TẠP CHÍ NGHIÊN CỨU Y HỌC Summary PM2.5 INCREASED RESPIRATORY MORTALITY IN HO CHI MINH CITY: A MULTI-SOURCE DATA STUDY The purpose of this study is to map PM2.5 spatial distribution and to assess the impact of PM2.5 on lung cancer mortality and respiratory cancer mortality in Ho Chi Minh City, Vietnam in 2018 The PM2.5 distribution was mapped based on 96 locations in the wet and dry seasons using the low-cost air pollution sensor (named AirBeam2) and the IDW algorithm Analysis of the PM2.5 impact on lung and respiratory cancer mortality was based on the A6/YTCS mortality data from the Ministry of Health Vietnam and the BenMAP model All districts showed that the annual average PM2.5 concentration was higher than the WHO health safety standard (PM2.5