Nghiên cứu giải pháp tối ưu hóa phèn sắt fecl3 và polymer trong xử lý nước tại công ty cp đầu tư và kinh doanh nước sạch sài gòn

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI Hồ Thanh Cường NGHIÊN CỨU GIẢI PHÁP TỐI ƯU HÓA PHÈN SẮT (FeCl3) VÀ POLYMER TRONG XỬ LÝ NƯỚC TẠI CÔNG TY CP ĐẦU TƯ VÀ KINH DOANH NƯỚC SẠCH SÀI GÒN LUẬN VĂN THẠC SĨ Chuyên ngành : Mã số : Kỹ thuật Cấp thoát nước 8580213 Người hướng dẫn khoa học : TS Võ Anh Tuấn TP.HỒ CHÍ MINH - 2019 i LỜI CAM ĐOAN Họ tên : HỒ THANH CƯỜNG Ngày sinh : 16/7/1968 Cơ quan công tác : Công ty CP Đầu tư Kinh doanh nước Sài Gòn Tác giả đề tài : Nghiên cứu giải pháp tối ưu hóa phèn sắt (FeCl3) Polymer xử lý nước Công ty CP Đầu tư Kinh doanh nước Sài Gòn Học viên lớp cao học : 25CTN11 – CS2 Chuyên ngành : Kỹ thuật Cấp nước Mã số : 8580213 Tơi xin cam đoan cơng trình cơng trình nghiên cứu thực cá nhân học viên hướng dẫn khoa học TS Võ Anh Tuấn Cơng trình chưa công bố lần Tất nội dung tham khảo trích dẫn nguồn đầy đủ theo quy định Tôi xin chịu trách nhiệm nội dung lời cam đoan Học viên thực Hồ Thanh Cường Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 ii LỜI CẢM ƠN Em xin chân thành gửi lời cảm ơn sâu sắc đến TS Võ Anh Tuấn người tận tình hướng dẫn giúp đỡ cho em hoàn thành đề tài này; xin chân thành cảm ơn Quý Thầy Cô Bộ mơn Cấp Thốt Nước Khoa Kỹ thuật – Tài nguyên nước , Trường Đại học Thủy lợi, truyền đạt kiến thức quý báu cho em suốt q trình học tập trường Bên cạnh đó, em xin chân thành gửi lời cảm ơn đến Ban Lãnh đạo tập thể cán nhân viên Công ty cổ phần Đầu tư Kinh doanh nước Sài Gòn cung cấp tài liệu tạo điều kiện cho em thu thập liệu số liệu quan trọng cần thiết, thông tin hữu ích để em hồn thành đề tài Với nỗ lực cố gắng thân trình thực đề tài chắn khơng tránh khỏi sai sót khuyết điểm trình thực luận văn Tuy nhiên, với ý kiến đóng góp chân thành từ quý Thầy Cô kiến thức trang bị trình học tập Trường tảng, hành trang quý báu giúp em hoàn thiện kỹ chun mơn nghiệp vụ, góp phần hồn thành nhiệm vụ sản xuất, cung cấp nước an toàn, liên tục, chất lượng, phục vụ ngày tốt cho nhân dân thành phố Hồ Chí Minh Trân trọng Tp.Hồ Chí Minh, ngày 06 tháng 01 năm 2019 Học viên thực Hồ Thanh Cường Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 MỤC LỤC Trang LỜI CAM ĐOAN i LỜI CẢM ƠN ii DANH MỤC CÁC HÌNH ẢNH DANH MỤC CÁC BẢNG BIỂU DANH MỤC CÁC TỪ VIẾT TẮT 10 LỜI MỞ ĐẦU 12 CHƯƠNG I TỔNG QUAN 13 1.1 Sơ lược Nhà máy nước trực thuộc Công ty CP Đầu tư Kinh doanh nước Sài Gòn 13 1.1.1 Quy mơ, cơng suất Nhà máy 13 1.1.2 Quy trình cơng nghệ Nhà máy 14 1.2 Quy trình cơng nghệ nhà máy nước tương tự So sánh phân tích ưu nhược điểm 25 1.2.1 Quy trình công nghệ nhà máy nước tương tự khu vực Thành phố Hồ Chí Minh 25 1.2.1.1 Nhà máy nước Thủ Đức trực thuộc Tổng Công ty Cấp nước Sài Gòn – Trách nhiệm hữu hạn thành viên (SAWACO) 25 1.2.1.2 Nhà máy nước trực thuộc Công ty CP BOO Thủ Đức 26 1.2.1.3 Nhà máy nước Tân Hiệp trực thuộc Tổng Công ty Cấp nước Sài Gòn – Trách nhiệm hữu hạn thành viên (SAWACO) 27 1.2.1.4 Nhà máy nước Tân Hiệp trực thuộc Công ty CP Đầu tư nước Tân Hiệp 28 1.2.2 Quy trinh công nghệ nhà máy nước giới 28 1.2.2.1 Nhà máy xử lý nước Auchel (CHLB Đức) 28 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 1.2.2.2 Nhà máy xử lý nước Biebesheim (CHLB Đức) 29 1.2.3 Phân tích ưu nhược điểm quy trình cơng nghệ mà Nhà máy áp dụng với quy trình cơng nghệ nhà máy nước tương tự khu vực Thành phố Hồ Chí Minh 30 1.2.3.1 Ưu điểm 30 1.2.3.2 Nhược điểm 30 1.3 Tính cấp thiết đề tài 30 1.4 Các bước tiếp cận vấn đề nghiên cứu 32 CHƯƠNG CƠ SỞ LÝ THUYẾT VÀ CƠ SỞ DỮ LIỆU NGHIÊN CỨU GIÚP TỐI ƯU HĨA QUY TRÌNH CHÂM HĨA CHẤT PHÈN SẮT (FECL3) VÀ POLYMER 33 2.1 Các quy định tiêu chuẩn chất lượng nước thô, chất lượng nước sau xử lý trước hòa vào mạng lưới cấp nước Thành phố 33 2.2 Cơ sở lý thuyết q trình keo tụ, tạo bơng lắng hệ thống xử lý nguồn nước mặt 33 2.2.1 Lý thuyết trình keo tụ 33 2.2.2 Các chất keo tụ thơng thường 35 2.2.3 Phương trình thủy phân muối kim loại vô 36 2.3 Hiện trạng trình vận hành nhà máy xử lý nước SWIC 39 2.4 Giới thiệu Phần mềm R bước thực phân tích liệu R 43 2.4.1 Giới thiệu sơ lược phần mềm R 43 2.4.2 Các bước phân tích liệu R 44 CHƯƠNG III NGHIÊN CỨU, ĐỀ XUẤT GIẢI PHÁP TỐI ƯU CHO QUY TRÌNH CHÂM HĨA CHẤT PHÈN SẮT (FeCl3) VÀ POLYMER TẠI NHÀ MÁY NƯỚC TRỰC THUỘC SWIC 50 3.1 Mơ tả quy trình, trang thiết bị sử dụng cho việc nghiên cứu 50 3.1.1 Thu thập biên tập liệu dùng cho việc nghiên cứu phân tích Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 50 3.1.2 Mơ tả quy trình trang thiết bị dùng cho việc nghiên cứu phân tích liệu 54 3.1.2.1 Mơ tả quy trình phân tích liệu 54 3.1.2.2 Các trang thiết bị dùng cho việc nghiên cứu 54 3.2 Phân tích liệu thu thập trình nghiên cứu 54 3.2.1 Phân tích diễn biến độ đục nước thơ sơng Đồng Nai so sánh độ đục nước thơ bình qn mùa khô mùa mưa qua năm 2016, 2017, 2018 54 3.2.2 Phân tích mối tương quan biến số 56 3.2.3 Xác định tầm quan trọng biến số 61 3.2.4 Phân tích phương sai biến số 65 3.2.4.1 Phân tích phương sai cho biến số “FeCl3” 3.2.4.2 Tương tự cho việc phân tích phương sai hậu định phương sai biến số độ đục nước thô “rw.ntu” qua năm 2016, 2017, 2018 3.2.4.3 Nhận xét 65 70 72 3.2.5 Phân tích mơ hình hồi quy tuyến tính cho biến số “FeCl3” có tính đến yếu tố ảnh hưởng tương tác biến số “rw.ntu” biến số “seasons” giải định 72 3.3 Kết thực mô hình nghiên cứu tối ưu 90 KẾT LUẬN VÀ KIẾN NGHỊ 99 TÀI LIỆU THAM KHẢO PHỤ LỤC SỐ CÁC KẾT QUẢ THỬ NGHIỆM JARTEST PHỤ LỤC SỐ CÁC CÂU LỆNH DÙNG TRONG PHÂN TÍCH DỮ LIỆU BẰNG R Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 DANH MỤC CÁC HÌNH ẢNH Trang Hình 1.1 Văn phịng Cơng ty CP Đầu tư Kinh doanh nước Sài Gịn 13 Hình 1.2 Sơ đồ khối quy trình cơng nghệ Nhà máy 14 Hình 1.3 Vị trí trạm bơm nước thơ họng thu nước Hóa An Google Map 15 Hình 1.4 Bên Trạm bơm nước thơ 15 Hình 1.5 Phối cảnh hạng mục khu xử lý nước nhà hành Cơng ty 16 Hình 1.6 Sơ đồ mặt khu nhà lắng 17 Hình 1.7 Minh họa sơ đồ cụm bể lắng 17 Hình 1.8 Sơ đồ mặt khu nhà lọc 20 Hình 1.9 Trạm bơm nước (Trạm cấp II) 22 Hình 1.10 Minh họa trạm châm hóa chất thiết bị bên trạm 23 Hình 1.11 Sơ đồ mặt quy trình thu hồi nước sau rửa lọc 24 Hình 1.12 Khu xử lý bùn thải 24 Hình 1.13 Minh họa hệ thống điều khiển SCADA 25 Hình 1.14 Quy trình cơng nghệ Nhà máy nước Thủ Đức 26 Hình 1.15 Quy trình xư lý nước nhà máy BOO Thủ Đức 27 Hình 1.16 Quy trình cơng nghệ Nhà máy nước Tân Hiệp 27 Hình 1.17 Quy trình xử lý nhà máy nước Tân Hiệp 28 Hình 1.18 Sơ đồ công nghệ xử lý nước Nhà máy nước Auchel 29 Hình 1.19 Cơng nghệ xử lý nước Nhà máy nước Biebesheim 29 Hình 2.1 Lược đồ phản ứng keo tụ 36 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 2.2 Tính chất keo tụ hợp chất polyme Fe(III) 38 Hình 2.3 Bể lắng bể lọc châm Polymer với liều lượng > 0.3ppm Ảnh chụp vào tháng 7/2017 42 Hình 2.4 Minh họa phần mềm phân tích liệu R 43 Hình 2.5 : Mơ tả mơ hình hồi quy tuyến tính đơn giản qua hai điểm 47 Hình 2.6 : Mơ tả phương pháp xác định tham số a b mơ hình hồi quy tuyến tính 48 Hình 3.1 Mẫu báo cáo số liệu ngày Ban Điều hành NMN 50 Hình 3.2 Chi tiết số liệu mẫu báo cáo Ban Điều hành NMN 51 Hình 3.3 Minh họa biến số liệu biên tập lại 52 Hình 3.4 Minh họa tập tin liệu dạng Excel 53 Hình 3.5 Biểu đồ diễn biến độ đục nước thô sông Đồng Nai qua năm 2016,2017,2018 55 Hình 3.6 Biểu đồ so sánh độ đục bình quân mùa năm 2016, 2017, 2018 55 Hình 3.7 Biểu đồ ma trận mối tương quan biến số tổng hợp năm: 2016, 2017, 2018 57 Hình 3.8 Biểu đồ ma trận mối tương quan biến số năm 2016 58 Hình 3.9 Biểu đồ ma trận mối tương quan biến số năm 2017 59 Hình 3.10 Biểu đồ ma trận mối tương quan biến số năm 2018 60 Hình 3.11 Các số liệu thể mức độ quan trọng biến số với biến số “FeCl3” năm 2016 62 Hình 3.12 Các số liệu thể mức độ quan trọng biến số với biến số “FeCl3” năm 2017 63 Hình 3.13 Các số liệu thể mức độ quan trọng biến số với biến số “FeCl3” năm 2018 64 Hình 3.14 Các số liệu thể mức độ quan trọng biến số với biến số “FeCl3” tổng hợp năm 2016, 2017, 2018 65 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.15 Biểu đồ lượng hóa chất FeCl3 bình qn qua năm 2016, 2017, 2018 66 Hình 3.16 Biểu đồ phân bố chuẩn biến số “FeCl3” 67 Hình 3.17 Kết phân tích phương sai biến số “FeCl3” 67 Hình 3.18 Kết phân tích hậu định phương sai theo phương pháp TuKey 68 Hình 3.19 Biếu đồ thể khác biệt biến số “FeCl3” 68 Hình 3.20 Kết phân tích hậu định phương sai biến số “FeCl3” theo Phương pháp Kruskal Wallis 69 Hình 3.21 Kết phân tích hậu định phương sai biến số “FeCl3” theo Phương pháp Kruskal Wallis 69 Hình 3.22 Kết phân tích phương sai biến số “FeCl3” theo mơ hình hóa 69 Hình 3.23 Kết phân tích phương sai biến số “FeCl3” theo mơ hình hóa 70 Hình 3.24 Biểu đồ độ đục nước thơ bình qn qua năm 2016, 2017, 2018 70 Hình 3.25 Biểu đồ phân bố chuẩn biến số “rw.ntu” 71 Hình 3.26 Kết phân tích phương sai biến số “rw.ntu” 71 Hình 3.27 Kết phân tích hậu định phương sai biến số “rw.ntu” theo phương pháp TuKey 71 Hình 3.28 Biếu đồ thể khác biệt biến số “rw.ntu” 72 Hình 3.29 Các ước số để xây dựng mơ hình hồi quy tuyến tính năm 2016 73 Hình 3.30 Các ước số để xây dựng mơ hình hồi quy tuyến tính có loại bớt biến số rw.ntu:seasons năm 2016 74 Hình 3.31 Biểu đồ mơ hình hồi quy tuyến tính năm 2016 75 Hình 3.32 Các ước số để xây dựng mơ hình hồi quy tuyến tính năm 2017 76 Hình 3.33 Các ước số để xây dựng mơ hình hồi quy tuyến tính có loại bớt biến số rw.ntu:seasons năm 2017 77 Hình 3.34 Biểu đồ mơ hình hồi quy tuyến tính năm 2017 78 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.35 Các ước số để xây dựng mơ hình hồi quy tuyến tính năm 2018 79 Hình 3.36 Các ước số để xây dựng mơ hình hồi quy tuyến tính có loại bớt biến số rw.ntu:seasons năm 2018 80 Hình 3.37.Biểu đồ mơ hình hồi quy tuyến tính năm 2018 81 Hình 3.38.Các ước số để xây dựng mơ hình hồi quy tuyến tính tổng hợp năm 82 Hình 3.39.Các ước số để xây dựng mơ hình hồ quy tuyến tính có loại bỏ biến số rw.ntu:seasons tổng hợp năm 83 Hình 3.40.Biểu đồ mơ hình tuyến tính tổng hợp năm 84 Hình 3.41 Biểu đồ mơ phần lượng dư mơ hình hồi quy tuyến tính 85 Hình 3.42 Biểu đồ phân bố chuẩn phần dư liệu tệp dat2018R.csv 86 Hình 3.43 Biểu đồ thể đường quan sát (màu hồng) “gần” với đường tuyến tính (màu xanh đứt khúc - đường kỳ vọng) 86 Hình 3.44 Biểu đồ xác định phương sai phần dư 87 Hình 3.45 Giá trị p = < 0.05 , nghĩa biến FeCl3 xem biến độc lập 87 Hình 3.46 Kết thấy khơng có giá trị ngoại vi đối tượng có id 222 có giá trị độ lệch chuẩn cao 88 Hình 3.47 Có đối tượng id: 221, 222, 223 có khả làm ảnh hưởng đến ước số mơ hình hồi quy tuyến tính 88 Hình 3.48 Các ước số xây dựng mơ hình hổi quy tuyến tính năm 2018 theo phương pháp BMA 89 Hình 3.49 Biếu đồ so sánh thực nghiệm mơ hình đề xuất 91 Hình 3.50 Các ước số mơ hình hồi quy tuyến tính thơng thường 92 Hình 3.51 Các ước số mơ hình hồi quy tuyến tính theo phương pháp BMA 92 Hình 3.52 Biểu đồ so sánh lần kết mơ hình hồi quy tuyến tính 94 Hình 3.53 Các ước số xây dựng mơ hình hồi quy tuyến tính thơng thường, lần 95 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 4/22/2017 11000 11.5 7.55 18 0.21 0.71 7.48 4/23/2017 12000 11.6 7.34 18 0.19 0.85 8.6 4/24/2017 12000 10.6 7.35 17 0.19 0.58 7.87 4/25/2017 11500 13.9 7.45 17 0.18 1.22 7.8 4/26/2017 10500 27.4 7.11 18 0.18 1.6 8.36 4/27/2017 11000 16 7.11 18 0.19 0.53 8.06 4/28/2017 12000 18.5 7.21 18 0.19 0.67 7.57 4/29/2017 11500 23.2 7.11 24 0.22 0.86 7.26 4/30/2017 11000 16.5 6.92 22 0.22 0.53 7.53 Phụ lục 1.2 Kết thử nghiệm Jartest tháng 6/2017 Nước thô Phèn sắt Ngày Q Polymer Kết sau làm Jartest NTU pH ppm ppm NTU pH (m3/h) 6/1/2017 12500 39.1 6.94 30 0.3 1.08 7.85 6/2/2017 12000 46.1 6.96 29 0.28 1.26 8.72 6/3/2017 11000 51.6 6.81 29 0.3 1.15 8.33 6/4/2017 12000 18.5 7.06 27 0.27 1.38 7.3 6/5/2017 10500 20.5 6.7 22 0.24 1.18 7.44 6/6/2017 11000 21.2 7.01 22 0.24 1.13 8.29 6/7/2017 11000 30.9 7.3 23 0.22 0.73 8.2 6/8/2017 10500 35.9 7.09 30 0.3 1.19 8.68 6/9/2017 12000 30.3 7.09 26 0.26 0.95 8.54 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 6/10/2017 13000 66.1 7.01 33 0.33 0.86 8.14 6/11/2017 10500 43.4 7.03 29 0.3 0.84 8.54 6/12/2017 10700 58.6 7.01 31 0.31 1.31 7.5 6/13/2017 12000 42.2 7.01 29 0.3 0.92 8.39 6/14/2017 11300 39.9 7.05 31 0.34 1.06 8.64 6/15/2017 9600 107 6.93 37 0.37 1.26 7.93 6/16/2017 10000 66 36 0.36 1.13 7.54 6/17/2017 12000 39.1 7.14 34 0.32 0.83 7.66 6/18/2017 11500 29.2 7.16 29 0.27 1.3 8.8 6/19/2017 11000 21 7.14 27 0.26 0.91 8.22 6/20/2017 10000 31 7.16 26 0.25 1.06 8.35 6/21/2017 12500 31.8 7.28 32 0.32 1.05 7.88 6/22/2017 9600 28.4 7.2 27 0.27 1.44 7.81 6/23/2017 10500 46.3 7.04 30 0.3 0.79 7.98 6/24/2017 12000 24.5 7.02 29 0.29 1.34 8.4 6/25/2017 11000 22.4 7.08 29 0.29 1.26 8.09 6/26/2017 10000 80.9 7.06 37 0.37 1.3 8.49 6/27/2017 10500 96.1 6.89 41 0.41 1.23 7.57 6/28/2017 11500 69 6.96 40 0.4 1.23 8.6 6/29/2017 10000 48.2 6.97 40 0.39 1.32 8.8 6/30/2017 11000 54.2 7.02 37 0.36 1.4 8.4 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Phụ lục 1.3 Kết thử nghiệm Jartest tháng 9/2017 Nước thô Ngày Q Phèn sắt Polymer Kết sau làm Jartest NTU pH ppm ppm NTU pH (m3/h) 9/1/2017 13300 87.1 6.95 42 0.4 1.26 7.48 9/2/2017 13000 147 7.26 44 0.44 1.02 7.5 9/3/2017 12700 86.3 7.08 42 0.42 1.19 7.62 9/4/2017 13000 65.4 7.12 36 0.35 1.26 9/5/2017 13300 47 7.04 28 0.28 0.86 8.1 9/6/2017 13000 48.5 6.95 31 0.3 0.71 7.06 9/7/2017 12500 35.4 7.06 26 0.38 1.12 7.87 9/8/2017 13000 36.4 6.92 33 0.32 1.07 7.79 9/9/2017 13000 36 6.89 35 0.34 0.93 7.69 9/10/2017 12600 39.5 6.97 36 0.35 1.08 7.57 9/11/2017 12700 42.4 6.77 38 0.38 0.76 7.53 9/12/2017 13000 53.1 6.84 39 0.38 1.36 7.89 9/13/2017 13300 102 6.78 40 0.4 1.44 7.8 9/14/2017 13200 101 7.01 44 0.42 1.26 7.65 9/15/2017 12600 122 7.11 43 0.4 1.39 7.53 9/16/2017 13200 76.5 7.02 39 0.36 0.97 7.61 9/17/2017 13300 57.6 6.84 35 0.33 1.3 7.5 9/18/2017 13100 47.3 6.92 34 0.31 1.02 7.5 9/19/2017 12800 47.8 6.86 35 0.33 1.13 8.08 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 9/20/2017 13000 74.6 7.05 38 0.36 1.05 7.59 9/21/2017 13300 81.8 6.86 43 0.41 1.33 7.53 9/22/2017 13300 85 44 0.41 0.82 7.8 9/23/2017 13000 90 7.03 44 0.43 1.18 7.62 9/24/2017 13000 70.6 6.98 43 0.41 1.09 7.56 9/25/2017 13000 97 6.77 39 0.37 1.17 7.68 9/26/2017 13200 63.4 6.83 40 0.39 1.5 7.8 9/27/2017 13200 76.9 7.09 45 0.42 1.35 7.81 9/28/2017 13300 64.6 7.02 41 0.39 1.07 7.58 9/29/2017 13300 92.6 6.82 42 0.4 1.16 8.1 9/30/2017 13300 64.4 6.81 38 0.36 1.19 7.42 Phụ lục 1.4 Kết thử nghiệm Jartest tháng 12/2017 Nước thô Ngày Phèn sắt Polymer Kết sau làm Jartest Q NTU pH ppm ppm NTU pH 12/1/2017 9019 16.4 21 0.18 1.35 7.98 12/2/2017 9019 25 6.92 23 0.21 8.85 12/4/2017 13680 27.7 6.98 20 0.18 1.18 8.12 12/5/2017 14349 34.9 6.9 20 0.15 0.82 8.49 12/6/2017 14174 29.7 6.82 16 0.21 1.21 7.96 12/7/2017 14192 39.9 6.98 24 0.15 1.05 8.32 12/8/2017 14116 37.9 6.69 23 0.21 0.95 8.86 12/9/2017 13099 36.9 6.91 21 0.15 1.34 8.2 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 12/10/2017 14072 36.4 6.92 17 0.15 1.31 8.13 12/11/2017 13310 29.2 20 0.21 1.23 8.45 12/12/2017 13680 24.1 7.01 19 0.24 1.04 8.02 12/13/2017 14005 31 20 0.16 1.05 8.62 12/14/2017 14069 30.5 6.97 14 0.21 1.13 7.93 12/15/2017 13619 28 6.87 18 0.15 0.96 8.4 12/16/2017 13702 28.9 7.02 17 0.24 1.17 8.11 12/17/2017 13680 23 19 0.15 1.16 7.65 12/18/2017 14319 22.5 7.02 17 0.18 1.36 8.28 12/19/2017 14422 39.5 7.01 20 0.12 1.12 8.5 12/20/2017 13622 26.9 7.04 19 0.15 1.02 8.16 12/21/2017 14231 29 7.02 18 0.15 1.03 9.14 12/22/2017 13000 24.4 6.8 15 0.24 1.34 8.36 12/23/2014 13570 33.5 6.84 16 0.15 1.55 8.46 12/24/2017 13685 25.3 6.94 17 0.24 1.14 8.76 12/25/2017 13684 27.6 6.86 16 0.21 1.19 9.76 12/26/2017 14288 16.9 7.23 15 0.12 0.93 8.8 12/27/2017 13971 32.5 7.2 16 0.12 0.96 8.52 12/28/2017 13700 17 7.03 10 0.15 1.24 7.85 12/29/2017 13544 18.7 16 0.21 0.96 8.43 12/30/2017 14291 14.4 7.1 18 0.18 0.65 7.77 12/31/2017 13705 17.1 7.07 15 0.15 1.6 8.4 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Phụ lục 1.5 Kết thử nghiệm Jartest tháng 3/2018 Nước thô Ngày Phèn sắt Polymer Kết sau làm Jartest Q NTU pH ppm ppm NTU pH 3/2/2018 12800 25.4 6.79 13 0.09 0.92 8.26 3/3/2018 12439 23.5 6.91 10 0.09 1.04 8.58 3/4/2018 12885 32.5 6.87 15 0.12 0.98 8.32 3/5/2018 13126 33.5 6.84 12 0.11 0.9 8.2 3/6/2018 12774 47.6 14 0.06 1.08 8.12 3/7/2018 12806 30 6.8 13 0.09 1.35 7.62 3/13/2018 12945 15 7.17 12 0.11 1.28 8.04 3/17/2018 12663 15.3 7.07 10 0.06 1.08 8.67 3/18/2018 13072 22.9 6.95 17 0.21 0.93 8.39 03/19/2018 12601 20 16 0.15 1.02 7.56 03/20/2018 12700 37.2 7.03 23 0.21 0.95 8.86 03/21/2018 12600 37.9 6.85 14 0.1 1.3 7.57 03/22/2018 13072 27.7 6.99 12 0.12 1.29 8.21 03/23/2018 12548 22.5 7.1 17 0.21 0.93 8.39 03/24/2018 12863 14.8 6.89 15 0.12 1.25 8.82 03/25/2018 13000 16.7 7.16 16 0.18 1.46 8.4 03/26/2018 13240 13.6 6.42 20 0.15 0.84 8.81 03/27/2018 12099 19.3 7.2 12 0.09 1.27 7.98 03/28/2018 11657 14.6 6.98 15 0.12 1.25 8.82 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 03/29/2018 13519 14.3 7.11 5.5 0.035 0.755 7.378 03/30/2018 12443 20.1 13 0.06 0.84 8.32 03/31/2018 12443 22 6.9 13 0.06 0.84 8.32 Phụ lục 1.6 Kết thử nghiệm Jartest tháng tháng 6/2018 Nước thô Ngày Q Phèn Polymer sắt Kết sau làm Jartest NTU pH ppm ppm NTU pH 5/11/2018 13769 116 7.15 28 0.26 1.83 7.38 5/12/2018 13000 109 7.30 25 0.2 1.02 6.51 5/13/2018 13703 124 7.23 29 0.23 0.72 5.91 5/14/2018 13000 140 7.45 30 0.24 1.45 7.13 5/15/2018 13904 156 7.34 31 0.25 1.04 7.77 5/16/2018 13000 102 7.34 23 0.21 1.5 6.13 5/17/2018 13934 167 7.23 34 0.21 0.72 9.23 5/18/2018 13000 110 7.30 21 0.18 1.27 9.14 5/19/2018 13425 196 7.37 43 0.3 1.11 9.77 5/20/2018 13256 121 7.48 29 0.24 6.74 5/21/2018 13419 136 7.23 31 0.20 1.07 8.56 5/22/2018 13000 158 7.09 30 0.23 1.82 6.42 5/23/2018 12415 136 6.57 30 0.20 1.32 7.87 5/24/2018 12952 160 7.43 31 0.15 1.2 8.62 5/25/2018 12854 147 7.47 29 0.22 1.33 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 5/26/2018 12361 129 7.33 28 0.13 0.94 9.25 5/27/2018 12319 177 7.1 27 0.15 1.37 9.57 5/28/2018 12693 170 7.2 32 0.2 1.29 7.42 5/29/2018 12488 200 6.2 44 0.31 0.99 9.20 5/30/2018 13500 210 7.1 33 0.29 1.29 7.41 5/31/2018 12500 91.9 7.1 20 0.23 1.67 8.45 6/1/2018 13090 198 7.2 33 0.28 0.9 6.32 6/2/2018 13059 165 7.2 30 0.24 1.55 6.14 6/3/2018 13500 114 6.8 21 0.18 1.15 6.24 6/4/2018 13000 99.7 6.8 25 0.19 1.6 6.64 6/5/2018 12177 118 6.9 34 0.24 1.37 8.62 6/6/2018 12202 188 7.1 31 0.3 1.61 6.92 6/7/2018 13000 195 7.1 31 0.29 1.93 8.23 6/8/2018 12741 216 7.2 46 0.41 1.09 8.25 6/9/2018 13050 105 7.1 22 0.18 1.14 6.91 6/10/2018 13000 243 7.4 44 0.43 1.36 9.34 6/11/2018 13000 253 7.1 44 0.43 1.12 9.65 6/12/2018 12598 210 7.3 29 0.36 1.39 8.54 6/13/2018 13000 202 7.4 27 0.38 1.31 9.91 6/14/2018 13404 183 7.4 44 0.32 0.93 8.31 6/15/2018 12700 228 7.4 29 0.39 1.16 8.85 6/16/2018 12800 82.1 18 0.15 6.76 6/17/2018 13099 94.5 7.1 20 0.15 1.34 8.26 Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 6/18/2018 13000 240 7.2 43 0.43 0.93 8.94 Nguồn : Sổ tay thí nghiệm Jartest Nhà máy nước trực thuộc SWIC Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 PHỤ LỤC CÁC CÂU LỆNH DÙNG TRONG PHÂN TÍCH DỮ LIỆU BẰNG R Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Nội dung câu lệnh > setwd("E:cuong_dataR/luanvanR") > dat2016 = read.csv("dat2016R.csv", header = T) > dat2017 = read.csv("dat2017R.csv", header = T) > dat2018 = read.csv("dat2018R.csv", header = T) > dat.th = read.csv("tonghopR.csv", header = T) > attach(dat.th) > library(sciplot) > lineplot.CI(months,rw.ntu,years, ylab="Do duc (NTU)", xlab="Thang", ylim=c(0,150), legend=T, col=c("red","blue", "green"), main="Bieu do duc nuoc tho tu nam 2016 den 2018") > library(ggplot2) > p = ggplot(data = dat.th, aes(x=years, y=rw.ntu, fill=years)) > p + stat_summary(fun.y="mean", geom="bar") + facet_wrap(~seasons) + theme_bw()+ xlab("NAM")+ ylab("DO DUC (NTU)") > d.th=cbind(rw.ntu, rw.pH, rw.cap, seasons, FeCl3, Plm, Lime) > library(psych) > pairs.panels(d.th) > detach(dat.th) > attach(dat2016) > d.2016=cbind(rw.ntu, rw.pH, rw.cap, seasons, FeCl3, Plm, Lime) > pairs.panels(d.2016) > detach(dat2016) > attach(dat2017) > d.2017=cbind(rw.ntu, rw.pH, rw.cap, seasons, FeCl3, Plm, Lime) > pairs.panels(d.2017) > detach(dat2017) > attach(dat2018) > d.2018=cbind(rw.ntu, rw.pH, rw.cap, seasons, FeCl3, Plm, Lime) > pairs.panels(d.2018) > library(relaimpo) > attach(dat2016) > model2016 = lm(FeCl3 ~ rw.ntu + seasons + Plm + Lime + rw.pH) > met = calc.relimp(model2016, type = c(“lmg”)) > met > library(relaimpo) > attach(dat2017) Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hiển thị kết Các câu lệnh khai báo trước thực phân tích liệu Hình 3.5 Hình 3.6 Hình 3.7 Hình 3.8 Hình 3.9 Hình 3.10 Hình 3.11 Hình 3.12 > model2017 = lm(FeCl3 ~ rw.ntu + seasons + Plm + Lime + rw.pH) > met = calc.relimp(model2017, type = c(“lmg”)) > met > library(relaimpo) > attach(dat2018) > model2018 = lm(FeCl3 ~ rw.ntu + seasons + Plm + Lime + rw.pH) > met.2018 = calc.relimp(model2018, type = c(“lmg”)) > met.2018 > library(relaimpo) > attach(dat.th) > model.th = lm(FeCl3 ~ rw.ntu + seasons + Plm + Lime + rw.pH) > met = calc.relimp(model.th, type = c(“lmg”)) > met > attach(dat2016) > fe.2016 = na.omit(FeCl3) > detach(dat2016) > attach(dat2017) > fe.2017 = na.omit(FeCl3) > detach(dat2017) > attach(dat2018) > fe.2018 = na.omit(FeCl3) > feth = c(fe.2016, fe.2017, fe.2018) > group = c(rep("nam2016",365), rep("nam2017",363), rep("nam2018",365)) > dat = data.frame(group, feth) > boxplot(feth ~ group, col="brown", ylim = c(0,200), ylab="g/m3", xlab="Nam") > qqnorm(feth) > qqline(feth, col = "red") > av = aov(feth ~ group) > summary(av) > TukeyHSD(av) > plot(TukeyHSD(av)) > library(pgirmess) > kruskal.test(fe.th ~ as.factor(group)) > kruskalmc(fe.th, group) > mod = lm(fe.th ~ group) > summary(mod) > anova(mod) > attach(dat2016) Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.13 Hình 3.14 Hình 3.15 Hình 3.16 Hình 3.17 Hình 3.18 Hình 3.19 Hình 3.20 Hình 3.21 Hình 3.22 Hình 3.23 Hình 3.24 > rwntu.2016 = na.omit(rw.ntu) > detach(dat2016) > attach(dat2017) > rwntu.2017 = na.omit(rw.ntu) > detach(dat2017) > attach(dat2018) > rwntu.2018 = na.omit(rw.ntu) > rwntu.th = c(rwntu.2016, rwntu.2017, rwntu.2018) > rwntu.group = c(rep("Nam.2016",304), rep("Nam.2017",362), rep("Nam.2018",294)) > dat.rwntu = data.frame(rwntu.group, rwntu.th) > boxplot(rwntu.th ~ rwntu.group, col="blue", ylim = c(0,200), ylab="Do duc (NTU)", xlab="Nam") > qqnorm(rwntu.th) > qqline(rwntu.th) > av.rwntu = aov(rwntu.th ~ rwntu.group) > summary(av.rwntu) > TukeyHSD(av.rwntu) > plot(TukeyHSD(av.rwntu), ordered=T) > summary(lm(FeCl3 ~ rw.ntu + seasons + rw.ntu:seasons, data = dat2016)) > summary(lm(FeCl3 ~ rw.ntu + seasons , data = dat2016)) > library(ggplot2) > qplot(x=rw.ntu, y=FeCl3, col=seasons) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ xlab("Do duc (NTU)")+ ylab("Lieu luong FeCl3(g/m3)")+ggtitle("Bieu phan bo luong FeCl3 nam 2016") > summary(lm(FeCl3 ~ rw.ntu + seasons + rw.ntu:seasons, data = dat2017)) > summary(lm(FeCl3 ~ rw.ntu + seasons , data = dat2017)) > detach(dat2016) > attach(dat2017) > qplot(x=rw.ntu, y=FeCl3, col=seasons) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ xlab("Do duc (NTU)")+ ylab("Lieu luong FeCl3(g/m3)")+ggtitle("Bieu phan bo luong FeCl3 nam 2017") > summary(lm(FeCl3 ~ rw.ntu + seasons + rw.ntu:seasons, data = dat2018)) > summary(lm(FeCl3 ~ rw.ntu + seasons , data = dat2018)) > detach(dat2017) > attach(dat2018) Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.25 Hình 3.26 Hình 3.27 Hình 3.28 Hình 3.29 Hình 3.30 Hình 3.31 Hình 3.32 Hình 3.33 Hình 3.34 Hình 3.35 Hình 3.36 Hình 3.37 > qplot(x=rw.ntu, y=FeCl3, col=seasons) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ xlab("Do duc (NTU)")+ ylab("Lieu luong FeCl3(g/m3)")+ggtitle("Bieu phan bo luong FeCl3 nam 2018") > summary(lm(FeCl3 ~ rw.ntu + seasons + rw.ntu:seasons, data = dat.th)) > summary(lm(FeCl3 ~ rw.ntu + seasons , data = dat2018)) > detach(dat2018) > attach(dat.th) > qplot(x=rw.ntu, y=FeCl3, col=seasons) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ xlab("Do duc (NTU)")+ ylab("Lieu luong FeCl3(g/m3)")+ggtitle("Bieu phan bo luong FeCl3 cua nam 2016,2017,2018") > detach(dat.th) > attach(dat2018) > h2 = lm(FeCl3~rw.ntu+seasons) > summary(h2) > res2 = resid(h2) > hist(res2) > require(car) > crPlots(h2) > ncvTest(h2) > spreadLevelPlot(h2) > spreadLevelPlot(h2, smooth = F) > durbinWatsonTest(h2) > outlierTest(h2) > influencePlot(h2,id.method="identify",main = "Influence Plot", sub = "Circle size is proportial to Cook's distance") > detach(dat.th) > attach(dat2018) > library(BMA) > X = cbind(rw.ntu, seasons) > Y = dat2018$FeCl3 > search = bicreg(X, Y, strict = F, OR=20) > summary(search) > setwd("E:/cuong_dataR/luanvanR") > complm = read.csv("compare_lmR.csv", header = T) > library(ggplot2) > attach(complm) > qplot(x=rw.ntu, y=FeCl3, col=Models) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.38 Hình 3.39 Hình 3.40 Hình 3.42 Hình 3.43 Hình 3.44 Hình 3.45 Hình 3.46 Hình 3.47 Hình 3.48 Hình 3.49 xlab("Do duc(NTU)")+ ylab("Luong FeCl3 (g/m3)")+ ggtitle("Bieu so sanh thuc nghiem va mo hinh") > models = lm(FeCl3~rw.ntu+Models) > summary(models) > library(BMA) > X = cbind(rw.ntu,Models) > Y = complm$FeCl3 > search = bicreg(X,Y,strict = F,OR = 20) > summary(search) > setwd("E:/cuong_dataR/luanvanR") > complm = read.csv("compare2_lmR.csv", header = T) > library(ggplot2) > attach(complm) > qplot(x=rw.ntu, y=FeCl3, col=Models) + stat_smooth(method="lm")+ theme_bw()+ theme_classic()+ xlab("Do duc(NTU)")+ ylab("Luong FeCl3 (g/m3)")+ ggtitle("Bieu so sanh thuc nghiem va mo hinh") > models = lm(FeCl3~rw.ntu+Models) > summary(models) > library(BMA) > X = cbind(rw.ntu,Models) > Y = complm$FeCl3 > search = bicreg(X,Y,strict = F,OR = 20) > summary(search) Học viên : Hồ Thanh Cường , Lớp : CH25CTN11-CS2 Hình 3.50 Hình 3.51 Hình 3.52 Hình 3.53 Hình 3.54 ... Cơ quan công tác : Công ty CP Đầu tư Kinh doanh nước Sài Gòn Tác giả đề tài : Nghiên cứu giải pháp tối ưu hóa phèn sắt (FeCl3) Polymer xử lý nước Công ty CP Đầu tư Kinh doanh nước Sài Gòn Học... châm hóa chất phèn sắt (FeCl3) Polymer áp dụng Công ty - Nghiên cứu, đề xuất giải pháp mơ hình tối ưu hóa chất phèn sắt (FeCl 3) Polymer xử lý nước Công ty - Ứng dụng công nghệ thông tin quản lý, ... tài liệu giới thiệu Cơng ty CP Đầu tư Kinh doanh nước Sài Gòn, quy mơ quy trình cơng nghệ Nhà máy nước trực thuộc Công ty sau : Công ty Cổ phần Đầu tư Kinh doanh nước Sài Gòn (gọi tắt SWIC) thành