NGHIÊN CỨU ÁP DỤNG MÔ HÌNH TOÁN MIKE 11 TÍNH TOÁN CHẤT LƯỢNG NƯỚC SÔNG NHUỆ - SÔNG ĐÁY

NGHIÊN CỨU ÁP DỤNG MÔ HÌNH TOÁN MIKE 11 TÍNH TOÁN CHẤT LƯỢNG NƯỚC SÔNG NHUỆ - SÔNG ĐÁY Lê Vũ Việt Phong, Trần Hồng Thái, Phạm Văn Hải Viện Khoa học Khí tượng Thủy văn và Môi trường 1. Giới thiệu Trong những năm gần đây, với sự tăng trưởng của dân số và các ngành kinh tế xã hội, lưu vực sông Nhuệ - Đáy đang đứng trước nhiều nguy cơ và thách thức, nổi bật trong số đó là vấn đề chất lượng môi trường nước của hai con sông huyết mạch trong lưu vực, sông Nhuệ và sông Đáy. Sông Đáy nguyên là một phân lưu lớn của sông H ồng, có chiều dài khoảng 247 km, bắt đầu từ cửa Hát Môn cho đến cửa Đáy trước khi đổ ra biển Đông. Do ảnh hưởng của đập Đáy nên đoạn thượng nguồn (từ sau đập Đáy đến Ba Thá dài 71 km) của sông Đáy gần như là một sông chết. Lượng nước chủ yếu cung cấp cho sông Đáy lấy từ các sông nhánh đổ vào, trong đó quan trọng nhất là sông Tích, sông Bôi, sông Đào và sông Nhuệ. Sông Nhuệ lấ y nước từ sông Hồng qua cống Liên Mạc, phục vụ nước tưới cho hệ thống thủy nông Đan Hoài. Bên cạnh đó, sông Nhuệ còn có nhiệm vụ tiêu nước cho thành phố Hà Nội, thị xã Hà Đông trước khi đổ vào sông Đáy tại thị xã Phủ Lý. Đây cũng là hai nguồn gây ô nhiễm chủ yếu cho sông Nhuệ, đặc biệt là nước thải từ thành phố Hà Nội qua đập Thanh Liệt. Lưu lượng nước trung bình trên sông Tô Lị ch qua đập Thanh Liệt đổ vào sông Nhuệ từ 11-17 (m 3 /s), cực đại là 30 (m 3 /s). Những nghiên cứu đo đạc gần đây cho thấy chất lượng nước tại một số điểm quan trắc trên sông Nhuệ và sông Đáy đều đang ở mức ô nhiễm nặng nề, thậm chí còn có những nơi chất lượng nước không đạt tiêu chuẩn nước mặt loại B (TCVN 5945 - 1995). Số liệu quan trắc nồng độ oxy hòa tan (Dissolved Oxygen - DO) và nhu cầu oxy sinh hóa (Biological Oxygen Demand - BOD) tại các điểm đo trong hệ thống lưu vực sông Nhuệ-Đáy được thu thập và cung cấp bởi Cục Bảo vệ Môi trường - Bộ Tài nguyên và Môi trường (Bảng 1 và Bảng 2). Bảng 1. Nồng độ Oxy hòa tan DO (mg/l) và nhu cầu ôxy sinh hoá BOD (mg/l) tại các điểm quan trắc trên sông Nhuệ tháng 11/2005 Trạm Liên Mạc Phúc La Cự Đà Đồng Quan Cống Thần Đò Kiều DO 6.3 1.64 0.99 1.77 1.4 1.26 BOD 22 82 95 85 75 80 (Nguồn: Cục Bảo vệ Môi trường) Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 269 Bảng 2. Nồng độ Oxy hòa tan DO (mg/l) và nhu cầu ôxy sinh hoá BOD (mg/l) tại các điểm quan trắc trên sông Đáy tháng 11/2005 Trạm Ba Thá Tế Tiêu Cầu Quế Hồng Phú Trung Hiếu Hạ Khánh Phú Độc Bộ Đò Mười DO 3.76 1.37 1.43 0.92 1.6 4.22 4.27 4.64 BOD 35 40 35 45 35 35 30 30 (Nguồn: Cục Bảo vệ Môi trường) Dựa trên các số liệu quan trắc cho thấy sông Nhuệ và sông Đáy đang bị ô nhiễm nghiêm trọng, đặc biệt là sông Nhuệ. Sự trong sạch và khỏe mạnh của dòng sông có ảnh hưởng rất lớn đến sinh hoạt và sản xuất của cả một vùng dân cư rộng lớn, chính vì vậy, việc nghiên cứu áp dụng các biện pháp khoa học góp phần cải thiện tình trạng ô nhiễm tại các dòng sông trên là vô cùng cấ p thiết. Ở đây, trong bài báo này, nhóm các tác giả đã nghiên cứu ứng dụng mô hình MIKE 11 trong tính toán chất lượng nước (DO, BOD và nhiệt độ) sông Nhuệ và sông Đáy bởi đây là mô hình có tích hợp nhiều các tính năng ưu việt và đồng bộ, thuận tiện cho việc nghiên cứu. Hình 1. Bản đồ lưu vực sông Nhuệ - Đáy 2. Nội dung nghiên cứu - cơ sở lý thuyết MIKE 11 là một phần mềm kỹ thuật chuyên dụng, được ứng dụng để mô phỏng lưu lượng, chất lượng nước và vận chuyển bùn cát ở cửa sông, sông, hệ thống tưới, kênh dẫn và các hệ thống dẫn nước khác. Mô hình giải quyết những bài toán liên quan thông qua các mô đun độc lập được tích hợp lại với nhau để thực hiện các quy trình tính toán c ần thiết. Trong nghiên cứu này, các mô đun được ứng dụng bao gồm:  Mô đun mưa - dòng chảy (Rainfall Runoff - RR);  Mô đun thuỷ lực (Hydrodynamic - HD);  Mô đun tải khuyếch tán (Advection Dispersion - AD);  Mô đun sinh thái (Ecolab). Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 270 Bài báo này chỉ đi sâu vào phần nội dung mô phỏng và tính toán thuỷ lực, chất lượng nước lưu vực sông Nhuệ và sông Đáy. L.V.V.Phong [7] đã cung cấp chi tiết nội dung nghiên cứu ứng dụng mô hình mưa dòng chảy MIKE - NAM khôi phục số liệu dòng chảy phục vụ cho mô hình thuỷ lực. a. Cơ sở lý thuyết mô hình thủy lực MIKE 11 Mô hình thủy lực (Hydrodynamic module - HD) là phần quan trọng nhất trong bộ mô hình MIKE 11, được xây dựng trên hệ 2 phương trình Saint - Venant vi ết cho dòng một chiều không ổn định với một số giả thiết đi kèm [2, 4, 5]. - Phương trình liên tục: QA q xt ∂∂ + = ∂∂ (1) - Phương trình động lượng: 2 2 0 AR Q gQ Q A Qh gA tx xC α ⎛⎞ ∂ ⎜⎟ ∂∂ ⎝⎠ = ++ ∂∂ ∂ = (2) Trong đó: Q: Lưu lượng (m 3 /s) A: Diện tích mặt cắt (m 2 ) q: Lưu lượng nhập lưu trên một đơn vị chiều dài dọc sông (m 2 /s) C: Hệ số Chezy [2]. α: Hệ số sửa chữa động lượng. R: Bán kính thuỷ lực (m) Hiện nay hệ phương trình Saint - Venant chỉ tìm được nghiệm giải tích trong trường hợp hệ thống đơn giản và phải kèm theo rất nhiều các giả thiết. Hệ phương trình này chủ yếu được giải bằng phương pháp gần đúng, cụ thể trong mô hình MIKE 11 đó là phương pháp sai phân hữu hạn [4, 5, 6]. Tài li ệu yêu cầu cho mô hình thủy lực tính toán cho bài toán dòng chảy một chiều không ổn định trong sông tự nhiên bao gồm: - Số liệu địa hình, thể hiện bằng mặt cắt ngang sông; - Số liệu lưu lượng và mực nước theo thời gian tại các biên trên và biên dưới, các điểm đo dùng để hiệu chỉnh và kiểm nghiệm mô hình; - Số liệu lưu lượng và mực nước tại thời điểm t=0 dùng làm điều kiện ban đầu để tính toán cho mô hình. Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 271 b. Cơ sở lý thuyết mô hình chất lượng nước MIKE 11 Ngoài phần mô hình thủy lực có nhiệm vụ tính toán về chế độ thủy lực của dòng chảy trong sông, bài toán chất lượng nước được mô hình MIKE 11 giải quyết thông qua mô đun truyền tải khuyếch tán và mô đun sinh hóa. Mô đun truyền tải khuyếch tán có chức năng tính toán sự lan truyền của các chất huyền phù hoặc hòa tan (phân hủy) dưới tác động của dòng chảy. Bài toán đặt ra trong nghiên cứu này là bài toán có liên quan đến các chu trình sinh hóa diễn ra trong dòng sông, do đó các tác giả không những ch ỉ sử dụng mô đun truyền tải khuyếch tán mà còn phải kết hợp với mô đun sinh hóa để giải quyết vấn đề. Phương trình cơ bản trong hai mô đun này là phương trình truyền tải khuyếch tán [6]: - Phương trình truyền tải - khuyếch tán 2 AC QC C AD AKC C q txx x ∂∂ ∂ ∂ ⎛⎞ +− =−+ ⎜⎟ ∂∂∂ ∂ ⎝⎠ (3) Trong đó: A: Diện tích mặt cắt (m 2 ); C: Nồng độ (kg/m 3 ); D: Hệ số khuyếch tán; q : Lưu lượng nhập lưu trên 1 đơn vị chiều dài dọc sông (m 2 /s); K : Hệ số phân huỷ sinh học, K chỉ được dùng khi các hiện tượng hay quá trình xem xét có liên quan đến các phản ứng sinh hoá. Hệ số phân huỷ sinh học K bao hàm trong đó rất nhiều các hiện tượng và phản ứng sinh hoá. Hệ số này không cần xem xét trong bài toán lan truyền chất thông thường. 3. Nghiên cứu áp dụng mô hình MIKE 11 tính toán chất lượng nước sông Nhuệ và sông Đáy a. Số liệu đầu vào phục vụ cho mô hình thuỷ lực và chất lượng nước Một trong những mục tiêu của nghiên cứu này đó là mô phỏng lại toàn bộ quá trình chất lượng nước diễn ra trong lưu vực sông Nhuệ, sông Đáy từ những số liệu cần thiết thu thập được tại những điểm rời rạc trên sông. Nội dung này được thực hiện qua công tác hiệu chỉnh và kiểm nghiệm mô hình, từ đó thu được những thông số mô hình phù hợp với thực tế. Trước hết, hệ thống sông nghiên cứu tính toán trong mô hình thuỷ lực và chất lượng nước có 6 sông với 219 mặt cắt (Hình 1) bao gồm: sông Hồng, sông Nhuệ, sông Đáy, sông Đuống, sông Đào và sông Hoàng Long, trong đó: - Thượng nguồn sông Nhuệ nối với sông Hồng tại cống Liên Mạc. - Hạ lưu sông Nhuệ nối với sông Đáy tại thị xã Phủ Lý. - Hạ lưu sông Hoàng Long nối với sông Đáy tại Gián Khẩu. - Sông Đào nối với sông Đáy tại Độc Bộ. Mô hình cần có các số liệu địa hình, thuỷ văn, thuỷ lực và chất lượng nước. Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 272 i. Số liệu mặt cắt các sông nghiên cứu (219 mặt cắt) được tiến hành đo đạc vào năm 2000, riêng sông Nhuệ số liệu đo năm 1990, chi tiết như sau:: - Sông Hồng: 89 mặt cắt (Từ trạm Sơn Tây - trạm Hưng Yên) - Sông Đuống: 3 mặt cắt (Từ ngã ba sông Đuống - trạm Thượng Cát) - Sông Nhuệ: 43 mặt cắt (Từ cống Liên Mạc - Phủ Lý) - Sông Đáy: 60 mặt cắt (Từ trạm Ba Thá - trạm Như Tân) - Sông Đào: 10 mặt cắt (Từ nhập lưu với sông Đáy - trạm Nam Định) - Sông Hoàng Long: 14 mặt cắt (Từ trạm Hưng Thi - trạm Gián Khẩu) ii. Số liệu lưu lượng và mực nước: Hình 2. Sơ đồ mạng sông tính toán Số liệu lưu lượng và mực nước thực đo dùng để làm điều kiện biên trên, biên dưới và hiệu chỉnh, kiểm nghiệm mô hình thuỷ lực. Các số liệu thuỷ văn (∆t =1 ngày) từ tháng 11/1998 đến tháng 5/1999 được dùng hiểu chỉnh mô hình thuỷ lực, số liệu thuỷ văn (∆t =1 ngày) từ 11/1999 đến tháng 5/2000 được dùng để kiểm nghiệm mô hình thuỷ lực. Danh sách các trạm thuỷ văn sử d ụng trong mô hình như Bảng 3. Bảng 3. Danh sách các trạm thủy văn dùng trong tính toán mô hình thuỷ lực TT Tên trạm Sông Số liệu Mục đích 1 Sơn Tây Hồng Q, H Biên trên 2 Ba Thá (*) Đáy Q Biên trên 3 Hưng Thi (**) Bôi Q Biên trên 4 Nam Định Đào H Biên trên 5 Phủ Lý Đáy H Hiệu chỉnh, kiểm nghiệm 6 Hà Đông Nhuệ H Hiệu chỉnh, kiểm nghiệm 7 Hà Nội Hồng Q, H Hiệu chỉnh, kiểm nghiệm 8 Gián Khẩu Hoàng Long H Hiệu chỉnh, kiểm nghiệm 9 Ninh Bình Đáy H Hiệu chỉnh, kiểm nghiệm 10 Hưng Yên Hồng H Biên dưới 11 Thượng Cát Đuống Q, H Biên dưới 12 Như Tân Đáy H Biên dưới (*, **) Lưu lượng Q được xác định từ mô hình mưa dòng chảy NAM Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 273 iii. Số liệu chất lượng nước: - Số liệu quan trắc nước mặt Được dùng để làm điều kiện biên trên, biên dưới và hiệu chỉnh, kiểm nghiệm trong mô hình chất lượng nước. Số liệu chất lượng nước được quan trắc tại các điểm trên trên với tần suất 12 lần/ năm. Các số liệu chất lượng nước quan trắc (DO, BOD, nhiệt độ) từ tháng 10 đến tháng 12 năm 2005 được dùng để hiệu chỉnh và ki ểm nghiệm mô hình. Danh sách các trạm quan trắc trong Bảng 4. - Số liệu nguồn thải và nguồn gây ô nhiễm Nguồn thải là những nguồn gây ô nhiễm chính trong các sông trên lưu vực, được hình thành từ các hoạt động dân sinh kinh tế của con người. Số liệu nguồn thải được thu thập cho hai nguồn chính: 9 Số liệu nước thải công nghiệp dựa trên các báo cáo quan trắc môi trường khu công nghiệp, quy mô và tần suất hoạt động. 9 Số liệu nước thải sinh hoạt ước tính dựa trên dân số và lượng nước cấp trong từng khu vực, phân bổ theo mạng lưới kênh và sông nhánh Ngoài ra còn có một số các nguồn phụ khác, được thu thập đầy đủ nhất có thể nhưng cũng không hoàn toàn chính xác do công tác quản lý các nguồn này còn nhiều hạn chế và bất cập như nước thải trực tiếp trên sông, nước thải nông nghiệp .v.v. Bảng 4. Danh sách các trạm quan trắc chất lượng nước mặt dùng trong tính toán mô hình chất lượng nước TT Tên trạm 1 Sông Nhuệ Sông Đáy 2 Liên Mạc Ba Thá 3 Phúc La Tế Tiêu 4 Thanh Liệt Cầu Quế 5 Cự Đà Hồng Phú 6 Đồng Quan Trung Hiếu Hạ 7 Cống Thần Khánh Phú 8 Đò Kiều Độc Bộ 9 Đò Mười 10 Kim Tân b. Một số kết quả Kết quả hiệu chỉnh và kiểm nghiệm cho mô hình thủy lực được nêu trong Hình 3, Hình 4 và Bảng 5. Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 274 Bảng 5. Phân tích hiệu quả và sai số của hiệu chỉnh mô hình thủy lực TT Trạm kiểm nghiệm Mức hiệu quả của mô hình (%) Sai số tuyệt đối (m) 1 Hà Nội 97% 0.15 2 Phủ Lý 88% 0.25 3 Ninh Bình 91% 0.20 4 Hà Đông 60% 0.35 5 Gián Khẩu 95% 0.20 Kết quả hiệu chỉnh và kiểm nghiệm cho mô hình chất lượng nước được nêu trong các hình từ Hình 5 đến Hình 8 21-11- 1998 11-1 2-1998 31-12-1998 20-1-1999 9-2-1999 1-3-19 99 21-3-1999 10-4-1999 30-4-1999 20-5-1999 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 [meter] Time Series Water Level Water Level SONGDAY 100575.00 SONG HONG 26100.00 SONG HONG 42335.00 External TS 1 HLienMac HHaNo i HPhuLy 11-12- 1998 31-12-1998 20-1-1999 9-2- 1999 1-3- 1999 21-3-1999 10-4-1999 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 [meter] Time Series Water Level Water Level SONG NHUE 15 7 00 .00 HOA NGLONG 50658.00 Ex t e r n a l T S 1 HHaDo ng HGianKhau Hình 3. So sánh giữa kết quả hiệu chỉnh mô hình diễn toán MIKE11 với số liệu mực nước thực đo, trạm Hà Nội, Liên Mạc, Phủ Lý, Hà Đông và Gián Khẩu tháng 11/1998 đến tháng 5/1999 21-11-1999 11-1 2-1999 31 -12-1999 20-1-2000 9-2- 2000 29-2-2000 20-3-200 0 9-4-2000 29-4-2000 19-5-2000 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 [meter] Time Series Water Level Water Level SONGDAY 100575.00 SONG HONG 26100.00 SONG HONG 42335.00 External TS 1 HPhuLy HLienMac HHaNo i 11-12-1999 31-12-1999 20-1-2000 9-2 -2000 29-2-2000 20-3- 2000 9-4- 2000 29-4-2000 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 [meter] Time Series Water Level Water L evel SONG NHUE 1570 0. 00 HOANGLONG 50658.00 Ex ternal TS 1 HHaDo ng HGia nKh au Hình 4. So sánh giữa kết quả kiểm nghiệm mô hình diễn toán MIKE11 với số liệu mực nước thực đo, trạm Hà Nội, Liên Mạc, Phủ Lý, Hà Đông và Gián Khẩu tháng 11/1998 đến tháng 5/1999 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 275 So sánh kết quả tính toán hiệu chỉnh nồng độ DO v ới số liệu thực đo, dọc theo sông Đáy , tháng 11/2005. 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 Tế Tiêu Cầu Quế Khánh Phú Độc BộĐò Mười Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B So sánh kết quả tính toán hiệu chỉnh nồng độ DO v ới số liệu thực đo, dọc theo sông Nhuệ , tháng 11/2005. 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Phúc La Cự Đà Đò Quan Cống Thần Đò Kiều Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B Hình 5. So sánh giữa kết quả hiệu chỉnh DO với số liệu thực đo dọc sông Đáy và sông Nhuệ tháng 11/2005 So sánh kết quả tính toán hiệu chỉnh nồng độ BOD với số liệu thực đo, dọc sông Đáy , tháng 11/2005. 0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 Tế Tiêu Trung Hiếu Hạ Khánh Phú Độc BộĐò Mười Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B So sánh kết quả tính toán hiệu chỉnh nồng độ BOD với số liệu thực đo, dọc sông Nhuệ , tháng 11/2005. 0.0 20.0 40.0 60.0 80.0 100.0 120.0 Phúc La Cự Đà Đò Quan Cống Thần Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B Hình 6. So sánh giữa kết quả hiệu chỉnh BOD với số liệu thực đo dọc sông Đáy và sông Nhuệ tháng 11/2005 So sánh kết quả tính toán hiệu chỉnh nồng độ DO v ới số liệu thực đo, dọc theo sông Đáy, tháng 12/2005. 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 Tế Tiêu Cầu Quế Khánh Phú Độc BộĐò Mười Vị trí quan tr ắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B So sánh kết quả tính toán hiệu chỉnh nồng độ DO với số liệu thực đo, dọc theo sông Nhuệ, tháng 12/2005. 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Phúc La Cự Đà Đò Quan Cống Thần Đò Kiều Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B Hình 7. So sánh giữa kết quả kiểm nghiệm DO với số liệu thực đo dọc sông Đáy và sông Nhuệ tháng 12/2005 Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 276 So sánh kết quả tính toán hiệu chỉnh nồng độ BOD với số liệu thực đo, dọc sông Đáy, tháng 12/2005. 0.0 10.0 20.0 30.0 40.0 50.0 60.0 Tế Tiêu Trung Hiếu Hạ Khánh Phú Độc BộĐò Mười Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B So sánh kết quả tính toán hiệu chỉnh nồng độ BOD với số liệu thực đo, dọc sông Nhuệ, tháng 12/2005. 0.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0 200.0 Phúc La Cự Đà Đò Quan Cống Thần Vị trí quan trắc Nồng độ (mg/l) Tính toán Thực đo TCVN 5942-B Hình 8. So sánh giữa kết quả kiểm nghiệm BOD với số liệu thực đo dọc sông Đáy và sông Nhuệ tháng 12/2005 Nhìn chung, các kết quả tính toán thu được còn có sự khác nhau giữa một số trạm thủy văn dùng để hiệu chỉnh mô hình, đặc biệt là kết quả tính toán tại trạm Hà Đông còn thấp, không được như mong muốn. Kết quả này có thể giải thích do trạm Hà Đông không nằm trong Mạng lưới trạm quốc gia của Bộ Tài Nguyên và Môi trường nên số liệu thu thập tại trạm không được chính xác. Hơn thế nữa, mặt cắ t sử dụng trên sông Nhuệ đo đạc cách đây đã lâu nên độ chính xác đã giảm đi rất nhiều. Kết quả hệ số nhám thu được trên các sông nghiên cứu có giá trị trung bình từ 0.03 đến 0.055, hệ số nhám này nói chung tương đối cao do thời gian tính toán trong mùa kiệt. Sau khi tính toán được các đặc trưng thuỷ lực của hệ thống sông đang nghiên cứu, nhóm tác giả đã áp dụng mô đun mô hình chất lượng nước để tính toán sự lan truy ền các chất ô nhiễm, cụ thể là sự lan truyền của 2 chỉ tiêu: BOD và DO. Sau khi hiệu chỉnh, mô hình đã cho kết quả với sai số của mô hình với số liệu đo đạc nhỏ hơn 10%. Điều đó chứng tỏ mô hình MIKE có thể sử dụng hiệu quả để mô phỏng sự lan truyền chất ô nhiễm trên sông. 4. Kết luận Mô hình MIKE 11 là một công cụ hữu ích dùng để tính toán và mô phỏng chất lượng nước một chiều trong sông tự nhiên với độ chính xác tương đối cao. Đây là một bộ phần mềm mang tính thương mại và phải có bản quyền sử dụng, tuy nhiên với sự chính xác, độ ổn định cao của các kết quả tính toán, chúng tôi đề xuất triển khai sử dụng bộ phần mềm này trên các sông thực tế để góp phần đánh giá hiệ n trạng ô nhiễm tại các sông ở Việt Nam. Căn cứ vào các kết quả tính toán cũng như quan trắc, khảo sát, có thể nhận thấy chất lượng nước trên sông Nhuệ, sông Đáy đã xấu hơn nhiều so với tiêu chuẩn nước mặt Việt Nam. Nếu các khu đô thị mới, khu công nghiệp mới được bố trí xây dựng không hợp lý thì chắc chắn tình trạng ô nhiễm môi trường nước sẽ ngày càng trầm trọng, dẫn đến kết quả cuối cùng là các dòng sông sẽ chết. Chính vì vậy, chúng tôi khuyến nghị áp dụng công cụ mô hình để tính toán dự báo chất lượng môi trường nước cho các sông ở Việt Nam với các kịch bản xả thải khác nhau, phù hợp với các kế hoạch phát triển kinh tế xã hội trong khu vực, nhằm lựachọn được phương án quy hoạch phát triển tối ưu, bảo đảm phát triển bền vững. Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 277 Tài liệu tham khảo 1 PGS TSKH Nguyễn Văn Cư - Báo cáo tổng kết đề tài cấp nhà nước: Xây dựng đề án tổng thể bảo vệ môi trường lưu vực sông Nhuệ và sông Đáy - Hà Nội, 2005. 2 V. Techow, D.R Maidment, L.W. Mays - Thuỷ văn ứng dụng - Đỗ Hữu Thành dịch - NXB Giáo dục -1994. 3 Trần Hồng Thái và nnk - Ứng dụng mô hình tính toán dự báo ô nhiễm môi trường nước cho các lưu vực sông: Cầu, Nhuệ - Đáy, Sài Gòn - Đồng Nai, Hà Nội, năm 2006. 4 TS Trần Hồng Thái - Numerical Methods for Parameter Estimation and Optimal Control of the Red River Network - Heidelberg, 2005. 5 P.G.Ciarlet and J.L.Lions - Finite Difference Methods- Elservier Science Publisher B.V - 1996. 6 DHI - MIKE 11 Reference Manual - DHI software 2004. 7 Lê Vũ Việt Phong - Đồ án tốt nghiệp - Nghiên cứu áp dụng mô hình MIKE11 tính toán chất lượng nước sông Nhuệ và sông Đáy - Hà Nội, 2006. Tuyển tập báo cáo Hội thảo khoa học lần thứ 10 - Viện KH KTTV & MT 278 . nghiệm cho mô hình chất lượng nước được nêu trong các hình từ Hình 5 đến Hình 8 21 -1 1- 1998 1 1- 1 2-1 998 3 1-1 2-1 998 2 0-1 -1 999 9-2 -1 999 1-3 -1 9 99 2 1-3 -1 999 1 0-4 -1 999 3 0-4 -1 999 2 0-5 -1 999 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 [meter] Time. bài toán lan truyền chất thông thường. 3. Nghiên cứu áp dụng mô hình MIKE 11 tính toán chất lượng nước sông Nhuệ và sông Đáy a. Số liệu đầu vào phục vụ cho mô hình thuỷ lực và chất lượng nước. NGHIÊN CỨU ÁP DỤNG MÔ HÌNH TOÁN MIKE 11 TÍNH TOÁN CHẤT LƯỢNG NƯỚC SÔNG NHUỆ - SÔNG ĐÁY Lê Vũ Việt Phong, Trần Hồng Thái, Phạm Văn Hải Viện Khoa học Khí tượng Thủy văn và Môi trường