8. Tình hình nghiên cứu trong nước và ngoài nước
1.4. Cơ sở lý luận về phương pháp hấp phụ
1.4.1 Các khái niệm cơ bản
Hấp phụ là sự tích lũy các chất trên bề mặt phân cách pha (khí – rắn, lỏng – rắn, khí – lỏng, lỏng – lỏng). Chất có bề mặt trên đó xảy ra sự hấp phụ gọi là chất hấp phụ, còn chất được tích lũy trên bề mặt chất hấp phụ gọi là chất bị hấp phụ [57, 58, 59, 60].
20
Hiện tượng hấp phụ xảy ra do lực tương tác giữa chất hấp phụ và chất bị hấp phụ. Tùy theo bản chất lực tương tác mà người ta có thể chia hấp phụ thành 2 loại: hấp phụ vật lý và hấp phụ hóa học.
Hấp phụ vật lý:
Các phân tử chất bị hấp phụ liên kết với những tiểu phân (nguyên tử, phân tử, các ion…) ở bề mặt phân chia pha bởi lực Van der Walls yếu. Đó là tổng hợp của nhiều loại lực khác nhau: tĩnh điện, tán xạ, cảm ứng và lực định hướng. Trong hấp phụ vật lý, các phân tử của chất bị hấp phụ và chất hấp phụ không tạo thành hợp chất hoá học (không tạo thành các liên kết hóa học) mà chất bị hấp phụ chỉ ngưng tụ trên bề mặt phân chia pha và bị giữ lại trên bề mặt chất hấp phụ. Do vậy, trong quá trình hấp phụ vật lý không có sự biến đổi đáng kể cấu trúc điện tử của cả chất hấp phụ và chất bị hấp phụ. Ở hấp phụ vật lý, nhiệt hấp phụ không lớn, năng lượng tương tác thường ít khi vượt quá 10 kcal/mol, phần nhiều từ 3 ÷ 5 kcal/mol và năng lượng hoạt hóa không vượt quá 1 kcal/mol [57,58].
Hấp phụ hóa học:
Xảy ra khi các phân tử chất hấp phụ tạo hợp chất hóa học với các phân tử chất bị hấp phụ. Lực hấp phụ hóa học khi đó là lực liên kết hóa học thông thường (liên kết ion, cộng hóa trị, liên kết phối trí…) Nhiệt hấp phụ hóa học tương đương với nhiệt phản ứng hóa học và có thể đạt tới giá trị 100 kcal/mol. Cấu trúc điện tử của cả chất hấp phụ và chất bị hấp phụ đều có sự biến đổi sâu sắc, tạo thành liên kết hóa học.
Trong thực tế, sự phân biệt hấp phụ vật lý và hấp phụ hóa học chỉ là tương đối vì ranh giới giữa chúng không rõ rệt. Trong một số quá trình hấp phụ xảy ra đồng thời cả hấp phụ vật lý và hấp phụ hóa học [57, 58, 59, 60].
Giải hấp phụ:
Giải hấp phụ là sự đi ra của chất bị hấp phụ khỏi bề mặt chất hấp phụ. Quá trình này dựa trên nguyên tắc sử dụng các yếu tố bất lợi đối với quá trình hấp phụ. Đây là phương pháp tái sinh vật liệu hấp phụ nên nó mang đặc trưng về hiệu quả kinh tế [57, 58, 59, 60].
21
Dung lượng hấp phụ (q) là lượng chất bị hấp phụ (độ hấp phụ) bởi 1 gam chất hấp phụ rắn [57, 58, 59, 60] được tính theo công thức:
q = (Co− C) × V
m
(1.1)
Trong đó:
q: Lượng chất bị hấp phụ (mg/g).
C0, C: Nồng độ ban đầu và nồng độ cân bằng của chất bị hấp phụ (mg/l). V: Thể tích dung dịch (l).
m: khối lượng chất hấp phụ (g).
Hiệu suất hấp phụ
Hiệu suất hấp phụ (H) là tỉ số giữa nồng độ dung dịch bị hấp phụ (C) và nồng độ dung dịch ban đầu C0, [61 - 65].
H = Co− C
Co × 100(%)
(1.2)
1.4.2. Phương trình hấp phụ đẳng nhiệt
Một hệ hấp phụ khi đạt đến trạng thái cân bằng, lượng chất bị hấp phụ là một hàm của nhiệt độ, áp suất hoặc nồng độ của chất bị hấp phụ:
q = f(T, PhoặcC) (1.3)
Ở nhiệt độ không đổi (T = const), đường biểu diễn q = fT(P hoặc C) được gọi là đường hấp phụ đẳng nhiệt. Đường hấp phụ đẳng nhiệt biểu diễn sự phụ thuộc của dung lượng hấp phụ tại một thời điểm vào nồng độ cân bằng hoặc áp suất của chất bị hấp phụ tại thời điểm đó ở một nhiệt độ xác định [66- 70].
Đối với chất hấp phụ là chất rắn, chất bị hấp phụ là chất lỏng, khí thì đường hấp phụ đẳng nhiệt được mô tả qua các phương trình hấp phụ đẳng nhiệt Henry, Freundlich, Langmuir…
22
mô tả cân bằng hấp phụ như: Dubinin, Frumkin, Tempkin tùy thuộc vào bản chất của hệ và các điều kiện tiến hành quá trình hấp phụ.
Luận văn này sẽ nghiên cứu cân bằng hấp phụ của vật liệu hấp phụ (VLHP) đối với ion kim loại Cu2+ trong môi trường nước theo mô hình đường hấp phụ đẳng nhiệt Langmuir và Freundlich.
1.4.3. Phương trình hấp phụ đẳng nhiệt Langmuir
Khi thiết lập phương trình hấp phụ [66- 70], Langmuir đã xuất phát từ các giả thuyết sau:
- Tiểu phân bị hấp phụ liên kết với bề mặt tại những trung tâm xác định. - Mỗi trung tâm chỉ hấp phụ một tiểu phân.
- Bề mặt chất hấp phụ là đồng nhất, nghĩa là năng lượng hấp phụ trên các trung tâm là như nhau và không phụ thuộc vào sự có mặt của các tiểu phân hấp phụ trên các trung tâm bên cạnh.
Phương trình Langmuir được xây dựng cho hệ hấp phụ khí rắn, nhưng cũng có thể áp dụng cho hấp phụ trong môi trường nước để phân tích các số liệu thực nghiệm. Trong pha lỏng phương trình phi tuyến tính có dạng:
q = qmax KL. C
1 + KL. C
(1.4)
Trong đó
KL: hằng số (cân bằng) hấp phụ Langmuir
q: dung lượng hấp phụ (lượng chất bị hấp phụ/1 đơn vị chất hấp phụ)
qmax: dung lượng hấp phụ tối đa của chất hấp phụ (lượng chất bị hấp phụ/1 đơn vị chất hấp phụ)
23 C/q (g/l) α M O C(mg/l) Hình 1.2. Đường hấp phụ đẳng nhiệt Langmuir [71]
Hình 1.3. Đồ thị sự phụ thuộc của C/q vào C [71]
Để xác định các hệ số trong phương trình hấp phụ đẳng nhiệt Langmuir, người ta chuyển phương trình (1.4) về dạng tuyến tính (1.5):
1 q= 1 qmax + 1 KL. qmax. C (1.5)
Từ đồ thị (hình 1.3) biểu diễn sự phụ thuộc của C/q vào C ta sẽ tính được KL và qmax: OM ̅̅̅̅̅ = 1 qmax. KL ; tgα = 1 qmax
Theo [72,73], từ giá trị KL có thể xác định được tham số cân bằng RL:
RL = 1 1 + KL. C0 (1.6) Trongđó: RL: tham số cân bằng C0: Nồng độ ban đầu (mg/l) KL: Hằng số Langmuir (l/mg)
24
Mối tương quan giữa các giá trị của KR và các dạng của mô hình hấp phụ đẳng nhiệt Langmuir thực nghiệm được thể hiện trong bảng 1.2.
Bảng 1.2. Mối tương quan của RL và dạng mô hình [73]
Giá trị RL Dạng mô hình
RL> 1 Không phù hợp
RL = 1 Tuyến tính
0 < RL< 1 Phù hợp
RL = 0 Không thuận nghịch
Phương trình Langmuir xác định được dung lượng hấp phụ cực đại và mối tương quan giữa quá trình hấp phụ và giải hấp phụ thông qua hằng số Langmuir KL, sự phù hợp của mô hình với thực nghiệm, do vậy đây là cơ sở để lựa chọn chất hấp phụ thích hợp cho hệ hấp phụ, [74-78].
1.4.4. Phương trình hấp phụ đẳng nhiệt Freundlich
Khi nghiên cứu về khả năng hấp phụ trong pha lỏng, trong trường hợp chất hấp phụ có lỗ xốp, Freundlich thiết lập được phương trình phi tuyến tính đẳng nhiệt trên cơ sở số liệu thực nghiệm [66- 70].
q = KF. C1/n (1.7)
Trong đó:
KF: là hằng số hấp phụ Freundlich.
Nếu C=1đơn vị thì a=KF tức là KF chính là dung lượng hấp phụ tại C=1, vậy nó là đại lượng có thể dùng để đặc trưng cho khả năng hấp phụ của hệ, giá trị KF lớn đồng nghĩa với hệ có khả năng hấp phụ cao.
1/ n (n > 1) là bậc mũ của C luôn nhỏ hơn 1, nó đặc trưng định tính cho bản chất lực tương tác của hệ, nếu 1/n nhỏ (n lớn) thì hấp phụ thiên về dạng hóa học và ngược lại, nếu 1/n lớn (n nhỏ) thì bản chất lực hấp phụ thiên về dạng vật lý, lực hấp phụ yếu.
25
Hình 1.4. Đường hấp phụ đẳng nhiệt hấp phụ Freundlich (bên trái), đồ thị để tìm các hằng số trong phương trình Frendlich (bên phải) [71].
Với hệ hấp phụ lỏng – rắn, n có giá trị nằm trong khoảng từ 1÷ 10 thể hiện sự thuận lợi của mô hình [78]. Như vậy, n cũng là một trong các giá trị đánh giá được sự phù hợp của mô hình với thực nghiệm.
Vì 1/n luôn nhỏ hơn 1 nên đường biểu diễn của phương trình (1.7) là một nhánh của đường parabol, và được gọi là đường hấp phụ đẳng nhiệt Freundlich (hình 1.4, bên trái).
Để xác định các hằng số trong phương trình Freundlich, người ta cũng sử dụng phương pháp đồ thị (hình 1.4, bên phải). Phương trình Freundlich có thể viết dưới dạng:
lg q = lg KF + 1/nlgC (1.8)
Như vậy lg a tỉ lệ bậc nhất với lg C. Đường biểu diễn trên hệ toạ độ lg q- lg C sẽ cắt trục tung tại N.
Tacó: ON̅̅̅̅ = lgKF; tgγ = 1/n
Mô hình hấp phụ Langmuir và Freundlich được ứng dụng nhiều trong nghiên cứu mô hình hấp phụ đối với hệ rắn – lỏng, đặc biệt trong các nghiên cứu hấp phụ chống ô nhiễm môi trường [74, 75, 77-81].
1.4.5. Động học hấp phụ
Đối với hệ hấp phụ lỏng – rắn, động học hấp phụ xảy ra theo một loạt các giai đoạn kế tiếp nhau [57, 58]:
26
- Chất bị hấp phụ chuyển động tới bề mặt chất hấp phụ. Đây là giai đoạn khuếch tán trong dung dịch
- Phân tử chất bị hấp phụ chuyển động đến bề mặt ngoài của chất hấp phụ chứa các hệ mao quản – giai đoạn khuếch tán màng.
- Chất bị hấp phụ khuếch tán vào bên trong hệ mao quản của chất hấp phụ - giai đoạn khuếch tán trong mao quản.
- Các phân tử chất bị hấp phụ chiếm chỗ các trung tâm hấp phụ - giai đoạn hấp phụ thực sự.
Trong tất cả các giai đoạn đó, giai đoạn nào có tốc độ chậm nhất sẽ quyết định hay khống chế chủ yếu toàn bộ quá trình động học hấp phụ. Với hệ hấp phụ trong môi trường nước, quá trình khuếch tán thường chậm và đóng vai trò quyết định. Tốc độ của một quá trình hấp phụ được xác định bởi sự thay đổi nồng độ của chất bị hấp phụ theo thời gian. Một vài mô hình động học hấp phụ đã được đưa ra để giải thích cơ chế hấp phụ.
1.4.5.1. Mô hình giả động học hấp phụ bậc 1
Theo đó, tốc độ của quá trình hấp phụ phụ thuộc bậc nhất vào dung lượng chất hấp phụ theo phương trình [72-82].
dqt
dt = k1x(qe− qt) (1.9)
Trong đó:
k1: hằng số tốc độ phản ứng theo mô hình động học bậc 1 (thời gian-1). qe, qt: dung lượng hấp phụ tại thời điểm cân bằng và thời điểm t (mg/g).
Áp dụng điều kiện biên tại thời điểm t = 0 và qt = 0, phương trình (1.9) trở thành:
ln qe
qe− qt = k1t
(1.10)
Và dạng phi tuyến tính bậc nhất:
qt = qe (1 – e –k1t) (1.11)
27
Hệ Số Xác Định (Hệ Số Tương Quan Bội Số R Bình Phương) (R2) [129]
R2 = 1 −∑NN=1(qt.exp.N− qt.model.N)2
∑NN=1(qt.exp.N − q̅̅̅̅̅̅̅̅̅)t.exp.N 2
(1.12)
Trong đó:
qt.exp: Dung lượng hấp phụ thực nghiệm (mg/g)
qt.model: Dung lượng hấp phụ tính theo mô hình (mg/g)
Lỗi bình phương gốc (Residual root mean square error: RMSE) [130]
RMSE = √ 1
n − 1∑(qe.exp.n− qe.model.n)2
n n−1
(1.13)
Kiểm định chi bình phương (chi-square test :χ 2 ) [130]
χ 2 = ∑(qe.exp.n− qe.model.n)2
qe.exp.n
n n−1
(1.14)
Phương trình (1.11) có thể chuyển về dạng tuyến tính bậc nhất
lg(qe − qt) = lgqe − k1/2,303 (1.15)
Từ (1.15) ta xác định được qe và hằng số k1. tgα = -k1/2,303; OM = lg qe
Phương trình (1.10) được gọi là phương trình giả động học bậc 1 [77-84]. Ngay từ khi công bố, phương trình đã sớm được áp dụng cho quá trình hấp phụ của triaxetat cellulozơ từ clorofom trên canxi silicat [73]. Trong suốt 4 thập kỉ tiếp theo cho đến nay, phương trình động học này đã được áp dụng phổ biến cho việc nghiên cứu động học hấp phụ với các chất ô nhiễm trong môi trường nước [84].
28
Hình 1.5. Đồ thị sự phụ thuộc của lg(qe – qt) vào t
1.4.5.1. Mô hình giả động học hấp phụ bậc 2
Theo mô hình, tốc độ của quá trình hấp phụ phụ thuộc bậc hai vào dung lượng của chất hấp phụ theo phương trình [72-75, 79-82, 84]:
dq1
dt = k2(qe− qt)
2 (1.16)
Trong đó:
k2: Hằng số tốc độ phản ứng theo mô hình giả động học bậc 2 (g/mg.thời gian). qe, qt: dung lượng hấp phụ tại thời điểm cân bằng và thời điểm t (mg/g).
Áp dụng điều kiện biên tại thời điểm t = 0 và qt = 0, phương trình (1.16) trở thành:
1 qe − qt = 1 qe + kt (1.17) Hoặc dạng tuyến tính 1 qt = 1 k2q2e+ 1 qe (1.18)
Có thể chuyển về dạng phi tuyến tính:
qt = k2qe
2t
1 + k2qet
29
Đặt h = k2qe2 phản ánh tốc độ hấp phụ ban đầu khi khi qt/t tiến dần đến 0, phương trình (1.17) và (1.18) trở thành: qt =1 1 h+ t qe (1.20) t qt = 1 h+ 1 qet (1.21)
Từ đồ thị sự phụ thuộc của t/qt vào t, ta xác định được qe và k2 [72, 73].
1.4.6. Động lực hấp phụ
Một quá trình hấp phụ có thể thực hiện theo hai dạng kỹ thuật: gián đoạn hoặc liên tục [66,67].
Kỹ thuật gián đoạn hay còn gọi là kỹ thuật theo mẻ với phương thức thực hiện nhận được sản phẩm sau một chu kỳ thời gian trước khi bắt đầu một chu kỳ mới tiếp theo. Kỹ thuật gián đoạn còn gọi là phương pháp hấp phụ tĩnh.
Kỹ thuật liên tục hay còn gọi là kỹ thuật dòng với phương thức thực hiện là nguyên liệu được liên tục đưa vào cột hấp phụ và sản phẩm cũng được lấy ra liên tục. Kỹ thuật liên tục còn gọi là phương pháp hấp phụ động.
Hấp phụ động là kỹ thuật có nhiều ưu điểm trong thực tiễn ứng dụng, tuy vậy việc thiết kế đúng một hệ hấp phụ khá phức tạp, trong nhiều trường hợp cần phải tiến hành nghiên cứu dạng pilot trước khi thiết kế hệ hoạt động sản xuất.
Một thiết bị chứa chất bị hấp phụ gọi là cột hấp phụ. Một dòng chất (lỏng, khí) chứa chất bị hấp phụ được đưa liên tục vào cột. Chất bị hấp phụ được giữ lại trong cột và chỉ xuất hiện ở phía đầu ra khỏi cột khi chất hấp phụ đã bão hòa dung lượng. Nếu theo dõi sự biến đổi của nồng độ chất hấp phụ theo thời gian, tức là theo dõi dải nồng độ theo thời gian (t =L/vC) sẽ nhận được đường cong thoát (hình 1.6).
Mục tiêu thực tiễn của nghiên cứu động lực hấp phụ là xác định (thiết kế) được thời gian hoạt động của một cột hấp phụ từ các số liệu thực nghiệm, đánh giá sự hao hụt dung lượng hấp phụ khi sử dụng phương pháp dòng chảy.
30
Hình 1.6. Đường cong thoát của cột hấp phụ [66,67]
Dựa vào mối quan hệ giữa nồng độ của chất bị hấp phụ tại thời điểm ban đầu (C0) và tại thời điểm t (Ct) vào thời gian, người ta đã đưa ra một số mô hình hấp phụ cho hệ hấp phụ động, như: Mô hình Clack, Thomas, Bohart – Adams; Yoon – Nelson, Wang, Wolborska. Dưới đây là một số mô hình động học hấp phụ của hệ hấp phụ động.
1.4.6.1. Mô hìnhThomas
Mô hình Thomas dựa trên giả thiết cho rằng quá trình hấp phụ tuân theo mô hình động học bậc hai và mô hình hấp phụ đẳng nhiệt Langmuir; đồng thời Thomas cũng bỏ qua quá trình chuyển khối bên trong và bên ngoài của chất hấp phụ. Do đó tốc độ quá trình khuếch tán được quyết định bởi phản ứng trên bề mặt giữa chất bị hấp phụ và dung lượng chưa bị sử dụng của chất hấp phụ [85]. Phương trình có dạng như sau [78, 86-89]:
Ce
C0 =
1
1 + exp [KT
31 Phương trình có dạng tuyến tính: ln (C0 Ce− 1) = KTq0M Q − KTCo Q V (1.23) Hoặc: (với t = V/Q) Ce C0 = 1 1 + exp [KT Q (qM − C0V)] (1.24) Trong đó:
C0, Ce: nồng độ đầu vào vào đầu ra của dung dịch hấp phụ (mg/l) q0: dung lượng hấp phụ cực đại (mg/g)
m: khối lượng chất hấp phụ (g) Q: Tốc độ dòng chảy (ml/phút)
V: lượng thể tích chảy qua cột hấp phụ (ml) KT: Hằng số tốc độ Thomas (ml/phút/mg).
Hình 1.7. Đồ thị sự phụ thuộc ln[(C0/Ce)-1] vào t
Xây dựng đồ thị biểu diễn sự phụ thuộc của ln (C0/Ct – 1) vào V (theo phương trình
1.23) hoặc ln (C0/Ct – 1) vào t (theo phương trình 1.24), ta sẽ xác định được các hệ số
trong phương trình.
32
1.4.6.2. Mô hình Yoon –Nelson
Mô hình này dựa trên giả thiết cho rằng độ giảm tốc độ của quá trình hấp phụ tỉ lệ với tỉ số giữa nồng độ sau khi hấp phụ và nồng độ ban đầu của chất bị hấp phụ [78]. Phương trình Yoon-Nelson [78, 85 - 89] có dạng: Ce C0 = 1 1 + exp[KYN(τ − t)] (1.25) Dạng tuyến tính: ln Ce C0− Ce = KYNt − τKYN (1.26) Trong đó:
C0, Ce: nồng độ đầu vào vào đầu ra của dung dịch hấp phụ (mg/l)