Một số vấn đề khoa học trong nghiên cứu mô phỏng và tính toán chất lượng nước

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ NGUYỄN HỒNG PHONG MỘT SỐ VẤN ĐỀ KHOA HỌC TRONG NGHIÊN CỨU MƠ PHỎNG VÀ TÍNH TỐN CHẤT LƯỢNG NƯỚC Ngành: Cơ học kỹ thuật Chuyên ngành: Cơ học kỹ thuật Mã số: 60 52 02 TÓM TẮT LUẬN VĂN THẠC SĨ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Trần Thu Hà Hà Nội - 2012 MỤC LỤC MỞ ĐẦU CHƯƠNG MƠ HÌNH TÍNH TỐN MÔ PHỎNG CHẤT LƯỢNG NƯỚC 1.1 Tổng quan kết ngồi nước mơ hình tính tốn mơ chất lượng nước 1.2 Mơ hình lan truyền ô nhiễm chiều 1.3 Thuật toán giải hệ phương trình dịng chảy chiều 1.4 Thuật tốn giải phương trình truyền tải khuyếch tán chiều 10 1.5 Phát triển mơ hình truyền tải đa chất 11 1.6 Hiệu chỉnh số tính tốn thuật tốn 13 CHƯƠNG CÁC TÍNH TỐN KIỂM ĐỊNH CHẤT LƯỢNG MƠ HÌNH 14 2.1 Kiểm tra thuật giải số cho toán thủy lực chiều 14 2.2 Kiểm tra thuật giải số cho toán lan truyền chất 19 3.1 Số liệu địa hình thơng số đo đạc 25 3.2 Các kết tính tốn mơ 29 KẾT LUẬN 36 TÀI LIỆU THAM KHẢO 38 PHỤ LỤC 40 DANH MỤC BẢNG BIỂU VÀ HÌNH VẼ Bảng 3.1: Số liệu đo đạc tiêu ô nhiễm sử dụng mô hình Bảng 3.2: So sánh kết tính tốn vị trí có điểm đo đạc Hình 1.1: Lưới khơng cấu trúc dạng tam giác Hình 1.2: Quan hệ phần tử miền Hình 2.1: Cấu trúc lưới tính sử dụng tốn mẫu Hình 2.2 : So sánh mực nước tính tốn nghiệm xác theo chiều dài kênh Hình 2.3: Sai số tuyệt đối mực nước nghiệm xác tính tốn theo chiều dài kênh Hình 2.4: Cao độ đáy kênh dẫn tốn mẫu thủy lực Hình 2.5: So sánh mực nước tính tốn nghiệm xác theo chiều dài kênh Hình 2.6: Sai số tuyệt đối mực nước nghiệm xác tính tốn theo chiều dài kênh Hình 2.7: Cao độ đáy kênh tốn mẫu thủy lực Hình 2.8: So sánh mực nước tính tốn nghiệm xác theo chiều dài kênh Hình 2.9: Sai số tuyệt đối mực nước nghiệm xác tính tốn theo chiều dài kênh Hình 2.10: Điều kiện đầu tốn mẫu Hình 2.11: So sánh nghiệm tính tốn nghiệm xác tốn mẫu Hình 2.12: Phân bố chất ô nhiễm dọc kênh thời điểm 1000s Hình 2.13: So sánh kết tính tốn nghiệm xác tốn mẫu Hình 2.14: Kết tính tốn so sánh với nghiệm xác toán mẫu tương ứng với hệ số k khác Hình 2.15: Sai số tương đối (%) toán mẫu tương ứng với giá trị khác hệ số k Hình 3.1: Bản đồ địa hình khu vực hồ Thanh Nhàn Hình 3.2: Xử lý số liệu địa hình lưới khơng cấu trúc Hình 3.3a: Kết tính trường vận tốc Hình 3.3b: Vị trí điểm đo nồng độ nhiễm Hình 3.4: Kết tính BOD Hình 3.5: Kết tính COD Hình 3.6: Kết tính NH3 Hình 3.7: Kết tính PO4 Hình 3.8: Kết tính SO4 Hình 3.9: Kết tính NO3 MỞ ĐẦU Sự phát triển nhanh chóng kinh tế sản xuất tiêu dùng gây sức ép trực tiếp gián tiếp lên mức độ ô nhiễm môi trường sống Đặc biệt thành phố lớn, nơi tập trung khu sản xuất tiêu dùng sản phẩm công nghiệp, nông nghiệp, dịch vụ phục vụ sống dân cư Tình trạng nhiễm mơi trường nói chung nhiễm môi trường nước từ trước đến xuất nhiều nơi diễn nghiêm trọng Bên cạnh việc nhiễm sơng ngịi nhiễm nguồn nước hồ chứa đô thị nguồn nước thải sinh hoạt sản xuất đáng lo ngại chúng ảnh hưởng trực tiếp đến sức khỏe phận không nhỏ dân cư sinh sống quanh hồ ảnh hưởng cảnh quan đô thị Do đặc điểm địa hình lịng hồ quy mơ diện tích lịng hồ nhỏ lịng sơng nhiều nên mơ hình chiều khơng cịn phù hợp để nghiên cứu đánh giá q trình nhiễm hồ Do mơ hình dịng chảy hai chiều lan truyền chất ô nhiễm hai chiều cần sử dụng để tính tốn q trình lan truyền nhiễm lịng hồ từ nguồn thải Mơ hình tính tốn truyền tải chất nhiễm PGS Hoàng Văn Lai cộng bắt đầu phát triển từ nghiên cứu chương trình nghiên cứu Viện học Từ kết nghiên cứu ban đầu mơ hình, học viên lựa chọn hướng nghiên cứu với hy vọng đạt bước tiến thêm việc tìm hiểu, phát triển đánh giá chất lượng mơ hình, sau bước đầu ứng dụng để mơ lan truyền chất ô nhiễm khu vực cụ thể với số tiêu nhiễm Bố cục luận văn gồm phần mở đầu, phần kết luận chương chính: Chương 1: MƠ HÌNH TÍNH TỐN MƠ PHỎNG CHẤT LƯỢNG NƯỚC Chương 2: CÁC TÍNH TỐN KIỂM ĐỊNH CHẤT LƯỢNG MƠ HÌNH Chương 3: THỬ NGHIỆM MÔ PHỎNG CHẤT LƯỢNG NƯỚC HỒ THANH NHÀN CHƯƠNG MƠ HÌNH TÍNH TỐN MƠ PHỎNG CHẤT LƯỢNG NƯỚC 1.1 Tổng quan kết nước mơ hình tính tốn mơ chất lượng nước 1.1.1 Các nghiên cứu giới Trên giới, việc nghiên cứu mô chất lượng nước quan tâm ys đến từ sớm đạt nhiều thành tựu Không dừng lại mức nghiên cứu, có nhiều mơ hình phát triển thành sản phẩm phần mềm thương mại ví dụ Mike viện thủy lực Đan Mạch (DHI) Một số mơ hình khác tiếng liệt kê sau: - Mơ hình WASP (Water Quality Analysis Simulation Program 7) phát triển từ mơ hình (WASP – Di Toro, 1983; Connolly vaf Winfield, 1984; Ambrose, R.B, 1988) Mơ hình cho phép tính toán 1, 2, chiều với nhiều thành phần chất nhiễm Mơ hình WASP liên kết với mơ hình thủy động lực vận chuyển trầm tích để tính tốn trường dịng chảy, nhiệt độ, độ muối thơng lượng trầm tích - Mơ hình QUAL2K (hay Q2K) (River and Stream Water Quality Model) phát triển từ mơ hình QUAL2E (hay Q2E (Brown Barnwell 1987)) Đây mơ hình mô chất lượng nước suối sông chiều Mơ hình tính tốn chu trình Nitơ thơng qua chu trình chuyển hóa nitơ để biểu diễn hợp chất cacbon, loại cacbon hữu không sống - Bộ phần mềm MIKE (DHI) phát triển môi trường đồ họa hệ điều hành Windows bao gồm nhiều modul tính tốn có tính tốn chất lượng nước MIKE tích hợp tốt với hệ thống thông tin đồ địa lý (GIS) nhiều đề tài nghiên cứu nước chọn để sử dụng cơng cụ hỗ trợ tính toán 1.1.2 Các nghiên cứu nước Ở nước gần có nhiều nghiên cứu quan tâm đến việc sử dụng mơ hình để mơ tính tốn chất lượng nước Các đề tài nghiên cứu chia làm hai hướng: - Hướng nghiên cứu thứ sử dụng mơ hình phần mềm nước ngồi để mơ tính tốn số yếu tố chất lượng nước khu vực cụ thể Hướng thiên nghiên cứu ứng dụng mơ hình tác động chất nhiễm đến khu vực nghiên cứu để phục vụ nhu cầu nghiên cứu, đánh giá tổn thất ô nhiễm hoạch định sách phát triển địa phương Ví dụ chương trình hợp tác với Cơ quan hợp tác Quốc tế Nhật Bản - JICA Viện Tài nguyên Môi trường biển – Viện Khoa học Công nghệ Việt Nam sử dụng phần mềm chuyên dụng CABARET LOICZ (Mỹ) để đánh giá mức độ tích tụ khuếch tán vật chất số điểm thuộc vịnh Hạ Long Một số nghiên cứu khác sử dụng mơ hình eco Lab MIKE - Hướng nghiên cứu thứ hai vào xây dựng mơ hình thuật toán riêng học tập giải pháp cơng bố tạp chí khoa học quốc tế để giải tốn lan truyền nhiễm với hay nhiều chất ô nhiễm Đánh giá hiệu tính tốn mơ hình tiến hành thử nghiệm với khu vực nghiên cứu cụ thể Ví dụ mơ hình tính tốn chất lượng nước cho vùng biển Quảng Ninh-Hải Phịng xây dựng mơ hình thủy lực hai chiều CHB-2D phòng Cơ học biển (Viện học)[5] Ngồi cịn kể đến mơ hình khác nghiên cứu phát triển HydroGIS TS Nguyễn Hữu Nhân, SALBOD GS.TS Nguyễn Tất Đắc vv Mơ hình tính tốn lan truyền chất nhiễm phịng thủy tin học (Viện Cơ học) theo hướng học tập xây dựng mơ hình để giải tốn lan truyền chất nhiễm với phạm vi ứng dụng lan truyền ô nhiễm đoạn sông ô nhiễm hồ chứa Mơ hình xây dựng mơ hình thủy lực IMECH_2DUSZ đề tài cấp Bộ: “Áp dụng tính tốn song song việc kết nối mơ hình thuỷ văn mơ hình thuỷ lực phục vụ dự báo kiểm sốt lũ lụt lưu vực sơng Hồng”[1] Đây mơ hình thuỷ lực sử dụng hệ phương trình Saint Venant 2D đầy đủ để mơ tả chuyển động nước miền chiều giải phương pháp FVM (Finite Volume Method) lưới khơng cấu trúc Mơ hình học viên lựa chọn để nghiên cứu đánh giá tiếp tục phát triển Phạm Tiến Đạt(2009), luận văn thạc sĩ "Sử dụng mơ hình eco lab đánh giá số đặc trưng môi trường khu vực nuôi trồng thủy sản" 1.2 Mơ hình lan truyền nhiễm chiều Mơ hình lan truyền chất gây nhiễm chiều bao gồm mơ hình thủy lực chiều mơ hình truyền tải khuyếch tán vật chất chiều Hệ phương trình S Venant 2D mơ tả định luật bảo toàn khối lượng động lượng nước chiều Hệ phương trình thiết lập từ hệ phương trình Navier - Stock cách trung bình hóa theo chiều sâu điều kiện áp suất tuân theo quy luật thuỷ tĩnh Nếu bỏ qua thành phần nhớt rối, ảnh hưởng gió ảnh hưởng lực Coriolis, hệ phương trình S Venant 2D viết dạng sau([1]-[2]) Phương trình liên tục:  z     uh   t x    y vh (1) Phương trình động lượng theo hướng x: u (2) u t Phương trình động lượng theo hướng y: (3) v u t Phương trình truyền tải khuyếch tán vật chất ([3],[4]): C i t  Trong phương trình trên: z – mực nước, h độ sâu dòng chảy u,v – vận tốc g – gia tốc trọng trường Kx, Ky – hệ số Strickler lực cản đáy Ci – nồng độ tiêu ô nhiễm thứ i fi – thành phần nguồn So,x,So,y – độ dốc đáy D – Hệ số khuyếch tán  x uC Theo lý thuyết đường đặc trưng, số điều kiện biên cho hệ phương trình S.venant (1)(3) phụ thuộc vào vận tốc un số Froude - Fr theo hướng pháp tuyến điểm biên un = unx + vny Fr  un / c, c  gh Có loại điều kiện biên sau: Biên có dịng chảy vào miền (inflow): un  0, dịng chảy có chiều ngược với pháp tuyến - Trạng thái chảy êm: un > c: cần cho điều kiện biên - Trạng thái chảy xiết: un  c: cần cho điều kiện biên Biên có dịng chảy khỏi miền (outflow): un  0, dịng chảy theo chiều pháp tuyến ngồi - Trạng thái chảy êm: un < c: cần cho điều kiện biên - Trạng thái chảy xiết: un  c : không cần cho điều kiện biên Đối với phương trình truyền tải khuyếch tán vật chất cần cho điều kiện biên vào điều kiện biên Ví dụ tốn thủy lực có dịng chảy vào êm, cần cho điều kiện biên vào điều kiện biên ra, để tốn ổn định giảm sai số điều kiện biên vào thường chọn vận tốc điều kiện biên chọn cao trình mực nước Với tốn truyền tải khuyếch tán điều kiện biên vào thường chọn nồng độ tiêu ô nhiễm điều kiện biên đạo hàm nồng độ tiêu nhiễm 1.3 Thuật tốn giải hệ phương trình dịng chảy chiều Hệ phương trình Saint – Venant 2D (1) - (3) giải số theo phương pháp khối hữu hạn (FVM - Finite Volume Method) lưới không cấu trúc dạng tam giác cho đoạn sông hồ chứa Lưới tam giác mô tả tốt dịng chảy miền có địa hình phức tạp (cả biên đáy) Trong phương pháp khối hữu hạn, miền tính tốn D chia nhỏ thành khối có hình học đơn giản Si dạng tam giác khối tạo thành lưới phủ kín miền tính tốn Các tham số để mơ tả lưới nút lưới liên kết nút lưới để tạo thành khối Trên hình phần lưới không cấu trúc Tính khơng cấu trúc lưới thể qua việc đánh số nút lưới đánh số khối lưới không cần theo quy luật Hình 1.1: Lưới khơng cấu trúc dạng tam giác Việc chia miền D thành phần tử Si cần đảm bảo điều kiện sau: N 1) D S i i 1 2) Si  Sj =  3) Sj= , Si không nằm cạnh Sj Si Mi,j, Si nằm cạnh Sj, Mi,j cạnh chung phần tử Hình 1.2: Quan hệ phần tử miền Ở ta dùng ký hiệu lý thuyết tập hợp  hợp,  giao,  tập rỗng S i phủ Si, đỉnh phần tử nằm hai đầu cạnh, khơng có đỉnh phần tử nằm cạnh phần tử khác Tính khơng cấu trúc lưới thể qua việc đánh số nút k đánh số phân tử phần tử không theo quy luật Những lưới đánh số nút (hoặc phần tử) theo số i, j, i theo chiều x, j theo chiều y thường phải liền kề có trật tự gọi lưới có cấu trúc Lưới không cấu trúc nút, phần tử đánh số số khơng có điều kiện ràng buộc cho số Trong chia lưới miền D, phần tử không bắt buộc có cấu trúc hình học tam giác tứ giác Thí dụ miền D chia lưới gồm tam giác tứ giác Tính chất hình học phần tử không tuân theo nguyên tắc Tuy nhiên đa giác phải đa giác lồi Để đảm bảo kết tính tốn tốt độ dài cạnh đa giác khơng khác q nhiều Vì vậy, đa giác chỗ miền nhỏ, chỗ khác to, phụ thuộc vào tính chất miền, từ chỗ nhỏ sang chỗ to có q trình thay đổi khơng q đột ngột Cần lưu ý mơ hình sử dụng nhiều thơng số hình học khoảng cách, diện tích 44 ir8=fvvolume(invo,1) else endif if(fvconect(invo,k).le.0.and.fvvolbty(invo).eq.-1) goto 560 if(fvconect(invo,k).le.0.and.fvvolbty(invo).eq.-2) goto 570 if(fvconect(invo,k).le.0.and.fvvolbty(invo).eq.-3) goto 580 ir3=fvconect(invo,k) k2=fvcfacek(invo,k) ufk=fvuface(ir3,k2) vfk=fvvface(ir3,k2) ! ufi=0.5*(ufi+ufk) ! vfi=0.5*(vfi+vfk) i1=1,nchat hfi=(concent(invo,i1)+concent(ir3,i1)+concen(invo,i1)+concen(ir3,i1))/4 r7=(concentx(invo,i1)+concentx(ir3,i1)+concentx1(invo,i1)+concentx1(ir3,i1))/4 r8=(concenty(invo,i1)+concenty(ir3,i1)+concenty1(invo,i1)+concenty1(ir3,i1))/4 hs2(i1)=normalxp(invo,k)*(ufi*hficoef_dis_mu(invo,i1)*r7)+normalyp(invo,k)*(vfi*hfi-coef_dis_mu(invo,i1)*r8) enddo goto 590 560 continue i1=1,nchat hfi=(concent(invo,i1)+concen(invo,i1))/2 r7=(concentx(invo,i1)+concentx1(invo,i1))/2 ! ??? Hinh chieu vecto len canh? r8= (concenty(invo,i1) +concenty1(invo,i1))/2 hs2(i1)=normalxp(invo,k)*(ufi*hficoef_dis_mu(invo,i1)*r7)+normalyp(invo,k)*(vfi*hfi-coef_dis_mu(invo,i1)*r8) 1/4/08 i1=1,nchat hfi=(concent(invo,i1)+concen(invo,i1)+2*fvbou2c(fvnb2vol(invo),i1))/4 !XL tam r7=(concentx(invo,i1)+concentx1(invo,i1))/2 r8=(concenty(invo,i1) +concenty1(invo,i1))/2 hs2(i1)=normalxp(invo,k)*(ufi*hficoef_dis_mu(invo,i1)*r7)+normalyp(invo,k)*(vfi*hfi-coef_dis_mu(invo,i1)*r8) hfi=(concent(invo,i1)+concen(invo,i1))/2 r7=(concentx(invo,i1)+concentx1(invo,i1))/2 r8=(concenty(invo,i1) +concenty1(invo,i1))/2 hs2(i1)=normalxp(invo,k)*(ufi*hficoef_dis_mu(invo,i1)*r7)+normalyp(invo,k)*(vfi*hfi-coef_dis_mu(invo,i1)*r8) enddo goto 590 !other bou.con.? 590 continue ! if(abs(rfk1).le.eps) rfk1=0 ! if(abs(rfk2).le.eps) rfk2=0 ! if(abs(rfk3).le.eps) rfk3=0 i1=1,nchat ! hs2(i1)=hs2(i1)*fvvollen(invo,k) hs1(i1)=hs1(i1)+hs2(i1) enddo 45 end !print *,invo,rf1,rf2,rf3 ! fvsource(invo,2),s0bottx(invo),r1 ! if(invo/40*40.eq.invo)pause ! Cong cac chat chuyen hoa subcen(1:nchat)=0 i1=1,nchat r1=0 k1=1,nkoef_sou(i1) r1=r1+concen(invo,conc_sou(i1,k1))*koef_sou(i1,k1,invo)*det if (i1.ne.conc_sou(i1,k1)) then subcen(conc_sou(i1,k1))=subcen(conc_sou(i1,k1))+ concen(invo,conc_sou(i1,k1))*koef_sou(i1,k1,invo)*det endif enddo concen(invo,i1)=r1+concent(invo,i1)- hs1(i1)*det/fvarvol(invo) !concen(invo,i1)=(concent(invo,i1)- hs1(i1)*det/fvarvol(invo))/(1koef_sou(1,1,invo)*det) enddo ! tru` cac chat da chuyen hoa i1=1,nchat concen(invo,i1)=concen(invo,i1)subcen(i1) if (concen(invo,i1)0) then ir1=1,nbienc i1=1,ncpoibou(ir1) concentp(cpoibou(i1,ir1),1:nchat)=valpoic(1:nchat,ir1) enddo enddo endif ! TINH LAI DAO HAM CONCENT call tinhlaidaohamconcent() err=concen-concent r1=max(err) if (r110) then write(*,*)'Mo hinh khuyech tan so' stop endif enddo ! =================================== =============== B Số liệu địa hình hồ Thanh Nhàn số liệu đo đạc số điểm lịng hồ Số liệu địa hình lịng hồ Tọa độ X 589233.31 589072.62 589000.12 589025.12 47 589093.38 589287 588993.12 589132.81 589071 589041.88 589133.69 589279 589059.69 589131.12 589011.19 589097.19 589067.81 589078.12 589266.81 589087.69 589167.38 589068.88 589151.81 589113.31 589132 589007.12 589242.38 589011.5 589273.38 589123.88 589129.81 589094.38 589015.31 589034.62 589060.62 589019.31 589049.12 589133.81 589106.81 589018.12 589134 589127.88 589131.31 589124.88 589132 589123 589041.88 589207.38 589057 589051.88 589014 589121.31 589090.5 589021.12 589128.88 589132.19 588994.19 589119.19 589117.81 589254.88 589044.38 589071.31 589134 589114.38 589081 589122 589187.38 589280 589017.81 589037 589046.62 48 589132 589132.69 589132 589061.38 589249.19 589100.81 589084.12 589063.88 589022.69 589162 589233 589084.62 589090.69 589133.12 589129.69 589098.5 589185.31 589220.12 589128.31 589119.5 589048.38 589126 588999.31 589204.88 589173.12 589126.81 589107.88 589261.62 589111.12 589054.12 589056.19 589195.19 589039.31 589152.62 589134 589123.12 589024.12 589125.62 589068.38 589075.19 589015 589112.81 589130.62 589240.88 589066.31 589015 589087.31 589032.5 589064.69 589226 589266.62 589109.19 589016 589058 589028.62 589004.12 589126.81 589037.12 589070 589254 589081.12 589052.5 588997.12 589137.5 589134 589063.12 589016 49 589131.31 589116 589124.31 589142.38 589133.38 589076.31 589231.31 589103.81 589003.88 589032.12 589103 589227.62 589038.19 589063.12 589257.81 589063.69 589072.69 589125.69 589065.69 589073 589148.62 589162.62 589122.38 589115 589085.81 589081.5 589195.5 589103.5 589013.12 589044.81 589075.62 589117.5 589111.19 589132.5 589104.38 589041.19 589096.5 589092.62 589068.19 589106.12 589086.88 589131.19 589107.81 589110.5 589058.31 589113.81 589080.12 589079.69 589084 589092.19 589067 589115.5 589077.19 589112.5 589059.88 589061.31 589087.19 589066.69 589097.69 589076.5 589032.38 589199.5 589065.12 589106.31 589153.62 589109.31 589106 50 589054.5 589098.38 589078.69 589059.5 589062.19 589057.69 589099.81 589137.31 589080.62 589066.38 589110.62 589115.81 589077.81 589123.69 589071.62 589094.62 589122.62 589129.19 589102 589082 589091.69 589088.88 589063.88 589080.5 589179.12 589024 589084.81 589099.81 589174.88 589141.81 589104.5 589065.19 589179.31 589049.38 589130.38 589183.62 589054.69 589132 589070.62 589089.38 589107 589110 589105.81 589083.81 589072 589088.19 589123.88 589254.81 589117.81 589123.31 589062.81 589094.5 589066.88 589086.12 589124.19 589071.81 589120.38 589114.88 589244.69 589208.81 589109.12 589077.69 589089.62 589104.81 589101.5 589097.12 589064.88 51 589129.88 589142.81 589049.69 589167.31 589201 589066 589026.19 589027.69 589077.12 589153.5 589085.81 589157.81 589043.38 589223.88 589200.88 589078.88 589091.88 589067 589094.31 589045.81 589058.88 589015.88 589106.81 589109.81 589112.69 589136 589123.12 589106.31 589037.69 589047 589225.81 589130.12 589092.62 589145.88 589088.38 589070.38 589076.12 589079.31 589100.88 589083.12 589071.81 589121.31 589133.62 589115.62 589058.5 589126.5 589074.88 589146.5 589100.88 589105.5 589049 589040.88 589224.81 589089.69 589052.31 589075.81 589101.81 589111.62 589178.81 589112.62 589114.38 589072.31 589140.38 589059.5 589111.31 589086.62 589099.38 52 589062.62 589115.31 589082.12 589107.62 589107.31 589073.19 589070 589096.31 589043.19 589231.62 589123.81 589121.19 589090.88 589086.12 589082 589051.81 589131.69 589071.69 589088.69 589221.38 589074.62 589119.62 589110.88 589118.69 589108.81 589233.38 589064.88 589092.88 589070.38 589115.88 589057.31 589122.88 589119.62 589116.38 589113.31 589124.88 589119.69 589107.81 589124.81 589144.5 589092.5 589128.38 589147 589096.19 589174.31 589109.81 589119.69 589075.12 589117.81 589022.81 589102 589115.38 589069.31 589076.38 589077.31 589101.69 589105.5 589083.62 589111.19 589027.19 589066.69 589056.69 589058.38 589087.81 589130.5 589121.19 589090.69 53 589077.5 589092.38 589148.62 589167.88 589166 589042.31 589071.12 589126.12 589074.5 589144.62 589109.12 589065.19 589110.19 589211.5 589078.81 589095.5 589150.81 589104.19 589029.31 589043.38 589095.19 589012.38 589054.5 589099.81 589127.31 589093.12 589130.81 589082.69 589079.81 589127.88 589081.31 589112.12 589131.62 589073.88 589118.5 589122.19 589071.31 589270.19 589109.81 589092.5 589095.31 589075.81 589060 589052.12 589069.31 589055.19 589012 589090.31 589103.88 589078.31 589142.69 589084.81 589105.69 589058.19 589046.69 589069.31 589065.62 589061.81 589123.12 589113 589116.5 589113.81 589097.62 589089.12 589093.31 589122.62 589108.81 54 589125.88 589051.19 589117.62 589124.81 589121 589112.5 589066.5 589102.69 589055.5 589156.88 589095.81 589211.62 589109.69 589192.88 589034.81 589087.5 589106.81 589103.38 589085.62 589121.88 589111.69 589124.5 589124.31 589029.5 589073.62 589114.62 589087 589082.62 589131.69 589075.81 589089.5 589108 589114.69 589129.88 589080.88 589102.19 589074.31 589094.19 589110.12 589115.88 589050.38 589070.19 589028 589240.5 589066.69 589171.69 589113.81 589056.5 589061.19 589119 589071.5 589068.31 589064.5 589099.12 589131 589122.5 Số liệu đo đạc thu thập HỒ THANH NHÀN (11-2001) : 55 Nằm sau bệnh viện Thanh Nhàn,hồ có diện tích 8,1 ha, độ sâu 1,5 - m, dung tích trung bình 162.0000 m3 Lượng nước thải đổ vào hồ khoảng 2100 m3/ngày-đêm Đặc tính thuỷ lý o o - Nhiệt độ nước hồ dao động từ 23,4 C đến 23,8 C mùa khô - Độ pH dao động mức kiềm ( 7,88 - 8, 57 ) không chênh lệch nhiều - Độ đục dao động từ - 34 mg/l - Độ dẫn không thay đổi trạm thu mẫu 0,04 S/m - Độ muối ( NaCL) không thay đổi điểm thu mẫu 0.01% - Hàm lượng xy hồ tan cao dao động từ 7,9 mg/l - 11,6 mg/l Đặc tính thuỷ hố : * Nhóm yếu tố hữu : - Hàm lượng yếu tố có nguồn gốc Nitơ dạng ionNO2,NO3, NH3 , có diễn biến sau : + Hàm lượng N02 dao động từ 0,055 - 0,062 mg/l cao giá trị tới hạn cho phép tiêu chuẩn nước mặt Việt Nam ( 0,01- 0,05mg/l) + Hàm lượng NO3 dao động từ 0,76 - 0,81 mg/l thấp nhiều so với giá trị giới hạn cho phép tiêu chuẩn nước mặt Việt nam ( 10 - 15 mg/l) + trị Hàm lượng NH3 dao động từ 0,51 đến 0,62 mg/l không cao so với giá giới hạn tiêu chuẩn nước mặt Việt nam ( mg/l) + Hàm lượng phốt dạng muối phốt phát PO4 hoà tan dao động từ 0,44 đến 46 mg/l + - + Hàm lượng muối sunphat ( S04 ) dao động từ 25 - 27 mg/l + Hàm lượng sunphit dạng H2S dao động từ 0,116 - 0,151 mg/l Nhu cầu ô xy sinh hoá (BOD5) dao động 15,5 đến 24,0 mg/l thấp giá trị giới hạn tiêu chuẩn Việt nam chất lượng nước mặt ( 25 mg/l) + Nhu cầu xy hố học (COD) dao động từ 345 - 400 mg/l Như nhu cầu ô xy hoá học cao khoảng 10 lần giá trị giới hạn cho phép theo tiêu chuẩn nước mặt Việt nam ( 35 mg/l) * Nhóm độc tố hoá học : Bao gồm yếu tố vi lượng kim loại nặng (cadmi, thuỷ ngân, chì măng gan ) Nhìn chung, hàm lượng trung bình yếu tố kim loại nặng thấp giá trị tới hạn cho phép chất lượng nước tầng mặt theo tiêu chuẩn Việt Nam (1995) 56 Bảng Chất lượng nước hồ Thanh nhàn tháng 11/2001 Điểm thu mẫu Nhiệt độ oC Oxy hoà tan (mg/l) PH Độ dẫn (S/m) NaCL(%) Độ đục (mg/l) NO2-N (mg/l) NO3-N (mg/l) NH3 (mg/l) SO4( mg/l) PO4 ( mg/l) COD( mg/l) BOD520oC (mg/l) Total Coliform (MPN/100ml) Fecalcoliform (MPN/100ml) Cu (mg/l) Mn ( mg/l) Pb (g/l) Cd (g /l) Hg (g /l) Sơ đồ ô nhiễm hệ số ô nhiễm sử dụng tính hồ Thanh Nhàn "config_onhiem.txt" BOD 0.0000017 (n he so koef_sou(i )+ chat onhiem lien quan) -1.75E-04 BOD (n he so koef_sou(i )+ chat onhiem lien quan) -1.6E-04 NH3 1.24E-05 BOD NO30.0000017 57 (n he so koef_sou(i )+ chat onhiem lien quan) -1.6E-04 NO3 1.6E-04 NH3 (n he so koef_sou(i )+ chat onhiem lien quan) 0.610E-06 BOD (n he so koef_sou(i )+ chat onhiem lien quan) -0.000045 COD (n he so koef_sou(i )+ chat onhiem lien quan) -0.000045 SO4 [so lieu Bien o nhiem] BOD 24 NH3 0.52 NO3 0.81 PO4 0.44 COD 400 SO4 27 ... CHẤT LƯỢNG MƠ HÌNH Chương 3: THỬ NGHIỆM MÔ PHỎNG CHẤT LƯỢNG NƯỚC HỒ THANH NHÀN 3 CHƯƠNG MƠ HÌNH TÍNH TỐN MƠ PHỎNG CHẤT LƯỢNG NƯỚC 1.1 Tổng quan kết ngồi nước mơ hình tính tốn mơ chất lượng nước. .. nhiều đề tài nghiên cứu nước chọn để sử dụng công cụ hỗ trợ tính tốn 1.1.2 Các nghiên cứu nước Ở nước gần có nhiều nghiên cứu quan tâm đến việc sử dụng mơ hình để mơ tính tốn chất lượng nước Các đề. .. nước Các đề tài nghiên cứu chia làm hai hướng: - Hướng nghiên cứu thứ sử dụng mơ hình phần mềm nước ngồi để mơ tính toán số yếu tố chất lượng nước khu vực cụ thể Hướng thiên nghiên cứu ứng dụng