MÔ HÌNH TÍNH TOÁN Ô NHIỄM KHÔNG KHÍ 1. Phương trình cơ bản để tính nồng độ chất ô nhiễm trong khí quyển: Khi mô tả quá trình khuyếch tán chất ô nhiễm trong không khí bằng mô hình toán học thì mức độ ô nhiễm không khí thường được đặc trưng bằng trị số nồng độ chất ô nhiễm phân bố trong không gian và biến đổi theo thời gian. Trong trường hợp tổng quát, trị số trung bình của nồng độ ô nhiễm trong không khí phân bố theo thời gian và không gian được mô tả từ phương trình chuyển tải vật chất (hay là phương trình truyền nhiệt) và biến đổi hoá học đầy đủ như sau: (1) Trong đó: C: Nồng độ chất ô nhiễm trong 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 của hệ số khuyếch tán rối theo các 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 của các chất ô nhiễm : Hệ số tính đến sự liên kết của chất ô nhiễm với các phần tử khác của môi trường không khí. : Hệ số tính đến sự biến đổi chất ô nhiễm thành các chất khác do những quá trình phản ứng hoá học xảy ra trên đường lan truyền. Tuy nhiên Phương trình (1) trên rất phức tạp và nó chỉ là một hình thức mô phỏng sự lan truyền ô nhiễm. Trên thực tế để giải phương trình này người ta phải tiến hành đơn giản hoá trên cơ sở thừa nhận 1 số điều kiện gần đúng bằng cách đưa ra các 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ì thành phần tốc độ gió chiếu lên trục Oy sẽ bằng 0, có nghĩa là v = 0. Tốc độ gió thẳng đứng thường nhỏ hơn rất nhiều so với tốc độ gió nên có thể bỏ qua, có nghĩa là w = 0. Trong nhiều trường hợp, nếu xét bụi nhẹ thì Ws = 0 (trong trường hợp bụi nặng thì lúc đó ta sẽ cho Ws 0). Nếu bỏ qua hiện tượng chuyển pha (biến đổi hoá học) của chất ô nhiễm cũng như không xét đến chất ô nhiễm được bổ sung trong quá trình khuyếch tán thì . Như vậy sau các giả thiết và chấp nhận 1 số điều kiện gần đúng thì phương trình ban đầu được viết dưới dạng là: (2) Nếu giả sử rằng các hệ số là không đổi thì pt (2) được viết lại là : (3) Trong trường hợp không tính đến thành phần phi tuyến thì phương trình (3) được viết là: (4) Ta thấy phương trình (4) là dạng phương trình truyền nhiệt 2 chiều. Tuỳ theo điều kiện ban đầu và điều kiện biên mà ta có các nghiệm giải tích khác nhau. Để tìm nghiệm giải tích phương trình (4), đầu tiên xét bài toán truyền nhiệt 1 chiều có dạng sau: , t = 0
1 MÔ HÌNH HÓA MÔI TRƯỜNG MÔ HÌNH TÍNH TOÁN Ô NHIỄM KHÔNG KHÍ Phương trình để tính nồng độ chất ô nhiễm khí quyển: Khi mô tả trình khuyếch tán chất ô nhiễm không khí mô hình toá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 = kx ÷+ k z ÷+ k y ÷+ α C − β C + wc ∂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 Phương trình (1) 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 1 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 hoá 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 = ∂ k ∂C ÷+ ∂ k ∂C ÷ (2) ∂t ∂x ∂y y ∂y ÷ ∂z z ∂z ÷ Nếu giả sử hệ số k y , kz không đổi pt (2) viết lại : ∂C + u ∂C = k ∂2C + k ∂2C z y ∂t ∂x ∂y ∂z (3) Trong trường hợp không tính đến thành phần phi tuyến phương trình (3) viết là: u ∂C ∂x ∂C = k ∂2C + k ∂2C z y ∂t ∂y ∂z (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 toán truyền nhiệt chiều có dạng sau: ∂u = a2 ∂ 2u −∞ < x < +∞ ∂t ∂x2 (5) Với điều kiện ban đầu : u ( x, t ) = ϕ ( x) ,t=0 −∞ < 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 hay X a 2T Từ suy : X ''+ λ X = T '+ a 2λ 2T = Nghiệm phương trình (7) X = C eiλ x (6) (7) (8) 1 X = C2e−iλ x (Xem cách giải phương trình trang 53 [7]) 2 Nghiệm phương trình (8) T = C3e−a λ t 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 (9) 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λ x d λ (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) Ngoài ta 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 λζ u( x,t ) = ϕ (ζ )e dζ ÷÷e−a2λ 2t +iλ xd λ ∫ ∫ 2π −∞ −∞ ÷ +∞ +∞ 2 = ∫ ∫ e−a λ t +iλ ( x−ζ )d λ ÷÷ϕ (ζ )dζ 2π −∞ −∞ 2 e Xét tích phân I = e−a λ t +iλ ( x−ζ )d λ = 2π π a2t Như u( x, t ) = +∞ −( x−ζ )2 4a2t −( x−ζ )2 e 4a2t ϕ (ζ )dζ ∫ −∞ π a 2t (13) e Đặt G( x,ζ , t ) = π a2t Ta có u( x, t ) = −( 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: − ( x − ζ ) x − ζ Gx = − e π 2(a2t )3/2 4a2t )2 e −( x − ζ )2 Gxx = − 21 3/2 + ( x −2ζ 5/2 π (a t ) 4a2t 4(a t ) ( x − ζ )2 −( x − ζ )2 a a e Gt = − 3/2 + 5/2 π 2(a t ) 4a2t 4(a t ) Vậy Gt = a 2Gxx 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 Đặt a2 = kx (15) nghiệm phương trình (15) viết lại là: − x2 Q C ( x, t ) = e 4tkx 1/2 π tk x (16) Đây nghiệm cảu toán lan truyền ô nhiễm chiều với nguồn thải Q Cùng với điều kiện biên x → ∞ C → (Nồng độ ô nhiễm điểm giảm điểm tiến xa khỏi châ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 4(π t)(k xk y )1/2 2 ÷ − x + y ÷÷ 4t kx k y ÷÷ (17) Đối với toán chiều ta có: C ( x, y, z, t ) = Q e 8(π t )3/2 (k xk y kz )1/2 2 ÷ − x + y + z ÷÷ 4t kx k y kz ÷ (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 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 2.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 đó, toán lan truyền chất ô nhiễm toán hai chiều ta chọn công thức (17) để áp dụng cho trường hợp này: c) d) a) u 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= M e 4π ut (k y k z )1/2 Đặt : 2 ÷ − y + z ÷÷ 4t k y kz ÷ (19) k y = 0.5σ y u x kz = 0.5σ z u x t=x u (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σ z ÷÷ = M e 2π uσ yσ z ÷ − y ÷÷ − z2 ÷÷ 2σ y2 ÷÷ 2σ z2 ÷÷ 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ở 2.2 Diễn giải công thức mô hình Gauss sở phương pháp phân tích thứ nguyên: Công thức (23) 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 phản ứng hoá 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 Từ ta viết C: M uyz (24) 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 (25) σ 2π 2σ Á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 e 2σ y e 2σ z σy (26) σ 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 3.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: C= M 2π uσ yσ z − y2 −( z − H )2 2 e 2σ y e 2σ z (27) 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 hoàn toà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), 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à: C= M 2π uσ yσ z − y2 −( z − H )2 2σ y 2σ z e e −( z + H )2 σ z +e (29) 10 Đây công thức tính toán khuyếch tán chất ô nhiễm từ nguồn điểm cao liên tục Khi tính toá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= M π uσ yσ z −( H )2 e 2σ z (31) 2.4 Hệ số khuyếch tán: Từ biểu thức (20) (21) ta có: 2k x σ y = y ÷÷ u ÷ 2 σ z = 2kz x ÷÷ u (32) 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 xuôi 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 toá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 xuôi theo chiều gió kể từ nguồn, tính km Các hệ số a, b, c, d cho bảng Bảng 2: Các hệ số a, b,c, d công thức (33) Cấp ổn x >1 km x ≤ km b c d b c d định 10 11 A B C D E F 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 22.8 0.678 -1.3 55.4 0.305 -34.0 14.35 0.740 -0.35 62.6 0.180 -48.6 34 2.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 - Chiều cao hiệu ống khói: 11 12 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 xuôi 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 trình hoà 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 ∆h = D ω ÷÷ u 1+ ∆T Tkhoi ÷, m ÷ ÷ (34) 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 12 13 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 - Sự lắng đọng bụi trình khuyếch tán khí thải từ nguồn điểm cao: Những công thức tính toá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 xuôi 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ụ hoàn toàn Ngoài ra, chất ô nhiễm thể khí không chịu ảnh hưởng lực trọng trường, 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 bụi rơi khoảng thời gian t Đoạn 13 14 đườ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: Cb = Mb e 2π uσ yσ z − y2 −( z − H + vur x )2 2σ y 2σ 2z (36) e 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 độ) = Diện tích (Diện tích) (Vận tốc) (Nồng độ) = vrCb(x) Gb(x) = = Gb( x) = • (38) Công thức tính vận tốc bụi: ρb gδ vr = 18µ Với δ µ • M bvr 2π uσ yσ z −( H − vur x )2 e 2σ z - (39) ρ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 µ : 14 15 3/2 µtoC = µ0o C 387 273 + t ÷ 387 + t 273 ÷ (40) với µ0 C = 17,17.10 Pa.s II Tính toá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 : −6 o 2C ∂ ∂C ∂ C ∂ C ∂ ÷+ α C (41) u + w = k y ÷÷+ kz ∂x ∂z ∂y ÷ ∂z ∂z ÷ 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ạ hoàn toà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 15 u H1+n − − y k1(1+n)2 x 4k0 x (42) 16 Nồng độ cực đại mặt đất: 2M k1 0.116(1 + n ) Cmax = k0u1 u1H 1.5(1+n) (43) Khoảng cách từ XM từ nguồn đến vị trí có nồng độ max: u1H 1+n X max = k1(1+ n)2 (44) 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 • Công thức tính vận tốc độ cao z: Z ÷ 10 uz = u10 (45) Với : uz - Vận tốc gió độ cao z u10 - Vận tốc gió độ cao 10m V.2 Đối với bụi nặng cỡ hạt đồng chất: Nồng độ bụI mặt đất: MH λ (1+n)u λ1 C( x, y, z =0) = e + λ + λ 2(1+ n) Γ(1+ λ )(k1x) π k0 x u H1+n − − y (1+n)2 k1x 4k0 x vr Trong λ = (1 + n)k (47) Nồng độ bụi cực đại mặt đất: Cmax = 0.063(1 + n) M u1 H 1.5(1+ n ) k1 (1.5 + λ )1.5+ λ k0u1 eλ Γ(1 + λ ) Khoảng cách từ nguồn đến vị trí có nồng độ max: 16 (48) (46) 17 u1H 1+n xM = (1+ n)2 (1.5 + λ )k1 (49) 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 toá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ụ: ∂ ( Rk ∂C ) + ∂ k ∂C + M δ (r )δ ( z − H ) = r R ∂r ∂r ∂z z ∂z (50) Với điều kiện biên sau đây: ∂C =0 ∂r - Khi z = 0: - Sự phân bố nồng độ ô nhiễm mặt nằm ngang có kz tính đối xứng qua tâm nguồn R2 + z2 → ∞ : C → R=0→ ∂C =0 ∂r Berlaind Kurebin giảI phương trình (50) vói nghiệm có dạng: C (R, z) = a( H 1+ n + z1+ n ) + R M 2π k1(1+ n) a( H (1+ n)/2 − z(1+ n)/2 )2 + R 3/2 (1+ n)/2 + z(1+ n)/2 )2 + R a( H (51) a= Trong 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 ) Cmax = M (1+ n)3 32π k13H 2(1+n) (54) 17 3/2 18 III Tính toán nồng độ trung bình chất ô nhiễm mặt đất nguồn thải gây ra: Khi tính toá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 toá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 toán với cấp ổn định Ngoà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 α Ngoài tần suất hướng gió khác có tần suất lặng gió P lặng Đó tỷ lệ thời gian 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) tần suất gió hướng α Plang tần suất lặng gió m - Số hướng gió thông thường Trong với Plang < Pα m ∑ Pα = α =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 toán nguồn thải thứ i gây sau : m Cxy (i) = Plang Clang (i) + ∑ kα Cα (i) (53) α =1 18 19 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 toá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 toá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 toán trị số trung bình (ngày đêm, tháng năm) IV 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… 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 toán nồng độ ô nhiễm yếu tố độc hại cách riêng rẽ mức độ ô 19 20 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: C0 = C1 C + + + Cn Ccf (1) Ccf (2) Ccf (n) (55) Trong đó: C1, C2,…, Cn nồng độ thực tế đo tính toá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) 20 21 Đạ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 hoá 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 21 22 Câu 8: Xác định nồng độ oxy bão hòa hồ nước mặn có độ cao so với mặt nước biển km, độ mặn 10 ppt nhiệt độ nước hồ 25oC? Lượng Ôxy bão hòa nước theo nhiệt độ độ mặn atm T DO (mg/L) ppm độ muối (0C) ppm độ muối 10 11,28 10,92 11 11,02 10,67 12 10,77 10,43 13 10,53 10,20 14 10,29 9,98 15 10,07 9,77 16 9,86 9,56 17 9,65 9,36 18 9,45 9,17 21 8,90 8,64 22 8,73 8,46 23 8,56 8,32 24 8,40 8,16 25 8,24 8,01 26 8,09 7,87 27 7,95 7,73 28 7,81 7,59 29 7,67 7,46 22 23 Bảng Nồng độ Oxy hòa tan nước (mg/l) phụ thuộc vào nhiệt độ, độ muối điều kiện không khí ẩm, áp suất 760 mmHg Nhiệt Độ muối (%0, g\l) độ 0C 10 15 20 25 30 15 10.07 9.77 9.47 9.19 8.91 8.64 8.38 16 9.86 9.56 9.28 9.00 8.73 8.47 8.21 17 9.65 9.36 9.09 8.82 8.55 9.30 8.05 18 9.45 9.17 8.90 8.64 8.38 8.14 7.90 19 9.26 8.99 8.73 8.47 8.22 7.98 7.75 20 9.08 8.81 8.56 8.31 8.06 7.83 7.60 21 8.90 8.64 8.39 8.15 7.91 7.68 7.46 22 8.73 8.48 8.23 8.00 7.77 7.54 7.33 23 8.56 8.32 8.08 7.85 7.63 7.41 7.20 24 8.40 8.16 7.93 7.71 7.49 7.28 7.07 25 8.24 8.01 7.79 7.57 7.36 7.15 6.95 26 8.09 7.87 7.65 7.44 7.23 7.03 6.83 27 7.95 7.73 7.51 7.31 7.10 6.91 6.72 28 7.81 7.59 7.38 7.18 6.98 6.79 6.61 29 7.67 7.46 7.26 7.06 6.87 6.68 6.50 30 7.54 7.33 7.14 6.94 6.75 6.57 6.39 (Em dựa vào bảng dựng đường chuẩn tính nhé) 23 35 8.13 7.97 7.81 7.66 7.52 7.38 7.25 7.12 6.99 6.87 6.75 6.64 6.53 6.42 6.32 6.22 24 Câu 13: Dòng thải từ trạm xử lý nước thải áp dụng công nghệ bùn hoạt tính có lưu lượng m3/s; BOD5 = 10 mg/L; hệ số f = 1,4, thải vào dòng sông có lưu lượng m3/s; U = 0,3 m/s BOD nước sông = Dòng sông có Ka (20oC) = 0,2/ngày Nhiệt độ nước sông T = 28oC a) Xác định nồng độ BOD vị trí điểm xả (điểm xáo trộn nước thải nước sông)? Qn = m3/s Ln = 10 mg/L f = 1,4 Qs = m3/s Us= 0,3 m/s Ls = mg/L Ka(20oC) = 0,2 / ngày Ts = 28oC Áp dụng công thức: La = (QnLn + QsLs)/(Qn+Qs) = (2*10+5*0)/(2+5) = 2,87 (mg/L) b) Xác định khoảng cách mà nồng độ BOD nước sông giảm xuống 5% nồng độ BOD ban đầu? Tính Ka(28oC) = Ka(20) (θ)T-20 Trong θ = 1,056 Nồng độ BOD vị cách điểm xả thải x (m) 5% lượng BOD điểm xáo trộn: BOD(x) = 5% La = 5%*2,87 = 0,1435 mg/L Mà BOD(x) = La*e-Ka*x/0,3 = 0,1435 Từ thay La, Ka(280C) vào để tìm x (m) 24 [...]... giá và so sánh chất lượng của môi trường giữa nơi này và nơi khác, giữa khu công nghiệp này và khu công nghiệp khác Cụ thể là nếu có hai khu vực khác nhau, ở khu vực 1 môi trường chỉ bị ô nhiễm bởi một chất duy nhất với nồng độ bằng nồng độ giới hạn cho phép, trong lúc đó ở khu vực 2 môi trường bị ô nhiễm đồng thời bởi nhiều chất độc hại khác nhau mà nồng độ mỗi chất đều đạt giới hạn cho phép thì theo... khác nhau : khí SO2, CO, CO2, NOx, H2Sv.v…và bụi Khi khí thải lan toả ra môi trường xung quanh các yếu tố độc hại nêu trên cùng song song tồn tại trong không khí và gây tác hại một cách tổng hợp đối với môi trường sống của con người cũng như của sinh vật nói chung Tuy nhiên, từ trước đến nay khi xem xét đánh giá chất lượng môi trường không khí của một khu vực đồ thị, nhà máy hay khu công nghiệp nào... khu vực nói trên đều được xem là đạt yêu cầu và đều được chấp nhận như nhau, trong khi trên thực tế, môi trường ở khu vực 2 chắc chắn là gây ảnh hưởng xấu đến sức khoẻ của con người nhiều hơn khu vực 1 Để kể đến tác dụng tổng hợp đồng thời của 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 của nhiều chất ô nhiễm cùng tồn tại đồng thời trong môi trường không khí... đó, người ta thường khảo sát đo đạc hoặc tính toán nồng độ ô nhiễm của các yếu tố độc hại một cách riêng rẽ và mức độ ô 19 20 nhiễm môi trường không khí cũng được đánh giá một cách độc lập đối với từng chất ô nhiễm khác nhau Mặc khác hiện nay trong các tiêu chuẩn về môi trường của Việt Nam cũng mới chỉ đưa ra giới hạn nồng độ cho phép của từng chất ô nhiễm riêng biệt trong không khí xung quanh hoặc... giới hạn nồng độ cho phép của từng chất 1, 2,…,n riêng biệt (xem phụ lục 3) 20 21 Đại lượng C0 có thể được xem như một chỉ số tổng hợp đánh giá chất lượng môi trường không khí về mặt ô nhiễm hoá học Nó cho phép ta đánh giá, so sánh chất lượng môi trường không khí theo một chỉ tiêu thống nhất 21 22 Câu 8: Xác định nồng độ oxy bão hòa của một hồ nước mặn có độ cao so với mặt nước biển là 1 km, độ mặn... mặt đất Dựa vào mô hình Gauss cơ sở - tức là mô hình Gauss chưa kể đến sự phản xạ của mặt đất đối với chất ô nhiễm được thể hiện bằng công thức (27) Đối với đa số chất ô nhiễm thể khí thì mặt đất không hấp thu mà phản xạ ngược trở lại vào khí quyển Còn đối với bụi ta có thể xem mặt đất là vật hấp thụ hoàn toàn Ngoài ra, chất ô nhiễm thể khí hầu như không chịu ảnh hưởng của lực trọng trường, còn bụi... theo chiều cao ta có các trường hợp đẳng nhiệt, đoạn nhiệt, siêu đoạn nhiệt hoặc 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 0 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 0 -10 -8 -6 -4 -2 0 2 4 6 8 Độ 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... ) (xem phụ lục 2) V.3 Khuyếch tán chất ô nhiễm từ nguồn điểm cao trong điều kiện không gió: Các phương pháp tính toán khuyếch tán chất ô nhiễm nêu ra trên đây đều áp dụng cho trường hợp trời có gió Để kể đến ảnh hưởng của khuyếch tán rối trong trường hợp trời không gió, Berlaind và Kurebin đã đưa ra phương trình sau đây trong hệ toạ độ trụ: 1 ∂ ( Rk ∂C ) + ∂ k ∂C + M δ (r )δ ( z − H ) = 0 r R ∂r ∂r... việc xác định nồng độ ô nhiễm tức thời, ta cần phải dự báo được sự phân bố nồng độ trung bình ngày đêm, trung bình tháng hoặc trung bình năm Trong trường hợp tính toán nồng độ trung bình cho thời gian ngắn, như trung bình ngày đêm ta có thể đơn giản hoá vấn đề bằng cách giả thiết rằng trong từng mùa nhất định, hè hoặc đông, cấp ổn định của khí quyển có thể thay đổI trong ngày đêm xung quanh cấp ổn định... nhiễm thể khí hầu như không chịu ảnh hưởng của lực trọng trường, còn bụi thì rơi trong khí quyển với vận tốc vr nhất định tuỳ thuộc vào kích thước hạt và khối lượng đơn vị của nó Do đó đại lượng H trong mô hình Gauss cần được hiệu chỉnh bằng cách trừ bớt đi đoạn đường mà hạt bụi rơi được trong khoảng thời gian t Đoạn 13 14 đường đó là vr t mà t = x/u vớI x là khoảng cách dọc theo trục gió tính từ nguồn