Nghiên cứu sản xuất sunphat amon từ chất thải gyps của nhà máy dap số 1 vinachem tại khu công nghiệp đình vũ hải phòng

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI - HỌ VÀ TÊN: LƯU NGỌC VĨNH ĐỀ TÀI: NGHIÊN CỨU SẢN XUẤT SUNPHÁT AMÔN TỪ CHẤT THẢI GYPS CỦA NHÀ MÁY DAP SỐ – VINACHEM TẠI KCN ĐÌNH VŨ HẢI PHỊNG CHUN NGÀNH: KỸ THUẬT HĨA HỌC LUẬN VĂN THẠC SĨ KỸ THUẬT NGƯỜI HƯỚNG DẪN 1: TS NGUYỄN HUY DŨNG NGƯỜI HƯỚNG DẪN 2: PGS.TS PHẠM THANH HUYỀN Hà Nội – Năm 2013 LỜI CẢM ƠN Để hồn thành chương trình cao học luận văn này, nhận hướng dẫn, giúp đỡ góp ý nhiệt tình q thầy trường Đại học Bách Khoa Hà Nội Trước hết, xin chân thành cảm ơn đến quý thầy cô trường Đại học Bách Khoa Hà Nội, đặc biệt thầy tận tình dạy bảo cho tơi suốt thời gian học tập trường Tôi xin gửi lời biết ơn sâu sắc đến TS Nguyễn Huy Dũng PGS.TS Phạm Thanh Huyền dành nhiều thời gian tâm huyết hướng dẫn nghiên cứu giúp hồn thành luận văn tốt nghiệp Mặc dù tơi cố gắng hoàn thiện luận văn tất nhiệt tình lực mình, nhiên khơng thể tránh khỏi thiếu sót, mong nhận đóng góp q báu thầy bạn Hà Nội, tháng năm 2013 Học viên Lưu Ngọc Vĩnh LỜI CAM ĐOAN Tôi xin cam đoan toàn nội dung số liệu luận văn tơi nhóm nghiên cứu, khảo sát thực Tôi xin cam đoan giúp đỡ cho việc thực luận văn cảm ơn thơng tin trích dẫn luận văn rõ nguồn gốc Học viên thực luận văn Lưu Ngọc Vĩnh Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” DANH MỤC CÁC KÝ HIỆU CHỮ VIẾT TẮT SA: Amôn Sunphát DAP: Di-amôn phốtphát MAP: Amôn phốtphát NPK: Nitơ – phốtpho - kali MMA: Methyl methacryalate CAN: Acrylonitrile PTN: Phịng thí nghiệm Gyps: Gypsum Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” DANH MỤC CÁC KÝ HIỆU TỐN HỌC Nồng độ mol cấu tử j (kmol/m3) D Hệ số khuếch tán (m2/h) Hệ số cấp khối (m/h) Dòng phát sinh (kmol/m3h) Vận tốc dòng đối lưu (m/h) Hệ số tỷ lượng Tốc độ phản ứng Hằng số tốc độ phản ứng (s-1) Bậc phản ứng m Khối lượng sản phẩm thu sau trình phản ứng yu Giá trị đo lần lặp thứ u S2 i Phương sai chung S2 ts Phương sai tái sinh chung fi Số bậc tự mẫu f ts Số bậc tự phương sai chung S ts Sai số tái sinh yi Giá trị thực nghiệm i t(p,f ) Giá trị tính từ hàm hồi quy lý thuyết Tiêu chuẩn Student tra bảng mức có nghĩa p bậc tự lặp f Sb Độ lệch chuẩn phân bố b y0 a Giá trị thực nghiệm lặp thứ a l Giá trị trung bình cộng thực nghiệm lặp Chuẩn số fisher (mức có nghĩa p, bậc tự dư f = N-l, bậc tự lặp f = m-1) Hệ số có nghĩa mơ tả thống kê T Nhiệt độ tiến hành phản ứng F(p,f ,f ) M1 Khối lượng mẫu Gyps nghiên cứu chứa 74% thạch cao; [g] M 1' Khối lượng thạch cao mẫu Gyps; [g] Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” M2 Khối lượng Amôn Sunphát thu theo phản ứng; [g] M 2' Khối lượng Amôn Sunphát thu thực tế; [g] ω Tốc độ khuấy trộn; [vòng/phút] t Thời gian tiến hành phản ứng (phút) η Hiệu suất phản ứng r Kích thước tinh thể (nm) Hàm mục tiêu (hàm hiệu suất) Z1 Biến thực tốc độ khuấy trộn cấu tử nhiệt độ phản ứng Z2 Biến thực thời gian phản ứng Z3 Biến thực nhiệt độ phản ứng x1 Biến mã hóa tốc độ khuấy trộn cấu tử nhiệt độ phản ứng x2 Biến mã hóa thời gian phản ứng x3 Biến mã hóa nhiệt độ phản ứng Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phòng” DANH MỤC CÁC BẢNG Bảng 1: So sánh thành phần Gyps nhà máy DAP giới DAP Đình Vũ Hải Phịng Bảng 2: Thành phần Gyps nhà máy DAP Bảng 3: Các đặc tính hóa lý Amôn Sunphát tinh thể Bảng 4: Tiêu thụ phân bón Amơn Sunphát theo khu vực từ 2002 đến 2011 Bảng 5: Sản xuất phân bón Amơn Sunphát theo khu vực từ 2002 đến 2011 Bảng 6: Tổng lượng phân bón Amơn Sunphát nhập qua năm Bảng 7: Ảnh hưởng nhiệt độ đến hiệu suất phản ứng Bảng 8: Ảnh hưởng tốc độ khuấy trộn đến hiệu suất phản ứng Bảng 9: Ảnh hưởng thời gian phản ứng đến hiệu suất phản ứng Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” DANH MỤC CÁC HÌNH Hình 1: Các nhà sản xuất Gyps giới năm 2011 Hình 2: Mẫu Gyps chưa qua xử lý Hình 3: Mẫu Gyps xử lý Hình 4: Sản lượng Amơn Sunphát phân theo nguồn gốc năm 2010 - 2011 Hình 5: Sơ đồ bố trí thí nghiệm chuyển hóa Gyps Hình 6: Khảo sát ảnh hưởng nhiệt độ đến hiệu suất phản ứng Hình 7: Khảo sát ảnh hưởng tốc độ khuấy đến hiệu suất phản ứng Hình 8: Khảo sát ảnh hưởng thời gian phản ứng đến hiệu suất phản ứng Hình 9: Ảnh hưởng nhiệt độ đến hiệu suất phản ứng Hình 10: Ảnh hưởng tốc độ khuấy đến hiệu suất phản ứng Hình 11: Ảnh hưởng thời gian phản ứng đến hiệu suất phản ứng Hình 12: Mẫu SA nghiên cứu PTN Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phòng” MỤC LỤC DANH MỤC CÁC KÝ HIỆU CHỮ VIẾT TẮT DANH MỤC CÁC KÝ HIỆU TOÁN HỌC DANH MỤC CÁC BẢNG MỤC LỤC .6 MỞ ĐẦU B.1 Tính cấp thiết đề tài B.2 Mục tiêu nghiên cứu luận văn: 10 B.3 Cách tiếp cận phương pháp nghiên cứu 10 B.4 Nội dung nghiên cứu 10 B5 Ý nghĩa khoa học, thực tiễn đóng góp luận văn 11 B6 Bố cục luận văn 11 CHƯƠNG 12 TỔNG QUAN VỀ CHẤT THẢI GYPS VÀ SẢN PHẨM AMÔN SUNPHÁT 12 I TỔNG QUAN VỀ CHẤT THẢI GYPS 12 1.1 Định nghĩa phân loại 12 1.1.1 Định nghĩa 12 1.1.2 Phân loại 12 1.2 Các ứng dụng quan trọng Gyps 13 1.2.1 Vật liệu xây dựng .13 1.2.2 Sản xuất phân bón Amơn Sunphát 13 1.2.3 Sản xuất xi măng Portland .14 1.3 Tình hình sản xuất tiêu thụ Gyps giới Việt Nam 14 1.3.1 Tình hình sản xuất Gyps giới 14 1.3.2 Tình hình sản xuất, tiêu thụ sản phẩm Gyps Việt Nam 16 1.4 Thành phần hợp chất lưu huỳnh chất thải Gyps – nhà máy DAP Đình Vũ – Hải Phòng .18 1.4.1 Sơ đồ lưu trình cơng nghệ q trình sản xuất DAP 18 1.4.2 Cơng đoạn sản xuất Axít Sunphuric từ lưu huỳnh 19 1.4.3 Cơng đoạn sản xuất Axít Phốtphoríc từ Axít Sunphuríc 20 1.4.4 Thành phần Gyps thải sau cơng đoạn sản xuất Axít Phốtphoric 21 Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” 1.5 Nhận xét 23 II TỔNG QUAN VỀ AMÔN SUNPHÁT 24 2.1 Đặc điểm Amôn Sunphát 24 2.2 Tình hình sản xuất, tiêu thụ Amơn Sunphát (SA) giới .25 2.3 Tình hình sản xuất, tiêu thụ SA thị trường Việt Nam 28 2.4 Công nghệ sản xuất sử dụng nguyên liệu từ Gyps 30 2.5 Quy trình sản xuất Amơn Sunphát từ bã thải Gyps giới 31 2.5.1 Công đoạn chuẩn bị nguyên liệu 33 2.5.2 Công đoạn Cácbonát hóa dung dịch Amơniắc .33 2.5.3 Công đoạn phản ứng trao đổi 34 2.5.4 Công đoạn bốc kết tinh 34 2.5.5 Công đoạn ly tâm, sấy, phân loại, đóng bao lưu kho .35 2.6 Nhận xét 35 CHƯƠNG - PHƯƠNG PHÁP NGHIÊN CỨU 37 CƠ SỞ LÝ THUYẾT .37 1.1 Quy hoạch thực nghiệm mơ tả q trình hóa học 37 1.1.1 Thiết lập mô tả thống kê 37 1.1.2 Các phương pháp kế hoạch hoá thực nghiệm cực trị chủ yếu 39 1.1.3 Xác định giá trị tối ưu hàm mục tiêu .40 1.2 Phương pháp phân tích thành phần Gyps .41 1.2.1 Phương pháp hóa học ướt 41 1.2.2 Phương pháp nhiệt trọng lượng .43 1.2.3 Phương pháp nhiễu xạ tia X .43 1.3 Phương pháp phân tích Amơni Sunphát 43 1.3.1 Hóa chất phân tích 43 1.3.2 Tiến hành 43 1.3.3 Tính tốn kết .43 1.3.4 Xác định hàm lượng cặn không tan nước .44 CƠ SỞ THỰC NGHIỆM 44 2.1 Quy trình sản xuất phịng thí nghiệm .44 2.2 Mô hình thực nghiệm quy mơ phịng thí nghiệm 46 2.2.1 Dụng cụ thiết bị thí nghiệm .46 Học viên: Lưu Ngọc Vĩnh Trang Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” 3.3.7 Triển vọng thị trường SA Việt Nam: Qua phân tích tình hình sản xuất, tiêu thụ dự báo cung cầu cho sản phẩm Gyps, phân bón SA giới Việt Nam, đồng thời có xét đến ảnh hưởng biến động cung cầu, giá nguyên liệu Amôniắc, báo cáo thấy việc đầu tư xây dựng nhà máy sản xuất phân bón SA, nhằm sử dụng hiệu hợp chất lưu huỳnh bã thải Gyps kèm theo dây chuyền tinh chế Gyps cần thiết vì: - Chủ động nguồn cung phân bón SA phục vụ cho nhu cầu nước - Giúp tiết kiệm lượng lớn ngoại tệ dành cho việc nhập phân bón SA - Góp phần xử lý bã thải Gyps nhà máy sản xuất phân bón DAP Đình Vũ – Hải Phịng, vừa giúp tăng nguồn thu cho nhà máy tạo thêm công ăn việc làm cho người dân khu vực 3.3.8 Địa điểm áp dụng Nếu dự án triển khai quy mô công nghiệp, địa điểm dự kiến xây dựng nhà máy lựa chọn khu vực lân cận mặt nhà máy sản xuất DAP Đình Vũ Hải Phịng lý sau: - Có thể tận dụng nguồn nguyên liệu chỗ, tiết kiệm chi phí vận chuyển có khả giới hóa nguồn nguyên liệu sản xuất + Gyps: Lấy bãi thải Nhà máy DAP gần + Khí NH : Cung cấp từ nhà máy sản xuất NH vận chuyển đường đường thủy kho chứa nhà máy + Khí CO : Nguyên liệu cung cấp từ nhiều nguồn như: nhà máy sản xuất đạm nhà máy tổng hợp etannol vận chuyển đường đường thủy kho chứa nhà máy Khí CO xem xét tận dụng từ nguồn khí thải lò đốt than đem xử lý tách lấy CO giống dây chuyền hoạt động nhà máy đạm Phú Mỹ - Có thể sử dụng số hệ thống phụ trợ nhà máy DAP nhằm tiết kiệm chi phí đầu tư - Phù hợp với quy hoạch tập trung khu cơng nghiệp sản xuất hóa chất, phân bón Dễ dàng quản lý vận hành sản xuất môi trường Học viên: Lưu Ngọc Vĩnh Trang 77 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” Trên thực tế việc sản xuất SA từ Gyps quy mơ cơng nghiệp lựa chọn với cơng suất 200.000 tấn/năm Tuy nhiên, trước sản xuất quy mô công nghiệp cần lập dự án sản xuất thử nghiệm (có thể lựa chọn cơng suất sản xuất SA 2.000kg/h) để đánh giá tính khả thi đề tài áp dụng sản xuất với điều kiện môi trường thực tế Học viên: Lưu Ngọc Vĩnh Trang 78 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phòng” KẾT LUẬN *********** Đề tài nghiên cứu thăm dò khả thu hồi sử dụng hiệu hợp chất lưu huỳnh Gyps phản ứng tạo thành Amơn Sunphát, từ xác định thơng số cơng nghệ ảnh hưởng đến hiệu suất phản ứng tạo thành Amơn Sunphát - Đề tài có tính thực tiễn cao sản phẩm Amôn Sunphát Việt Nam gần hồn tồn phải nhập khẩu, việc nghiên cứu cơng nghệ sử dụng nguyên liệu chất thải Gyps thải để sản xuất Amôn Sunphát thành công bước đột phá để xây dựng nhà máy sản xuất Amơn Sunphát cho tập đồn Hóa chất Việt Nam - Mặt khác, chuyển sang sản xuất quy mơ cơng nghiệp có ý nghĩa cao sử dụng hiệu nguồn thải Gyps để sản xuất sản phẩm có ích Tuy nhiên, tiến hành áp dụng từ nghiên cứu phịng thí nghiệm sang cơng nghiệp cần phải qua mơ hình pilot để kiểm nghiệm lại kết nghiên cứu đạt trước đưa vào sản xuất quy mô công nghiệp Do đề tài cần có bước nghiên cứu tiếp quy mơ pilot để đánh giá xác thơng số ảnh hưởng đến q trình phản ứng tạo thành SA trước chuyển sang quy mô công nghiệp  Phạm vi nghiên cứu đề tài: - Trên sở thu thập tài liệu tham khảo nước kết hợp với khảo sát sơ trình nghiên cứu thực nghiệm, hướng nghiên cứu đề tài tập vào nội dung sau: - Khảo sát nhu cầu sản xuất tiêu thụ Gyps, SA giới Việt Nam - Khảo sát, xác định thành phần Gyps thải nhà máy DAP Đình Vũ – Hải Phòng, đánh giá thành phần nguyên liệu Gyps cho việc nghiên cứu sản xuất SA - Đưa quy trình sản xuất SA từ Gyps thải Amơn Cácbonát - Xây dựng mơ hình thí nghiệm khảo sát thông số công nghệ nhiệt độ, tốc độ khuấy trộn, thời gian phản ứng ảnh hưởng đến hiệu suất tạo thành SA - Đưa phương trình hồi quy thông số công nghệ ảnh hưởng đến hiệu suất phản ứng, xác định thông số công nghệ cho hiệu suất cao, hiệu kinh tế cao áp dụng vào quy mô công nghiệp - Đánh giá khả áp dụng vào sản xuất thực tế quy mô công nghiệp Học viên: Lưu Ngọc Vĩnh Trang 79 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng”  Kết đề tài thực hiện: Lựa chọn đưa phương pháp thích hợp, tiên tiến, cho hiệu suất cao, phù hợp với điều kiện nghiên cứu nước Đã xác định hàm lượng hợp chất lưu huỳnh hữu ích Gyps thải nhà máy DAP Hải Phòng dạng CaSO 2H O chiếm 75%, tạp chất khác 25%, phù hợp với sản xuất Amôn Sunphát Đã xây dựng phương trình hồi quy thực nghiệm ảnh hưởng thông số công nghệ đến hiệu suất phản ứng tạo SA Đã xác định quy trình cơng nghệ phịng thí nghiệm để điều chế Amơn Sunphát từ Gyps thải với nguyên liệu Gyps thải Nhà máy DAP Đình Vũ tinh chế qua khâu : nghiền, rửa, sấy khô để tách nước, nguyên liệu thứ Amôn Cácbonát (25%); khâu phản ứng nghiên cứu hiệu suất phản ứng theo số thơng số q trình nhiệt độ, tốc độ khuấy, thời gian phản ứng, có kết Hiệu suất phản ứng tạo thành Amôn Sunphát 93% điều kiện: + Nhiệt độ phản ứng : 45 ± 0C + Tốc độ khuấy : 45 ± vòng/phút + Thời gian phản ứng : 2,5 ± 0,1 (giờ)  Hạn chế đề tài: - Chưa tiến hành nghiên cứu ảnh hưởng thành phần tạp chất đến hiệu suất phản ứng  Kiến nghị hướng mở rộng: - Nghiên cứu q trình phản ứng tạo thành Amơn Sunphát chịu ảnh hưởng nhiều trình phức tạp Từ trình phản ứng Amơniắc với Cácboníc tạo thành Amơn Cácbonát đến q trình phản ứng Amơn Cácbonát với bã thải Gyps tạo thành Amơn Sunphát - Do kinh phí cịn hạn hẹp, thời gian có hạn, việc nghiên cứu cịn nhiều hạn chế, nghiên cứu quy mô phịng thí nghiệm cần cấp thêm kinh phí thời gian để nghiên cứu thông số ảnh hưởng rộng hơn, phản ánh đầy đủ chuyển mơ hình từ phịng thí nghiệm sang quy mơ cơng nghiệp Học viên: Lưu Ngọc Vĩnh Trang 80 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” TÀI LIỆU THAM KHẢO Tiếng Việt [1] Nguyễn Minh Tuyển, Phạm Văn Thiêm Kĩ thuật hệ thống cơng nghệ hóa học, tập Nhà xuất khoa học kĩ thuật, Hà Nội, 1997 [2] Nguyễn Minh Tuyển – Quy hoạch thực nghiệm- Nhà xuất khoa học kỹ thuật Hà Nội, 2005 [3] Nguyễn Bin, Các q trình, thiết bị cơng nghệ hóa chất thực phẩm, Tập 5, Nhà xuất Khoa học & Kỹ thuật Hà nội [4] Đào Trần Cao, Giáo trình Vật lý Chất rắn sở, Viện Khoa Học Vật liệu, TTKHTN & CNQG, Hà nội, 10/2001, trang II-1 đến II-17 [5] Chu Thị Thơm, Kỹ thuật sản xuất chế biến phân bón – Nhà xuất lao động 2006 [6] Thuyết minh dự án đầu tư xây dựng nhà máy sản xuất DAP Đình Vũ - Hải Phịng cơng suất 330.000 tấn/năm Tiếng Anh [7] Arab Gulf, Scient Res Math.Phys A6 (l) PP77-83(1988),"Chemical Processing of Phốtphátghypsum for Production of Ammonium Sulfate Fertilizer by Selim F Estefan [8] Wikkiam C Burnett, Michael S., Carter H." Behavior of Radionuclides During Ammonocarbon of PhốtphátGyps", Florida State Univ [9] Gopinath, N.D., 1968 Disposal or use of Gyps in production of ammonium sulfate In: Phosphoric Acid, vol l, Part II, (ed A.V Slack), Marcel Dekker, New York, 541-566 Chen, Jing-Yu, (1990) The trial production results from demonstration plant for production of ammonium sulfate from phosphogypsum Phosphogypsum Conference Proceedings, Orlando, [10] Florida Department of Environmental Protection (FDEP) management requirements for phosphogypsum stack systems (17-673.100-900, Florida Administrative Code) [11] CPCB under the publication series HAZWAMS/17/2000-01 of Hazardous Waste Management i.e criteria for hazardous waste landfills [12] Guidelines for Development of Basic Infrastructure of Gypsum Handling and Storage on Railway Siding and Stack Yard of State Pollution Control Board, Orissa issued in August 2010 [13] Mineral Fertilizer Production and the Environment Part I The Fertilizer Industry's Manufacturing Processes and Environmental Issues of United Nations Environmental Programme: ISBN: 92-807-1640-9 [14] Handbook on Fertilizer Technology: FAI HO, New Delhi Học viên: Lưu Ngọc Vĩnh Trang 81 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” [15] Fertilizer Annual Report 2009 & 2010 published by the ‘Fertilizer Association of India’ [16] ‘Study on The Production of Ammonium Sulfate Fertilizer from Phosphogypsum’ published by Khalid K Abbas, Material Engineering Department, University of Technology, Baghdad in ‘Eng & Tech Journal, Vol.29, No.4, 2011 [17] Eco- balance features and significance of Hemihydrate Phosphogypsum Reprocessing into Gypsum Binding Materials’ S Gaiducis et al published in JOURNAL OF CIVIL ENGINEERING AND MANAGEMENT, 2009, 15(2): 205– 213 [18] IFA Technical Conference held at Marrakech, Morocco (during September-1 October 1998); Beijing, China (during 20-23 April 2004) [19] Potential groundwater contamination by fluoride from two South African phosphogypsum by MP Motalane* and CA Strydom, Department of Chemistry, University of Pretoria, Pretoria 0002, South Africa published in Water SA Vol 30 No October 2004, 465-468 [20] Phosphate Fertilizer Plant :Pollution Prevention and Abatement Handbook WORLD BANK GROUP Effective July 1998 [21] ‘Phosphogypsum Stack Closure’ by Dean Kleinschmidt, Senior Civil Engineer/Land Manager, Cargil Fertilizer, Inc (USA.) [22] Atomic Energy Regulatory Board Directive No 01/09-Use of Phosphogypsum in Building & Construction Materials & in Agriculture (Ref No No CH/AERB/IPSD/78/2009 dated March 19, 2009) downloaded from AERB website i.e http://www.aerb.gov.in/cgi-bin/constitution/directives.asp#dir3 [23] Indian Standard: Specification for By-product Gypsum IS:10170-1982 and Indian Standard: Specification for gypsum for agricultural use (first revision) IS:6046-1982 [24] The Fertilizer (Control) Order, 1985 published in the Gazette of India vide No GSR 758 (E) dated the 25 th September 1985 and latest amended vide No S.O 2886 (E) dated rd December 2010 [25] Study on The Production of Ammonium Sulfate Fertlizer From Phosphogypsum Khalid K Material Engineering Department, University of Technology / Baghdad 2011, 814-821 Học viên: Lưu Ngọc Vĩnh Trang 82 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phịng” PHỤ LỤC Phần mềm quy hoạch tuyến tính (phi tuyến) tối ưu hóa hai mức tối ưu clear;clc;clf p=0.05;%Muc co nghia 0.01 0.05 0.2 % -Quy hoach tuyen tinh % % y=[91 90.3 93.7 93.4 90.85 89.3 93.95 94.5]%ke hoach tai nhan 2^k y0=[95.2 95.5 95.4]%Thi nghiem lap tai tam % -Quy hoach phi tuyen -% y=[ %ke hoach tai nhan 2^k % % Canh tay don 2k % ]';%Tam n0=1 % y0=[ ]';%Thi nghiem lap tai tam % [hang cot]=size(y); switch hang case 2^3 qh31(p,y,y0); shg; case 2^3+2*3+1 qh32(p,y,y0); shg; case 2^4 qh41(p,y,y0);shg; case 2^4+2*4+1 qh42(p,y,y0);shg; end function qh31(p,yy,yy0) % So cac yeu to k=3;N=2^k;shs=k+2+(k-1)*k/2; % -% Học viên: Lưu Ngọc Vĩnh Trang 83 Đề tài: “Nghiên cứu sản xuất Sunphát amôn từ chất thải Gyps nhà máy DAP số – VINACHEM KCN Đình Vũ – Hải Phòng” % Toa khong thu nguyen xx0=[1 1 1 1 ]'; xx1=[-1 +1 -1 +1 -1 +1 -1 +1 ]'; xx2=[-1 -1 +1 +1 -1 -1 +1 +1 ]'; xx3=[-1 -1 -1 -1 1 1]'; X=[xx0 xx1 xx2 xx3 xx1.*xx2 xx1.*xx3 xx2.*xx3 xx1.*xx2.*xx3]; % he so pt hoi quy -B=(X'*X)^(-1)*X'*yy; % Danh gia tinh co nghia he so Student f2=length(yy0)-1; S11binh=sum((yy0-sum(yy0)/length(yy0)).^2)/f2 Sb=sqrt(S11binh/N) t=abs(B)/Sb hsvn=0; for i=1:shs if t(i)