TÌM HIỂU PHẦN MỀM ENVIMAP - TÍNH TOÁN MÔ PHỎNG Ô NHIỄM KHÔNG KHÍ CHO NHIỀU NGUỒN ĐIỂM
TRƯỜNG ĐẠI HỌC THỦ DẦU MỘT KHOA TÀI NGUYÊN MÔI TRƯỜNG - HỌC PHẦN: MÔ HÌNH HÓA ỨNG DỤNG ĐỀ TÀI: TÌM HIỂU PHẦN MỀM ENVIMAP - TÍNH TOÁN MÔ PHỎNG Ô NHIỄM KHÔNG KHÍ CHO NHIỀU NGUỒN ĐIỂM GVHD: ThS Nguyễn Thị Khánh Tuyền SVTH: Nhóm 6 Nguyễn Thị Thanh Nguyễn Thị Cẩm Tiên Nguyễn Hữu Duy Mai Thanh Điền Nguyễn Ngọc Hữu Nguyễn Ngọc Sơn Mai Thế Tâm 1220510147 1220510157 1220510201 1220510195 1220510204 1220510223 1220510225 Năm học 2015-2016 MỤC LỤC DANH MỤC BẢNG DANH MỤC HÌNH I ĐẶT VẤN ĐỀ Hiện nay, tình hình ô nhiễm môi trường, đặc biệt ô nhiễm không khí vấn đề cấp bách toàn giới khó kiểm soát chất thải không khí Mà nước ta đà phát triển theo hướng công nghiệp hóa, nên lượng khí thải từ khu công nghiệp, khu sản xuất thải ngày nhiều Vì vậy, cần phải có phương pháp tính toán mô hình để mô phỏng, kiểm soát ô nhiễm không khí Hơn nữa, biết mô hình hóa toán học trình môi trường sản phẩm túy khoa học mà coi phương pháp tiếp cận để hiểu biết sâu sắc tượng thiên nhiên mục tiêu cuối để nhận thông tin giới thực Thông tin thúc đẩy phát triển vấn đề khoa học phương pháp giải chúng, làm sở để thông qua định tiến hành dự án cụ thể Trong năm gần mối quan tâm xây dựng mô hình toán ô nhiễm không khí, nước, đất, dự báo đánh giá khía cạnh kinh tế ô nhiễm dựa phương pháp mô tăng lên Việc xây dựng mô hình toán cho hệ thống kiểm soát quản lý ô nhiễm không khí, luận chứng phương pháp dự báo dài hạn phục vụ cho công tác qui hoạch không ngừng tăng lên Để ứng dụng mô hình giải vấn đề thực tiễn cần thiết xây dựng phần mềm chuyên dụng Tuy nhiên để sử dụng phần mềm cần thiết phải hiể biết rõ phần mềm để mô trạng ô nhiễm khu vực xét Vì vậy, nhóm chúng em tìm hiểu phần mềm ENVIMAP để giúp trình mô ô nhiễm nhanh xác II NỘI DUNG NGHIÊN CỨU 2.1 Tổng quan phần mềm ENVIMAP phiên 3.0 Phần mềm ENVIMAP phần mềm tích hợp GIS, CSDL môi trường mô hình toán học mẫu lan truyền chất ô nhiễm môi trường không khí môi trường nước Các CSDL khí tượng, điểm lấy mẫu chất lượng không khí, nguồn thải điểm, phát thải phải gắn với vị trí địa lí nhằm thể hiện trạng môi trường Bên cạnh để mô trình lan truyền chất ô nhiễm thường xuất cần thiết phải thực toán chuẩn như: lấy CSDL cần thiết cho mô từ liệu khí tượng từ ống khói (hay phát thải), biểu diễn kết tính toán mô dạng khác nhau, chồng lớp thông tin lớp kết tính toán mô với lớp không gian gắn với địa phương cụ thể Quá trình làm việc, ENVIMAP diễn phối hợp thông tin chặt chẽ với GIS Các mô hình nhận liệu từ GIS, kết làm việc thông báo ngược trở lại vào GIS hình thức thông tin đồ Như thấy diễn trình phân công chức tự nhiên: GIS đảm nhận việc biểu diễn thông tin, mô hình đảm nhận thực xử lý thuật toán xử lý thông tin Phiên ENVIMAP đời vào năm 2002 Tới tháng 11/2005 phần mềm ENVIMAP 1.0 nâng cấp thành phiên 2.0 Tháng 9/2006, phiên 3.0 đời Từ năm 2008 trở ENVIMAP đặt tên theo năm tác giả cập nhật thường xuyên Phần mềm ENVIMENVIMAP nâng cấp nhằm mục tiêu sau đây: - Quản lý nguồn thải cố định (cụ thể ống khói) - Cho phép tính toán ảnh hưởng nguồn thải lên tranh ô nhiễm chung - Thực báo cáo nguồn thải kết tính toán - Cho phép nhập lưu trữ liệu liên quan tới khí tượng - Cho phép nhập lưu trữ liệu liên quan tới chất lượng không khí xung quanh - Tích hợp văn pháp lý liên quan tới quản lý chất lượng không khí 2.2 Cấu trúc chức phần mềm ENVIMAP 2.2.1 Cấu trúc Hình 2.1 Sơ đồ cấu trúc phần mềm ENVIMAP ENVIMAP 2.0 gồm khối liên kết với - Khối CSDL môi trường (liên quan tới môi trường không khí) - Khối mô hình (trong phần mềm ENVIMAP mô hình Berliand) - Khối GIS – quản lý đối tượng cách trực diện đồ - Khối thực báo cáo thống kê - Khối hỗ trợ văn pháp qui Hình 2.2: Sơ đồ cấu trúc CSDL môi trường ENVIMAP MÔ HÌNH BERLIAND MÔ HÌNH MÔ HÌNH GAUSS MÔ HÌNH ISC3 Hình 2.3: Mô hình Berliand tích hợp ENVIMAP 2.2.2 Chức Tạo đối tượng cần quản lý Tạo Cơ sở sản xuất (CSSX): Đối tượng nguồn thải ENVIMAP quản lý thuộc CSSX ENVIMAP cho phép người dùng chức tạo CSSX đưa vào CSDL - CSSX nằm khu công nghiệp (KCN) CSSX thuộc cấp - CSSX nằm KCN Khi KCN đối tượng cấp CSSX thuộc cấp Các CSDL liên quan tới CSSX gồm: - Tên CSSX - Thuộc Ban ngành chức (ví dụ thuộc Sở Công nghiệp,…) - Thành phần kinh tế: (Cơ quan nhà nước, Công ty liên doanh, Tư nhân, …) - Thuộc Ngành công nghiệp - Mã ngành công nghiệp - Danh mục nguồn thải điểm thuộc CSSX - Tên giám đốc (điện thoại, fax, e-mail, web site) - Tên người phụ trách môi trường Các thông tin ENVIMAP sử dụng để xuất báo cáo cần thiết phục vụ cho mục tiêu quản lý Tạo ống khói đồ số: thông tin thuộc tính ống khói trình bàytrong bảng đây: Thông tin lưu trữ Id Mã ống khói Tên ống khói Kinh độ Vĩ độ X Y Chiều cao Đường kính Kiểu liệu Int Char Nvarchar Nvarchar Nvarchar float float float float Kích thước tối đa (byte) 50 53 53 53 8 8 Tạo vị trí lấy mẫu không khí cho mục tiêu quan trắc: Theo chương trình khác nhau, hàng năm vùng có tiến hành lấy mẫu phân tích chất lượng không khí vùng Các liệu có ý nghĩa pháp lý quan trọng nhằm đưa định quản lý hành Hiện nhiều tỉnh thành đất nước, số liệu nằm rải rác nhiều tài liệu khác khó khai thác Bởi cần thiết phải xây dựng công cụ tin học trợ giúp công tác quản lý tổng hợp thống số liệu quan trắc Phần mềm ENVIMAP 2.0 hỗ trợ chức tạo vị trí lấy mẫu đồ nhập kết quan trắc vào ENVIMAP Tạo điểm kiểm soát chất lượng không khí: Trong vùng cần quản lý, có số điểm nhạy cảm cần giám sát chặt chẽ chất lượng môi trường như: vị trí gần trường học, bệnh viện, khu dân cư, khu du lịch,… Phần mềm ENVIMAP cho phép người dùng tạo vị trí đồ Các vị trí đưa vào CSDL phục vụ cho mục tiêu quản lý Ví dụ tính toán theo mô hình Berliand lan truyền ô nhiễm, phần mềm ENVIMAP 2.0 xuất kết tính toán điểm kiểm soát chất lượng không khí Quản lý thông tin quan trắc Các số liệu quan trắc thay đổi theo thời gian quản lý phần mềm ENVIMAP 2.0 thông qua giao diện thân thiện xây dựng riêng Như biết thông tin môi trường thường xuyên thay đổi theo thời gian, để quản lý nhóm thông tin này, ENVIMAP sử dụng phần mềm quản trị CSDL mạnh dễ sử dụng MS Access Nhóm thông tin thay đổi theo thời gian quản lý ENVIMAP gồm: số liệu quan trắc chất lượng không khí, số liệu khí tượng, tải lượng ô nhiễm ống khói (cùng thông tin khác kèm nhiệt độ khí thoát ra, vận tốc khí phụt,…) Để nhập số liệu vào ENVIMAP, người sử dụng cần chọn vị trí quan trắc (trạm), ngày, tháng, năm lấy liệu Cấu trúc liệu chất lượng không khí xây dựng dựa thực tiễn quan trắc số tỉnh thành Việt Nam Phần trình bày số cấu trúc liệu thay đổi theo thời gian ENVIMAP quản lý Báo cáo thống kê ENVIMAP 2.0 cho phép thực dạng báo cáo dựa số liệu thay đổi theo thời gian sau đây: - Kết quan trắc chất lượng không khí điểm quan trắc (một hay nhiều trạm người dùng tùy chọn) - Kết tính toán nồng độ chất ô nhiễm điểm kiểm soát theo mô hình toán Berliand; Kết tính toán nồng độ cực đại chất ô nhiễm không khí tính nhóm ống khói (tính theo mô hình Berliand) ống khói ENVIMAP hỗ trợ người dùng lựa chọn thêm số chức như: thể dạng đường đồng mức, có lưới, đường đồng mức tô đặc hay không tô đặc, đường đồ thị Cho phép in % phát thải đóng góp nguồn thải ô nhiễm không khí - Cho phép đưa Text lẫn đồ thị vào Báo cáo tổng hợp - Lưu kết thành file theo định dạng doc, excel, pdf,… Cho phép gửi báo cáo qua Email Tính toán phát tán ô nhiễm theo mô hình Nghiên cứu ứng dụng mô hình phát tán ô nhiễm không khí thực tế chủ đề riêng rẽ Nhóm tác giả thực đề tài có nhiều năm nghiên cứu mô hình không khí ứng dụng cho mức độ địa phương Mô hình cài đặt ENVIMAP mô hình Berliand, mà theo phân loại Tổ chức khí tượng giới thuộc nhóm Mô hình thống kê thủy động, lý thuyết nửa thứ nguyên (còn gọi mô hình K) Mô hình Berliand (Nga) hoàn thiện áp dụng Liên Xô Ở Việt Nam, KS Nguyễn Cung người áp dụng mô hình cho số công trình, dự án phía Bắc Trong thời gian qua nhóm nghiên cứu GS Lê Đình Quang GS Phạm Ngọc Hồ ứng dụng mô hình cho nhiều tỉnh phía Bắc Hà Nội, Hải Phòng,… Giả thiết vận tốc gió theo hướng gió V x hệ số khuếch tán rối theo phương đứng kz phương ngang ky cho dạng sau: Trong đó: U1, k1: Là vận tốc gió hệ số rối đo đạc chỉnh lý độ cao z1 = 1mét n: Là tham số không thứ nguyên chỉnh lý tính toán từ số liệu đo đạc tầng không khí sát đất khu công nghiệp (thường người ta lấy xấp xỉ n ≈ 0.15, z = 1m k0 xác định sở giải toán ngược khuếch tán rối Đối với nồng độ chất khí hay bụi nhẹ mặt đất, xuất phát từ giả thiết xem xét lời giải toán biên lan truyền, khuếch tán chất ô nhiễm, Berliand đưa công thức sau: Trong đó: M: Là công suất nguồn thải H: Là độ cao hữu dụng (nghĩa chiều cao vật lýcủa ống khói cộng với vệt nâng ống khói) Từ công thức Berliand nhận hai đặc trưng quan trọng nồng độ cực đại khoảng cách (cách ống khói theo hướng gió) đạt giá trị cực đại đó: 2.3 Nguyên tắc mô Để mô tả trình lan truyền khuếch tán chất ô nhiễm không khí theo không gian thời gian phương trình toán học người ta xem xét trị số trung bình nồng độ chất ô nhiễm 10 4.2.6 Xây dựng kịch Berliand 4.2.6.1 Xây dựng kịch tính nồng độ bụi nhẹ cho ống khói - Vào menu kịch -> kịch Berliand -> Xuất hộp thoại -> Thêm -> nhập thông số (hình 4.13) Hình 4.13: Hộp thoại kịch Berliand – Trang thông tin - Sau chọn vận tốc – tần suất gió (Mô hình tự cập nhật thông tin khí tượng) Hình 4.14: Kịch Berliand – Trang Vận tốc – Tần suất gió 29 - Các thông số kịch ta để mặc định Hình 4.15: Hộp thoại Kịch Berliand – Trang Thông số kịch Berliand Hình 4.16: Hộp thoại Kịch Berliand – Trang Số liệu phát thải kịch - Chọn lưu thoát khỏi hộp thoại 4.2.6.2 Xây dựng kịch tính phát tán ô nhiễm cho nguồn OK1, OK2, OK3, OK4, OK5 OK6 theo trung bình ngày điểm nhạy cảm 30 4.2.6.2.1 Tạo điểm nhạy cảm - Vào menu thông tin -> chọn Điểm nhạy cảm -> Xuất hộp thoại điểm nhạy cảm -> chọn thêm -> Xuất hộp thoại -> nhập thông tin điểm nhạy cảm -> Chọn lưu Hình 4.17: Hộp thoại điểm nhạy cảm Hình 4.18: Thông tin điểm nhạy cảm 4.2.6.2.2 Tạo kịch - Tạo kịch cho bụi nhẹ, SO2, CO, NO2 - Vào menu kịch -> kịch Berliand -> Xuất hộp thoại -> Thêm -> nhập thông số (hình 4.19) 31 Hình 4.19: Hộp thoại kịch Berliand – Trang thông tin - Sau chọn vận tốc – tần suất gió (Mô hình tự cập nhật thông tin khí tượng) Hình 4.20: Kịch Berliand – Trang Vận tốc – Tần suất gió - Các thông số kịch ta để mặc định 32 Hình 4.21: Hộp thoại Kịch Berliand – Trang Thông số kịch Berliand Hình 4.22: Hộp thoại Kịch Berliand – Trang Số liệu phát thải kịch - Chọn lưu thoát khỏi hộp thoại - Làm tương tự cho chất SO2, CO, NO2 33 Hình 4.23: Thông tin kịch Berliand cho chất 4.2.7 Chạy kịch - Click vào menu Mô hình -> Chọn loại kịch -> chọn kịch -> Click chạy mô hình -> Xuất hộp thoại (Hình 4.24) Hình 4.24: Chọn kịch chạy mô hình - Click tiếp -> Xuất hộp thoại (hình 4.25) 34 Hình 4.25: Nhập thông số cho lưới tính - Click Hoàn tất để chạy mô hình 4.3 Kết chạy mộ hình - Sau chạy mô hình ta thu số kết sau: 35 4.3.1 Kịch câu a Hình 4.26: Bản đồ mô – kịch câu a Hình 4.27: Kết chạy mô hình – Trang nồng độ cực đại - Theo yêu cầu toán tính nồng độ cực đại bụi nhẹ ống khói doanh nghiệp gây vào thời điểm 13 trưa ngày 20 tháng Dựa vào kết chạy mô hình (hình 4.28) ta thấy nồng độ cực đại bụi nhẹ OK1; OK2; OK3; OK4; OK5;OK6 36 0.530542; 0.455872; 0.426508; 0.533135; 0.385427; 0.42397 mg/m theo TCVN 5937 – 2005 nồng độ cực đại bụi nhẹ ống khói vượt tiêu chuẩn 4.3.2 Kịch câu b 4.3.2.1 Kịch bụi nhẹ Hình 4.28: Bản đồ mô kịch câu b – Bụi nhẹ 37 Hình 4.29: Kết chạy mô hình – Trang nồng độ bụi nhẹ điểm nhạy cảm - Dựa vào kết chạy mô hình, nồng độ điểm nhạy cảm 0.219487 0.291238 mg/m3 theo TCVN 5937 – 2005 – Trung bình ngày bụi nhẹ 0,2mg/m3 -> Nồng độ bụi nhẹ điểm nhạy cảm vượt tiêu chuẩn 4.3.2.2 Kịch SO2 Hình 4.30: Kịch mô câu b – SO2 38 Hình 4.31: Kết chạy mô hình – Trang nồng độ SO2 điểm nhạy cảm - Dựa vào kết chạy mô hình, nồng độ điểm nhạy cảm 0.305869 0.401537 mg/m3 tương ứng 0.305869 0.401537 µg/l theo TCVN 5937 – 2005 – Trung bình ngày SO 300 µg/l -> Nồng độ SO điểm nhạy cảm thấp tiêu chuẩn 4.3.2.3 Kịch CO Hình 4.32: Kịch mô câu b – CO 39 Hình 4.33: Kết chạy mô hình – Trang nồng độ CO điểm nhạy cảm - Dựa vào kết chạy mô hình, nồng độ điểm nhạy cảm 0.001095 0.001606 mg/m3 tương ứng 0.001095 0.001606 µg/l theo TCVN 5937 – 2005 – Trung bình ngày CO 5000 µg/l -> Nồng độ CO điểm nhạy cảm thấp tiêu chuẩn 4.3.2.4 Kịch NO2 Hình 4.34: Kịch mô câu b – NO2 40 Hình 4.35: Kết chạy mô hình – Trang nồng độ NO2 điểm nhạy cảm - Dựa vào kết chạy mô hình, nồng độ điểm nhạy cảm 0.067826 0.101634 mg/m3 tương ứng 0.067826 0.101634 µg/l theo TCVN 5937 – 2005 – Trung bình ngày NO 400 µg/l -> Nồng độ NO2 điểm nhạy cảm thấp tiêu chuẩn 41 V KẾT LUẬN Qua tìm hiểu chạy thử mô hình nhóm chúng em thấy số tính phần mềm ENVIMAP sau: - Trợ giúp công tác quản lý tổng hợp thống CSDL môi trường CSSX thuộc khu công nghiệp - Cho phép vẽ đồ ô nhiễm theo kịch khác - Tính toán mô lan truyền ô nhiễm không khí - Biểu diễn kết chạy mô hình dạng bảng, đồ thị, đồ số 42 TÀI LIỆU THAM KHẢO http://luanvan.net.vn/luan-van/luan-van-ung-dung-cong-nghe-thong-tin-giam-sat-chatluong-khong-khi-khu-cong-nghiep-amata-loteco-dong-nai-36636/ http://www.kilobooks.com/ung-dung-phan-mem-envimap-giam-sat-o-nhiem-khong- khi-cho-khu-kinh-te-dung-quat-quang-ngai-57760 http://www.zbook.vn/ebook/xay-dung-mo-hinh-giam-sat-moi-truong-khong-khi-khu- cong-nghiep-bien-hoa-1-bang-ky-thuat-tin-hoc-38505/ http://blogv.p-gis.com/2011/02/giam-sat-o-nhiem-khong-khi.html Bùi Tá Long, Các phần mềm môi trường, Đại học Quốc gia Hồ Chí Minh, năm 2010 Bùi Tá Long, Bài giảng Mô hình hóa môi trường 43 [...]... khu công nghiệp Mở phần mềm ENVIMAP -> vào menu Thông tin -> Khu công nghiệp -> Thêm mới -> Xuất hiện hộp thoại như hình 4.2 -> nhập các thông tin cần thiết -> Chọn lưu 22 Hình 4.2: Hộp thoại thông tin khu công nghiệp Hình 4.3: Cửa sổ thông tin khu công nghiệp 4.2.3 Dữ liệu cơ sở sản xuất Trên thanh công cụ Envimap chọn “Thông tin” -> Cơ sở sản xuất -> Thêm mới -> Xuất hiện hộp thoại như hình 4.4 -> ... tính sự phát tán ô nhiễm cho các nguồn OK1, OK2, OK3, OK4, OK5 OK6 theo trung bình ngày tại 2 điểm nhạy cảm 30 4.2.6.2.1 Tạo điểm nhạy cảm - Vào menu thông tin -> chọn Điểm nhạy cảm -> Xuất hiện hộp thoại điểm nhạy cảm -> chọn thêm mới -> Xuất hiện hộp thoại -> nhập các thông tin điểm nhạy cảm -> Chọn lưu Hình 4.17: Hộp thoại điểm nhạy cảm Hình 4.18: Thông tin điểm nhạy cảm 4.2.6.2.2 Tạo kịch bản -. .. nặng thì thành phần ở phương trình 2.1 sẽ bằng tốc độ rơi của hạt (dấu âm), còn đối với chất ô nhiễm khí và bụi nhẹ thì =0 Nếu bỏ qua hiện tượng chuyển “pha” 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 khuếch tán thì ==0 Như vậy có thể sử dụng phương trình mô tả sự phân tán các chất ô nhiễm từ nguồn điểm sau đây vào mục đích tính toán sự nhiễm bẩn không khí: 12 (2.5)... Chọn kịch bản chạy mô hình - Click tiếp -> Xuất hiện hộp thoại (hình 4.25) 34 Hình 4.25: Nhập thông số cho lưới tính - Click Hoàn tất để chạy mô hình 4.3 Kết quả chạy mộ hình - Sau khi chạy mô hình ta thu được một số kết quả như sau: 35 4.3.1 Kịch bản câu a Hình 4.26: Bản đồ mô phỏng – kịch bản câu a Hình 4.27: Kết quả chạy mô hình – Trang nồng độ cực đại - Theo yêu cầu bài toán là tính nồng độ cực đại... điểm tính toán - Clặng(i): Nồng độ tức thời do nguồn thải i gây ra tại điểm tính toán khi lặng gió (u=0) - Ca(i): Nồng độ tức thời do nguồn thải i gây ra tại điểm tính toán khi có gió thỗi theo hướng a ứng với vận tốc trung bình trên hướng gió và tốc độ ổn định trung bình của khí quyển trong suốt khoảng thời gian tính toán trị số trung bình ( ngày đêm, tháng hoặc năm) Trị số Ca(i)=0 khi điểm tính toán. .. liệu phát thải Trên thanh công cụ Envimap chọn menu Số liệu -> Phát thải tại ống khói -> Thêm mới -> Xuất hiện hộp thoại như hình 4.8 -> Chọn ngày tính toán -> chọn lưu 25 Hình 4.8: Tạo thông tin phát thải cho 1 ngày cụ thể - Click 2 lần vào số 1 xuất hiện hộp thoại (hình 4.9) -> Click tiếp -> chọn các thông số đo (hình 4.10) -> Click tiếp -> Xuất hiện hộp thoại (hình 4.11) -> Nhập số liệu phát thải... xung quanh - Phương pháp Gauss - Phương pháp ISC3 - Phương pháp Berliand Trong đề tài này chúng em chủ yếu tìm hiểu về phương pháp Berliand Công thức Berliand trong trường hợp chất khí và bụi nặng Để giải bài toán (2.5), (2.6), (2.8), (2.10) với điều kiện (2.11), Berliand và các học trò của mình đã sử dụng nhiều công cụ toán học khác nhau Để áp dụng vào thực tế tính toán ô nhiễm không khí, Berliand... Oy, Oz 11 ,,: Các thành phần của tốc độ trung bình theo 3 trục Ox, Oy,Oz : Hệ dố tính đến sự liên kết của chất ô nhiễm với các thành 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 quá trình phản ứng chất hóa học xảy ra trên đường lan truyền Phương trình 2.1 rất phức tạp (mặc dù vậy nó chỉ mô phỏng sự lan truyền chất ô nhiễm) trên thực tế để giải... Xuất hiện hộp thoại như hình 4.4 -> Nhập thông tin các công ty -> Chọn lưu 23 Hình 4.4: Hộp thoại thông tin cơ sở sản xuất Hình 4.5: Cửa sổ thông tin cơ sở sản xuất 4.2.4 Dữ liệu ống khói Trên thanh công cụ Envimap chọn “Thông tin” -> Ống khói -> Thêm mới -> Xuất hiện hộp thoại như hình 4.6 -> Nhập thông tin các ống khói -> Chọn lưu 24 Hình 4.6: Hộp thoại thông tin các ống khói Hình 4.7: Cửa sổ thộng... chất ô nhiễm (mg/m3) M: Công suất nguồn thải , : Là các hàm Dirac Nếu như nguồn thải không phải ống khói mà thải ra mặt đất thì một số tác giả cho rằng, tại thời điểm t=0 chất phát thải chưa hoạt động, khi đó giả sử nguồn đặt ở gốc tọa độ thì: (2.7) 2.3.2 Điều kiện biên Trong lớp không khí khảo sát thường giới hạn bởi mặt đất, còn độ cao thường là vô hạn hoặc hữu hạn tùy theo sự phân lớp của khí quyển