Mô hình hóa ô nhiễm không khí

MƠ HÌNH HĨA Ơ NHIỄM KHƠNG KHÍ PGS.TSKH Bùi Tá Long , Khoa Môi trường Tài nguyên, Đại học Bách khoa Tp Hồ Chí Minh Giám sát nhiễm khơng khí C C   C   C u  Ky  Kz  t x y  y  z  z  Lan truyền ô nhiễm không diện rộng Toshihiko Takemura  Thực tế tình trạng nhiễm khơng khí xun biên giới phụ thuộc chủ yếu vào điều kiện thời tiết Ví dụ, khơng khí nhiễm thường lan từ Trung Quốc sang Nhật Bản vào mùa Xn  Ơ nhiễm khơng khí châu Á thực vấn đề xã hội nghiêm trọng Cần phải tính đến lộ trình tối ưu nhằm giảm thiểuđến mức tối đa nhiễm khơng khí tình trạng trái đất ấm lên  Ơ nhiễm khơng khi, va đặc biệt việc giảm ô nhiễm không đến mức chấp nhận được, vấn đê mơi trường quan trọng;  Ơ nhiễm khơng khơng giới hạn vùng có nguồn thải lớn  Các chất ô nhiễm không lắng đọng bê mặt khơng khu vực có nguồn thải lớn mà khu vực xung quanh  http://www.ft.com/cms Vai trị mơ hình tốn Vấn đê cần quan tâm Cần xác định xác loại cần giảm va nơi cần giảm Nồng đô va mức lắng đọng chấp nhận (hoặc tới hạn) cần xác định xác va nhiễm khơng cần giảm đến mức đô không nên giảm hơn, (do chi phí thêm có thê đắt, va vậy, có thê gây nên kho khăn lớn vê mặt kinh tê)  SẤM CHỚP MÂY NÚ MƯA Nhà I KHU CÔNG NGHIỆP Cây xanh  N Ơ N G Thành phố đẹp Cây xanh N G H I Ệ P Cây Cơng viênxanh SƠNG ĐẤT Ô TÔ NÔNG NGHIỆP $ Cây I NÚ RỪNG GIĨ NẮNG Người Các mơ hình tốn biểu diễn tượng ô nhiễm không công cụ không thê thiếu q trình nơ lực giải toán phác họa Để giảm tối ưu mức ô nhiễm không đến mức độ chấp nhận có thê giải thành cơng sư dụng mơ hình tốn học nhiễm không đáng tin cậy PHÁT THẢI TẢI + PHÂN TÁN LẮNG ĐỌNG NGUỒN NGUỒN TIẾP NHẬN ĐẤT Nhà thờ Xe hàng Nhà 10 Vì cần phải sử dụng mơ hình Lợi ích dùng mơ hình 11 Mơ hình khí tượng, phát thải phát tán Các bước mơ hình hố nhiễm khơng khí 14 Sự phân bố nhiệt độ theo độ cao 2000 m Khí tiêu chuẩn = 0,65 K/100 m Trạng thái khí Khí đoạn nhiệt khơ = = K/100 m Giữa buổi chiều Bình minh + h Bình minh + h 300 m Lúc sáng sớm Khí siêu đoạn nhiệt > 15 20 ·C 35 ·C Nhiệt độ 15 Phân bố nhiệt độ theo chiều cao vào thời gian khác ngày thời tiết nắng ráo, khơng mây gió nhẹ 16 Trạng thái khí  Với trạng thái biến thiên nhiệt độ theo chiều cao mà lực tác động vật lý khí làm cho chất nhiễm khuếch tán theo chiều cao khơng dễ dàng gọi khí trạng thái ổn định, tức khí ổn định cản trở khuếch tán pha lỗng chất nhiễm;  Ngược lại với profile biến thiên nhiệt độ theo chiều cao mà hịa trộn khơng khí theo chiều cao dễ dàng làm cho chất nhiễm khơng khí khuếch tán (pha lỗng) khí dễ dàng gọi khí trạng thái không ổn định  Ở hai trạng thái trung tính 17 CÁC HÌNH DẠNG KHÁC NHAU CỦA LUỒNG KHĨI 18 Hình dạng khác luồng khói       Looping – luồng khói uốn lượn; Coning – luồng khói hình cơn; Fanning – Luồng khói hình quạt; Lofting – luồng khói khuếch tán mạnh biên Fumigation – luồng khói “xơng khói”; Trapping – luồng khói mắc kẹt; 19 20 Uốn lượn (looping) Độ cao      Khi không ổn định mạnh (khi phân bô nhiệt đô theo chiều cao siêu đoạn nhiệt); Thường xảy vào ban ngày mặt trời đốt nóng mặt đất với cường xa lớn; Luồng khói uốn lượn – khơng ổn định mạnh Chiều cao Hình (coning)   Hình thành điều kiện trung tính gần trung tính trời có mây che phu làm cho mặt trời vào ban ngày bị che phu hay xa hồng ngoại tư mặt đất vào ban đêm bị giảm; Góc mơ khoảng 20 21 Hình quạt (fanning) Sự biến thiên nhiệt độ khối khí ( Hình quạt   ,   Độ cao ổn định Vồng lên (lofting) Xảy điều kiện ổn định với phân bô nhiệt theo chiều cao nghịch nhiệt; Rối theo chiều đứng bị triệt tiêu, có phát triển theo chiều ngang; Luồng khói hình quạt- lớp nghịch nhiệt tư bê mặt đất đến đô cao ống khói Chiều cao Sự biến thiên nhiệt độ môi trường (β) Nhiệt độ 22    0 Phía sát mặt đất phạm vi ống khói có nghịch nhiệt, cịn phía khuếch tán mạnh (siêu đoạn nhiệt); Chất nhiễm tích tụ gần mép lớp nghịch nhiệt; Chất ô nhiễm khuếch tán thuận lợi phía trên;  Thường xảy mặt trời xuống 23 24 Lớp nghịch nhiệt Xơng khói (fumigation)  0   Ngược lại với trên:lớp nghịch nhiệt bên lớp siêu đoạn nhiệt; Luồng khói khuếch tán mạnh phía dưới; Thường xảy vào buổi sáng sớm; Luồng khói khuếch tán mạnh biên 25 … Tiếp theo 26 … Tiếp theo: nghịch nhiệt 27 28 Key words: tư khóa cần nhơ Mắc bẫy (trapping)  Adiabatic (đoạn nhiệt): tính chất mơi trường khơng khí sư trao đổi nhiệt với môi trường không đu lớn đê diễn sư cân nhiệt nên có thê bỏ qua;  Gradian nhiệt đô đoạn nhiệt: đô hạ tăng nhiệt đô khối khi lên cao xuống thấp;  Đô ổn định quyển: khơng ổn định, trung tính, ổn định; Nghịch nghiệt bên bên ống khói Luồng khói bị hạn chê biên lẫn biên Chất ô nhiễm kho khuếch tán lên lẫn xuống 29 30 MƠ HÌNH GAUSS - PASQUILL 31 Q trình lan truyền chất nhiễm Đối tượng thật Độ cao hữu dụng 34 Các bước mơ hình hóa ô nhiễm không khí        Chọn hệ tọa độ Chọn hệ tọa độ Xây dựng mơ hình ý niệm Xác định biến tham số Xây dựng phương trình, thiết lập tốn Tin học hóa thành phần mềm Hiệu chỉnh, kiểm định Ứng dụng vào thực tế 35 Xây dựng mơ hình ý niệm Xác định biến tham số 37 Vệt khói tức thời trung bình Xây dựng phương trình, thiết lập tốn 40 10 Mơ hình Gauss sở (Basic Gaussian formula) cho nguồn thải cao Giải toán biên U  M  y z  H  UC (,x ,y )z x0 C  x, y , z  y  Kz u u Ky 0.5y2 ,Kz 0.5z2 x x C   C    C   K   K  , C  x, y , z  x y  y y  z  z z  C z z 0  C  x, y , z  z  0 0   z  H   U    y 1   C  x, y, z   C0 x exp   Kz  4x    K y   C  x, y , z    zH     exp     2 z   49     50 Gauss biến đổi (Gaussian plume formula) Mơ hình Gauss sở cho nguồn thải mặt đất u u Ky 0.5y2 ,Kz 0.5z2 x x   y2 M z2 C  x, y , z   exp     2 u  y  z    y  z  y2   z2  M exp    exp     2   2  u  y  z y   z    y2 M exp    2 u y  z y   y2   M  C(,x ,y )exp z  exp     2   2uyz y           z  H 2  exp   22z      z  H 2      2z    Phản xạ hoàn toàn mặt phẳng z =  y2 M C (,x ,y )exp z  exp   2 u y  z y  y , z - gọi hệ số phạm vi khuếch tán theo phương ngang va phương đứng, có thứ nguyên độ dài (do Ky , Kz – có thứ nguyên m2/s)       zH 2     exp  22z     z  H        22z     Hấp thụ hoàn toàn mặt phẳng z = 51 52 13 Tính nồng độ chất ô nhiễm mặt đất Gauss tức thời (Gaussian puff formula)  xut2 y2   zH2   zH2    exp   exp  exp 2 2 xyz  2x 2y   2z   2z2  Q C(x, y,z,t) C ( x, y)  3/2  y2   H2  exp   exp     u y z  2 z   2 y  M Nồng đô mặt đất dọc theo trục gio (trục x) Phản xạ hoàn toàn mặt phẳng z =  xut2 y2   zH2   zH2   exp     C(x, y,z,t) 3/2 exp    exp 2 xyz  2x2 2y2   2z2   2z2  Q C( x)   H2  exp   u y z  2 z  M Hấp thụ hoàn toàn mặt phẳng z = 53 Ý nghĩa hệ số phạm vi khuếch tán Hệ số phạm vi khuếch tán ngang đứng 1/  2Ky x    u   y   1/  2K z x    u  z   54 y , z - phụ thuộc vào khoảng cách x, độ ổn định khí vận tốc gió ()X, Y, Z Z X (X, 0, 0) H (X, Y, Z) h Y 55 X 56 14 Cơng thức tính y(x), z(x) cho vùng thành thị Cơng thức tính y(x), z(x) cho vùng nông thôn Loại ổn định y(x) z(x) A 0.22x(1+0.0001x)-0.5 0.20x B 0.16x(1+0.0001x)-0.5 0.12x C 0.11x(1+0.0001x)-0.5 0.08x(1+0.0002x)-0.5 D 0.08x(1+0.0001x)-0.5 0.06x(1+0.0015x)-0.5 E 0.06x(1+0.0001x)-0.5 0.03x(1+0.0003x)-1 F 0.04x(1+0.0001x)-0.5 0.16x(1+0.0003x)-1 Độ ổn định y(x) A- B 0.32x(1+0.0004x)-0.5 0.24x(1+0.0001x)0 C 0.22x(1+0.0004x)-0.5 0.12x D 0.16x(1+0.0004x)-0.5 0.14x(1+0.0003x)-0.5 E-F 0.11x(1+0.0004x)-0.5 0.08x(1+0.00015x)-0.5 z(x) 57 Mức ôn đinh cua quyên       58 Phân loại độ bền vững theo Pasquill A -rất không bền vững B - không bền vững loại trung bình C - khơng bền vững loại yếu D - trung hòa E - bền vững yếu F - bền vững loại trung bình Vận tốc gió độ cao 10 m Điều kiện thời tiết ban ngày Điều kiện thời tiết ban đêm Bức xạ mặt trời ban ngày Độ che phủ ban đêm (hệ số mây) Mạnh (biên độ > 600) Trung bình (Biên độ 35600) Yếu (Biên độ 15 – 350) Lớn 50% Nhỏ 50% F 6 C C–D D D D Lưu ý А – không ổn định; В – không ổn định vừa phải; С – không ổn định yếu; D – điều kiện trung tính; E – điều kiện ổn định yếu; F – điều kiện ổn định vừa phải 59 60 15 Công thức Davidson   h  D   u  Một sơ cơng thức tính tốn độ nâng vệt khói    hv  D   u 1/ ,m 1/  T   T khoi    , m  1/    T ht  D  ,m  u  Tkhoi D – đường kính ống khói; - vận tốc ban đầu luồng khói miệng khói, m/s; u – vận tốc gió, m/s; Tkhoi – nhiệt độ khói miệng ống khói, 0K ΔT – chênh lệch nhiệt độ khói khơng khí xung quanh, 0K 61 Công thức Holland Δh = 62 Hạn chế mơ hình Gauss Tkhoi - Txung_quanh  ωDa  -3 1.5 + 2.68.10 P.D   u  TKhoi       D – đường kính; P – áp suất khí quyển, milibar (1 atm = 1013 mbar); Chỉ ứng dụng trường hợp bề mặt tương đối phẳng; Rất khó lưu ý tới yếu tố cản; Các điều kiện khí tượng khơng đổi diện rộng; Chỉ ứng dụng tốt khí có mật độ gần với khơng khí; Vận tốc gió trung bình u > m/s Tkhoi, T xung quanh – nhiệt độ khí khơng khí xung quanh, K  a – hệ số hiệu chỉnh Cấp A , B C – nhận với hệ số 1.1 – 1.2; Cấp D, E, F – nhân hệ số 0.8 – 0.9 63 64 16 Công thức y(x), z(x) cho vùng nông thôn Step by step mơ hình Gauss C( x, y)  C ( x)   y2   H  exp   exp    2  u y z y   2 z   M  H2   exp    u  y z  2 z  M Loại ổn định u  uH H  h  h M – Lượng thải (công suất nguồn thải) (mg/s); y(x) z(x) A 0.22x(1+0.0001x)-0.5 0.20x B 0.16x(1+0.0001x)-0.5 0.12x C 0.11x(1+0.0001x)-0.5 0.08x(1+0.0002x)-0.5 D 0.08x(1+0.0001x)-0.5 0.06x(1+0.0015x)-0.5 E 0.06x(1+0.0001x)-0.5 0.03x(1+0.0003x)-1 F 0.04x(1+0.0001x)-0.5 0.016x(1+0.0003x)-1 65 Cơng thức tính vận tốc gió Cơng thức y(x), z(x) cho vùng thành thị Loại ổn định y(x) A- B 0.32x(1+0.0004x)-0.5 0.24x(1+0.001x)- 0.5 C 0.22x(1+0.0004x)-0.5 0.12x 66 z(x) D 0.16x(1+0.0004x)-0.5 0.14x(1+0.0003x)-0.5 E-F 0.11x(1+0.0004x)-0.5 0.08x(1+0.0005x)-0.5 Loại tầng kết p   z  U , z  200 m U  z    10 m  10  U  200 m , z  200 m 67 p Thành phố Nông thôn A 0.15 0.07 B 0.15 0.07 C 0.20 0.10 D 0.25 0.15 E 0.30 0.35 F 0.30 0.55 68 17 Tự động hóa tính tốn vệt nâng Tự động hóa tính tốn phát tán h  Tkhoi Txung_quanh  3 1.52.68.10 PD  u  TKhoi  Da 69 70 Tính tốn theo mơ hình Gauss C  y2 M exp    2 2 u  y z y  2         exp    z  H    exp    z  H    2  2 z  2 z      MƠ HÌNH BERLIAND 71 72 18 Tóm tắt Sự phân loại theo WMO  Mơ hình thống kê kinh nghiệm;  Mơ hình thống kê thủy động;  Mơ hình sơ trị,  Các phương trình mơ hình Berliand;  Bài tập, ứng dụng phần mềm CAP, ENVIMAP; 73 74 Mơ hình thống kê thủy động Mơ hình Gauss  Dưa sơ ly thuyêt toán học Gauss  Các nha toán học co cơng phat triên mơ hình Taylor (1915), Sutton (1925 – 1953), Turner (1961 – 1964), Pasquill (1962 – 1971), Seifeld (1975)  Gân đươc nha khoa học môi trương cua nươc My, Anh, Phap, ứng dụng hoan thiên mơ hình tính theo điêu kiên cua mơi nươc 75  Cịn gọi ly thuyết nửa thư ngun (cịn gọi mơ hình K) Mơ hình Berliand (Nga) hồn thiện va áp dụng Liên Xô;  Ở Việt Nam, sơ nha khoa học đa áp dụng mơ hình cho sơ cơng trình, dư án 76 19 Phương trình xuất phát Mơ hình sơ trị C C C C   C    C    C    Kz Vx Vy Vz  Kx   K y   αC t x y z x  x  y  y  z  z   Tưc giải phương trình vi phân băng phương phap sơ;  Viêc triên khai mơ hình Việt Nam đoi hỏi thơi gian sơ liêu cho mơ hình thiêu phương tiên tính tốn chưa đu mạnh C – nồng trung bình chất nhiễm (mg/m3 ); x,y,z – thành phần tọa đô theo trục Ox, Oy, Oz; t – thời gian; Kx, Ky, Kz – thành phần sô khuếch tán rối theo trục Ox, Oy, Oz;  Vx, Vy , Vz – thành phần tốc đô trung bình theo ba trục Ox, Oy, Oz ;  α - sơ tính đến sư biến đổi chất nhiễm thành chất khác trình phản ứng hóa học xảy đường lan truyền     77 78 Các giả thiết làm đơn giản toán Các giả thiết làm đơn giản toán Trên thực tê thành phần khuếch tán rối theo chiều gio nho nhiều so với thành phần khuếch tán rối theo phương vng góc với chiều gio, đo: Công suất nguồn điểm phát thải liên tục va coi trình ổn định, nghĩa C 0 t  C (Kx )0 x x Nếu hướng trục Ox trùng với hướng gio thi thành phần vận tốc gio chiếu lên trục Oy  Vx  V  u  Vy  79 80 20 Sau phép đơn giản Điều kiện ban đầu  Tốc đô thẳng đứng thường nho so với tốc gio nên có thê bo qua, trục z thường lấy chiều dương hướng lên trên, đo bụi nặng thi thành phần Vz tốc rơi hạt (dấu âm), cịn chất ô nhiễm va bụi nhe thi Vz =  x    y   z  H u.C  M.  y. z  H  Berliand Vx C   C    C    Kz  Ky  x y  y  z  z  81 82 Giả thiết kèm Điều kiện biên n  z Vx u1 z  ;  1 Điều kiện xa vô x    z    y   Điều kiện phản xạ hoàn toàn Kz C z0 0 z C 0 m  z kz k1 z    1 k y  k 0V x  Trong đo u1 , k1 – vận tốc gio va sô rối đo đạc va chỉnh ly đô cao z1 = m; n va m tham sô không thư ngun chỉnh ly tính tốn tư sơ liệu đo đạc tầng không sát đất khu công nghiệp (thường thi người ta lấy xấp xỉ m ≈ 1, n ≈ 0,15, z1 = 1m), Điều kiện hấp thụ hoàn toàn C z0   k0 xác định sơ giải toán ngược khuếch tán rối (kết nhận cho thấy k0 0.1 – m phu thuộc vào mức ổn định tầng kết) 83 84 21 Nồng độ mặt đất Nghiệm giải tích xấp xỉ u C   C    C    Ky    Kz  x y  y  z  z  n  z Vx uu1  ; z1 Đối với độ cao z vài m với x > 10 – 20 m đặt z = 0; Giới hạn trường hợp m ≈1 m  z Kz k1  K y  k u z1 δ y δ z - H u.C x=0,y=0,z=H = M     Kz C z z=0 Cx, y,0  =0   mn u z mn Mz12nm exp  1 z2nm  H 2nm  4yk x  k  m  n x     C  x, y, z   mn 1n 1m 2nm 2nm 1n k0 x   n  m   2nm  k1 x   2u1  2nm  u H1n y2  exp   1 n2 k x 4k0x  21 nk1 k0 x3   M 85 Các biểu thức mơ hình Nồng độ cực đại mặt đất Cx, y,0  1n M 21 nk1 k0 x 32 Cx, y,0   uH y  exp   1 n2 k x 4k0x    Cm  u1H 1.51n  u H1n y2  exp   1 n2 k x 4k0x  21 nk1 k0 x3   M n  z Uxz=u1  ; z1 C C  0 x y 0.1161  n2 M 86 ky  z kz =k1  , k0 = U z1  M – Lượng thải (công suất nguồn thải) (mg/s); k1 k 0u1  k1 – sô khuếch tán rối đứng đô cao z1 = 1m (m2/s); u1 H 1 n xm  k1 1  n 2  n – sô mũ hàm biến thiên tốc đô gio n = 0.14 – 0.23 Thường lấy n=0.14 điều kiện bất ổn định va n = 0.2 khi ổn định;  k0 – Kích thước khuếch tán rối ngang (m): k =ky/U;  u1 – Tốc đô gio đô cao z1 = m 87 88 22 Xác định hệ số khuếch tán rối đứng Xác định hệ số khuếch tán rối đứng (tt) Nếu khơng có số liệu quan trắc nhiệt độ độ cao m 0.5 m ta sử dụng Profile thẳng đứng có qui luật luỹ thừa   T  k1  0,104.U 1  1,38   U     z K z  K1    z1  m T  Tk  Tk 850 1000 Nếu khơng có giá trị nhiệt độ quan trắc độ cao 1500m chấp nhận U = U2 - U0,5 - hiệu tốc độ gió mực 0,5 mét T = T0,5 - T2 - hiệu nhiệt độ khơng khí mực 0,5 mét T  1,5 * 0,63.10 2  0,945.10 2 89 Xác định hệ số khuếch tán rối ngang K0 Xác định hệ số khuếch tán rối đứng (tt) Chỉ số n giá trị Z0 theo tháng năm Hà Nội Độ dày lớp biên Tháng 10 Năm n 0.19 0.19 0.17 0.23 0.20 Z0 (mét) 0.0085 0.016 0.0083 0.017 0.016 K h  0.05 ln Z  ln Z ln10  ln z0 Z  Vx  z   V10    10  K 12 z12  z h  0.05 K1 z1 z z 7,29.105 grad1 Vận tốc gió ngồi lớp biên Tốc độ gió độ cao 2m, m 0.5 m tính theo cơng thức: Vx  z   V10 90 n Vh  V10 hệ số K0 ln h  ln Z ln 10  ln Z h Vh  V10    10  n K  K h Vh 91 92 23 Đối với tạp chất nặng Cx, y,0  MH  1 n u1 1 2 21  n  1 1    k0 x k1 x  Nồng độ cực đại mặt đất  u H 1 n y   exp    1  n  k x 4k x    Cm  w 1.3.102.prp2  w k1 1  n     w - tốc rơi hạt có dạng hình cầu, đo ρp mật hạt bụi, rp - bán kính chúng Trong cơng thức w xác định cm/s, ρp va rp cho g/cm3 va m tương ứng   0.0631  n 2 M u1H 1.51n  k1 1.5   1.5 u1 H 1 n , xm   k 0u1 1   e (1  n ) (1   ) k1 M – Lượng thải (công suất nguồn thải) (mg/s); k1 – sô khuếch tán rối đô cao z1 = 1m (m2/s); n – sô mũ hàm biến thiên tốc đô gio n = 0.14 – 0.2 Thường lấy n=0.14 điều kiện bất ổn định va n = 0.2 khi ổn định (GS Lê Đình Quang tính cho Hà Nội); k0 – Kích thước khuếch tán rối ngang (m): k =ky/U; u1 – Tốc đô gio đô cao z1= m 93 Trường hợp lặng gió 94 Độ cao hữu dụng theo Berliand H  M Cx, y,0  2k1 1  n   H 1n   x2  y2    1  n k1  Cm  MK  1  n  H 1,5 W R V10  3,3 g R  T  ,5    T k  273 ,1.V10      T Tr Tk  V10 – Vận tốc gio đô cao 10 m;  ΔT – Hiệu nhiệt đô thải khỏi miệng ống khói (Tr , 0C) va nhiệt mơi trường không xung quanh (Tk , 0C);  W0 – Tốc đô luồng thải (m/s); 1  n   R0 – Bán kính miệng ống khói (m);  g – Gia tốc trọng trường (m/s2);   2k1 95 96 24 Quá trình tự động hóa H  1,5.W0 R0 V10  3,3.g R0 T  2,5  Tk  273,1.V102  h 1n  uH y  M  Cx, y,0  exp   21 nk1 k0 x3  1 n2 k1x 4k0x   Tính tốn nồng độ trung bình; Tính tốn trường nồng độ ảnh hưởng nhiều nguồn thải Hhh 97 Tóm tắt 98 Nồng độ trung bình  Tính tốn nồng độ trung bình;  Tính tốn trường nồng độ ảnh hưởng nhiều nguồn thải; Trung bình: Ngày đêm, tháng, năm 99 100 25 Nồng độ trung bình ngày Các số liệu khí tượng cho vào thời điểm ngày: 1g, 7g, 13g, 19g Kịch tính tốn Khí tượng (thay đổi ngày lần): - Vận tốc gió - Hướng gió - Nhiệt độ C ( x , y , ),i   C  ,i , ( x , y , )  1 Nguồn thải (trung bình ngày) - Nhiệt độ khí - Vận tốc khí - Tải lượng chất nhiễm 101 Hệ số trung bình tần suất gió 102 Tính hệ số trung bình k gió - Thời gian có gio tất hướng; - Tần suất lặng gio: Plặng = lặng - Thời gian lặng gio; o = gió+ lặng - Tổng thời gian quan trắc (ngày đêm, tháng năm);  lang 0 - Tần suất gio hướng  : P  Plang   - thời gian có gio hướng ; m – sơ hướng gio, thông thường m = 16 , m P   1   0 Plang  m P   1 1 - Hê sô trung bình kα hướng gio  trường hợp k  P 103 104 26 Công thức tính nồng độ trung bình tháng C xy ( i )  Plang C lang C xy ( tong )  n  [P i 1 lang (i)  C lang m   1 (i)  Cxy(tổng) - nồng tổng cộng trung bình n nguồn thải gây điểm tính tốn; k  C  (i) m   1 Các ký hiệu Clặng(i) -nồng đô tức thời nguồn thải thư i gây điểm tính tốn lặng gio (u=0); k  C  (i) ] C(i) - nồng đô tức thời nguồn thải thứ i gây điểm tính tốn có gio thổi theo hướng  ứng với vận tốc gio trung bình hướng đo va ổn định trung bình suốt khoảng thời gian tính tốn trị sơ trung bình (ngày đêm, tháng năm) Cxy(i) – nồng trung bình vị trí có tọa x,y nguồn thư i gây ra; 105 106 27