Các mô hình tính toán lan truyền ô nhiễm không khí

MƠ HÌNH TÍNH TỐN Ơ NHIỄM KHƠNG KHÍ I Phương trình để tính nồng độ chất nhiễm khí quyển: Khi mơ tả q trình khuyếch tán chất nhiễm khơng khí mơ hình tốn học mức độ nhiễm khơng khí thường đặc trưng trị số nồng độ chất ô nhiễm phân bố không gian biến đổi theo thời gian Trong trường hợp tổng quát, trị số trung bình nồng độ nhiễm khơng khí phân bố theo thời gian không gian mô tả từ phương trình chuyển tải vật chất (hay phương trình truyền nhiệt) biến đổi hoá học đầy đủ sau: C C C C   C    C    C  C u v w  k x   k y   k z  C  C  w c t x y z x  x  y  y  z  z  z (1) Trong đó: C : Nồng độ chất nhiễm khơng khí x,y,z: Các thành phần toạ độ theo trục Ox, Oy, Oz t : Thời gian Kx, Ky, Kz : Các thành phần hệ số khuyếch tán rối theo trục Ox, Oy Oz u,v,w : Các thành phần vận tốc gió theo trục Ox, Oy, Oz Wc : Vận tốc lắng đọng chất ô nhiễm  : Hệ số tính đến liên kết chất nhiễm với phần tử khác mơi trường khơng khí  : Hệ số tính đến biến đổi chất ô nhiễm thành chất khác trình phản ứng hoá học xảy đường lan truyền Tuy nhiên pt (20) phức tạp hình thức mơ lan truyền ô nhiễm Trên thực tế để giải phương trình người ta phải tiến hành đơn giản hoá sở thừa nhận số điều kiện gần cách đưa giả thuyết phù hợp với điều kiện cụ thể sau: - Nếu hướng gió trùng với trục Ox thành phần tốc độ gió chiếu lên trục Oy 0, có nghĩa v = - Tốc độ gió thẳng đứng thường nhỏ nhiều so với tốc độ gió nên bỏ qua, có nghĩa w = Trong nhiều trường hợp, xét bụi nhẹ Ws = (trong trường hợp bụi nặng lúc ta cho Ws 0) - Nếu bỏ qua tượng chuyển pha (biến đổi hố học) chất nhiễm không xét đến chất ô nhiễm bổ sung trình khuyếch tán     Như sau giả thiết chấp nhận số điều kiện gần phương trình ban đầu viết dạng là: � C u � C  �� � C� ��k � C� k � � �z � (2) y � � � � � � t � x � y� � y z z � � � � Nếu giả sử hệ số k y , k z không đổi pt (2) viết lại : 2C 2C � C u � C k � �  k z y � t � x � y2 � z2 Trong trường hợp khơng tính đến thành phần phi tuyến u 2C 2C � C k � �  kz y � t � y2 � z (3) � C phương trình (3) viết là: � x (4) Ta thấy phương trình (4) dạng phương trình truyền nhiệt chiều Tuỳ theo điều kiện ban đầu điều kiện biên mà ta có nghiệm giải tích khác Để tìm nghiệm giải tích phương trình (4), xét tốn truyền nhiệt chiều có dạng sau: � u  a2 �2u � x  � � t � x2 ,t=0 (5) � x  � Với điều kiện ban đầu : u ( x, t )   ( x)  ( x) : hàm liên tục Đặt u(x, t) = X(x)T(t) vào phương trình truyền nhiệt ta XT '  a X ''T X ''  T '    const X a2T hay ''  X  T ' a2 2T  Từ suy : X Nghiệm phương trình (7) X1  C1ei x X  C2ei x (6) (7) (8) (Xem cách giải phương trình trang 53 [7]) Nghiệm phương trình (8) T  C3ea2 2t Xem cách giải phương trình (8) trang 262 [6] Khi nghiệm phương trình vi phân (5) có dạng u ( x,t )  A( )e 2a2t �i x  số thực  �   � Vì ta chọn dấu dương phương trình (9) lập hàm số � u( x,t )  A( )ea2 2t i xd  � � (9) (10) Nếu đạo hàm phương trình (5) tính cách vi phân thành phần dấu tích phân (10) có nghĩa phương trình (10) thoả mãn phương trình (5) hay phương trình (10) nngiệm phương trình (5) Ngồi ta cịn phải thoả mãn điều kiện ban đầu t = Khi ta có:  ( x, t )  � i x d  A( )e � � (11) Sử dụng công thức tính tích phân Fourier ngược ta � A( )  �  ( )ei d 2 � (12) thay (12) vào (10) ta � � � �  i  � � a2 2t i x d  u( x,t )   (  ) e d  e � � 2 �� � ��� � � �� � a2 2t i ( x ) �  �� e d   ( )d � � � 2 ��� � � Xét tích phân Như 2 e I  ea  t i ( x  )d   2  a2t � u( x,t )  � e �2  a2t e Đặt G( x, ,t )   a2t Ta có u( x,t )  ( x )2 4a2t ( x )2 4a2t  ( )d (13) ( x )2 4a2t � G( x, , t ) ( )d � � (14) Hàm số G ( x,  , t ) gọi nghiệm sở phương trình truyền nhiệt Hàm số thoả mãn phương trình truyền nhiệt theo biến (x,t) kiểm tra trực tiếp cách lấy đạo hàm:  e ( x   )2 Gx   x   2(a2t )3/2 4a2t � � ( x   )2 ( x   ) 1 � Gxx    5/2 � e � 3/2 2� 4a2t ( a t ) 4( a t ) � � � ( x   )2 � ( x   )2 a a � � Gt   3/2  e � 5/2 2� 4a2t 2( a t ) 4( a t ) � � Vậy Gt  a Gxx Trở lại với phương trình lan truyền nhiễm chiều () viết lại với nguồn thải Q x = � C  k �2C x � t � x2 (15) Đặt a  k x nghiệm phương trình (15) viết lại là: C ( x,t )  Q e  tkx1/2  x2 4tkx (16) Đây nghiệm cảu tốn lan truyền nhiễm chiều với nguồn thải Q Cùng với điều kiện x � � C � (Nồng độ ô nhiễm điểm giảm điểm tiến xa khỏi chân biên nguồn thải ) Đối với toán hai chiều ta có phương trình tương tự C ( x, y,t)  Q e 1/2 4( t )(kxk y ) � � � � � � � � � � � � � � � � � � 2 �� � 1 x y � � 4t kx k y � � � � � (17) � � Đối với tốn chiều ta có: C ( x, y, z, t )  Q e 8( t )3/2 (kx k y kz )1/2 � � � � � � � � � � � � � � � � x  y2  z � � 1� � � 4t � k k k � x z �� y � (18) Trong công thức Q – lương phát thải chất ô nhiễm nguồn điểm tức thời, g kg II Công thức xác định phân bố nồng độ chất ô nhiễm theo luật phân phối chuẩn Gauss II.1 Công thức sở: Lượng chất ô nhiễm luồng khói xem tổng hợp vộ số khói tức thời, khối gió mang nở rộng khí xa ống khói giống ổ bánh mì cắt thành nhiều lát mỏng xếp chồng kề mép lên (hình 1) Lượng chất ô nhiễm lát mỏng luồng khói xem nhau, tức bỏ qua trao đổi chất từ lát sang lát kề bên trục x Từ cách lập luận đó, tốn lan truyền chất nhiễm tốn hai chiều ta chọn công thức (17) để áp dụng cho tr ường hợp này: a) c) u d) b) Hình 1:Biểu đồ luồng khói khối tức thời liên tục Nếu ta thiết lập cân vật chất “lát” khói có bề dày 1m theo chiều x vá chiều y, z vô cực lát khói chuyển động với vận tốc gió u thời gian để lát qua khỏi ống khói m/u lượng chất nhiễm chứa “lát” khói Q = M x 1/u Ngoài ra, cấn lưu ý toán hai chiều chiều y z thay cho chiều x y cơng thức (17) Khi cơng thức (17) trở thành : C Đặt : M e 4 ut (k y k z )1/2 � � � � � � � � � � � � � � � 2 �� �  y z � � � 4t k y kz � � � k y  0.5 y u x u kz  0.5 z x x t u (19) � � (20) (21) (22) Trong  y  z - gọi hệ số khuyếch tán theo phương ngang phương đứng, có thứ nguyên độ dài m Thay (20), (21), (22) vào (19) ta được: C M e 2 u y z � � � � � � � � � � � � � � � � � � 2 ��  y  z � � � � 2 y2 2 z2 � � �  � � � M e 2 u y z � � � � � � � � � � � � z2 �  y � � � � � 2� 2 y2 � � �2 z � e � � � (23) Đây công thức sở mơ hình lan truyền chất nhiễm theo luật phân phối chuẩn Gauss hay gọi “mơ hình Gauss” sở II.2 Diễn giải cơng thức mơ hình Gauss sở phương pháp phân tích thứ ngun: Cơng thức (23) cịn diễn giải phương pháp phân tích thứ nguyên sau: Từ miệng ống khói chất nhiễm gío mang theo trục x trùng với hướng gió với vận tốc vận tốc gió u, m/s Nếu lượng phát thải chất ô nhiễm M, g/s không đổi theo thời gian mật độ chất nhiễm tất mặt cắt trực giao với trục gió (cũng trục luồng khói) M/u, g/m Cường độ phát thải M = đơn vị/s u = m/s u = m/s 4 Khoảng cách dọc theo trục gió (x), m Hình 2:Sơ đồ minh hoạ ảnh hưởng vận tốc gió đến nồng độ chất nhiễm nguồn phát thài liên tục số gây Nếu giả thiết chất nhiễm khơng có phản ứng hố học với khơng khí xung quanh tức khơng sản sinh khơng phân huỷ đi, mật độ chất ô nhiễm tất mặt cắt trực giao với trục gió khoảng cách x hình Nhưng nồng độ chất nhiễm luồng khói giảm dần khoảng cách x tăng có tượng khuyếch tán theo phương ngang (trục y) theo phương đứng (trục z) mà luồng khói lan rộng xung quanh trục luồng Càng xa khỏi trục luồng theo phương y z theo phương y z nồng độ giảm nhỏ, tức nồng độ nghịch biến với khoảng cách y z M Từ ta viết C : (24) uyz Bằng nghiên cứu lý thuyết thực nghiệm người ta thấy phân bố nồng độ mặt cắt trực giao với trục luồng theo chiều ngang y theo chiều đứng z tuân theo dạng hình chuông luật phân phối chuẩn Gauss với sai phương chuẩn  Khi biểu thức phân phối chuẩn Gauss có dạng   ( )  e  2 2 (25) Áp dụng biểu thức (25) vào trường hợp cụ thể  y z hàm  ( y ) ,  ( z ) nghịch biến với |y| |z| Do biểu thức (24) ta viết thành: M C  M  ( y)  ( z )  u 2 u y z  y2  z2 2 y 2 2z e e (26)  y  z hệ số khuyếch tán theo phương ngang y phương đứng z hàm số khoảng cách x kể từ nguồn đến mặt cắt xem xét Các hệ số xác định thực nghiệm phụ thuộc vào khoảng cách x với điều kiện khác Chính dấu tỉ lệ biểu thức (24) thay dấu = đằng thức (26) Biểu thức (26) nghiệm cách giải ph ương trình vi phân đạo hàm riêng trình khuyếch tán II.3 Sự biến dạng mơ hình Gauss sơ: Điều cần lưu ý trước tiên công thức (19),(23), (26) toạ độ y z tính từ trục luồng khói Khi chuyển hệ trục x, y, z mà gốc O trùng với chân ống khói mặt đất y khơng thay đổi z phải thay z - H H – z (hình 3), công thức (26) trở thành:  y2 ( z  H )2 2 y M 2 z (27) C 2 u y z e e Ngoài tuỳ thuộc theo độ xa x luồng khói nở rộng chạm mặt đất mặt đất cản trở không cho luồng tiếp tục phát triển, ngược lại chiều hướng khuyếch tán bị mặt đất phản xạ ngược trở lên thể có nguồn ảo hồn tồn đối xứng qua mặt đất mặt đất xem gương phản chiếu (hình 3) Để kể đến ảnh hưởng mặt đất phản xạ khuyếch tán, nồng độ điểm A, B giả thiết hai nguồn giống hệt gây ra, có nguồn thực nguồn ảo hoàn toàn đối xứng với qua mặt đất Nồng độ điểm xem xét (A B) nguồn thực gây tín cơng thức (27), cịn nguồn ảo gây tính biểu thức : C M 2 u y z  y2 ( z  H )2 2 e 2 y e 2 z (28) Nồng độ tổng cộng tính từ (27), (28) là:  y2 � ( z  H )2 ( z  H )2 � � � 2 y � 2 2z M �  z C e e e � � (29) 2 u y z � � � � � Đây cơng thức tính tốn khuyếch tán chất nhiễm từ nguồn điểm cao liên tục Khi tính tốn nồng độ ô nhiễm mặt đất z = công thức (29) trở thành: C M  u y z  y2 ( H )2 2 e 2 y e 2 z (30) Trường hợp tính phân bố nồng độ mặt đất dọc theo trục gió (trục x) , ta cho y = thu C II.4 M  u y z ( H )2 e 2 z (31) Hệ số khuyếch tán : Từ biểu thức (20) (21) ta có: 1 �2k x � �2k x � y � � z (32) � � y � z � � � u u � � � � � � Như  y  z phụ thuộc vào khoảng cách x, độ rối khí vận tốc gió Pasquill Gifford thực nghiệm thiết lập mối quan hệ hệ số  y ,  z phụ thuộc vào khoảng cách x xi theo chiều gió ứng với mức độ ổn định khí khác A, B, C, D, E F Mối quan hệ đuợc cho dạng biểu đồ (hình 3.9 hình 3.10 trang 83, 84 sách [1]) Tuy nhiên để thuận tiện tính tốn lập trình D.O.Martin đưa cơng thức tính  y ,  z sau:  y  ax0.894  z  bxc  d (33) Trong x – khoảng cách xi theo chiều gió kể từ nguồn, tính km Các hệ số a, b, c, d cho bảng Cấp định Bảng 2: Các hệ số a, b,c, d công thức (33) ổn x >1 km x  km a b c d b A B C D E F c d 213 440.8 1.941 9.27 459.7 2.094 -9.6 156 106.6 1.149 3.3 108.2 1.098 2.0 104 61 0.911 61 0.911 68 33.2 0.725 -1.7 44.5 0.516 -13.0 50.5 22.8 0.678 -1.3 55.4 0.305 -34.0 34 14.35 0.740 -0.35 62.6 0.180 -48.6 I.5 Các cấp ổn định khí quyển: Theo Pasquill Gifford cấp ổn định khí có liên quan chặt chẽ tới biến thiên nhiệt độ khơng khí theo chiều cao Tuỳ theo chiều hướng mức độ thay đổi theo chiều cao ta có trường hợp đẳng nhiệt, đoạn nhiệt, siêu đoạn nhiệt nghịch nhiệt (hình 3) CÁC CẤP KHÍ HẬU Khơng ổn định trung tính Ổn định Rất ổn định 200 180 160 140 Đoạn nhiệt (-1oC/100m) Đẳng nhiệt 120 100 80 60 Siêu đoạn nhiệt (Ví dụ -6oC/100m) Nghịch nhiệt 40 20 -10 -8 -6 -4 -2 Độ gia tăng nhiệt độ Nhiệt độ giảm Nhiệt độ tăng 10 Hình 3: Các trường hợp biến thiên nhiệt độ khơng khí theo chiều cao mặt đất III Chiều cao hiệu ống khói: x x+x Tại miệng ống khói, nhờ vận tốc phụt, luồng khói có động ban đầu làm cho có xu hướng bốc thẳng lên Mặc khác, nhiệt độ khói cao nhiệt độ khơng khí xung quanh, luồng khói chịu tác dụng “lực nổi” chênh lệch nhiệt độ gây Cùng với lực nâng, luồng khói chịu tác động lực gió nằm ngang, đỉnh cao luồng khói nằm cách ống khói khoảng cách định xi theo chiều gió Khi đạt độ cao tức lúc động ban đầu luồng khói bị triệt tiêu nhiệt độ khói trở nên cân với nhiệt độ khí kết q trình hồ trộn với khơng khí xung quanh, luồng khói giữ phương nằm ngang song song với chiều gió Dựa kềt thực nghiệm Davidson W.F đưa công thức sau –cịn gọi cơng thức Brayant – Davidson 1.4 � � � � �   T � h  D � 1 ,m (34) � �u � � � T � � � � khoi � � Cơng thức phân biệt thành thành phần: Trong cơng thức trên: D - đường kính miệng ống khói, m;  - vận tốc ban đầu luồng khói miệng ống khói, m/s u - vận tốc gió , m/s; Tkhoi - Nhiệt độ tuyệt đối khói miệng ống khói, K; T - Chênh lệch nhiệt độ khói khơng khí xung quanh, oC K *Cơng thức tính  : L S (35) Trong L – lưu lượng khí thải (m3/s) S - Diện tích miệng ống thải (m2) Lúc chiều cao hiệu ống khói H = h + h ,m với h chiều cao thực ống khói IV Sự lắng đọng bụi q trình khuyếch tán khí thải từ nguồn điểm cao: Những cơng thức tính tốn khuyếch tán nêu áp dụng cho chất khí Đối với bụi nhẹ lơ lửng, cách gần xem vận tốc rơi chúng dướI tác dụng trọng lực không đáng kể mức độ khuyếch tán chúng gần khí, lúc ta áp dụng cơng thức để xác định nồng độ bụi mặt đất Tuy nhiên, khí thải có chứa bụi với thành phần cỡ hạt khác nhau, vận tốc rơi cỡ hạt bụi thơ nặng kích thước   20  m đáng kể, chúng lắng đọng nhanh xuống mặt đất vùng gần chân khói xi theo chiều gió Như có khác biệt đáng kể nồng độ bụi nồng độ khí mặt đất Dựa vào mơ hình Gauss sở - tức mơ hình Gauss chưa kể đến phản xạ mặt đất chất ô nhiễm thể công thức (27) Đối với đa số chất ô nhiễm thể khí mặt đất khơng hấp thu mà phản xạ ngược trở lại vào khí Cịn bụi ta xem mặt đất vật hấp thụ hồn tồn Ngồi ra, chất nhiễm thể khí khơng chịu ảnh hưởng lực trọng trường, cịn bụi rơi khí với vận tốc vr định tuỳ thuộc vào kích thước hạt khối lượng đơn vị Do đại lượng H mơ hình Gauss cần hiệu chỉnh cách trừ bớt đoạn đường mà hạt 10 bụi rơi khoảng thời gian t Đoạn đường vr t mà t = x/u vớI x khoảng cách dọc theo trục gió tính từ nguồn u – vận tốc gió Do đó, cơng thức (27) hiệu chỉnh thành:  y2 ( z  H  vur x )2 Mb 2 (36) C  e y e 2 z b 2 u y z Nồng độ bụi mặt đất dọc theo trục gió Cb( x)  Mb 2 u y z ( H  vur x )2 e 2 z (37) Trong đó: Cb : Nồng độ bụi tính theo g/m3 Mb: Lượng phát thảI bụi thuộc nhóm cỡ hạt cầt xem xét, g/s vr : Vận tốc rơi tới hạn trung bình nhóm cỡ hạt bụi xem xét m/s x : Khoảng cách dọc theo trục gió kể từ nguồn, m Cường độ lắng đọng bụi mặt đất dọc theo trục gió Gb(x) Khối lượng vận chuyển (Lưu lượng)(Nồng độ) Gb(x) = = Diện tích (Diện tích) = (Vận tốc) (Nồng độ) = vrCb(x) Gb( x)   M bvr 2 u y z ( H vur x )2 e 2 z (38) Công thức tính vận tốc bụi: b g vr  18 (39) Với b - Khối lượng đơn vị bụI, kg/m3  - Đường kính hạt bụi, m  - Hệ số tính theo cơng thức (40), Pa.s Cơng thức tính  : 3/2 � � toC  0o C 387 �273  t � 387  t � 273 � � � (40) 6 với  0o C  17,17.10 Pa.s V Tính tốn khuyếch tán chất ô nhiễm từ nguồn điểm cao theo ph ương pháp BERLIAND M.E: Ở Liên Xơ cũ, cơng trình nghiên cứu lý thuyết khuyếch tán tiến hành Đài Địa lý - Thiên văn Trung ương dựa phương trình vi phân trình khuyếch tán từ nguồn điểm cao viết dạng : 11 � 2C � �� � � C � C � C� � � �  C u  w  �k y �  k z � � � � x � z � � z � z � y � � � � � (41) C - Nồng độ chất nhiểm u - Vận tốc gió  - Vận tốc theo phương thẳng đứng chất ô nhiễm kz ky - hệ thống trao đổI theo phương đứng phương ngang  - Hệ số xác định thay đổi nồng độ chất ô nhiễm phân huỷ hoá học gội mưa, sương Các giả thiết ban đầu thừa nhận là: nguồn điểm đặt độ cao z = H (khi x = 0); có phản xạ hồn tồn chất nhiễm từ mặt đất ; nồng độ chất ô nhiễm tiến dần đến triệt tiêu khoảng cách tương đối xa so với nguồn V.1 Đối với khí bụi nhẹ: Với giả thiết trên, Berlaind chứng minh đưa kết sau: Nồng độ ô nhiễm mặt đất khí bụi nhẹ: C( x, y,z 0)  M e 2(1 n)k1  k0 x � � � � � � � � � u H1n   y � k1(1n)2 x 4k0 x � � � � Nồng độ cực đại mặt đất: 2M k1 0.116(1  n ) Cmax  k0u1 u1H 1.5(1n) (43) u1H 1n X max  k1(1 n)2 (44) Khoảng cách từ XM từ nguồn đến vị trí có nồng độ max: Trong công thức : u1 – vận tốc gió độ cao z = 1m, m/s n = 0.15  0.2 k1= 0.1  0.2 m/s k0= 0.5  m 12 (42)  V.2 Cơng thức tính vận tốc độ cao z: �Z � uz  u10 � � 10 � � Với : uz - Vận tốc gió độ cao z u10 - Vận tốc gió độ cao 10m Đối với bụi nặng cỡ hạt đồng chất: Nồng độ bụI mặt đất: (45) MH  (1n)u 1 C( x, y, z 0)  e     2(1 n) (1  )(k1x)  k0 x (46) Trong  vr (1  n)k1 � � � � � � � � � u H1n   y � � (1n)2 k1x 4k0 x � � � (47) Nồng độ bụi cực đại mặt đất: 0.063(1  n) M k1 (1.5   )1.5  Cmax  u1 H 1.5(1 n ) k0u1 e (1   ) Khoảng cách từ nguồn đến vị trí có nồng độ max: (48) u1H 1n xM  (1 n)2 (1.5   )k1 (49) � � C � � C ( Rkr )  kz  M  (r ) ( z  H )  R� r � r � z � z (50) Trong (1   ) hàm số Gamma (1   ) (xem phụ lục 2) V.3 Khuyếch tán chất ô nhiễm từ nguồn điểm cao điều kiện khơng gió: Các phương pháp tính tốn khuyếch tán chất nhiễm nêu áp dụng cho trường hợp trời có gió Để kể đến ảnh hưởng khuyếch tán rối trường hợp trời khơng gió, Berlaind Kurebin đưa phương trình sau hệ toạ độ trụ: Với điều kiện biên sau đây: � C  R  z � �: C � - Khi z = 0: k z � r - Sự phân bố nồng độ ô nhiễm mặt nằm ngang có tính đối xứng qua tâm nguồn � C 0 R  � � r Berlaind Kurebin giảI phương trình (50) vói nghiệm có dạng: 13 C ( R, z)  a( H 1n  z1n )  R M 3/2 3/2 2 k1(1 n) � (1n)/2 (1n)/2 � � � (1n)/2 (1 n)/2 a( H z )  R � �a( H z )  R� � � (51) Trong a � � 4k1 (1  n)2 � (52) Nồng độ mặt đất: CR  M 2 k1(1 n)(aH 1n  R)2 (53) Nồng độ cực đại mặt đất (tại chân ống khói ) M (1  n ) Cmax  32 k13H 2(1n) (54) VI Tính tốn nồng độ trung bình chất ô nhiễm mặt đất nguồn thải gây ra: Khi tính tốn dự báo mức độ nhiễm địa điểm nguồn thải khác gây ra, việc xác định nồng độ ô nhiễm tức thời, ta cần phải dự báo phân bố nồng độ trung bình ngày đêm, trung bình tháng trung bình năm Trong trường hợp tính tốn nồng độ trung bình cho thời gian ngắn, trung bình ngày đêm ta đơn giản hoá vấn đề cách giả thiết mùa định, hè đông, cấp ổn định khí thay đổI ngày đêm xung quanh cấp ổn định trung bình ta tính tốn với cấp ổn định Ngồi ra, cấp vận tốc gió thay trị số vận tốc gió trung bình uTB ( ) theo hướng  vớI tần suất xuất gió P( ) hướng  Ngồi tần suất hướng gió khác cịn có tần suất lặng gió P lặng Đó tỷ lệ thời gian khơng có gió hướng (u = 0)  Hệ số trung bình ứng với số liệu tần suất tần suất lặng gió: k  P � 1 Plang � (52) � � � � Trong P tần suất gió hướng  Plang tần suất lặng gió m - Số hướng gió thơng thường với Plang < m P  1  1 Công thức xác định nồng độ trung bình theo tần suất gió: Từ lập luận nêu đây, ta viết biểu thức xác định nồng độ trung bình ngày đêm chất ô nhiễm mặt đất vị trí tính tốn nguồn thải thứ i gây sau :  14 m Cxy (i)  Plang Clang (i )  �k C (i) (53)  1 Nồng độ tổng cộng trung bình điểm có toạ độ x, y n nguồn thải gây : n m i 1  1 Cxy (tong )  �(Plang Clang (i)  (1 Plang ) �k C (i ) ) (54) Trong cơng thức trên: Cxy(i) - nồng độ trung bình vị trí có toạ độ x, y nguồn thứ i gây Cxy(tong) - nồng độ tổng cộng trung bình n nguồn thải gây điểm tính toán Clang(i) - Nồng độ tức thời nguồn thải thứ i gây điểm tính tốn lặng gió (u = 0) C(i) - Nồng độ tức thời nguồn thải thứ i gây điểm tính tốn có gió thổi theo hướng  ứng với vận tốc gió trung bình hướng độ ổn định trung bình khí suốt khoảng thời gian tính tốn trị số trung bình (ngày đêm, tháng năm) VII Xác định nồng độ tương đối tổng cộng mặt đất nhiều nguồn điểm cao gây ra: Trong khí thải cơng nghiệp bao gồm khói thải từ lò nung, lò đốt, lò luyện thép v.v… ln có mặt đồng thời nhiều yếu tố độc hại khác : khí SO2, CO, CO2, NOx, H2Sv.v…và bụi Khi khí thải lan toả mơi trường xung quanh yếu tố độc hại nêu song song tồn khơng khí gây tác hại cách tổng hợp môi trường sống người sinh vật nói chung Tuy nhiên, từ trước đến xem xét đánh giá chất lượng mơi trường khơng khí khu vực đồ thị, nhà máy hay khu cơng nghiệp đó, người ta thường khảo sát đo đạc tính tốn nồng độ ô nhiễm yếu tố độc hại cách riêng rẽ mức độ ô nhiễm môi trường khơng khí đánh giá cách độc lập chất ô nhiễm khác Mặc khác tiêu chuẩn môi trường Việt Nam đưa giới hạn nồng độ cho phép chất ô nhiễm riêng biệt khơng khí xung quanh khơng khí bên nhà mà chưa có qui định giới hạn cho phép nhiều chất ô nhiễm đồng thời có mặt khơng khí Điều dẫn đến thiếu xác bất hợp lý việc xem xét, đánh giá so sánh chất lượng môi trường nơi nơi khác, khu công nghiệp khu công nghiệp khác Cụ thể có hai khu vực khác nhau, khu vực môi trường bị ô nhiễm chất với nồng độ nồng độ giới hạn cho phép, lúc khu vực mơi trường bị ô nhiễm đồng thời nhiều chất độc hại khác mà nồng độ chất đạt giới hạn cho phép theo tiêu chuẩn hành, hai khu vực nói xem đạt yêu cầu chấp nhận nhau, thực tế, môi trường khu vực chắn gây ảnh hưởng xấu đến sức khoẻ người nhiều khu vực Để kể đến tác dụng tổng hợp đồng thời nhiều chất gây ô nhiễm Người ta áp dụng khái niệm nồng độ tương đối tổng cộng C0 nhiều chất ô nhiễm tồn đồng thời môi trường khơng khí Nồng độ C0 - đại lượng khơng thứ nguyên biểu diễn công thức: 15 C0  C1 C    Cn Ccf (1) Ccf (2) Ccf (n) (55) Trong đó: C1, C2,…, Cn nồng độ thực tế đo tính tốn chất ô nhiễm 1, 2, 3, …,n, mg/m3 Ccf(1), Ccf(2) ,…, Ccf(3) : giới hạn nồng độ cho phép chất 1, 2,…,n riêng biệt (xem phụ lục 3) Đại lượng C0 xem số tổng hợp đánh giá chất lượng môi trường không khí mặt nhiễm hố học Nó cho phép ta đánh giá, so sánh chất lượng môi trường không khí theo tiêu thống 16 ... lượng môi trường nơi nơi khác, khu công nghiệp khu công nghiệp khác Cụ thể có hai khu vực khác nhau, khu vực môi trường bị ô nhiễm chất với nồng độ nồng độ giới hạn cho phép, lúc khu vực môi trường... máy hay khu công nghiệp đó, người ta thường khảo sát đo đạc tính tốn nồng độ nhiễm yếu tố độc hại cách riêng rẽ mức độ ô nhiễm mơi trường khơng khí đánh giá cách độc lập chất ô nhiễm khác Mặc... Trong công thức Q – lương phát thải chất ô nhiễm nguồn điểm tức thời, g kg II Công thức xác định phân bố nồng độ chất ô nhiễm theo luật phân phối chuẩn Gauss II.1 Công thức sở: Lượng chất nhiễm