TỔNG LIÊN ĐOÀN LAO ĐỘNG VIỆT NAM TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA MÔI TRƯỜNG VÀ BẢO HỘ LAO ĐỘNG NGÀNH KỸ THUẬT MƠI TRƯỜNG KHĨA LUẬN TỐT NGHIỆP ỨNG DỤNG MƠ HÌNH CORMIX MƠ PHỎNG Ơ NHIỄM NƯỚC SÔNG THỊ VẢI – XÉT NGUỒN THẢI VEDAN SVTH : NGUYỄN DUY HIẾU MSSV : 911086B Lớp : 09MT1N GVHD : Th.s Trần Thị Vân TP Hồ Chí Minh, tháng 12 năm 2009 TỔNG LIÊN ĐOÀN LAO ĐỘNG VIỆT NAM TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA MÔI TRƯỜNG VÀ BẢO HỘ LAO ĐỘNG NGÀNH KỸ THUẬT MƠI TRƯỜNG KHĨA LUẬN TỐT NGHIỆP ỨNG DỤNG MƠ HÌNH CORMIX MƠ PHỎNG Ô NHIỄM NƯỚC SÔNG THỊ VẢI – XÉT NGUỒN THẢI VEDAN Ngày giao nhiệm vụ luận văn: 16/09/2009 Ngày hoàn thành luận văn : 16/12/2009 Xác nhận GVHD TP Hồ Chí Minh, tháng 12 năm 2009 Lời cảm ơn Em xin gửi lời cảm ơn chân thành sâu sắc đến PGS.TSKH Bùi Tá Long, người tận tình giúp đỡ tạo điều kiện cho em nghiên cứu thực đề tài Xin cám ơn thầy lời động viên dạy bảo giúp em có thêm động lực để hồn thành luận văn Em xin gửi lời cảm ơn chân thành tới Th.S Trần Thị Vân, người hướng dẫn nhiệt tình, ủng hộ em gặp khó khăn tạo điều kiện thuận lợi giúp em thực luận văn tốt nghiệp thời gian vừa qua Em xin bày tỏ lịng lịng biết ơn đến thầy Khoa Môi Trường & Bảo Hộ Lao Động, thầy cô ân c ần dạy dỗ truyền đạt cho em kiến thức suốt nhũng năm học vừa qua Em xin gửi lời cảm ơn đến KS Nguyễn Thái Hòa, KS Lê Thị Hiền Th.S Cao Duy Trường phòng Tin học Môi trường, Viện Môi Trường Tài Nguyên, Đại học Quốc gia TP.HCM t ận tình giúp đ ỡ bảo em thời gian thực luận văn tốt nghiệp Em xin gửi lời cám ơn đến công ty quốc tế Mixzon t ạo điều kiện cho em sử dụng mơ hình CORMIX Xin gửi lời cám ơn đến bạn lớp 09MT1N, người bạn gắn bó với thời gian học tập, người bạn đem l ại cho niềm vui ngồi ghế giảng đường Cuối xin gửi lời cám ơn đến ba mẹ, ba mẹ hết lịng đ ộng viên khích lệ, tạo điều kiện cho hồn cảnh khó khăn Con xin cám ơn ba mẹ quản khó nhọc ni dạy nên người Một lần em xin chân thành cảm ơn tất Tp.Hồ Chí Minh, ngày 16 tháng 12 năm 2009 Nguyễn Duy Hiếu TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA MÔI TRƯỜNG VÀ BẢO HỘ LAO ĐỘNG Đề tài khóa luận tốt nghiệp: ỨNG DỤNG MƠ HÌNH CORMIX MƠ PHỎNG Ơ NHIỄM NƯỚC SÔNG THỊ VẢI – XÉT NGUỒN THẢI VEDAN Sinh viên thực : Nguyễn Duy Hiếu Giáo viên hướng dẫn : Th.S Trần Thị Vân NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN -TP Hồ Chí Minh, ngày tháng Giáo viên hướng dẫn năm TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA MÔI TRƯỜNG VÀ BẢO HỘ LAO ĐỘNG Đề tài khóa luận tốt nghiệp: ỨNG DỤNG MƠ HÌNH CORMIX MƠ PHỎNG Ơ NHIỄM NƯỚC SƠNG THỊ VẢI – XÉT NGUỒN THẢI VEDAN Sinh viên thực : Nguyễn Duy Hiếu Giáo viên hướng dẫn : Th.S Trần Thị Vân Giáo viên phản biện : NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN -TP Hồ Chí Minh, ngày tháng Giáo viên phản biện năm Mục lục DANH MỤC CHỮ VIẾT TẮT iii DANH MỤC THUẬT NGỮ TIẾNG ANH iv DANH MỤC KÍ HIỆU TOÁN HỌC vi DANH MỤC BẢNG vii DANH MỤC HÌNH ẢNH viii TÓM TẮT LUẬN VĂN x MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ MÔI TRƯỜNG LƯU VỰC SÔNG THỊ VẢI 1.1 ĐIỀU KIỆN TỰ NHIÊN VÀ ĐẶC ĐIỂM KHÍ TƯỢNG THỦY VĂN SƠNG THỊ VẢI 1.1.1 Đặc điểm địa hình 1.1.2 Đặc điểm khí hậu .4 1.2 HIỆN TRẠNG Ô NHIỄM SÔNG THỊ VẢI 1.2.1 Độ màu 11 1.2.2 Độ pH .11 1.2.3 Oxy hòa tan (DO) 11 1.2.4 Nhu cầu oxy sinh hóa (BOD) 12 1.2.5 Nhu cầu oxy hóa học (COD) 13 1.2.6 Amoni (N-NH4+) 14 1.2.7 Tổng Nitơ 15 1.2.8 Tổng Photpho 16 1.3 TỔNG QUAN CÔNG TY VEDAN VIỆT NAM 17 1.3.1 Hiện trạng xả thải Vedan 18 1.3.2 Những sai phạm Vedan .20 CHƯƠNG TỔNG QUAN MƠ HÌNH ĐƯỢC SỬ DỤNG 22 2.1 MỞ ĐẦU 22 2.1.1 Định nghĩa mơ hình mơ 22 2.1.2 Những lợi ích mơ hình 22 2.1.3 Vai trò ý nghĩa mơ hình hóa quản lý 22 2.2 TỔNG QUAN MƠ HÌNH CORMIX .23 i 2.2.1 Mơ hình mô thủy động lực học 24 2.2.2 Các công cụ hỗ trợ 26 2.2.3 Trình tự xử lý cấu trúc hệ thống 27 2.3 PHƯƠNG PHÁP TÍNH TỐN TRONG CORMIX .28 2.3.1 Q trình hịa trộn dịng phun có xu hướng 29 2.3.2 Tổng quát hình dạng loại dịng thải 32 2.3.3 Tỉ lệ độ dài .34 2.3.4 Kiểu dịng chảy điển hình 36 2.3.5 Quá trình lan tỏa 38 2.3.6 Khuếch tán thụ động 38 CHƯƠNG TÍNH TỐN MƠ PHỎNG PHẠM VI VÀ MỨC ĐỘ ẢNH HƯỞNG CỦA VEDAN…… 40 3.1 SƠ ĐỒ TIẾN HÀNH MÔ PHỎNG .40 3.2 THÔNG SỐ ĐẦU VÀO VÀ CÁC THIẾT LẬP CHO MƠ HÌNH .40 3.2.1 Dữ liệu nguồn thải 40 3.2.2 Dữ liệu sông Thị Vải .42 3.3 KỊCH BẢN MÔ PHỎNG 48 3.4 KẾT QUẢ MÔ PHỎNG 50 3.4.1 Vị trí xả 50 3.4.2 Vị trí xả 56 3.5 Thảo luận 62 TÀI LIỆU THAM KHẢO 66 ii DANH MỤC CHỮ VIẾT TẮT BTNMT Bộ Tài nguyên Mơi trường BOD Nhu cầu oxy sinh hóa COD Nhu cầu oxy hóa học DO Oxy hịa tan GUI Giao diện sử dụng KCN Khu công nghiệp QCVN Quy chuẩn Việt Nam US EPA Cơ quan Bảo vệ Môi trường Mỹ VTX Vị trí xà iii DANH MỤC THUẬT NGỮ TIẾNG ANH Ambient crossflow Dòng chảy mặt cắt Buoyancy flux Tính dịng thải Buoyant jet Dịng phun Density current Mật độ dòng Density stratification Mật độ phân tầng Discharge volumn flux Lưu lượng thải Discharge length scale Phạm vi thải Far-field Khu vực tiếp nhận nước thải nơi trình lan rộng khuyếch tán thụ động kiểm soát hướng pha loảng luồng thải Jet/Plume length scale Tỉ lệ phụt/luồng Jet/Stratification length scale Tỉ lệ phụt/phân tầng Jet/Crossflow Length Scale Tỉ lệ phụt/dòng chảy ngang Momentum flux Động dòng thải Near-field Vùng tiếp nhận nước thải ban đầu, nơi đặc trưng động tính đặc điểm hình học cống ảnh hưởng đến hướng phun hòa trộn dòng thải Pure jet Nơi dòng thải ban đầu có động dạng tốc độ cao gây nên hịa trộn rối Nó gọi momentum jet nonbuoyancy jet Dòng phun túy Pure plume Nơi dòng thải ban đầu có khuynh hướng lên dẫn tới dâng lên theo phương đứng sau dẫn tới hịa trơn rối Dịng luồng túy Plume/crossflow length scale Tỉ lệ luồng/dòng chảy ngang Plume/stratification length scale Tỉ lệ luồng/phân tầng Pycnocline Lớp nằm ngang vùng tiếp nhận nước thải iv nơi có thay đổi đột ngột khối lượng riêng Passive diffusion Khuếch tán thụ động Stagnant Vùng tù đọng Weakly deflect jet Dòng phun lệch yếu dòng chảy cắt ngang v 29.29 -9.21 3.86 62.9 0.840E+03 2.22 0.043 73141E+02 29.59 -9.22 3.88 63.3 0.834E+03 2.23 0.042 74020E+02 29.99 -9.23 3.91 63.8 0.827E+03 2.24 0.041 75194E+02 30.39 -9.24 3.93 64.4 0.820E+03 2.25 0.041 76370E+02 30.79 -9.26 3.96 64.9 0.814E+03 2.26 0.040 77548E+02 31.19 -9.27 3.98 65.3 0.808E+03 2.27 0.039 78727E+02 31.59 -9.28 4.00 65.8 0.803E+03 2.27 0.039 79908E+02 31.99 -9.30 4.01 66.2 0.797E+03 2.28 0.038 81091E+02 32.39 -9.31 4.03 66.6 0.792E+03 2.29 0.038 82275E+02 32.79 -9.32 4.04 67.0 0.788E+03 2.30 0.038 83461E+02 33.19 -9.34 4.05 67.4 0.783E+03 2.30 0.037 84647E+02 33.59 -9.35 4.06 67.8 0.779E+03 2.31 0.037 85835E+02 33.99 -9.36 4.07 68.1 0.775E+03 2.32 0.037 87024E+02 34.39 -9.37 4.08 68.4 0.772E+03 2.32 0.036 88213E+02 34.78 -9.39 4.08 68.8 0.768E+03 2.33 0.036 89404E+02 35.18 -9.40 4.08 69.1 0.765E+03 2.33 0.036 90595E+02 Maximum jet height has been reached 35.48 -9.41 4.08 69.3 0.762E+03 2.34 0.036 91488E+02 Terminal level in stratified ambient has been reached Cumulative travel time = 91.4882 sec END OF CORJET (MOD110): JET/PLUME NEAR-FIELD MIXING REGION -BEGIN MOD131: LAYER BOUNDARY/TERMINAL LAYER APPROACH Control volume inflow: X 35.48 Y -9.41 Z S 4.08 C B TT 68.5 0.771E+03 2.34 91488E+02 Profile definitions: BV = top-hat thickness, measured vertically BH = top-hat half-width, measured horizontally in Y-direction ZU = upper plume boundary (Z-coordinate) ZL = lower plume boundary (Z-coordinate) S = hydrodynamic average (bulk) dilution C = average (bulk) concentration (includes reaction effects, if any) TT = Cumulative travel time X Y Z S C BV BH ZU ZL TT 33.15 -9.37 10.00 69.3 0.762E+03 0.00 0.00 10.00 10.00 91488E+02 33.85 -9.38 10.00 69.3 0.762E+03 2.95 1.49 10.00 7.05 91488E+02 34.55 -9.39 10.00 69.3 0.762E+03 3.50 2.10 10.00 6.50 91488E+02 35.25 -9.41 10.00 69.3 0.762E+03 3.85 2.57 10.00 6.15 91488E+02 35.95 -9.42 10.00 71.2 0.742E+03 4.11 2.97 10.00 5.89 93050E+02 36.65 -9.43 10.00 80.0 0.660E+03 4.31 3.32 10.00 5.69 95388E+02 37.36 -9.44 10.00 92.1 0.573E+03 4.45 3.64 10.00 5.55 97727E+02 38.06 -9.45 10.00 103.2 0.512E+03 4.56 3.93 10.00 5.44 10007E+03 38.76 -9.46 10.00 110.8 0.476E+03 4.64 4.20 10.00 5.36 10240E+03 39.46 -9.47 10.00 114.9 0.459E+03 4.68 4.46 10.00 5.32 10474E+03 40.16 -9.48 10.00 117.6 0.449E+03 4.70 4.70 10.00 5.30 10708E+03 Cumulative travel time = 107.0805 sec END OF MOD131: LAYER BOUNDARY/TERMINAL LAYER APPROACH -** End of NEAR-FIELD REGION (NFR) ** -BEGIN MOD141: BUOYANT AMBIENT SPREADING Profile definitions: BV = top-hat thickness, measured vertically BH = top-hat half-width, measured horizontally in Y-direction ZU = upper plume boundary (Z-coordinate) ZL = lower plume boundary (Z-coordinate) S = hydrodynamic average (bulk) dilution C = average (bulk) concentration (includes reaction effects, if any) TT = Cumulative travel time Plume Stage (not bank attached): X Y Z S C BV BH ZU ZL TT 40.16 -9.48 10.00 117.6 0.449E+03 4.70 4.70 10.00 5.30 10708E+03 41.76 -9.48 10.00 119.4 0.442E+03 4.49 4.98 10.00 5.51 11241E+03 43.36 -9.48 10.00 121.1 0.436E+03 4.32 5.26 10.00 5.68 11775E+03 44.96 -9.48 10.00 122.6 0.431E+03 4.16 5.54 10.00 5.84 12308E+03 46.56 -9.48 10.00 124.1 0.425E+03 4.02 5.80 10.00 5.98 12841E+03 48.16 -9.48 10.00 125.5 0.421E+03 3.89 6.06 10.00 6.11 13375E+03 49.76 -9.48 10.00 126.9 0.416E+03 3.78 6.32 10.00 6.22 13908E+03 51.36 -9.48 10.00 128.2 0.412E+03 3.67 6.57 10.00 6.33 14441E+03 52.96 -9.48 10.00 129.5 0.408E+03 3.58 6.81 10.00 6.42 14975E+03 54.56 -9.48 10.00 130.7 0.404E+03 3.49 7.05 10.00 6.51 15508E+03 56.16 -9.48 10.00 131.9 0.400E+03 3.41 7.29 10.00 6.59 16041E+03 57.76 -9.48 10.00 133.0 0.397E+03 3.33 7.52 10.00 6.67 16575E+03 59.36 -9.48 10.00 134.1 0.394E+03 3.26 7.75 10.00 6.74 17108E+03 60.96 -9.48 10.00 135.2 0.390E+03 3.20 7.97 10.00 6.80 17641E+03 62.56 -9.48 10.00 136.3 0.387E+03 3.14 8.20 10.00 6.86 18175E+03 64.16 -9.48 10.00 137.3 0.384E+03 3.08 8.42 10.00 6.92 18708E+03 65.76 -9.48 10.00 138.4 0.381E+03 3.03 8.63 10.00 6.97 19241E+03 67.36 -9.48 10.00 139.4 0.379E+03 2.98 8.85 10.00 7.02 19775E+03 68.96 -9.48 10.00 140.4 0.376E+03 2.93 9.06 10.00 7.07 20308E+03 70.56 -9.48 10.00 141.4 0.373E+03 2.89 9.27 10.00 7.11 20841E+03 72.16 -9.48 10.00 142.4 0.371E+03 2.84 9.47 10.00 7.16 21375E+03 73.76 -9.48 10.00 143.4 0.368E+03 2.80 9.68 10.00 7.20 21908E+03 75.36 -9.48 10.00 144.3 0.366E+03 2.77 9.88 10.00 7.23 22441E+03 76.96 -9.48 10.00 145.3 0.363E+03 2.73 10.08 10.00 7.27 22975E+03 78.56 -9.48 10.00 146.3 0.361E+03 2.70 10.28 10.00 7.30 23508E+03 80.16 -9.48 10.00 147.2 0.358E+03 2.66 10.48 10.00 7.34 24041E+03 81.76 -9.48 10.00 148.2 0.356E+03 2.63 10.67 10.00 7.37 24575E+03 83.36 -9.48 10.00 149.1 0.354E+03 2.60 10.86 10.00 7.40 25108E+03 84.96 -9.48 10.00 150.1 0.352E+03 2.58 11.05 10.00 7.42 25641E+03 86.56 -9.48 10.00 151.0 0.349E+03 2.55 11.24 10.00 7.45 26175E+03 88.16 -9.48 10.00 152.0 0.347E+03 2.53 11.43 10.00 7.47 26708E+03 89.76 -9.48 10.00 152.9 0.345E+03 2.50 11.62 10.00 7.50 27241E+03 91.36 -9.48 10.00 153.9 0.343E+03 2.48 11.80 10.00 7.52 27775E+03 92.96 -9.48 10.00 154.8 0.341E+03 2.46 11.99 10.00 7.54 28308E+03 94.56 -9.48 10.00 155.8 0.339E+03 2.44 12.17 10.00 7.56 28841E+03 96.16 -9.48 10.00 156.7 0.337E+03 2.42 12.35 10.00 7.58 29375E+03 97.76 -9.48 10.00 157.7 0.334E+03 2.40 12.53 10.00 7.60 29908E+03 99.36 -9.48 10.00 158.7 0.332E+03 2.38 12.71 10.00 7.62 30441E+03 100.96 -9.48 10.00 159.6 0.330E+03 2.36 12.88 10.00 7.64 30974E+03 102.56 -9.48 10.00 160.6 0.328E+03 2.35 13.06 10.00 7.65 31508E+03 104.16 -9.48 10.00 161.6 0.326E+03 2.33 13.23 10.00 7.67 32041E+03 105.76 -9.48 10.00 162.6 0.324E+03 2.32 13.41 10.00 7.68 32574E+03 107.36 -9.48 10.00 163.6 0.322E+03 2.30 13.58 10.00 7.70 33108E+03 108.96 -9.48 10.00 164.6 0.320E+03 2.29 13.75 10.00 7.71 33641E+03 110.56 -9.48 10.00 165.6 0.318E+03 2.27 13.92 10.00 7.73 34174E+03 112.16 -9.48 10.00 166.6 0.316E+03 2.26 14.09 10.00 7.74 34708E+03 113.76 -9.48 10.00 167.6 0.315E+03 2.25 14.26 10.00 7.75 35241E+03 115.36 -9.48 10.00 168.7 0.313E+03 2.24 14.42 10.00 7.76 35774E+03 116.96 -9.48 10.00 169.7 0.311E+03 2.23 14.59 10.00 7.77 36308E+03 118.56 -9.48 10.00 170.8 0.309E+03 2.22 14.75 10.00 7.78 36841E+03 120.16 -9.48 10.00 171.8 0.307E+03 2.21 14.92 10.00 7.79 37374E+03 121.76 -9.48 10.00 172.9 0.305E+03 2.20 15.08 10.00 7.80 37908E+03 123.36 -9.48 10.00 174.0 0.303E+03 2.19 15.24 10.00 7.81 38441E+03 124.96 -9.48 10.00 175.0 0.301E+03 2.18 15.40 10.00 7.82 38974E+03 126.56 -9.48 10.00 176.1 0.299E+03 2.17 15.56 10.00 7.83 39508E+03 128.16 -9.48 10.00 177.3 0.297E+03 2.16 15.72 10.00 7.84 40041E+03 129.76 -9.48 10.00 178.4 0.296E+03 2.16 15.88 10.00 7.84 40574E+03 131.36 -9.48 10.00 179.5 0.294E+03 2.15 16.04 10.00 7.85 41108E+03 132.96 -9.48 10.00 180.6 0.292E+03 2.14 16.20 10.00 7.86 41641E+03 134.56 -9.48 10.00 181.8 0.290E+03 2.14 16.36 10.00 7.86 42174E+03 136.16 -9.48 10.00 182.9 0.288E+03 2.13 16.51 10.00 7.87 42708E+03 137.76 -9.48 10.00 184.1 0.286E+03 2.13 16.67 10.00 7.87 43241E+03 139.36 -9.48 10.00 185.3 0.284E+03 2.12 16.82 10.00 7.88 43774E+03 140.96 -9.48 10.00 186.5 0.283E+03 2.12 16.97 10.00 7.88 44308E+03 142.56 -9.48 10.00 187.7 0.281E+03 2.11 17.13 10.00 7.89 44841E+03 144.16 -9.48 10.00 188.9 0.279E+03 2.11 17.28 10.00 7.89 45374E+03 145.76 -9.48 10.00 190.1 0.277E+03 2.10 17.43 10.00 7.90 45908E+03 147.36 -9.48 10.00 191.4 0.275E+03 2.10 17.58 10.00 7.90 46441E+03 148.96 -9.48 10.00 192.6 0.274E+03 2.10 17.73 10.00 7.90 46974E+03 150.56 -9.48 10.00 193.9 0.272E+03 2.09 17.88 10.00 7.91 47508E+03 152.16 -9.48 10.00 195.2 0.270E+03 2.09 18.03 10.00 7.91 48041E+03 153.76 -9.48 10.00 196.5 0.268E+03 2.09 18.18 10.00 7.91 48574E+03 155.36 -9.48 10.00 197.8 0.266E+03 2.09 18.33 10.00 7.91 49108E+03 156.96 -9.48 10.00 199.1 0.265E+03 2.08 18.47 10.00 7.92 49641E+03 158.56 -9.48 10.00 200.4 0.263E+03 2.08 18.62 10.00 7.92 50174E+03 160.16 -9.48 10.00 201.8 0.261E+03 2.08 18.77 10.00 7.92 50708E+03 161.76 -9.48 10.00 203.1 0.259E+03 2.08 18.91 10.00 7.92 51241E+03 163.36 -9.48 10.00 204.5 0.258E+03 2.08 19.06 10.00 7.92 51774E+03 164.96 -9.48 10.00 205.9 0.256E+03 2.08 19.20 10.00 7.92 52308E+03 166.56 -9.48 10.00 207.3 0.254E+03 2.08 19.35 10.00 7.92 52841E+03 168.16 -9.48 10.00 208.7 0.252E+03 2.08 19.49 10.00 7.92 53374E+03 Cumulative travel time = 533.7428 sec -Plume is ATTACHED to LEFT bank/shore Plume width is now determined from LEFT bank/shore Plume Stage (bank attached): X Y Z S C BV BH ZU ZL TT 168.16 10.00 10.00 208.7 0.252E+03 2.08 38.98 10.00 7.92 53374E+03 179.10 10.00 10.00 217.3 0.242E+03 2.12 39.95 10.00 7.88 57021E+03 190.04 10.00 10.00 226.2 0.233E+03 2.16 40.90 10.00 7.84 60667E+03 200.98 10.00 10.00 235.4 0.224E+03 2.20 41.85 10.00 7.80 64313E+03 211.92 10.00 10.00 245.1 0.215E+03 2.25 42.79 10.00 7.75 67960E+03 222.86 10.00 10.00 255.1 0.206E+03 2.29 43.71 10.00 7.71 71606E+03 233.80 10.00 10.00 265.5 0.198E+03 2.34 44.63 10.00 7.66 75252E+03 244.73 10.00 10.00 276.2 0.190E+03 2.40 45.53 10.00 7.60 78899E+03 255.67 10.00 10.00 287.4 0.183E+03 2.45 46.43 10.00 7.55 82545E+03 266.61 10.00 10.00 298.9 0.176E+03 2.51 47.32 10.00 7.49 86191E+03 277.55 10.00 10.00 310.8 0.169E+03 2.56 48.20 10.00 7.44 89838E+03 288.49 10.00 10.00 323.0 0.163E+03 2.62 49.07 10.00 7.38 93484E+03 299.43 10.00 10.00 335.7 0.157E+03 2.69 49.94 10.00 7.31 97130E+03 310.37 10.00 10.00 348.8 0.151E+03 2.75 50.79 10.00 7.25 10078E+04 321.31 10.00 10.00 362.2 0.145E+03 2.81 51.64 10.00 7.19 10442E+04 332.25 10.00 10.00 376.1 0.140E+03 2.88 52.49 10.00 7.12 10807E+04 343.19 10.00 10.00 390.4 0.134E+03 2.95 53.33 10.00 7.05 11172E+04 354.12 10.00 10.00 405.0 0.130E+03 3.02 54.16 10.00 6.98 11536E+04 365.06 10.00 10.00 420.1 0.125E+03 3.09 54.98 10.00 6.91 11901E+04 376.00 10.00 10.00 435.6 0.120E+03 3.17 55.80 10.00 6.83 12265E+04 386.94 10.00 10.00 451.4 0.116E+03 3.24 56.61 10.00 6.76 12630E+04 397.88 10.00 10.00 467.7 0.112E+03 3.32 57.42 10.00 6.68 12995E+04 408.82 10.00 10.00 484.5 0.108E+03 3.39 58.22 10.00 6.61 13359E+04 419.76 10.00 10.00 501.6 0.105E+03 3.47 59.01 10.00 6.53 13724E+04 430.70 10.00 10.00 519.1 0.101E+03 3.55 59.80 10.00 6.45 14089E+04 441.64 10.00 10.00 537.1 0.976E+02 3.64 60.59 10.00 6.36 14453E+04 452.58 10.00 10.00 555.5 0.943E+02 3.72 61.37 10.00 6.28 14818E+04 463.51 10.00 10.00 574.3 0.912E+02 3.81 62.14 10.00 6.19 15183E+04 474.45 10.00 10.00 593.6 0.882E+02 3.89 62.91 10.00 6.11 15547E+04 485.39 10.00 10.00 613.3 0.854E+02 3.98 63.68 10.00 6.02 15912E+04 496.33 10.00 10.00 633.4 0.826E+02 4.07 64.44 10.00 5.93 16276E+04 507.27 10.00 10.00 653.9 0.800E+02 4.16 65.20 10.00 5.84 16641E+04 518.21 10.00 10.00 674.9 0.775E+02 4.25 65.95 10.00 5.75 17006E+04 529.15 10.00 10.00 696.3 0.751E+02 4.34 66.70 10.00 5.66 17370E+04 540.09 10.00 10.00 718.2 0.728E+02 4.44 67.44 10.00 5.56 17735E+04 551.03 10.00 10.00 740.5 0.706E+02 4.53 68.18 10.00 5.47 18100E+04 561.97 10.00 10.00 763.3 0.685E+02 4.63 68.92 10.00 5.37 18464E+04 572.90 10.00 10.00 786.5 0.665E+02 4.73 69.65 10.00 5.27 18829E+04 583.84 10.00 10.00 810.1 0.645E+02 4.83 70.38 10.00 5.17 19194E+04 594.78 10.00 10.00 834.2 0.626E+02 4.93 71.10 10.00 5.07 19558E+04 605.72 10.00 10.00 858.8 0.608E+02 5.03 71.82 10.00 4.97 19923E+04 616.66 10.00 10.00 883.8 0.591E+02 5.13 72.54 10.00 4.87 20287E+04 627.60 10.00 10.00 909.3 0.574E+02 5.24 73.25 10.00 4.76 20652E+04 638.54 10.00 10.00 935.2 0.558E+02 5.34 73.96 10.00 4.66 21017E+04 649.48 10.00 10.00 961.6 0.543E+02 5.45 74.67 10.00 4.55 21381E+04 660.42 10.00 10.00 988.5 0.528E+02 5.56 75.37 10.00 4.44 21746E+04 671.36 10.00 10.00 1015.8 0.514E+02 5.66 76.07 10.00 4.34 22111E+04 682.30 10.00 10.00 1043.6 0.500E+02 5.77 76.77 10.00 4.23 22475E+04 693.23 10.00 10.00 1071.9 0.486E+02 5.88 77.47 10.00 4.12 22840E+04 704.17 10.00 10.00 1100.6 0.474E+02 6.00 78.16 10.00 4.00 23205E+04 715.11 10.00 10.00 1129.8 0.461E+02 6.11 78.85 10.00 3.89 23569E+04 726.05 10.00 10.00 1159.5 0.449E+02 6.22 79.53 10.00 3.78 23934E+04 736.99 10.00 10.00 1189.7 0.438E+02 6.34 80.21 10.00 3.66 24298E+04 747.93 10.00 10.00 1220.3 0.427E+02 6.46 80.89 10.00 3.54 24663E+04 758.87 10.00 10.00 1251.5 0.416E+02 6.57 81.57 10.00 3.43 25028E+04 769.81 10.00 10.00 1283.1 0.406E+02 6.69 82.24 10.00 3.31 25392E+04 780.75 10.00 10.00 1315.2 0.396E+02 6.81 82.91 10.00 3.19 25757E+04 791.69 10.00 10.00 1347.8 0.386E+02 6.93 83.58 10.00 3.07 26122E+04 802.62 10.00 10.00 1380.9 0.377E+02 7.06 84.25 10.00 2.94 26486E+04 813.56 10.00 10.00 1414.4 0.368E+02 7.18 84.91 10.00 2.82 26851E+04 824.50 10.00 10.00 1448.5 0.359E+02 7.30 85.57 10.00 2.70 27215E+04 835.44 10.00 10.00 1483.1 0.351E+02 7.43 86.23 10.00 2.57 27580E+04 846.38 10.00 10.00 1518.1 0.342E+02 7.55 86.89 10.00 2.45 27945E+04 857.32 10.00 10.00 1553.7 0.335E+02 7.68 87.54 10.00 2.32 28309E+04 868.26 10.00 10.00 1589.8 0.327E+02 7.81 88.19 10.00 2.19 28674E+04 879.20 10.00 10.00 1626.3 0.319E+02 7.94 88.84 10.00 2.06 29039E+04 890.14 10.00 10.00 1663.4 0.312E+02 8.07 89.48 10.00 1.93 29403E+04 901.08 10.00 10.00 1701.0 0.305E+02 8.20 90.13 10.00 1.80 29768E+04 912.01 10.00 10.00 1739.1 0.299E+02 8.33 90.77 10.00 1.67 30133E+04 922.95 10.00 10.00 1777.7 0.292E+02 8.47 91.41 10.00 1.53 30497E+04 933.89 10.00 10.00 1816.9 0.286E+02 8.60 92.05 10.00 1.40 30862E+04 944.83 10.00 10.00 1856.5 0.280E+02 8.74 92.68 10.00 1.26 31226E+04 955.77 10.00 10.00 1896.7 0.274E+02 8.87 93.31 10.00 1.13 31591E+04 966.71 10.00 10.00 1937.4 0.268E+02 9.01 93.94 10.00 0.99 31956E+04 977.65 10.00 10.00 1978.6 0.262E+02 9.15 94.57 10.00 0.85 32320E+04 988.59 10.00 10.00 2020.4 0.257E+02 9.29 95.20 10.00 0.71 32685E+04 999.53 10.00 10.00 2062.6 0.251E+02 9.43 95.82 10.00 0.57 33050E+04 1010.47 10.00 10.00 2105.4 0.246E+02 9.57 96.45 10.00 0.43 33414E+04 1021.40 10.00 10.00 2148.8 0.241E+02 9.71 97.07 10.00 0.29 33779E+04 1032.34 10.00 10.00 2192.6 0.236E+02 9.85 97.69 10.00 0.15 34144E+04 1043.28 10.00 10.00 2237.0 0.232E+02 10.00 98.30 10.00 0.00 34508E+04 Cumulative travel time = 3450.8184 sec END OF MOD141: BUOYANT AMBIENT SPREADING -Due to the attachment or proximity of the plume to the bottom, the bottom coordinate for the FAR-FIELD differs from the ambient depth, ZFB = m In a subsequent analysis set "depth at discharge" equal to "ambient depth" -BEGIN MOD162: PASSIVE AMBIENT MIXING IN STRATIFIED AMBIENT Vertical diffusivity (initial value) = 0.513E-04 m^2/s Horizontal diffusivity (initial value) = 0.681E+00 m^2/s Profile definitions: BV = Gaussian s.d.*sqrt(pi/2) (46%) thickness, measured vertically = or equal to layer depth, if fully mixed BH = Gaussian s.d.*sqrt(pi/2) (46%) half-width, measured horizontally in Y-direction ZU = upper plume boundary (Z-coordinate) ZL = lower plume boundary (Z-coordinate) S = hydrodynamic centerline dilution C = centerline concentration (includes reaction effects, if any) TT = Cumulative travel time Plume Stage (bank attached): X Y Z S C BV BH ZU ZL TT 1043.28 10.00 10.00 2237.0 0.232E+02 10.00 98.30 10.00 0.00 34508E+04 1055.24 10.00 10.00 2245.1 0.231E+02 10.00 98.74 10.00 0.00 34907E+04 1067.20 10.00 10.00 2253.3 0.230E+02 10.00 99.17 10.00 0.00 35305E+04 1079.16 10.00 10.00 2261.5 0.229E+02 10.00 99.61 10.00 0.00 35704E+04 Plume interacts with BOTTOM The passive diffusion plume becomes VERTICALLY FULLY MIXED within this prediction interval 1091.12 10.00 10.00 2269.7 0.228E+02 10.00 100.04 10.00 0.00 36103E+04 1103.08 10.00 10.00 2277.8 0.227E+02 10.00 100.48 10.00 0.00 36501E+04 1115.04 10.00 10.00 2286.1 0.226E+02 10.00 100.91 10.00 0.00 36900E+04 1127.00 10.00 10.00 2294.3 0.225E+02 10.00 101.35 10.00 0.00 37299E+04 1138.95 10.00 10.00 2302.6 0.225E+02 10.00 101.79 10.00 0.00 37697E+04 1150.91 10.00 10.00 2310.9 0.224E+02 10.00 102.23 10.00 0.00 38096E+04 1162.87 10.00 10.00 2319.3 0.223E+02 10.00 102.67 10.00 0.00 38494E+04 1174.83 10.00 10.00 2327.7 0.222E+02 10.00 103.11 10.00 0.00 38893E+04 1186.79 10.00 10.00 2336.1 0.221E+02 10.00 103.55 10.00 0.00 39292E+04 1198.75 10.00 10.00 2344.6 0.220E+02 10.00 103.99 10.00 0.00 39690E+04 1210.71 10.00 10.00 2353.1 0.219E+02 10.00 104.43 10.00 0.00 40089E+04 1222.67 10.00 10.00 2361.6 0.219E+02 10.00 104.88 10.00 0.00 40488E+04 1234.63 10.00 10.00 2370.2 0.218E+02 10.00 105.32 10.00 0.00 40886E+04 1246.59 10.00 10.00 2378.8 0.217E+02 10.00 105.76 10.00 0.00 41285E+04 1258.54 10.00 10.00 2387.5 0.216E+02 10.00 106.21 10.00 0.00 41684E+04 1270.50 10.00 10.00 2396.2 0.215E+02 10.00 106.65 10.00 0.00 42082E+04 1282.46 10.00 10.00 2404.9 0.214E+02 10.00 107.10 10.00 0.00 42481E+04 1294.42 10.00 10.00 2413.7 0.214E+02 10.00 107.54 10.00 0.00 42879E+04 1306.38 10.00 10.00 2422.4 0.213E+02 10.00 107.99 10.00 0.00 43278E+04 1318.34 10.00 10.00 2431.3 0.212E+02 10.00 108.44 10.00 0.00 43677E+04 1330.30 10.00 10.00 2440.1 0.211E+02 10.00 108.89 10.00 0.00 44075E+04 1342.26 10.00 10.00 2449.0 0.210E+02 10.00 109.34 10.00 0.00 44474E+04 1354.22 10.00 10.00 2458.0 0.209E+02 10.00 109.79 10.00 0.00 44873E+04 1366.18 10.00 10.00 2467.0 0.209E+02 10.00 110.24 10.00 0.00 45271E+04 1378.13 10.00 10.00 2476.0 0.208E+02 10.00 110.69 10.00 0.00 45670E+04 1390.09 10.00 10.00 2485.0 0.207E+02 10.00 111.14 10.00 0.00 46069E+04 1402.05 10.00 10.00 2494.1 0.206E+02 10.00 111.59 10.00 0.00 46467E+04 1414.01 10.00 10.00 2503.2 0.205E+02 10.00 112.04 10.00 0.00 46866E+04 1425.97 10.00 10.00 2512.4 0.205E+02 10.00 112.49 10.00 0.00 47264E+04 1428.26 10.00 10.00 2514.1 0.204E+02 10.00 112.58 10.00 0.00 Cumulative travel time = 4734.0815 sec CORMIX prediction has been TERMINATED at last prediction interval Limiting distance due to TIDAL REVERSAL has been reached END OF MOD162: PASSIVE AMBIENT MIXING IN STRATIFIED AMBIENT -CORMIX1: Single Port Discharges End of Prediction File 1111111111111111111111111111111111111111111111 PHỤ LỤC 4: Mô 1: Mô 2: Mô 3: Mô 4: ... suốt nhũng năm học vừa qua Em xin gửi lời cảm ơn đến KS Nguyễn Thái Hòa, KS Lê Thị Hiền Th.S Cao Duy Trường phịng Tin học Mơi trường, Viện Môi Trường Tài Nguyên, Đại học Quốc gia TP.HCM t ận tình... nên người Một lần em xin chân thành cảm ơn tất Tp.Hồ Chí Minh, ngày 16 tháng 12 năm 2009 Nguyễn Duy Hiếu TRƯỜNG ĐẠI HỌC TÔN ĐỨC THẮNG KHOA MÔI TRƯỜNG VÀ BẢO HỘ LAO ĐỘNG Đề tài khóa luận tốt nghiệp:... MƠ HÌNH CORMIX MƠ PHỎNG Ơ NHIỄM NƯỚC SƠNG THỊ VẢI – XÉT NGUỒN THẢI VEDAN Sinh viên thực : Nguyễn Duy Hiếu Giáo viên hướng dẫn : Th.S Trần Thị Vân NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN