ỨNG DỤNG MÔ HÌNH MIKE 11 TÍNH TOÁN THỦY LỰC, CHẤT LƯỢNG NƯỚC CHO LƯU VỰC SÔNG SÀI GÒN- ĐỒNG NAI Trần Hồng Thái, Hoàng Thị Thu Trang, Nguyễn Văn Thao, Lê Vũ Việt phong Viện Khoa học Khí tượng Thủy văn Môi trường Đặt vấn đề Trong năm gần đây, với phát triển công nghệ thông tin khoa học kỹ thuật nói chung, mô hình toán ứng dụng ngày phát triển nhiều Các mô hình toán với ưu điểm cho kết tính toán nhanh, giá thành rẻ, dễ dàng thay đổi kịch toán,vv trở thành công cụ mạnh, phục vụ đắc lực nhiều lĩnh vực, có lĩnh vực quản lý tài nguyên môi trường Lựa chọn mô hình khâu quan trọng phương pháp mô hình toán, phụ thuộc vào yêu cầu công việc, điều kiện tài liệu tiềm tài nguồn nhân lực sẵn có Trên giới có nhiều mô hình toán sử dụng Trong nghiên cứu này, với mục tiêu mô tính toán chế độ thủy văn, thủy lực chất lượng nước cho hệ thống sông Sài Gòn – Đồng Nai, lựa chọn áp dụng phần mềm MIKE 11, đáp ứng tiêu chí sau: - Là phần mềm tích hợp đa tính năng; - Là phần mềm kiểm nghiệm thực tế; - Cho phép tính toán thủy lực chất lượng nước với độ xác cao; - Giao diện thân thiện, dễ sử dụng; - Có ứng dụng kỹ thuật GIS, kỹ thuật với tính hiệu cao Giới thiệu mô hình MIKE 11 Mô hình MIKE 11 phần mềm kỹ thuật chuyên dụng Viện Thuỷ lực Đan Mạch (DHI) xây dựng phát triển khoảng 20 năm trở lại đây, ứng dụng để mô chế độ thủy lực, chất lượng nước vận chuyển bùn cát vùng cửa sông, sông, hệ thống tưới, kênh dẫn hệ thống dẫn nước khác MIKE 11 bao gồm nhiều mô đun có khả nhiệm vụ khác như: mô đun mưa dòng chảy (RR), mô đun thuỷ động lực (HD), mô đun tải - khuếch tán (AD), mô đun sinh thái (Ecolab) số mô đun khác Trong đó, mô đun thuỷ lực (HD) coi phần trung tâm mô hình, nhiên tuỳ theo mục đích tính toán mà kết hợp sử dụng với mô đun khác cách hợp lý khoa học Trong nghiên cứu này, sử dụng mô đun HD, AD Ecolab 2.1 Mô đun thủy lực HD Mô đun thủy lực xây dựng sở hệ phương trình Saint Venant chiều cho trường hợp dòng không ổn định, gồm hai phương trình sau: 342 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT Phương trình liên tục ∂Q ∂A + =q ∂x ∂t (1) Phương trình động lượng ⎛ Q2 ⎞ ∂ ⎜α ⎟ (2) A⎠ ∂Q ∂h gQ Q ⎝ = + gA + =0 ∂t ∂x ∂x C AR Là hệ hai phương trình vi phân đạo hàm riêng phi tuyến bậc nhất, hệ phương trình có dạng nói chung không giải phương pháp giải tích Do đó, người ta giải hệ phương trình phương pháp số với lược đồ sai phân hữu hạn điểm ẩn (Abbott-Inoescu) 2.2 Các mô đun chất lượng nước Để giải vấn đề chất lượng nước có liên quan đến phản ứng sinh hóa, mô hình MIKE 11 sử dụng đồng thời hai mô đun mô đun tải - khuếch tán (AD) mô đun sinh thái (Ecolab) tính toán a Mô đun truyền tải khuếch tán Mô đun truyền tải khuếch tán dùng để mô vận chuyển chiều chất huyền phù hoà tan (phân huỷ) lòng dẫn hở dựa phương trình để trữ tích luỹ với giả thiết chất hoà tan trộn lẫn Quá trình biểu diễn qua phương trình sau: ∂AC ∂QC ∂ ⎛ ∂C ⎞ + − ⎜ AD ⎟ = − AKC + C2 q ∂t ∂x ∂x ⎝ ∂x ⎠ (3) Trong đó, hệ số phân huỷ sinh học K dùng tượng hay trình xem xét có liên quan đến phản ứng sinh hoá Phương trình () thể hai chế truyền tải, truyền tải đối lưu tác dụng dòng chảy truyền tải khuếch tán Gradien nồng độ gây Phương trình giải theo phương pháp số với sơ đồ sai phân ẩn trung tâm b Mô đun sinh thái (Ecolab) Mô đun sinh thái mô hình MIKE 11 giải khía cạnh chất lượng nước sông vùng bị ảnh hưởng hoạt động dân sinh kinh tế.v.v Mô đun phải kèm với mô đun tải - khuếch tán, điều có nghĩa mô đun chất lượng nước giải trình biến đổi sinh học hợp chất sông mô đun tải - khuếch tán dùng để mô trình truyền tải khuếch tán hợp chất Phương pháp nghiên cứu Phương pháp nghiên cứu sử dụng nghiên cứu bao gồm phương pháp thu thập, phân tích tài liệu từ nguồn có, phương pháp khảo sát thực địa, lấy mẫu phân tích thực trạng chất lượng nước vị trí quan trắc chất lượng nước dọc sông, nhằm đánh giá trạng môi trường khu vực nghiên cứu thu thập Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 343 tài liệu cần thiết cho tính toán Tiếp đó, ứng dụng phương pháp mô hình toán để tính toán, mô chế độ thủy văn thủy lực diễn biến chất lượng nước hệ thống sông Sài Gòn – Đồng Nai Giới thiệu lưu vực nghiên cứu Lưu vực sông Đồng Nai lưu vực sông lớn Việt Nam giữ vai trò vô quan trọng phát triển kinh tế - xã hội đất nước Lưu vực nằm vùng khí hậu nhiệt đới gió mùa, có lưu vực tích thuỷ từ vùng cao nguyên Tây Nguyên đến hết đồng miền Đông Nam Bộ, giới hạn toạ độ từ 105o30'21'' đến 109o01'20" kinh độ Đông từ 10o19'55" đến 12o20'38" vĩ độ Bắc Tổng diện tích tự nhiên khoảng 43.450 km2 Lưu vực có hình nan quạt, kéo dài từ cuối sườn Tây dãy Trường Sơn thuộc Nam Trung Bộ, qua hết vùng Đông Nam Bộ đến giáp vùng Đồng Tháp Mười thuộc Đồng sông Cửu Long Dòng sông Đồng Nai chảy theo hướng Đông Bắc - Tây Nam, nhánh sông lớn quan trọng đổ vào dòng sông La Ngà (nằm bên trái dòng theo hướng từ thượng nguồn cửa sông), sông Bé, sông Sài Gòn sông Vàm Cỏ (nằm bên phải) Toàn hệ thống sông suối lưu vực tập trung cửa Gành Rái Xoài Rạp Trong xu phát triển kinh tế - xã hội nước, khu vực nảy sinh hàng loạt vấn đề ô nhiễm môi trường, đặc biệt ô nhiễm nguồn nước hai sông quan trọng sông Đồng Nai sông Sài Gòn Các khu đô thị, khu dân cư khu công nghiệp tập trung hình thành phát triển mạnh mẽ dọc theo sông từ thượng nguồn (Lâm Đồng, Bình Thuận, Bình Phước, Tây Ninh) đến trung lưu (Đồng Nai, Bình Dương) hạ lưu (Thành phố Hồ Chí Minh, Bà Rịa - Vũng Tàu) đã, nguồn gây ô nhiễm nước cho sông Các tài liệu sử dụng Với mục tiêu tính toán chất lượng nước, tài liệu, số liệu phục vụ cho mô hình tính toán cần nhiều, bao gồm số liệu thuỷ văn, thuỷ lực chất lượng nước, cụ thể sau: 5.1 Tài liệu mặt cắt Toàn hệ thống sông khu vực nghiên cứu gồm có sông với 52 mặt cắt, đó: Sông Đồng Nai: 26 mặt cắt ( Từ trạm Hoá An - trạm Nhà Bè) Sông Sài Gòn: 16 mặt cắt (Từ trạm Bình Phước – nam Cát Lái) Sông Nhà Bè: mặt cắt (Từ trạm Nhà Bè - trạm Vàm Sát) Sông Lòng Tàu: mặt cắt (Từ Phú Xuân – Nhà Bè - trạm Tam Thôn Hiệp) 344 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT Số liệu mặt cắt dùng tính toán chủ yếu Tổng cục Khí tượng Thuỷ văn cũ (nay Bộ Tài nguyên Môi trường) đo đạc năm 2003 Nói chung số liệu mặt cắt sử dụng mô hình có độ tin cậy tương đối cao Sơ đồ hóa mạng lưới sông lưu vực thể Hình Hình Sơ đồ mặt cắt hệ thống sông Sài Gòn - Đồng Nai 5.2 Tài liệu khí tượng thủy văn Tài liệu khí tượng thủy văn bao gồm tài liệu mực nước lưu lượng, đo đặc theo năm 2005, dùng để làm điều kiện biên (Q~t ), biên (H~t) để hiệu chỉnh kiểm định mô hình Danh sách trạm thủy văn sử dụng liệt kê bảng Bảng Danh sách trạm thuỷ văn sử dụng STT Tên trạm Sông Trạm thuỷ văn Hoá An Đồng Nai Trạm thuỷ văn Cát Lái Đồng Nai Trạm thuỷ văn Nhà Bè Đồng Nai Trạm thuỷ văn Bình Phước Sài Gòn Trạm thuỷ văn Phú An Sài Gòn Trạm thuỷ văn Tam Thôn Hiệp Trạm thuỷ văn Vàm Sát Lòng Tàu Nhà Bè 5.3 Tài liệu chất lượng nước − − Số liệu chất lượng nước mặt quan trắc điểm lấy mẫu vào năm 2003, 2004, 2005 tháng đầu năm 2006, số liệu quan trắc năm 2005 dùng làm điều kiện biên để hiệu chỉnh kiểm định mô hình, chuỗi số liệu từ tháng 10/2005 đến tháng 12/2005 sử dụng để hiệu chỉnh mô hình chất lượng nước Các tiêu chất lượng nước tính toán thu thập bao gồm: DO, BOD, TN, TP, tổng Coliforms Tài liệu nguồn gây ô nhiễm: Trong trình phát triển, công nghiệp hóa, đại hóa, địa phương lưu vực sông Sài Gòn – Đồng Nai Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 345 tiếp tục đối mặt với vấn đề ô nhiễm nguồn nước với xu hướng ngày gia tăng, đặc biệt khu vực hạ lưu hệ thống sông Các nguồn thải gây ô nhiễm hệ thống sông Sài Gòn – Đồng Nai nhận diện bao gồm: + Nguồn thải từ khu đô thị + Nguồn thải từ khu công nghiệp tập trung + Nguồn thải từ sở công nghiệp phân tán + Nguồn thải từ bãi rác Các số liệu thuỷ văn, thuỷ lực trước tiên sử dụng để diễn toán chế độ thuỷ lực hệ thống sông nghiên cứu, làm đầu vào cho toán chất lượng nước, sau kết hợp với số liệu chất lượng nước để tính toán mô chất lượng nước sông Kết ứng dụng mô hình MIKE 11 6.1 Kết tính toán thủy lực hệ thống sông Sài Gòn - Đồng Nai Mục tiêu tính toán thủy lực để xác định thông số thuỷ lực phù hợp lưu vực nghiên cứu, phục vụ cho mô hình chất lượng nước Căn vào nhiệm vụ tính toán, tài liệu mặt cắt tài liệu thu thập được, sơ đồ tính toán thủy lực cho mạng sông Sài Gòn – Đồng Nai thu hẹp lại Hình lý sau: − Do thiếu số liệu mặt cắt vùng thượng lưu sông Sài Gòn − Do thiếu số liệu lưu lượng, mực nước thực đo mặt cắt vùng Cửa Soài Rạp Lòng Tàu – Ngã Bảy a Các bước ứng dụng Cài đặt mô hình Các Thiết lập điều kiện biên Xây dựng điều kiện ban đầu Hiệu chỉnh mô hình Kiểm định mô hình Biên toán xác định sau: − Biên trình lưu lượng thực đo từ 01/01/2005 đến 31/12/2005 trạm Bình Phước, Hóa An − Biên trình mực nước thực đo từ 01/01/2005 đến 31/12/2005 trạm Vàm Sát, Tam Thôn Hiệp − Hiệu chỉnh, xác định thông số mô hình Việc hiệu chỉnh thông số mô hình thủy lực thực chủ yếu qua việc thay đổi hệ số nhám Manning, bước thời gian tính toán giá trị ban đầu Sau hiệu chỉnh thông số, mô hình đảm bảo độ xác cần thiết với bước thời gian tính toán giờ, hệ số nhám khoảng 0.02-0.03 Kết hiệu chỉnh mô hình thuỷ lực trạm Cát Lái, Phú An Nhà Bè thể hình 2, 3, sau: 346 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT Phu An Cat Lai 1.0 1.4 H CatLai Obs [meter] H Cat Lai Cal [meter] 0.9 1.3 H Phu An Obs [meter] H Phu An Cal [meter] 1.2 0.8 0.7 1.1 0.6 1.0 0.5 0.9 0.4 0.8 0.3 0.7 0.2 0.6 0.1 0.5 0.0 0.4 -0.1 0.3 -0.2 0.2 -0.3 0.1 -0.4 -0.5 0.0 -0.6 -0.1 -0.7 -0.2 -0.8 -0.3 -0.9 -0.4 -1.0 -0.5 -1.1 -0.6 -1.2 -0.7 -1.3 -0.8 -1.4 -0.9 -1.5 -1.0 -1.6 -1.1 -1.7 -1.8 -1.2 -1.9 -1.3 -2.0 09:00 2005-01-15 09:00 01-16 09:00 01-17 09:00 01-20 09:00 01-19 09:00 01-18 09:00 01-22 09:00 01-21 09:00 01-23 09:00 2005-03-02 09:00 01-24 Hình So sánh mực nước kết tính toán với số liệu thực đo trạm Cát Lái 09:00 03-04 09:00 03-03 09:00 03-05 09:00 03-06 09:00 03-07 09:00 03-08 09:00 03-09 Hình So sánh mực nước kết tính toán với số liệu thực đo trạm Phú An Phu An 1.4 H Phu An Obs [meter] H Phu An Cal [meter] 1.3 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 -0.5 -0.6 -0.7 -0.8 -0.9 -1.0 -1.1 -1.2 -1.3 09:00 2005-03-02 09:00 03-03 09:00 03-04 09:00 03-05 09:00 03-10 09:00 03-06 09:00 03-07 09:00 03-08 09:00 03-09 09:00 03-10 Hình So sánh mực nước kết tính toán với số liệu thực đo trạm Nhà Bè Sai số số liệu thực đo tính toán trạm đánh giá theo số Nash-Sutcliffe, cụ thể bảng Bảng Phân tích hiệu sai số hiệu chỉnh mô hình TT Trạm kiểm định Mức hiệu mô hình (%) Cát Lái Nhà Bè Phú An 90% 94% 89% b Kiểm định mô hình Mục đích công tác kiểm định mô hình nhằm đánh giá mức độ phù hợp thông số xác định phần hiệu chỉnh mô hình Ở sử dụng số liệu thực mực nước thực đo mùa kiệt từ tháng 1/2004 đến tháng 12/2004 để kiểm định Kết kiểm định mực nước trạm Nhà Bè trình bày Hình với mức độ sai lệch đạt 92% Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 347 [meter] Time Series Water Level 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 -0.2 -0.4 -0.6 -0.8 -1.0 -1.2 -1.4 -1.6 -1.8 -2.0 -2.2 4-2-2004 6-2-2004 8-2-2004 10-2-2004 12-2-2004 14-2-2004 Hình So sánh mực nước kết tính toán với số liệu thực đo trạm Nhà Bè năm 2004 6.2 Ứng dụng mô hình MIKE 11 tính toán chất lượng nước lưu vực sông Sài Gòn - Đồng Nai a Mục tiêu Mục tiêu để xác định thông số chất lượng nước phù hợp cho khu vực nghiên cứu Do hạn chế thời gian tài liệu chất lượng nước nên nghiên cứu tập trung vào tính toán số tiêu chất lượng nước quan tâm sông theo thời gian không gian DO, BOD, tổng Nitơ, tổng Photpho, tổng Coliforms tương ứng với điều kiện biên thủy lực nguồn thải b Thiết lập toán Dựa kết mô hình thủy lực tính toán phần Bài toán chất lượng nước cho lưu vực sông Sài Gòn – Đồng Nai xây dựng tiếp tục sau: − Thiết lập điều kiện biên chất lượng nước nguồn thải đổ vào sông: (i) thông số chất lượng nước biên (các vị trí quan trắc), (ii) thông số chất lượng nước biên dưới, (iii) vị trí nồng độ, lưu lượng nguồn xả thải; − Hiệu chỉnh mô hình để xác định thông số phù hợp; − Sử dụng mô hình với thông số xác định để mô thủy lực chất lượng nước theo không gian, thời gian với điều kiện biên thích hợp c Sơ đồ tính toán Lưới tính toán mô hình chất lượng nước Hình thu gọn lại lý sau: - Chủ yếu tập trung tính toán chất lượng nước sông Sài Gòn Đồng Nai - Mạng lưới quan trắc nước mặt chủ yếu tập trung sông Sài Gòn sông Đồng Nai - 348 Số liệu đầu vào cho mô hình chất lượng nước bao gồm: Các điều kiện thủy lực biên biên hệ thống sông sau thu gọn sơ đồ tính; Các điều kiện biên: + Điều kiện biên nồng độ chất lượng nước gồm tiêu chất lượng nước đo từ 01/10/2005 đến 31/12/2005 trạm Bình Phước, Phú An + Điều kiện biên nồng độ chất lượng nước gồm tiêu chất lượng nước đo từ 01/10/2005 đến 31/12/2005 trạm Tam Thôn Hiệp, Vàm Sát Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT - Các nguồn ô nhiễm: (i) vị trí nguồn thải đổ sông, (ii) nồng độ chất nước thải (iii) lưu lượng thải - Số liệu chất lượng nước số điểm quan trắc dọc sông sử dụng cho hiệu chỉnh mô hình Các số liệu đo đạc đồng với số liệu dùng làm biên biên mô hình thuỷ lực mô hình chất lượng nước d Hiệu chỉnh mô hình Kết hiệu chỉnh mô hình chất lượng nước trạm Phú An Nhà Bè với hai tiêu DO BOD mô tả hình 7, 8, 9,10 So sánh kết hiệu chỉnh nồng độ DO với số liệu thực đo trạm Nhà Bè năm 2005 So sánh kết hiệu chỉnh nồng độ DO với số liệu thực đo trạm Phú An năm 2005 DO Phú An Cal DO Phú An Obs DO Nhà Bè Cal TCVN 5942-B 4.5 3.5 TCVN 5942-B mg/l mg/l DO Nhà Bè Obs 2.5 1.5 0.5 0 1/10 8/10 15/10 22/10 1/11 8/11 15/11 22/11 1/12 8/12 15/12 1/10 8/10 15/10 22/10 1/11 8/11 15/11 22/11 1/12 8/12 15/12 Hình So sánh nồng độ DO kết tính Hình So sánh nồng độ DO kết tính toán với số liệu thực đo trạm Nhà Bè toán với số liệu thực đo trạm Phú An So sánh kết hiệu chỉnh nồng độ BOD với số liệu thực đo trạm Nhà Bè năm 2005 So sánh kết hiệu chỉnh nồng độ BOD với số liệu thực đo trạm Phú An năm 2005 BOD Phú An Obs BOD Nhà Bè Cal TCVN 5942-B 10 mg/l mg/l BOD Phú An Cal 12 BOD Nhà Bè Obs TCVN 5942-B 2 0 1/10 15/10 1/11 15/11 1/12 1/10 15/12 15/10 1/11 15/11 1/12 15/12 Hình So sánh nồng độ BOD kết tính Hình 10 So sánh nồng độ BOD kết tính toán với số liệu thực đo trạm Nhà Bè toán với số liệu thực đo trạm Phú An e Kiểm định mô hình Kiểm định mô hình nhằm đánh giá mức độ phù hợp thông số xác định phần hiệu chỉnh mô hình Số liệu thực đo từ tháng 1/2006 đến tháng 3/2006 dùng để kiểm định Kết kiểm định mô hình trạm Phú An Nhà Bè tiêu DO BOD trình bày tương ứng Hình 10 11 Kết kiểm định mô hình chất lượng nước 2006 Kết kiểm định mô hình chất lượng nước 2006 DO Cal 8-2006 BOD Cal 8-2006 DO Obs 8-2006 BOD Obs 8-2006 7.2 6.8 6.6 6.4 mg/l mg/l 6.2 5.8 5.6 5.4 5.2 Nhà Bè Phú An Hình 10 Kết kiểm định DO Nhà Bè Phú An Hình 11 Kết kiểm định BOD Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 349 Theo kết tính toán biểu diễn trên, ta thấy sai số kết tính toán với kết thực đo nhỏ 20%, mô hình chất lượng nước đạt độ xác yêu cầu áp dụng để dự báo chất lượng nước lưu vực tương lai Kết luận Trong nghiên cứu này, phương pháp mô hình toán sử dụng việc mô chế độ thủy văn, thủy lực chất lượng nước cho hệ thống sông Sài Gòn – Đồng Nai – lưu vực sông lớn giữ vai trò quan trọng phát triển kinh tế đất nước Kết tính toán, mô thủy văn, thủy lực chất lượng nước mô hình MIKE 11 tốt, cho thấy khả ứng dụng hiệu mô hình Tuy nhiên để sử dụng mô hình tốt tương lai, đòi hỏi số liệu quan trắc thủy văn, thủy lực chất lượng nước cần đồng bộ, dày đặc xác Tài liệu tham khảo Ven T Chow, Larry W Mays, David R Maidment – Applied Hydrology – Mc 350 Graw - Hill International Editions - 1988 PGS.TS Đặng Văn Bảng - Đại học Thuỷ Lợi – Mô hình toán thuỷ văn - Hà Nội, 2005 Lê Trình, Lê Quốc Hùng - Môi trường lưu vực sông Đồng Nai-Sài Gòn Nhà xuất khoa học kỹ thuật, Hà Nội 2004 Trần Hồng Thái nnk - Báo cáo tổng kết dự án “Ứng dụng mô hình tính toán dự báo ô nhiễm môi trường nước cho lưu vực sông: Cầu, Nhuệ – Đáy, Sài Gòn – Đồng Nai” – Viện Khoa h ọc Khí tượng thủy văn Môi trường, Hà Nội - 2006 Viện Khoa học Thủy lợi Miền nam - Nghiên cứu đề xuất giải pháp khoa học công nghệ để ổn định lòng dẫn hạ du hệ thống sông Đồng Nai – Sài Gòn phục vụ phát triển kinh tế - xã hội vùng Đông Nam Bộ TP Hồ Chí Minh, 01-2005 DHI software - MIKE software 2004 User Guide DHI software - MIKE 11 Reference Manual – 2004 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT ... mg/l DO Nhà Bè Obs 2.5 1.5 0.5 0 1/10 8/10 15/10 22/10 1 /11 8 /11 15 /11 22 /11 1/12 8/12 15/12 1/10 8/10 15/10 22/10 1 /11 8 /11 15 /11 22 /11 1/12 8/12 15/12 Hình So sánh nồng độ DO kết tính Hình... 5942-B 10 mg/l mg/l BOD Phú An Cal 12 BOD Nhà Bè Obs TCVN 5942-B 2 0 1/10 15/10 1 /11 15 /11 1/12 1/10 15/12 15/10 1 /11 15 /11 1/12 15/12 Hình So sánh nồng độ BOD kết tính Hình 10 So sánh nồng độ BOD... hệ thống sông Đồng Nai – Sài Gòn phục vụ phát triển kinh tế - xã hội vùng Đông Nam Bộ TP Hồ Chí Minh, 01-2005 DHI software - MIKE software 2004 User Guide DHI software - MIKE 11 Reference Manual