Đánh giá chất lượng môi trường không khí thông qua số liệu quan trắc tự động trên địa bàn tỉnh bắc ninh năm 2019

HỌC VIỆN NƠNG NGHIỆP VIỆT NAM NGUYỄN BÍCH DIỆP ĐÁNH GIÁ CHẤT LƯỢNG MƠI TRƯỜNG KHƠNG KHÍ THƠNG QUA SỐ LIỆU TRẠM QUAN TRẮC TỰ ĐỘNG TRÊN ĐỊA BÀN TỈNH BẮC NINH NĂM 2019 Ngành: Khoa học môi trường Mã ngành: 8.44.03.01 Người hướng dẫn khoa học: PGS.TS Ngô Thế Ân NHÀ XUẤT BẢN HỌC VIỆN NÔNG NGHIỆP - 2020 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi, kết nghiên cứu trình bày luận văn trung thực, khách quan chưa dùng để bảo vệ lấy học vị Tôi xin cam đoan giúp đỡ cho việc thực luận văn cám ơn, thơng tin trích dẫn luận văn rõ nguồn gốc Hà Nội, ngày 22 tháng 05 năm 2020 Tác giả luận văn Nguyễn Bích Diệp i LỜI CẢM ƠN Trong suốt thời gian học tập, nghiên cứu hoàn thành luận văn, nhận hướng dẫn, bảo tận tình thầy giáo, giúp đỡ, động viên bạn bè, đồng nghiệp gia đình Nhân dịp hồn thành luận văn, cho phép tơi bày tỏ lịng kính trọng biết ơn sâu sắc PGS.TS Ngơ Thế Ân tận tình hướng dẫn, dành nhiều công sức, thời gian tạo điều kiện cho tơi suốt q trình học tập thực đề tài Tơi xin bày tỏ lịng biết ơn chân thành tới Ban Giám đốc, Ban Quản lý đào tạo, Bộ môn Sinh thái Nông nghiệp, Khoa Môi trường - Học viện Nơng nghiệp Việt Nam tận tình giúp đỡ tơi q trình học tập, thực đề tài hồn thành luận văn Tơi xin chân thành cảm ơn tập thể lãnh đạo, cán viên chức Trung tâm quan trắc Tài nguyên môi trường Bắc Ninh giúp đỡ tạo điều kiện cho suốt trình thực đề tài Xin chân thành cảm ơn gia đình, người thân, bạn bè, đồng nghiệp tạo điều kiện thuận lợi giúp đỡ tơi mặt, động viên khuyến khích tơi hoàn thành luận văn./ Hà Nội, ngày 22 tháng 05 năm 2020 Tác giả luận văn Nguyễn Bích Diệp ii MỤC LỤC Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục từ viết tắt v Danh mục bảng vi Danh mục hình vii Trích yếu luận văn viii Thesis abstract x Phần Mở đầu 1.1 Tính cấp thiết đề tài 1.2 Giả thuyết khoa học 1.3 Mục tiêu nghiên cứu 1.3.1 Mục tiêu chung 1.3.2 Mục tiêu cụ thể 1.4 Phạm vi nghiên cứu 1.5 Những đóng góp, ý nghĩa khoa học thực tiễn Phần Tổng quan tài liệu 2.1 Các số đánh giá mơi trường khơng khí 2.1.1 Chỉ số chất lượng khơng khí (AQI) 2.1.2 Chỉ số nhiễm khơng khí tương đối (RAPI) 10 2.2 Các phương pháp đánh giá chất lượng khơng khí 13 2.2.1 Các phương pháp đánh giá chất lượng khơng khí giới 13 2.2.2 Các phương pháp đánh giá chất lượng khơng khí Việt Nam 14 2.3 Hiện trạng mơi trường khơng khí 2.3.1 Hiện trạng môi trường khơng khí giới 2.3.2 Hiện trạng môi trường không khí Việt Nam 2.4 Hệ thống quan trắc mơi trường khơng khí tỉnh Bắc Ninh 15 2.5 Công tác quản lý mơi trường khơng khí 19 2.5.1 Công tác quản lý mơi trường khơng khí giới 19 2.5.2 Công tác quản lý môi trường khơng khí Việt Nam 23 iii Phần Vật liệu phương pháp nghiên cứu 27 3.1 Địa điểm nghiên cứu 27 3.2 Thời gian nghiên cứu 27 3.3 Đối tượng/vật liệu nghiên cứu 27 3.4 Nội dung nghiên cứu 27 3.5 Phương pháp nghiên cứu 28 3.5.1 Nguồn số liệu sử dụng đề tài 28 3.5.2 Phương pháp đánh giá số chất lượng mơi trường khơng khí (AQI) 29 3.5.3 Phương pháp đồ 32 3.5.4 Phương pháp xử lý số liệu 33 Phần Kết thảo luận 34 4.1 Điều kiện tự nhiên, kinh tế, xã hội tỉnh Bắc Ninh 34 4.1.1 Vị trí địa lý 34 4.1.2 Tình hình phát triển kinh tế - xã hội 35 4.1.3 Điều kiện tự nhiên 43 4.2 Hệ thống trạm quan trắc môi trường khơng khí tự động tỉnh Bắc Ninh 51 4.3 Chất lượng mơi trường khơng khí bắc ninh năm 2019 58 4.4 Sự biến động chất lượng khơng khí tỉnh Bắc Ninh 61 4.4.1 Biến động chất lượng khơng khí theo khơng gian 61 4.4.2 Biến động chất lượng khơng khí theo tháng năm 65 4.5 Giải pháp giám sát chất lượng mơi trường khơng khí địa bàn tỉnh Bắc Ninh 67 Phần Kết luận kiến nghị 70 5.1 Kết luận 70 5.2 Kiến nghị 71 Tài liệu tham khảo 72 Phụ lục 75 iv DANH MỤC CÁC TỪ VIẾT TẮT Chữ viết tắt Nghĩa Tiếng Việt AQI Chỉ số chất lượng khơng khí BTNMT Bộ Tài nguyên Môi trường CLKK Chất lượng không khí CNN Cụm cơng nghiệp HEI Viện tác động đến sức khỏe KCN Khu công nghiệp MTKK Môi trường khơng khí ONKK Ơ nhiễm khơng khí QCCP Quy chuẩn cho phép QCVN Quy chuẩn Việt Nam QĐ-TCMT Quyết định-tiêu chuẩn môi trường QĐ-TTg Quyết định- thủ tướng TB Thông báo THCS Trung học sở TN&MT Tài nguyên Môi trường UBND Ủy ban nhân dân WHO Tổ chức Y tế giới v DANH MỤC BẢNG Bảng 2.1 Các phương pháp tính AQI giới (1) Bảng 2.2 Thang phân cấp đánh giá mức độ nhiễm khơng khí RAPId/RAPI*d có trọng số ứng với n thơng số khảo sát (2) 12 Bảng 3.1 Danh sách trạm quan trắc trực tuyến tỉnh Bắc Ninh (3) 28 Bảng 3.2 Khoảng giá trị AQI đánh giá chất lượng khơng khí (4) 29 Bảng 3.3 Các mức VN_AQI tương ứng ảnh hưởng tới sức khỏe (5) 31 Bảng 3.4 Các mức VN_AQI số hoạt động khuyến nghị (6) 32 Bảng 4.1 Tỷ trọng ngành địa bàn tỉnh Bắc Ninh năm 2015 – 2019 (7) 36 Bảng 4.2 Chỉ số sản xuất công nghiệp tỉnh Bắc Ninh năm 2015 – 2019 (8) 38 Bảng 4.3 Tỷ trọng đóng góp vào giá trị tổng sản phẩm địa bàn tỉnh (9) 39 Bảng 4.4 Dân số trung bình tỉnh Bắc Ninh năm 2015 – 2018 (10) 40 Bảng 4.5 Nhiệt độ khơng khí trung bình tháng, trạm Bắc Ninh (oC) (11) 43 Bảng 4.6 Số nắng tháng, trạm Bắc Ninh (h) (12) 44 Bảng 4.7 Lượng mưa trung bình tháng, trạm Bắc Ninh (mm) (13) 44 Bảng 4.8 Độ ẩm tương đối trung bình tháng, trạm Bắc Ninh (%) (14) 45 Bảng 4.9 Lượng bốc trung bình tháng, trạm Bắc Ninh (mm) (15) 45 Bảng 4.10 Hiện trạng sử dụng đất tỉnh Bắc Ninh năm 2019 (16) 50 Bảng 4.11 Mạng lưới quan trắc tự động tỉnh Bắc Ninh (17) 54 Bảng 4.12 Chỉ số chất lượng không khí AQI trạm quan trắc năm 2019 (18) 59 Bảng 4.13 Kết phân nhóm trạm nghiên cứu theo số AQI (19) 63 Bảng 4.14 Kết kiểm định khác biệt thống kê thông số quan trắc trạm tự động (20) 64 Bảng 4.15 Kết Phân tích tương quan AQI thơng số khí tượng (21) 66 vi DANH MỤC HÌNH Hình 3.1 Sơ đồ khu vực nghiên cứu vị trí trạm quan trắc khơng khí tự động liên tục Bắc Ninh 27 Hình 4.1 Bản đồ tỉnh Bắc Ninh 35 Hình 4.2 Bản đồ trạng hệ thống đô thị tỉnh Bắc Ninh 42 Hình 4.3 Biểu đồ biểu diễn thay đổi yếu tố khí tượng tỉnh Bắc Ninh 45 Hình 4.4 Các yếu tố khí tượng từ trạm quan trắc tự động tỉnh Bắc Ninh năm 2019 46 Hình 4.5 Mơ hình trạm quan trắc chất lượng mơi trường khơng khí xung quanh tự động 52 Hình 4.6 Trạm quan trắc chất lượng mơi trường khơng khí xung quanh tự động cố định 53 Hình 4.7 Bản đồ số khơng khí AQI tỉnh Bắc Ninh năm 2019 62 Hình 4.8 Đồ thị AQI tỉnh Bắc Ninh năm qua tháng năm 2019 65 vii TRÍCH YẾU LUẬN VĂN Tên tác giả: Nguyễn Bích Diệp Tên luận văn: Đánh giá chất lượng mơi trường khơng khí thơng qua số liệu trạm quan trắc tự động địa bàn tỉnh Bắc Ninh năm 2019 Ngành: Khoa học môi trường Mã số: 8.44.03.01 Tên sở đào tạo: Học viện Nông nghiệp Việt Nam Mục đích nghiên cứu - Đánh giá chất lượng mơi trường khơng khí tỉnh Bắc Ninh thơng qua số liệu trạm quan trắc khơng khí tự động liên tục - Đề xuất giải pháp khai thác, sử dụng nguồn số liệu từ trạm quan trắc tự động để cảnh báo giám sát chất lượng môi trường khơng khí địa bàn tỉnh Bắc Ninh Phương pháp nghiên cứu Đề tài sử dụng phương pháp thu thập tài liệu thứ cấp từ Trung tâm quan trắc Tài nguyên môi trường Bắc Ninh nhằm xác định thơng tin địa điểm nghiên cứu Tính AQI theo Quyết định số 1459/QĐ-TCMT phần mềm Excel; Lập đồ phương pháp nội suy (Kriging) phần mềm ArcGIS 10.3 Sau có kết tiến hành phương pháp kiểm chứng thống kê (determinant analysis) đề xuất định hướng khai thác số liệu quản lý chất lượng khơng khí Kết kết luận Bắc Ninh tỉnh có diện tích nhỏ Việt Nam, nằm vùng kinh tế trọng điểm Bắc Hiện nay, Bắc Ninh có 62 làng nghề 11.000 doanh nghiệp hoạt động Bắc Ninh nằm vùng khí hậu cận nhiệt đới ẩm Có chênh lệch rõ ràng nhiệt độ mùa hè nóng ẩm mùa đơng khơ lạnh (15-16 °C) Hàng năm có hai mùa gió chính: gió mùa Đơng Bắc gió mùa Đơng Nam Gió mùa Đơng Bắc thịnh hành từ tháng 10 năm trước đến tháng năm sau, gió mùa Đơng Nam thịnh hành từ tháng đến tháng mang theo ẩm gây mưa rào Hầu hết số chất lượng khơng khí AQI Bắc Ninh cao mức trung bình cho phép Đặc biệt trạm quan trắc online tự động Phong Khê tháng 11 năm 2019 tính AQI 1081 gấp 10 lần khoảng giá trị số AQI trung bình Chỉ số AQI có khác biệt trạm quan trắc tháng năm Xu thấy rõ ràng AQI tăng cao tháng thu - đông giảm thấp tháng mùa hè Giải pháp giám sát chất lượng khơng khí dựa vào số liệu quan trắc tự động gồm có: tăng cường lắp đặt thiết bị quan trắc khí tự động doanh nghiệp có nguồn khí thải lớn; Huy động viii nguồn lực cho quản lý chất lượng khơng khí; Đẩy mạnh công tác tra, kiểm tra việc thực quy định pháp luật bảo vệ môi trường khơng khí doanh nghiệp địa bàn tỉnh Đồng thời, huy động tham gia cộng đồng việc giám sát khí thải; Tăng cường cơng tác tuyên truyền, giáo dục, nâng cao nhận thức cộng đồng tác hại nhiễm khơng khí ix chất ô nhiễm thứ cấp sinh tương tác chất ô nhiễm NOx, HC, VOC với xạ tử ngoại mặt trời Tại trạm quan trắc tự động gần đường giao thông, nồng độ O3 vượt giới hạn QCVN nhiều ngày năm Một số khu vực gần nhà máy sản xuất công nghiệp nằm đô thị, nồng độ chất CO, NO2, SO2 số thời điểm vượt giới hạn QCVN 05:2013 4.4.2 Biến động chất lượng khơng khí theo tháng năm Sự khác biệt chất lượng khơng khí theo tháng năm tương đối rõ ràng số AQI (log10AQI) thể đồ thị (Hình 4.8) Theo đồ thị này, xu biến động AQI có giá trị cao vào tháng đầu cuối năm (mùa thu - đông) hạ thấp vào tháng mùa hè Hình 4.8 Đồ thị AQI tỉnh Bắc Ninh năm qua tháng năm 2019 (9) Sự khác biệt số AQI qua tháng kiểm chứng thơng qua phân tích thống kê với hệ số Wilks' Lambda đạt 0,659 (p < 0,000) Như vậy, quy luật thay đổi chất lượng khơng khí theo tháng năm 2019 có độ tin cậy cao kết tương đồng với kết nghiên cứu gần 65 nhóm tác giả khác Hà Nội (Tran et al, 2020) Sự giao động lớn AQI tháng 5, thay đổi nhu cầu sản xuất khu vực làng nghề nên lượng khí thải khu vực khác Nhu cầu tiêu dùng thấp hoạt động sản xuất tạm trì hỗn để đảm bảo trì sản xuất Nguyên nhân tượng ô nhiễm không khí tăng cao tháng đầu cuối năm hộ, sở sản xuất hoạt động với công suất lớn, thời gian nhiều vào khoảng thời gian Ngoài ra, nguyên nhân biến động khí tượng đóng vai trị quan trọng dẫn tới thay đổi chất lượng khơng khí Kết nghiên cứu từ số liệu trạm quan trắc tự động Hà Nội (Tran et al 2020) chứng minh yếu tố khí tượng có tương quan chặt chẽ với hàm lượng bụi nhẹ PM0.5 Ở Bắc Ninh nói riêng khu vực phía bắc nói chung mùa hè có nhiệt cao, mưa nhiều, ánh sáng nhiều, gió mạnh, chất nhiễm nhanh chóng phát tán giúp cho nồng độ bụi khơng khí thấp Ngược lại, vào mùa đơng với nhiệt thấp, ánh sáng, mưa, gió mùa đơng bắc tràn về, đem theo chất nhiễm từ phía Bắc khiến cho khơng khí nhiễm Đặc biệt, ngày xảy tượng nghịch nhiệt, chất lượng khơng khí thường nhiễm nặng chất ô nhiễm không phát tán Hiện tượng xảy lên cao, nhiệt độ khơng khí cao (trái với quy luật thơng thường lên cao, nhiệt độ thấp) Lớp nghịch nhiệt giống mũ, ngăn chất ô nhiễm phát tán lên cao Vì vậy, chất ô nhiễm bị giữ lại môi trường không khí Q trình làm nồng độ chất nhiễm tăng cao khiến mơi trường khơng khí bị nhiễm nặng nề Ảnh hưởng yếu tố khí tượng, bao gồm nhiệt độ, áp suất khí quyển, độ ẩm khơng khí, tốc độ gió, lượng mưa xạ mặt trời đến thay đổi số AQI phân tích thơng qua hệ số tương quan Spearman (Spearman’s correlation) bảng sau Bảng 4.15 Kết Phân tích tương quan AQI thơng số khí tượng (21) Các thơng số Temp Press Humid Windspeed Precip Solar Hệ số tương quan Spearman -0,187** 0,331** -0,317** -0,083 -0,205** -0,002 Mức ý nghĩa (p) N 0,008 200 0,000 201 0,000 199 0,240 201 0,004 191 0,977 207 Ghi chú: ** Hệ số tương quan đạt mức ý nghĩa 0.01 (2-tailed) 66 Theo kết phân tích có tương quan yếu tố nhiệt độ, áp suất khí quyển, độ ẩm khơng khí lượng mưa có ý nghĩa thống kê (p < 0,01) (xem thêm phần phụ lục 3) Trong đó, nhiệt độ, độ ẩm khơng khí, tốc độ gió lượng mưa có mối tương quan âm – nghĩa giá trị yếu tố tăng số AQI giảm xuống, mức độ nhiễm Kết chứng minh số liệu quan trắc có quy luật logic với nguyên lý trình bày 4.5 GIẢI PHÁP GIÁM SÁT CHẤT LƯỢNG MÔI TRƯỜNG KHƠNG KHÍ TRÊN ĐỊA BÀN TỈNH BẮC NINH Theo Sở tài nguyên môi trường tỉnh Bắc Ninh, mục tiêu quản lý môi trường đến năm 2020 phấn đấu 80% sở sản xuất thép, hóa chất, phân bón hóa học xử lý bụi khí thải SO2, NOx, CO đạt Quy chuẩn kỹ thuật môi trường; 90% sở sản xuất nhiệt điện, 80% sở sản xuất xi măng, 70% sở sản xuất thép hóa chất đầu tư lắp đặt thiết bị quan trắc khí thải tự động liên tục với thông số theo Quy chuẩn kỹ thuật môi trường Đối với sở thuộc đối tượng phải lắp đặt hệ thống quan trắc khí thải tự động, tồn tỉnh có sở Cơng ty TNHH cơng nghệ hóa chất mơi trường Vũ Hồng, lắp đặt hệ thống quan trắc khí thải tự động, cố định Tồn tỉnh triển khai nhiều chương trình nâng cao chất lượng khơng khí, giảm thiểu nhiễm bụi như: trồng xanh, giảm việc đốt rơm rạ, than tổ ong, hạn chế đốt vàng mã; khuyến khích người dân sử dụng phương tiện công cộng, hạn chế sử dụng phương tiện cá nhân tham gia giao thông, tăng cường sử dụng nhiên liệu Tăng cường cơng tác giám sát chất lượng khơng khí xung quanh thông qua việc tăng số lượng Trạm quan trắc khơng khí xung quanh tự động liên tục đô thị, giám sát thường xuyên thông số theo Quy chuẩn kỹ thuật môi trường thông số VOCs, HC Để giảm thiểu nguồn nhiễm khơng khí lâu dài, tỉnh lựa chọn dự án đầu tư bảo đảm công nghệ thân thiện môi trường, công nghệ cao Nghiên cứu ứng dụng công nghệ xử lý chất thải phù hợp cho quy mô, loại hình sản xuất, đặc biệt bụi khí thải sở luyện kim, xây dựng, hóa chất, xi măng, nhiệt điện Thực chương trình hỗ trợ sản xuất hơn, đặc biệt sở công nghiệp vừa nhỏ Kiên ngăn chặn việc chuyển giao công nghệ, phương tiện, thiết bị cũ, lạc hậu, phế liệu nhập khẩu, chất thải gây ô nhiễm môi trường Tích cực tuyên truyền nâng cao nhận thức 67 người dân doanh nghiệp bảo vệ môi trường sử dụng hợp lý, tiết kiệm, hiệu nguồn tài nguyên, giảm phát thải khí, sức khỏe cộng đồng Bắc Ninh tiến hành lập quy hoạch, xây dựng lộ trình bước di dời sở sản xuất làng nghề vào khu, cụm công nghiệp tập trung để đầu tư xây dựng hệ thống xử lý khí thải; Yêu cầu hộ gia đình sản xuất làng nghề bắt buộc phải đầu tư xây dựng hệ thống xử lý khí thải phát sinh; Dừng sản xuất sở sản xuất gây ô nhiễm môi trường nghiêm trọng biện pháp cắt điện Cùng với đó, tăng cường công tác kiểm tra, xử lý nghiêm vi phạm pháp luật bảo vệ môi trường; Triển khai thực Quy chế bảo vệ môi trường làng nghề, khu dân cư, cụm công nghiệp tỉnh Bắc Ninh; Đẩy mạnh công tác điều tra bản, quan trắc, dự báo, cảnh báo kiểm sốt nhiễm, đặc biệt, làng nghề, khu vực có nguy gây nhiễm khơng khí; Xây dựng Đề án mở rộng hệ thống mạng lưới giám sát chất lượng môi trường làng nghề truyền thống, khu dân cư đông đúc quy định bắt buộc sở cơng nghiệp tự báo cáo tình hình bảo vệ môi trường… Mạng lưới quan trắc môi trường tỉnh Bắc Ninh giai đoạn 2016-2020, đó, mơi trường khơng khí theo dõi chất lượng 52 điểm (có điểm theo dõi chất lượng trạng môi trường huyện, thị xã thành phố, 21 điểm theo dõi chất lượng khu vực chịu tác động từ hoạt động sản xuất kinh doanh) Lắp đặt 18 Trạm quan trắc mơi trường khơng khí tự động truyền liệu sở Tài nguyên Môi trường để giám sát Vị trí số lượng trạm quan trắc địa bàn tỉnh hợp lý Do phản ánh thực trạng mơi trường khơng khí tỉnh Trong tương lai cần phát triển trạm quan trắc online, tự động điểm phát thải khơng khí lớn cụ thể hơn: nhà máy, sở sản xuất… để tăng cường công tác quản lý đảm bảo việc bảo vệ môi trường Công tác quản lý mơi trường nói chung, quản lý chất lượng khơng khí nói riêng thực đồng sở tăng cường lực quan trắc môi trường (trực tiếp định kỳ) Từ năm 2016, mạng lưới quan trắc môi trường điều chỉnh theo hướng bổ sung điểm quan trắc mơi trường tác động, có quan trắc chất lượng mơi trường khơng khí Hệ thống tiếp nhận liệu quan trắc tự động đưa vào vận hành, theo dõi, giám sát thường xuyên sở liệu truyền Đẩy mạnh công tác thanh, kiểm tra giám sát việc chấp hành quy định pháp luật doanh nghiệp, đặc biệt việc xây dựng, vận hành hệ thống xử lý chất thải kiểm soát ô nhiễm định kỳ 68 Để thực Kế hoạch hành động quốc gia quản lý chất lượng không khí đến năm 2020, tầm nhìn đến năm 2025 Thủ tướng Chính phủ, tỉnh có chế sách nhằm khuyến khích sở sản xuất sử dụng nhiên liệu sạch, đổi công nghệ để giảm phát thải khí nhà kính, sử dụng khí nén tự nhiên làm nguyên liệu cho phương tiện giao thông; Tăng cường lắp đặt thiết bị quan trắc khí tự động doanh nghiệp có nguồn khí thải lớn; Nghiên cứu, thực Đề án Kiểm soát khí thải xe mơ tơ, xe gắn máy tham gia giao thông Mặt khác, huy động nguồn lực cho quản lý chất lượng khơng khí; Đẩy mạnh cơng tác tra, kiểm tra việc thực quy định pháp luật bảo vệ mơi trường khơng khí doanh nghiệp địa bàn tỉnh Đồng thời, huy động tham gia cộng đồng việc giám sát khí thải; Tăng cường cơng tác tun truyền, giáo dục, nâng cao nhận thức cộng đồng tác hại nhiễm khơng khí 69 PHẦN KẾT LUẬN VÀ KIẾN NGHỊ 5.1 KẾT LUẬN Bắc Ninh tỉnh có diện tích nhỏ Việt Nam, nằm vùng kinh tế trọng điểm Bắc Hiện nay, Bắc Ninh có 62 làng nghề, 11.000 doanh nghiệp hoạt động Bắc Ninh nằm vùng khí hậu cận nhiệt đới ẩm Có chênh lệch rõ ràng nhiệt độ mùa hè nóng ẩm mùa đông khô lạnh (15-16 °C) Mùa mưa kéo dài từ tháng năm đến tháng mười hàng năm Hàng năm có hai mùa gió chính: gió mùa Đơng Bắc gió mùa Đơng Nam Gió mùa Đơng Bắc thịnh hành từ tháng 10 năm trước đến tháng năm sau, gió mùa Đông Nam thịnh hành từ tháng đến tháng mang theo ẩm gây mưa rào Chỉ số chất lượng khơng khí AQI Bắc Ninh có khác biệt rõ ràng trạm quan trắc tháng năm Chỉ số AQI đặc biệt cao trạm quan trắc online tự động Phong Khê tháng 11 năm 2019, đạt mực 1081, gấp 10 lần khoảng giá trị số AQI trung bình Nguyên nhân khác biệt liên quan tới phân bổ nguồn gây ô nhiễm tập trung cao khu vực thành phố Mức độ ảnh hưởng phân loại theo màu sắc số AQI 0-50: Tốt (màu Xanh) 51-100: Trung bình (màu Vàng) 101-150: Xấu (Ảnh hưởng xấu tới nhóm sức khỏe nhạy cảm) (màu )Da cam 151-200: Kém (Ảnh hướng xấu tới sức khỏe) (màu Đỏ) 201-300: Rất (màu Tím) 301-500: Nguy hại (Ảnh hướng xấu tới sức khỏe) (màu Nâu) Nhóm nhạy cảm bao gồm: trẻ em, người già người mắc bệnh hô hấp, tim mạch… Xu biến động số chất lượng khơng khí AQI theo tháng năm 2019 khẳng định kết phân tích thống kê Nhìn chung, số AQI cao mùa thu – đông (trong tháng đầu cuối năm) thấp mùa hè Nguyên nhân khác biệt thời điểm tập trung sản xuất hàng hóa tập trung cao điểm vào dịp đầu năm cuối năm Ngoài ra, yếu tố khí tượng, bao gồm nhiệt độ, áp suất khí quyển, độ ẩm khơng khí lượng mưa có ý nghĩa chi phối biến động số AQI Bắc Ninh phát triển chế sách nhằm khuyến khích sở sản xuất sử dụng nhiên liệu sạch, đổi cơng nghệ để giảm phát thải khí nhà kính, sử dụng khí nén tự nhiên làm nguyên liệu cho phương tiện giao thông; Tăng 70 cường lắp đặt thiết bị quan trắc khí tự động doanh nghiệp có nguồn khí thải lớn; Mặt khác, huy động nguồn lực cho quản lý chất lượng khơng khí; Đẩy mạnh cơng tác tra, kiểm tra việc thực quy định pháp luật bảo vệ mơi trường khơng khí doanh nghiệp địa bàn tỉnh Đồng thời, huy động tham gia cộng đồng việc giám sát khí thải; Tăng cường công tác tuyên truyền, giáo dục, nâng cao nhận thức cộng đồng tác hại ô nhiễm khơng khí 5.2 KIẾN NGHỊ Tiếp tục phát triển nghiên cứu đề tài theo hướng: sử dụng lượng mặt trời để vận hành trạm quan trắc tự động, hạn chế sai sót q trình vận hành máy… để dự báo thông số AQI tốt phục vụ đời sống người bảo vệ môi trường Tăng cường công tác tuyên truyền, giáo dục cho người dân bảo vệ môi trường, bảo vệ sức khỏe, tích cực thực hành động xanh: trồng xanh, không đốt rơm rạ… Sử dụng hiệu thông số từ trạm quan trắc online tự động để có biện pháp bảo vệ mơi trường phù hợp tiết kiệm địa phương 71 TÀI LIỆU THAM KHẢO Al-Odaini, A.N., P.Zakaria, M.Zali, A.M.Juahir,H.Yaziz, I.M Surif (2012) Application of chemometrics in understanding the spatial distribution of human pharmaceuticals in surface water, Environ.Monit.Assess 184:6735-6748 DOI 10.1007/s10661-011-2454-3 An Bình (2019) Tồn cảnh rừng tự nhiên Việt Nam kể từ năm 1945 Tạp chí mơi trường http://tapchimoitruong.vn/ pages/article.aspx?item= To%C3%A0nc%E1%BA%A3nh-v%E1%BB%81-r%E1%BB%ABng-t%E1%BB%B1nhi%C3%AAn-c%E1%BB%A7a-Vi%E1%BB%87t-Nam-k%E1%BB%83t%E1%BB%AB-n%C4%83m-1945-45490 Bộ Tài nguyên môi trường (2016) Chương trình tổng thể quan trắc quốc gia mơi trường khơng khí, tiếng ồn độ rung giai đoạn 2017 – 2020 Tr 25 – 45 Bộ Tài nguyên môi trường (2019) Báo cáo Hiện trạng Môi trường quốc gia 2019: Nhận diện nguồn gây ô nhiễm Cục thống kê tỉnh Bắc Ninh (2019) Niên giám thống kê Bắc Ninh Nhà xuất thống kê Chi cục bảo vệ môi trường tỉnh Bắc Ninh (2011) Đánh giá chất lượng khơng khí khu cơng nghiệp tỉnh Bắc Ninh Chinh C Tran, Thi D Ta, An T Duong, Oanh T K Phan, Dung A Nguyen (2020) Analysis on Temporal Pattern of Fine Particulate Matter (PM2.5) in Hanoi, Vietnam and the Impact of Meteorological Conditions Journal of Environmental Protection, 2020, 11, 241-256 ISSN Online: 2152-2219 Farmaki, G.E., S.N.Thomaidis, V.Simeonov, E.C.Efstathiou (2012) A comparative chemometric study for water quality expertise of the Athenian water reservoirs,Environ.Monit.Assess, 184:7635-7652.DOI 10.1007/s10661-012-2524 Giáng Hương (2019) Hiểu để hành động nhiễm khơng khí Tạp chí mơi trường http://tapchimoitruong.vn/pages/article.aspx?item=Hi%E1%BB%83u%C4%91%C3%BAng-%C4%91%E1%BB%83-h%C3%A0nh%C4%91%E1%BB%99ng-%C4%91%C3%BAng-v%E1%BB%81-%C3%B4nhi%E1%BB%85m-kh%C3%B4ng-kh%C3%AD-50854 72 Hà Minh Họa (2014) Hiện trạng giải pháp cải thiện chất lượng môi trường làng nghề tỉnh Bắc Ninh Tạp chí mơi trường http://tapchimoitruong.vn/pages/article.aspx?item=Hi%E1%BB%87ntr%E1%BA%A1ng-v%C3%A0-gi%E1%BA%A3i-ph%C3%A1pc%E1%BA%A3i-thi%E1%BB%87n-ch%E1%BA%A5tl%C6%B0%E1%BB%A3ng-m%C3%B4i-tr%C6%B0%E1%BB%9Dngl%C3%A0ng-ngh%E1%BB%81-t%E1%BB%89nh-B%E1%BA%AFc-Ninh38643 http://tapchimoitruong.vn/pages/article.aspx?item=Qu%E1%BA%A3n-l%C3%BDch%E1%BA%A5t-l%C6%B0%E1%BB%A3ng-kh%C3%B4ng-kh%C3%ADt%E1%BA%A1i-Vi%E1%BB%87t-Nam:-Th%E1%BB%B1ctr%E1%BA%A1ng-v%C3%A0-gi%E1%BA%A3i-ph%C3%A1p-38525 Huy Hoàng (2019) Hà Nội: Thực nhiều giải pháp nhằm giảm thiểu nhiễm khơng khí Tạp chí môi trường http://tapchimoitruong.vn/pages/article.aspx?item=H%C3%A0N%E1%BB%99i:-Th%E1%BB%B1c-hi%E1%BB%87n-nhi%E1%BB%81ugi%E1%BA%A3i-ph%C3%A1p-nh%E1%BA%B1m-gi%E1%BA%A3mthi%E1%BB%83u-%C3%B4-nhi%E1%BB%85m-kh%C3%B4ng-kh%C3%AD50546 International Network for Environmental Compliance and Enforcement, Principles of Environmental Compliance and Enforcement Handbook, (2009) Phạm Ngọc Đăng (2016) Ðánh giá mức độ nhiễm khơng khí Hà Nội đề xuất số giải pháp Tạp chí mơi trường http://tapchimoitruong.vn/pages/article.aspx ?item=%C3%90%C3%A1nh-gi%C3%A1-m%E1%BB%A9c%C4%91%E1%BB%99-%C3%B4-nhi%E1%BB%85m-kh%C3%B4ngkh%C3%AD-%E1%BB%9F-H%C3%A0-N%E1%BB%99i-v%C3%A0%C4%91%E1%BB%81-xu%E1%BA%A5t-m%E1%BB%99t-s%E1%BB%91gi%E1%BA%A3i-ph%C3%A1p-48633 Pham Ngoc Ho (2011) Weighted and Standardized Total Environmental Quality Index (TEQI) Approach in Assessing Environmental Components (Air, Soil and Water) Tạp chí ĐHQG (tiếng Anh), số 3, tập 27, chuyên san Khoa học Trái đất, 10/2011 Pham Ngoc Ho (2017) Relative Air Pollution Index (RAPI) - A New Method for Aggregate Air Pollution Assessment Research Center for Environmental Monitoring and Modeling (CEMM), VNU University of Science 73 Sở Tài nguyên môi trường tỉnh Bắc Ninh (2019) Báo cáo trạng môi trường tỉnh Bắc Ninh 05 năm (2015 - 2019) Sở Y tế Bắc Ninh (2018) Báo cáo hàng năm Sở Y tế Bắc Ninh Tổng cục môi trường (2017) Báo cáo môi trường quốc gia Tổng cục môi trường (2019) Quyết định 1459/QĐ-TCMT 2019 Hướng dẫn kỹ thuật tính tốn cơng bố số chất lượng khơng khí TS Đặng Văn Lợi - Phó Cục trưởng ThS Nguyễn Hồng Đức (2014) Cục Kiểm sốt nhiễm Quản lý chất lượng khơng khí Việt Nam: Thực trạng giải pháp Tạp chí mơi trường Trần Sơn Tùng (2012) Đánh giá chất lượng khơng khí thơng qua số liệu trạm quan trắc tự động Hà Nội, Đà Nẵng đề xuất giải pháp giảm thiểu Trường đại học Khoa học Tự nhiên Trung tâm quan trắc Tài nguyên môi trường Bắc Ninh (2018) Báo cáo trạng môi trường tỉnh Bắc Ninh Trung tâm quan trắc Tài nguyên môi trường Bắc Ninh (2018) Đề án Quản lý vận hành trạm, xe quan trắc môi trường tự động, bảng thông tin môi trường Trung tâm quan trắc Tài nguyên môi trường Bắc Ninh (2019) Báo cáo trạng môi trường tỉnh Bắc Ninh Trung tâm Quan trắc Tài nguyên Môi trường tỉnh Bắc Ninh (2019) Báo cáo tổng hợp kết thực Mạng lưới quan trắc tài nguyên môi trường tỉnh Bắc Ninh năm 2019 74 PHỤ LỤC Kết chuyển đổi số liệu 75 Kết phân tích thống kê Discriminant Analysis Tests of Equality of Group Means Wilks' Lambda PM25log PM10log O3log SO2log NO2log COlog F 901 635 966 962 900 773 df1 df2 11.041 57.508 3.517 3.944 11.062 29.301 2 2 2 Sig 200 200 200 200 200 200 000 000 032 021 000 000 Pooled Within-Groups Matrices PM25log PM10log Correlation O3log SO2log NO2log COlog PM25log 1.000 149 -.189 083 547 -.024 PM10log 149 1.000 -.119 001 219 264 O3log -.189 -.119 1.000 028 -.100 -.001 SO2log 083 001 028 1.000 048 149 NO2log 547 219 -.100 048 1.000 276 COlog -.024 264 -.001 149 276 1.000 Analysis Box's Test of Equality of Covariance Matrices Log Determinants Wards Rank Pooled within-groups Log Determinant 6 6 -21.333 -23.197 -12.841 -15.429 The ranks and natural logarithms of determinants printed are those of the group covariance matrices Test Results Box's M F 1.346E3 27.394 Approx df1 42 df2 2.411E3 Sig .000 Tests null hypothesis of equal population covariance matrices Summary of Canonical Discriminant Functions Eigenvalues Function Eigenvalue 894a 134a % of Variance Cumulative % 87.0 13.0 87.0 100.0 a First canonical discriminant functions were used in the analysis 76 Canonical Correlation 687 344 Wilks' Lambda Test of Function(s) Wilks' Lambda through 2 Chi-square 466 882 150.944 24.822 df Sig 12 000 000 Standardized Canonical Discriminant Function Coefficients Function PM25log -.175 282 PM10log 776 -.004 O3log 230 -.060 SO2log -.224 -.371 NO2log -.195 743 COlog 445 072 Pooled within-groups correlations between discriminating variables and standardized canonical discriminant functions Variables ordered by absolute size of correlation within function * Largest absolute correlation between each variable and any discriminant function Canonical Discriminant Function Coefficients Function PM25log PM10log O3log SO2log NO2log COlog (Constant) -.724 3.493 718 -.512 -.981 1.062 -7.270 1.171 -.020 -.187 -.850 3.732 172 -5.801 Unstandardized coefficients Functions at Group Centroids Function Wards 1 -.077 346 -3.841 -.581 220 304 Unstandardized canonical discriminant functions evaluated at group means Classification Statistics Classification Processing Summary Processed Excluded Missing or out-of-range group codes At least one missing discriminating variable Used in Output 77 216 13 203 Classification Resultsb,c Predicted Group Membership Wards Original Count % Cross-validateda Count % Total 36 21 57 44 91 135 2 11 63.2 36.8 100.0 32.6 67.4 100.0 18.2 18.2 63.6 100.0 35 21 57 48 87 135 11 61.4 36.8 1.8 100.0 35.6 64.4 100.0 36.4 18.2 45.5 100.0 a Cross validation is done only for those cases in the analysis In cross validation, each case is classified by the functions derived from all cases other than that case b 66.0% of original grouped cases correctly classified c 62.6% of cross-validated grouped cases correctly classified 78 Phân tích thống kê Spearman’s correlation Correlations AQIlog Spearman's rho AQIlog Correlation Coefficient Humid Solar Windspeed Precip Solar 331** -.317** -.083 -.205** -.002 008 000 000 240 004 977 216 200 201 199 201 191 207 -.187** 1.000 -.828** -.084 -.033 -.270** 632** Sig (2-tailed) 008 000 239 648 000 000 N 200 200 197 198 198 186 196 331** -.828** 1.000 -.069 -.015 010 -.427** Sig (2-tailed) 000 000 335 828 897 000 N 201 197 201 198 200 187 199 -.317** -.084 -.069 1.000 282** 640** -.202** Sig (2-tailed) 000 239 335 000 000 004 N 199 198 198 199 199 185 197 -.083 -.033 -.015 282** 1.000 313** 372** Sig (2-tailed) 240 648 828 000 000 000 N 201 198 200 199 201 187 199 -.205** -.270** 010 640** 313** 1.000 -.208** Sig (2-tailed) 004 000 897 000 000 004 N 191 186 187 185 187 191 186 -.002 632** -.427** -.202** 372** -.208** 1.000 Sig (2-tailed) 977 000 000 004 000 004 N 207 196 199 197 199 186 207 Correlation Coefficient Correlation Coefficient Correlation Coefficient Windspeed Correlation Coefficient Precip Humid -.187** N Press Press 1.000 Sig (2-tailed) Temp Temp Correlation Coefficient Correlation Coefficient ** Correlation is significant at the 0.01 level (2-tailed) 79