Tổng quan luận án
Công tác quản lý chất thải rắn (CTR) tại Việt Nam đang đối mặt với nhiều áp lực do phần lớn chất thải rắn sinh hoạt (CTRSH) chưa được phân loại tại nguồn mà được thu gom hỗn hợp và vận chuyển thẳng tới các bãi chôn lấp (BCL). Theo thống kê từ Bộ Xây dựng và Báo cáo Hiện trạng Môi trường Quốc gia, cả nước phát sinh lượng CTRSH lớn với tốc độ tăng trưởng dự báo 8,4%/năm tại khu vực đô thị và khoảng 5%/năm trên phạm vi toàn quốc; tổng khối lượng chất thải ước tính sẽ đạt 54 triệu tấn vào năm 2030. Hiện có khoảng 71% tổng khối lượng CTRSH được thu gom (tương đương 34.000 tấn/ngày) được xử lý bằng biện pháp chôn lấp tại 904 bãi chôn lấp. Tuy nhiên, chỉ khoảng 20% trong số đó là bãi chôn lấp hợp vệ sinh, 2/3 còn lại là bãi chôn lấp hở không hợp vệ sinh hoặc các bãi rác lộ thiên cấp xã thiếu hệ thống đầm nén, thu gom khí sinh học, xử lý nước rỉ rác và quan trắc môi trường.
Ngay cả tại các BCL được thiết kế hợp vệ sinh, việc thi công, vận hành và quản lý sau khi đóng cửa vẫn tồn tại nhiều bất cập kỹ thuật. Nước rỉ rác hình thành từ quá trình phân hủy yếm khí kết hợp với sự xâm nhập của nước mưa chứa hàm lượng cao các hợp chất hữu cơ khó phân hủy và kim loại nặng (KLN). Khi lớp lót đáy bị hư hại hoặc rò rỉ, các ion KLN độc hại sẽ di chuyển theo dòng thấm vào môi trường đất và tầng chứa nước dưới đất, đe dọa trực tiếp đến hệ sinh thái và sức khỏe dân cư khu vực xung quanh. Trong khi đó, các nghiên cứu chuyên sâu về cơ chế và dự báo định lượng phạm vi lan truyền KLN từ các ô chôn lấp đã đóng tại Việt Nam còn rất hạn chế.
Mục tiêu nghiên cứu:
- Đánh giá đặc điểm và nguồn phát sinh ô nhiễm từ nước rỉ rác của bãi chôn lấp CTR hợp vệ sinh, xác định sự ảnh hưởng đến quá trình lan truyền chất ô nhiễm trong môi trường đất.
- Đánh giá và dự báo mức độ lan truyền ô nhiễm KLN bằng phương pháp định lượng (mô hình hóa toán học) nhằm thiết lập cơ sở dữ liệu phục vụ quản lý môi trường bãi chôn lấp.
Nội dung nghiên cứu cụ thể:
- Nghiên cứu các yếu tố ảnh hưởng đến việc hình thành lượng và tính chất của nước rỉ rác; khảo sát sự biến thiên nồng độ chất ô nhiễm trong nước rỉ rác theo thời gian.
- Xây dựng mô hình tính toán lượng nước rỉ rác phát sinh phù hợp với đối tượng và điều kiện khí hậu, vận hành tại Việt Nam.
- Nghiên cứu cơ sở lý thuyết về lan truyền chất ô nhiễm bằng mô hình toán học; áp dụng phương pháp phần tử hữu hạn (PTHH) để giải hệ phương trình vi phân đạo hàm riêng nhằm dự báo nồng độ kim loại nặng lan truyền trong đất.
- Đánh giá các nguy cơ ô nhiễm môi trường đất phát sinh từ bãi chôn lấp chất thải rắn hợp vệ sinh và xác định mức độ ô nhiễm KLN thực tế tại khu vực nghiên cứu.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Ô nhiễm kim loại nặng (As, Cd, Cr, Cu, Pb, Zn...) trong môi trường đất phát sinh từ ô chôn lấp đã đóng của bãi chôn lấp chất thải rắn hợp vệ sinh.
- Phạm vi không gian: Bãi chôn lấp chất thải rắn Kiêu Kỵ, xã Kiêu Kỵ, huyện Gia Lâm, thành phố Hà Nội (tổng diện tích 5,4 ha; khu vực chôn lấp gồm 5 ô chôn lấp với diện tích 31.190 m²).
- Phạm vi thời gian và tài liệu: Số liệu quan trắc thành phần nước rỉ rác trong giai đoạn 2008–2017; các ô chôn lấp số 1 và số 2 vận hành từ năm 2012 và kết thúc chôn lấp, đóng bãi vào năm 2017; khảo sát, khoan lấy mẫu đất tại hiện trường vào tháng 4.
Tổng quan tài liệu và vị trí của luận án
Luận án đã tổng thuật hệ thống các công trình nghiên cứu trong và ngoài nước về thành phần chất thải, sự hình thành nước rỉ rác và cơ chế di chuyển của các chất ô nhiễm:
- Nguồn gốc và đặc trưng kim loại nặng trong CTR: Tác giả dẫn dụng nghiên cứu của Vitalii Ishchenko (2019) về hàm lượng kim loại nặng phân bố trong các thành phần rác thải (nhựa, cao su, gỗ, giấy, dệt, chất thải nguy hại); nghiên cứu của Kjeldsen và cộng sự (2002) phân loại 4 nhóm chất ô nhiễm chính trong nước rỉ rác (hợp chất hữu cơ hòa tan COD/BOD/TOC, thành phần vô cơ Ca/Mg/Na/NH4+, kim loại nặng Cd/Cr/Cu/Pb/Ni/Zn, và hợp chất hữu cơ xenobiotic - XOCs).
- Biến thiên nồng độ ô nhiễm theo tuổi bãi chôn lấp và khí hậu: Nghiên cứu của Farquhar (1989) về đồ thị phát sinh nước rỉ rác theo thời gian; Ehrig (1989) và Krumpelbeck & Ehrig (1999) về sự ổn định của nồng độ amoni và kim loại nặng qua các pha axit và pha metan hóa tại các bãi chôn lấp của Đức; Akesson & Nilsson (1997) tại Thụy Điển và Chu cùng cộng sự (1994) tại Hồng Kông về tác động của mùa mưa đến việc pha loãng nồng độ nước rỉ rác; các nghiên cứu so sánh tính chất nước rỉ rác tại Châu Á, Châu Âu, Châu Mỹ và Châu Phi của Akgul và cộng sự (2013), Chou và cộng sự (2013), Guo và cộng sự (2012), Zainol và cộng sự (2012), Bazimenyera và cộng sự (2012), Boumechhour (2013), Joshua và cộng sự (2009).
- Hành vi và độ linh động của kim loại nặng trong đất: Dẫn theo Kabata-Pendias (1984) về sự phân loại độ linh động của các nguyên tố (Se, Hg, As, Cd, Pb, Bi, Ti...) trong các điều kiện môi trường đất khác nhau (oxy hóa, trung tính - kiềm, axit, khử).
- Các phương pháp tính toán lượng nước rỉ rác hiện hành:
| Tiêu chí so sánh |
Phương pháp Cân bằng nước cổ điển (WBM) |
Mô hình Định lượng Thủy văn BCL (HELP - US EPA) |
Phương pháp Cân bằng nước theo chuỗi (SWB - Mexico) |
| Cơ sở tính toán |
Phương trình cân bằng nước tổng quát |
Mô hình mô phỏng số trên máy tính |
Phương trình cân bằng nước theo trình tự thời gian |
| Yêu cầu dữ liệu khí hậu |
Lượng mưa, lượng bốc hơi theo tháng |
Số liệu khí tượng chi tiết theo từng ngày |
Lượng mưa, nhiệt độ theo tháng |
| Yêu cầu dữ liệu bãi rác |
Số liệu địa chất thủy văn chung |
Số liệu thiết kế chi tiết, tính chất các lớp đất |
Tính chất vật lý, độ ẩm của khối chất thải |
| Công cụ thực hiện |
Bảng tính thủ công đơn giản |
Phần mềm chuyên dụng (môi trường DOS/Windows) |
Bảng tính Excel |
| Giai đoạn áp dụng |
Mọi loại bãi và các giai đoạn |
BCL đang vận hành, chuẩn bị hoặc đã đóng |
BCL đang trong giai đoạn vận hành |
| Độ chuẩn xác dự báo |
Sai số lớn, chưa tính chi tiết đầm nén |
Thường dự báo lượng nước cao hơn thực tế |
Thường dự báo thấp hơn phương pháp WBM |
Khoảng trống nghiên cứu được luận án giải quyết:
Tại Việt Nam, việc dự báo lượng nước rỉ rác chủ yếu áp dụng công thức kinh nghiệm đơn giản, không phản ánh đúng chế độ vận hành thực tế (bãi thường quá tải, việc phủ đất hàng ngày bị bỏ qua, kỹ thuật đầm nén không đồng đều). Đồng thời, các mô hình lan truyền chất ô nhiễm thương mại 1D hoặc các giải pháp giải tích chưa mô tả được đầy đủ cấu trúc địa chất phân lớp bất đồng nhất và điều kiện biên thực tế trong không gian 3 chiều (3D). Luận án định vị giải quyết khoảng trống này bằng cách đề xuất công thức cân bằng nước cải tiến và ứng dụng phương pháp phần tử hữu hạn để mô phỏng 3D trường nồng độ kim loại nặng.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và khung phân tích
- Động học chất lưu và truyền chất trong môi trường rỗng: Cơ chế di chuyển của nước trong tầng đất chưa bão hòa và bão hòa dưới tác động của lực trọng trường và lực hút mao dẫn.
- Cơ chế đối lưu - khuếch tán - phân tán (Advection - Dispersion): Chất ô nhiễm trong nước rỉ rác lan truyền do dòng chất lỏng chuyển động cuốn theo (đối lưu) kết hợp với hiện tượng khuếch tán phân tử do chênh lệch nồng độ và phân tán thủy lực do cấu trúc lỗ rỗng của đất.
- Cơ chế hấp phụ và làm chậm: Sự tương tác giữa các ion kim loại nặng với bề mặt hạt khoáng sét và chất hữu cơ trong đất được mô tả thông qua hệ số làm chậm $R$ ($R \ge 1$).
- Toán tử trường không gian: Sử dụng khái niệm trường vectơ gradient nồng độ $\nabla C(x, y, z, t)$ cùng các phép toán vi phân phân kỳ $(\nabla \cdot)$ và tích vô hướng $\odot$ trong hệ tọa độ Descartes $Oxyz$ để biểu diễn phương trình truyền chất tổng quát dưới dạng phương trình vi phân đạo hàm riêng (PDE).
Phương pháp nghiên cứu
- Phương pháp thực nghiệm hiện trường và phòng thí nghiệm:
- Khảo sát thực địa bãi chôn lấp Kiêu Kỵ, thu thập mẫu nước rỉ rác tại hồ sinh học và các ô chôn lấp theo các mùa trong năm giai đoạn 2008–2017.
- Khoan địa chất lấy mẫu đất tại các lỗ khoan L4 (LK4) và L5 (LK5) ở ô chôn lấp đã đóng, chia cắt mẫu đất theo từng tầng độ sâu để đánh giá mức độ tích tụ kim loại nặng.
- Phân tích các thông số lý hóa của nước rỉ rác ($pH$, $BOD_5$, $COD$, $TOC$, $NH_4^+$, $NO_3^-$, $NO_2^-$, Tổng $N$, Tổng $P$, Coliform, Dầu mỡ) và kim loại nặng ($As, Cd, Cr, Cu, Pb, Zn$) tại phòng thí nghiệm của Viện Hàn lâm Khoa học và Công nghệ Việt Nam và Viện Khoa học và Kỹ thuật Môi trường (Trường Đại học Xây dựng Hà Nội).
- Phương pháp tính toán và mô hình hóa:
- Ứng dụng phương pháp phần tử hữu hạn (PTHH) để giải số phương trình vi phân đạo hàm riêng mô tả quá trình lan truyền ô nhiễm trong môi trường đất phân lớp với hình học và điều kiện biên phức tạp.
- Xây dựng thuật toán và lập trình tính toán trên phần mềm MATLAB để mô phỏng sự biến thiên nồng độ theo thời gian $t$ và không gian 3 chiều ($x, y, z$).
Nội dung chính theo từng chương
Chương 1: Tổng quan các vấn đề nghiên cứu
Chương 1 hệ thống hóa cơ sở pháp lý và tiêu chuẩn kỹ thuật về bãi chôn lấp tại Việt Nam (TCXDVN 261:2001, TCVN 6696:2009, QCVN 25:2009/BTNMT, Nghị định 38/2015/NĐ-CP). Đồng thời, tác giả trình bày bức tranh tổng thể về hiện trạng xử lý CTRSH trên 6 vùng kinh tế - xã hội.
| Vùng lãnh thổ |
Tổng số BCL |
BCL hợp vệ sinh |
BCL không hợp vệ sinh |
Tổng khối lượng tiếp nhận (tấn/năm) |
Quy mô $\ge 20\text{ ha}$ |
Quy mô $1 - 20\text{ ha}$ |
Quy mô $< 1\text{ ha}$ |
| Tây Bắc |
9 |
2 |
7 |
30.166 |
0 |
7 |
2 |
| Đông Bắc Bộ |
10 |
10 |
0 |
54.490 |
0 |
4 |
12 |
| KTTĐ Bắc Bộ |
96 |
12 |
84 |
366.041 |
0 |
12 |
83 |
| Bắc Trung Bộ & ĐB Sông Hồng |
31 |
14 |
17 |
172.666 |
2 |
31 |
21 |
| KTTĐ Miền Trung |
9 |
7 |
2 |
55.279 |
1 |
10 |
1 |
| Nam Trung Bộ & Tây Nguyên |
9 |
1 |
8 |
69.830 |
0 |
11 |
2 |
| KTTĐ Phía Nam |
9 |
3 |
6 |
214.198 |
3 |
5 |
4 |
| Đồng bằng Sông Cửu Long |
65 |
4 |
61 |
268.297 |
0 |
19 |
2 |
| Tổng cộng |
657 |
203 |
454 |
5.153 (tấn/ngày)* |
36 |
384 |
227 |
* Nguồn trích xuất từ Báo cáo Bộ Xây dựng (2015) trong tài liệu luận án.
Đặc điểm khu vực nghiên cứu Kiêu Kỵ (Gia Lâm, Hà Nội):
- Khí hậu: Nhiệt độ trung bình năm $23,5^\circ\text{C}$; độ ẩm tương đối trung bình $83% - 84,5%$; lượng mưa trung bình năm đạt $1.612\text{ mm}$ (giai đoạn 2006–2019 dao động từ $1.240\text{ mm}$ đến $2.268\text{ mm}$); lượng bốc hơi trung bình năm đạt $936 - 989\text{ mm}$.
- Cấu tạo địa chất và địa chất thủy văn: Khu đất có diện tích 5,4 ha nguyên là moong khai thác sét của xí nghiệp gạch Kiêu Kỵ. Mực nước dưới đất nông cách mặt đất tự nhiên khoảng $2\text{ m}$. Cấu trúc địa tầng từ trên xuống dưới gồm 6 lớp:
| STT |
Tên lớp đất |
Chiều dày trung bình (m) |
Dung tích tự nhiên $\gamma\text{ (g/cm}^3\text{)}$ |
Sức chịu tải $R_0\text{ (kg/cm}^2\text{)}$ |
Mô đun biến dạng $E_0\text{ (kg/cm}^2\text{)}$ |
| 1 |
Đất bề mặt lẫn sét pha nâu hồng rời rạc lẫn gạch vỡ |
0,8 |
- |
- |
- |
| 2 |
Sét pha nâu vàng đốm đen đáy hồ |
1,3 |
- |
- |
120 |
| 3 |
Sét pha xám xanh hồng tro |
0,7 |
- |
- |
75 |
| 4 |
Sét pha xám nâu |
0,6 |
- |
- |
100 |
| 5 |
Cát hạt trung xám nâu (độ ẩm 23%) |
4,0 |
2,0 |
- |
150 |
| 6 |
Cát hạt nhỏ xám tro (độ ẩm 22%) |
3,0 - 4,0 |
- |
- |
- |
-
Thiết kế kỹ thuật ô chôn lấp Kiêu Kỵ: Khu chôn lấp có diện tích $31.190\text{ m}^2$ chia làm 5 ô chôn lấp. Cấu tạo đáy ô từ trên xuống gồm: lớp sỏi $2\times 4\text{ cm}$ dày $0,05\text{ m}$; lớp đất sét đầm kỹ dày $0,2\text{ m}$; màng chống thấm HDPE dày $1,5\text{ mm}$ ($0,0015\text{ m}$); và lớp nền á sét đầm chặt. Rác được đầm nén từng lớp dày $0,8\text{ m}$ và phủ lớp đất hàng ngày dày $0,15 - 0,2\text{ m}$. Nước rỉ rác được thu gom qua hệ thống mương, ống dẫn có độ dốc ngang $i = 0,01$ và bơm về hồ sinh học xử lý trước khi xả ra sông Cầu Bây.
-
Hàm lượng kim loại nặng trong các thành phần rác thải (theo Ishchenko, 2019):
| Thành phần chất thải |
Chì - Pb (g/kg) |
Cadimi - Cd (g/kg) |
Niken - Ni (g/kg) |
Crom - Cr (g/kg) |
| Nhựa |
40,1 |
0,8 |
- |
38,6 |
| Vải dệt |
10,2 |
0,3 |
- |
14,0 |
| Cao su |
170,0 |
6,4 |
- |
1.228,0 |
| Gỗ |
10,6 |
0,2 |
- |
167,0 |
| Giấy |
7,7 |
0,1 |
- |
9,4 |
| Vật liệu xây dựng |
6,9 |
0,1 |
- |
10,5 |
| Thực phẩm |
4,0 |
- |
- |
7,0 |
| Chất thải nguy hại |
13.950,0 |
228.220,0 |
661.220,0 |
50.550,0 |
| Phần còn lại |
19,9 |
0,2 |
- |
49,3 |
Chương 2: Lý thuyết lan truyền ô nhiễm và phương pháp thực hiện nghiên cứu
Chương 2 tập trung thiết lập hệ thống phương trình toán - lý mô tả quá trình vận chuyển chất ô nhiễm trong môi trường xốp bão hòa và không bão hòa nước:
- Cơ chế dòng thấm và mao dẫn: Sự chuyển động của nước dưới tác động của trọng lực và hiện tượng mao dẫn trong các mao quản đất; xác định lưu lượng dòng thấm qua hệ số thấm $k$ và gradient thủy lực.
- Phương trình vi phân lan truyền chất ô nhiễm: Thiết lập phương trình truyền khối có tính đến quá trình phân tán cơ học, khuếch tán phân tử, dòng đối lưu và tương tác hóa lý (hấp phụ/nhả hấp phụ):
$$\frac{\partial (\theta C)}{\partial t} + \rho_b \frac{\partial S}{\partial t} = \nabla \cdot (\theta \mathbf{D} \nabla C) - \nabla \cdot (\mathbf{q} C)$$
Trong đó $C$ là nồng độ chất tan trong pha lỏng; $\theta$ là độ ẩm thể tích; $\rho_b$ là dung trọng khô của đất; $S$ là hàm lượng chất bị hấp phụ trên pha rắn; $\mathbf{D}$ là tensor hệ số khuếch tán - phân tán; $\mathbf{q}$ là vectơ vận tốc dòng thấm Darcy. Khi giả định cân bằng hấp phụ tuyến tính ($S = K_d C$), phương trình được rút gọn qua hệ số làm chậm $R = 1 + \frac{\rho_b K_d}{\theta}$.
- Rời rạc hóa bằng phương pháp phần tử hữu hạn (PTHH): Không gian nghiên cứu 3D được chia lưới thành các phần tử tứ diện/lục diện hữu hạn. Tác giả xây dựng phương trình ma trận cho từng phần tử và ghép nối thành hệ phương trình đại số tổng thể dạng:
$$[M] \left{\frac{dC}{dt}\right} + [K] {C} = {F}$$
Trong đó $[M]$ là ma trận khối lượng, $[K]$ là ma trận độ cứng/đối lưu - phân tán, và ${F}$ là vectơ tải trọng/nguồn biên. Thuật toán giải được lập trình và thực thi trong môi trường MATLAB.
Chương 3: Kết quả và thảo luận
Chương 3 trình bày các kết quả thực nghiệm và mô phỏng số:
1. Đề xuất phương pháp tính toán lượng nước rỉ rác cải tiến:
Tác giả xây dựng mô hình cân bằng nước tổng quát áp dụng cho bãi chôn lấp Kiêu Kỵ và điều kiện Việt Nam:
$$LLC = IC + WGC + WD + WEM$$
Trong đó:
- $IC$: Lượng nước rỉ rác sinh ra do sự xâm nhập của nước mưa qua bề mặt ô chôn lấp (xác định theo diện tích hở, lớp phủ hàng ngày và lượng mưa thực tế).
- $WGC$: Lượng nước rỉ rác phát sinh từ độ ẩm tự nhiên sẵn có trong CTR khi bị đầm nén vượt quá khả năng giữ nước bão hòa.
- $WD$: Lượng ẩm thoát ra từ quá trình phân hủy sinh học chất hữu cơ.
- $WEM$: Lượng chất lỏng bổ sung từ các chế phẩm sinh học hoặc nước tuần hoàn đưa vào ô rác.
Công thức này cho phép tính toán linh hoạt cho cả giai đoạn vận hành mở và giai đoạn sau khi đóng bãi (có lớp phủ đỉnh trồng cây hoặc không trồng cây, đồng thời xét đến sự suy thoái vật lý của lớp che phủ theo thời gian).
2. Khảo sát chất lượng nước rỉ rác tại bãi Kiêu Kỵ (2008–2017):
- Chất hữu cơ: Nồng độ $COD$ và $BOD_5$ rất cao trong những năm đầu sau khi đóng bãi và giảm dần theo thời gian. Tỷ lệ $BOD_5/COD$ giảm từ mức $> 0,5$ (giai đoạn mới chôn lấp, giàu axit béo dễ bay hơi VFA) xuống dưới $0,1 - 0,2$ (giai đoạn metan hóa ổn định, chủ yếu còn lại axit humic và axit fulvic khó phân hủy).
- Dinh dưỡng và vi sinh: Nồng độ $NH_4^+$ duy trì ở mức cao và kéo dài nhiều năm; nồng độ tổng photpho ($P$), Coliform và dầu mỡ biến thiên theo mùa mưa - khô.
- Kim loại nặng trong nước rỉ rác: Phát hiện sự có mặt của $As, Cd, Cr, Cu, Pb, Zn, Ni, Fe, Mn$. Trong đó, nồng độ một số kim loại như chì ($Pb$), asen ($As$) và crom ($Cr$) ở mức đáng kể do sự hòa tan từ rác thải công nghiệp, đồ điện tử và pin lẫn trong CTRSH.
3. Khảo sát ô nhiễm kim loại nặng trong đất bãi chôn lấp Kiêu Kỵ:
Kết quả phân tích mẫu đất khoan tại hai lỗ khoan L4 và L5 xuyên qua các tầng địa chất cho thấy sự tích lũy kim loại nặng giảm dần theo chiều sâu. Tại các lớp đất sát đáy ô chôn lấp, hàm lượng $As$, $Cr$ và $Pb$ ghi nhận giá trị vượt quy chuẩn kỹ thuật quốc gia về chất lượng đất, chứng minh có sự rò rỉ và thẩm thấu chất ô nhiễm qua các vị trí màng chống thấm bị khiếm khuyết.
4. Kết quả mô phỏng lan truyền kim loại nặng 1D và 3D bằng phương pháp PTHH:
- Mô hình 1D: Dự báo sự di chuyển theo phương thẳng đứng ($Z$) của Asen ($As$) và Cadimi ($Cd$) qua các lớp đất sét pha và cát.
- Mô hình 3D: Mô phỏng toàn diện trường nồng độ theo các phương $X, Y, Z$ tại các thời điểm (ví dụ: tháng thứ 10, tháng thứ 24, và các năm tiếp theo).
- Ảnh hưởng của hệ số làm chậm ($R$): Khi không xét đến quá trình hấp phụ ($R = 1$), chất ô nhiễm lan truyền nhanh và sâu vào các tầng cát bên dưới. Khi xét đến tương tác hấp phụ ($R = 2$ hoặc lớn hơn tùy loại đất), tốc độ di chuyển của đỉnh nồng độ $As$ và $Cd$ giảm rõ rệt, chất ô nhiễm có xu hướng lưu giữ lại ở các tầng sét pha trên cùng.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học
- Thiết lập công thức cân bằng nước cải tiến ($LLC = IC + WGC + WD + WEM$): Đã bổ sung và lượng hóa chi tiết các thành phần phát sinh nước rỉ rác phù hợp với đặc thù vận hành thực tế tại Việt Nam (chất thải hữu cơ cao, độ ẩm lớn, đầm nén không triệt để, lớp phủ hàng ngày bị gián đoạn, và sự suy thoái của lớp che phủ sau đóng bãi).
- Xây dựng mô hình toán - số 3D bằng phương pháp phần tử hữu hạn: Ứng dụng thành công toán tử gradient $\nabla C(x, y, z, t)$ và thuật toán PTHH trên nền tảng MATLAB để giải bài toán lan truyền chất ô nhiễm trong môi trường địa chất phân lớp, giải quyết được bài toán biên phức tạp mà phương pháp giải tích 1D không thể đáp ứng.
Đóng góp mới về mặt thực tiễn
- Cung cấp bộ dữ liệu thực nghiệm chi tiết: Xác định hiện trạng ô nhiễm kim loại nặng ($As, Cr, Pb$) trong nước rỉ rác và các lớp đất tại BCL Kiêu Kỵ (Gia Lâm, Hà Nội) giai đoạn 2008–2017.
- Dự báo không gian - thời gian vành đai lan truyền ô nhiễm: Mô phỏng định lượng phạm vi di chuyển của kim loại nặng trong lòng đất, hỗ trợ các nhà quản lý và kỹ sư môi trường xác định khoảng cách an toàn vệ sinh, lập kế hoạch quan trắc định kỳ và thiết kế phương án cải tạo, phục hồi môi trường cho các bãi chôn lấp đã đóng cửa.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Luận án tập trung khảo sát thực nghiệm tại một địa bàn cụ thể là bãi chôn lấp Kiêu Kỵ (Gia Lâm, Hà Nội) với đặc thù địa tầng phù sa sông Hồng (các lớp sét pha xen kẹp cát). Các thông số mô hình hóa (như hệ số phân tán, hệ số làm chậm $R$) được cố định theo một số kịch bản giả định tuyến tính.
- Hướng nghiên cứu tiếp: Mở rộng kiểm chứng mô hình toán 3D trên các bãi chôn lấp có cấu trúc địa chất, thủy văn khác nhau (vùng đồi núi, karst, hoặc vùng cát ven biển); tích hợp các mô hình hấp phụ phi tuyến (Langmuir, Freundlich) và xét thêm ảnh hưởng của phản ứng sinh - hóa phức tạp trong nước rỉ rác.
Giá trị tham khảo
- Đối với nghiên cứu sinh và giảng viên ngành Kỹ thuật Môi trường, Địa kỹ thuật: Cung cấp tài liệu tham khảo về phương pháp luận kết hợp giữa quan trắc thực địa, phân tích hóa học môi trường và mô hình hóa số bằng phương pháp phần tử hữu hạn (MATLAB).
- Đối với kỹ sư thiết kế công trình hạ tầng kỹ thuật môi trường: Hướng dẫn phương pháp tính toán dung tích nước rỉ rác theo công thức $LLC$ cải tiến, phục vụ thiết kế hồ chứa và trạm xử lý nước rỉ rác sát với thực tế vận hành.
- Đối với cơ quan quản lý nhà nước về môi trường (Sở TN&MT, Chi cục BVMT): Cung cấp công cụ khoa học và cơ sở dữ liệu định lượng để đánh giá rủi ro, quy hoạch khoảng cách ly an toàn và xây dựng phương án đóng cửa bãi chôn lấp hợp quy chuẩn.
Câu hỏi thường gặp
1. Luận án lựa chọn đối tượng và địa bàn nghiên cứu cụ thể nào?
Luận án nghiên cứu đặc điểm ô nhiễm và lan truyền kim loại nặng trong môi trường đất từ ô chôn lấp đã đóng của bãi chôn lấp chất thải rắn hợp vệ sinh Kiêu Kỵ, nằm trên tuyến đường liên xã Kiêu Kỵ - Đa Tốn, huyện Gia Lâm, thành phố Hà Nội.
2. Công thức tính toán lượng nước rỉ rác do luận án đề xuất có dạng như thế nào?
Luận án đề xuất công thức cân bằng nước: $LLC = IC + WGC + WD + WEM$, trong đó tính toán chi tiết lượng nước mưa xâm nhập ($IC$), lượng nước ép ra từ rác khi đầm nén ($WGC$), lượng nước phát sinh từ phân hủy sinh học ($WD$), và lượng chất lỏng sinh học bổ sung ($WEM$).
3. Phương pháp toán học và công cụ tính toán nào được sử dụng để dự báo lan truyền ô nhiễm?
Luận án sử dụng phương pháp phần tử hữu hạn (PTHH) để giải hệ phương trình vi phân đạo hàm riêng mô tả cơ chế đối lưu - khuếch tán trong không gian 3D, với thuật toán mô phỏng được lập trình và giải trên phần mềm MATLAB.
4. Những kim loại nặng nào trong đất tại bãi chôn lấp Kiêu Kỵ được phát hiện có nồng độ vượt quy chuẩn?
Kết quả phân tích mẫu đất tại các lỗ khoan khảo sát (L4 và L5) cho thấy hàm lượng asen ($As$), crom ($Cr$) và chì ($Pb$) ở các tầng đất sát đáy ô chôn lấp có giá trị vượt quá quy chuẩn kỹ thuật cho phép.
5. Tỷ lệ bãi chôn lấp hợp vệ sinh tại Việt Nam được ghi nhận trong luận án là bao nhiêu?
Theo số liệu Báo cáo Hiện trạng Môi trường Quốc gia được trích dẫn trong luận án, trong tổng số 904 bãi chôn lấp CTRSH trên cả nước, chỉ có khoảng 20% là bãi chôn lấp hợp vệ sinh, còn lại 80% là các bãi chôn lấp không hợp vệ sinh hoặc bãi rác lộ thiên cấp xã.
Kết luận
Luận án Tiến sĩ của tác giả Hoàng Ngọc Hà đã giải quyết bài toán đánh giá và dự báo lan truyền kim loại nặng từ bãi chôn lấp chất thải rắn hợp vệ sinh Kiêu Kỵ vào môi trường đất thông qua sự kết hợp giữa khảo sát thực nghiệm và mô hình hóa toán học. Công trình đã đề xuất công thức cân bằng nước cải tiến $LLC$ phù hợp với điều kiện vận hành tại Việt Nam và xây dựng thành công mô hình phần tử hữu hạn 3D trên nền tảng MATLAB để mô phỏng trường nồng độ các kim loại độc hại ($As, Cd, Cr, Pb$). Kết quả nghiên cứu là nguồn tư liệu khoa học có giá trị thực tiễn cao, đóng góp cơ sở dữ liệu quan trọng cho công tác quản lý, quan trắc và xử lý môi trường tại các bãi chôn lấp đã đóng cửa.