Tổng quan về luận án

Sự gia tăng nhanh chóng của chất thải rắn sinh hoạt (CTRSH) tại các quốc gia đang phát triển đang đặt ra những thách thức nghiêm trọng đối với an ninh môi trường và sức khỏe cộng đồng. Tại Việt Nam, theo Báo cáo Hiện trạng Môi trường Quốc gia, tốc độ phát sinh CTRSH đô thị gia tăng bình quân 8,4%/năm, ước tính tổng lượng rác thải toàn quốc sẽ đạt 54 triệu tấn vào năm 2030. Hiện nay, khoảng "71% (tương đương 34.000 tấn/ngày) được xử lý bằng phương pháp chôn lấp", tuy nhiên "trong số các bãi chôn lấp hiện nay chỉ có khoảng 20% là bãi chôn lấp hợp vệ sinh". Tình trạng rác thải không được phân loại tại nguồn dẫn đến nước rỉ rác chứa nồng độ kim loại nặng (KLN) phức tạp, độc tính cao và khả năng tích lũy sinh học lâu dài.

Khoảng trống nghiên cứu (research gap) trọng yếu hiện nay nằm ở sự thiếu hụt các mô hình định lượng 3D mô phỏng động thái lan truyền KLN trong điều kiện địa chất thủy văn và khí hậu nhiệt đới gió mùa tại các bãi chôn lấp (BCL) đã đóng cửa ở Việt Nam. Luận án tiến sĩ kỹ thuật môi trường của tác giả Hoàng Ngọc Hà với đề tài "Nghiên cứu đặc điểm 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 tại Kiêu Kỵ, Gia Lâm, Hà Nội" (thực hiện dưới sự hướng dẫn khoa học của GS.TS. Nguyễn Thị Kim Thái và PGS. Nguyễn Văn Hoàng tại Trường Đại học Xây dựng Hà Nội) đã giải quyết triệt để bài toán này.

Hệ thống câu hỏi và giả thuyết nghiên cứu được thiết lập chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Các yếu tố khí hậu, cơ chế vận hành và thành phần CTR tác động như thế nào đến động thái hình thành lưu lượng và tải lượng ô nhiễm trong nước rỉ rác theo thời gian?
    • Giả thuyết 1 (H1): Nước rỉ rác chuyển biến rõ rệt qua các pha phân hủy sinh hóa; tỷ lệ $BOD_5/COD$ suy giảm mạnh từ pha axit sang pha metan hóa kéo theo sự biến đổi khả năng hòa tan và giải phóng KLN.
  2. Câu hỏi nghiên cứu 2 (RQ2): Quy luật lan truyền không gian - thời gian của các ion kim loại nặng ($As, Cd, Cr, Pb, Cu, Zn, Ni$) qua các tầng địa chất đất á sét và cát ngầm diễn ra theo cơ chế thủy lực và hấp phụ nào?
    • Giả thuyết 2 (H2): Sự kết hợp giữa trường vectơ gradient nồng độ $\nabla C(x,y,z,t)$ với phương trình vi phân đạo hàm riêng đối lưu - phân tán - hấp phụ giải bằng phương pháp phần tử hữu hạn (PTHH) trên nền tảng MATLAB sẽ dự báo chính xác biên độ xâm nhập của KLN qua các lớp địa tầng bão hòa và chưa bão hòa.

Khung lý thuyết của luận án tích hợp Lý thuyết Vận chuyển Thủy lực ngầm Darcy, Lý thuyết Đối lưu - Khuếch tán phân tán thủy động lực học (Advection-Dispersion Theory) và Động học Hấp phụ đẳng nhiệt (Sorption Isotherm Kinetics). Nghiên cứu thực hiện trên quy mô bãi chôn lấp Kiêu Kỵ với diện tích 5,4 ha (khu vực 5 ô chôn lấp chiếm $31.190\ m^2$), phân tích chuỗi số liệu quan trắc địa kỹ thuật, thủy văn và mẫu nước rỉ rác xuyên suốt giai đoạn 2008–2017, kết hợp hai lỗ khoan chuyên sâu LK4 và LK5. Đóng góp đột phá của luận án là xác lập mô hình giải tích - số 3D dự báo chính xác vành đai lan truyền ô nhiễm kim loại nặng, phục vụ thiết lập hành lang cách ly an toàn môi trường.

Literature Review và Positioning

Nghiên cứu về nước rỉ rác và ô nhiễm bãi chôn lấp trên thế giới đã trải qua nhiều thập kỷ phát triển với các trường phái lý thuyết chuyên sâu:

[Mô hình Cân bằng Nước Cổ điển] (WBM / Farquhar, 1989)
[Định lượng Thủy văn Chuyên sâu] (HELP - US EPA / Schroeder et al., 1994)
[Mô hình Lan truyền Đối lưu - Phân tán 1D/2D] (Kjeldsen et al., 2002; Christensen et al., 2001)
[MÔ HÌNH TOÁN - SỐ 3D ĐA CHIỀU CỦA LUẬN ÁN] 
(Tích hợp Vector Gradient ∇C + PTHH + Điều kiện Khí hậu Nhiệt đới Gió mùa Việt Nam)

Trường phái thủy văn bãi chôn lấp khởi đầu với các công trình kinh điển của Farquhar (1989) về đường cong thủy văn nước rỉ rác và phương pháp cân bằng nước nối tiếp (Serial Water Balance - SWB). Tiếp đó, Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA) phát triển mô hình HELP (Hydrologic Evaluation of Landfill Performance) nhằm mô phỏng cân bằng ẩm một chiều. Tuy nhiên, các tranh biện học thuật lớn đã nảy sinh giữa hai quan điểm đối nghịch:

  • Quan điểm 1 (Schroeder et al., 1994): Cho rằng mô hình cân bằng ẩm HELP phản ánh chuẩn xác lưu lượng rò rỉ dựa trên các thông số khí tượng định lượng hàng ngày.
  • Quan điểm 2 (Blight et al., 1999; Ehrig, 1989): Phản biện rằng HELP thường ước tính thiên cao (overestimate) lượng nước rỉ rác tại các vùng khí hậu nhiệt đới hoặc bán khô hạn do bỏ qua khả năng đầm nén dị hướng của rác thải thực tế và sự suy thoái phi tuyến tính của lớp che phủ bề mặt.

Về mặt hóa - sinh môi trường, các nghiên cứu dài hạn của Ehrig (1989) và Krumpelbeck & Ehrig (1999) trên 50 bãi chôn lấp tại CHLB Đức chỉ ra rằng nồng độ amoni ($NH_4^+$) duy trì ở mức cao trung bình $740\ mg/L$ và hầu như không suy giảm ngay cả sau 30 năm đóng bãi do cơ chế rửa trôi chậm, trong khi nồng độ KLN giảm mạnh khi bãi chuyển từ pha axit ($pH < 6,5$) sang pha metan hóa ($pH = 7,5 - 8,5$) do hiện tượng kết tủa hydroxit và cacbonat. Nghiên cứu của Åkesson & Nilsson (1997) tại Thụy Điển và Chu et al. (1994) tại Hồng Kông chứng minh tính biến thiên mùa rõ rệt: nồng độ chất ô nhiễm loãng hơn trong mùa mưa nhưng tổng tải lượng phát thải ra môi trường xung quanh lại đạt đỉnh.

So sánh với các nghiên cứu quốc tế tại các nước đang phát triển (Joshua et al. tại BCL Thohoyandou - Nam Phi, Semur - Ấn Độ, Jebel Chakir - Tunisia), nước rỉ rác tại BCL Kiêu Kỵ có hàm lượng chất hữu cơ và kim loại nặng rất cao: $COD$ đạt đỉnh lên tới hàng chục nghìn $mg/L$ ở giai đoạn axit hóa trước khi giảm dần về sau năm 2012. Định vị của luận án nằm ở điểm giao thoa giữa cơ học chất lỏng trong môi trường xốp và hóa học môi trường đất: chuyển đổi các lý thuyết 1D truyền thống thành hệ phương trình đạo hàm riêng 3D, giải quyết triệt để tính dị hướng của địa tầng đồng bằng sông Hồng (sét pha dày 6m phủ trên cát hạt mịn - thô) dưới tác động của lượng mưa cực đoan lên tới $2.268\ mm/năm$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng Lý thuyết Lan truyền Thủy động lực học (Hydrodynamic Dispersion Theory) của Bear (1972) và Lý thuyết Cân bằng Nước Bãi chôn lấp bằng các đóng góp mang tính nền tảng:

  1. Thiết lập công thức cân bằng nước mở rộng ($LLC$): Luận án hoàn thiện công thức cân bằng nước phù hợp với điều kiện Việt Nam: $$LLC = I_C + W_{GC} + W_D + W_{EM}$$ Trong đó định lượng chi tiết lượng nước mưa xâm nhập ($I_C$), độ ẩm vật liệu che phủ ($W_{GC}$), lượng ẩm giải phóng do đầm nén cơ học ($W_D$) và lượng dung dịch sinh học bổ sung ($W_{EM}$), đồng thời tính đến hệ số suy thoái của lớp che phủ đỉnh bãi sau khi đóng cửa.
  2. Khái quát hóa toán học qua trường Gradient nồng độ: Ứng dụng toán tử vi phân không gian $\nabla C(x, y, z, t)$ cùng các toán tử phân kỳ $(\nabla \cdot)$ và nhân vô hướng tích hợp để mô tả liên tục thông lượng đối lưu ($\mathbf{J}{adv} = \mathbf{v} \cdot C$) và thông lượng phân tán ($\mathbf{J}{disp} = -\mathbf{D} \cdot \nabla C$).
  3. Chứng minh hiệu ứng cản dị pha (Retardation Dynamics): Chứng minh trên thực nghiệm vai trò của hệ số làm chậm $R = 1 + (\rho_b / \theta) K_d$, giải thích rõ sự phân tách vận tốc lan truyền giữa ion không bị hấp phụ ($R=1$) và ion kim loại nặng chịu tương tác tĩnh điện/trao đổi cation ($R=2$ đối với $As, Cd$).

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp ba lý thuyết trụ cột:

  • Lý thuyết Dòng thấm trong môi trường xốp chưa bão hòa và bão hòa (Richards & Darcy Laws): Xác định trường vận tốc thấm $\mathbf{v}(u, v, w)$ trong các lớp đất nền có hệ số thấm $k$ biến thiên từ $10^{-7}$ đến $10^{-4}\ cm/s$.
  • Lý thuyết Lan truyền Chất tan (Solute Transport Theory): Kết hợp tensor hệ số khuếch tán phân tán đa hướng ($D_{xx}, D_{yy}, D_{zz}$).
  • Lý thuyết Địa hóa Môi trường (Environmental Geochemistry): Phân loại trạng thái linh động của kim loại nặng theo Kabata-Pendias (1984), chỉ rõ tính linh động của $As, Cd, Cr, Pb$ dưới môi trường oxy hóa khử thay đổi.

Điều kiện biên và giới hạn áp dụng (boundary conditions) được xác định rõ ràng: biên loại 1 (Dirichlet) tại đáy ô chôn lấp với nồng độ rò rỉ $C_0(t)$ biến thiên theo tuổi bãi; biên loại 2 (Neumann) tại các tầng cản sét sâu; và miền tính toán không gian 3D giới hạn trong phạm vi tầng chứa nước Holocen ($qh_2, qh_1$) và Pleistocen ($qp_2, qp_1$).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án vận hành theo hệ hình thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), đi từ quan trắc thực nghiệm thực địa, phân tích hóa nghiệm tiêu chuẩn đến mô hình hóa toán học số trị nâng cao. Thiết kế nghiên cứu đa tầng (Multi-level experimental design) kết hợp giữa:

  • Tầng 1: Động thái phát sinh chất ô nhiễm tại nguồn (nước rỉ rác 5 ô chôn lấp, diện tích $31.190\ m^2$).
  • Tầng 2: Sự tích lũy và chuyển hóa địa hóa học trong tầng đất đệm và đất tự nhiên (chiều sâu khoan thăm dò đến 6–10m).
  • Tầng 3: Mô phỏng động học lan truyền 3D trong lòng đất và tầng ngậm nước lân cận sông Cầu Bây.
       CẤU TRÚC ĐỊA TẦNG VÀ LỚP CHỐNG THẤM TẠI BÃI CHÔN LẤP KIÊU KỴ

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích mẫu tuân thủ nghiêm ngặt các tiêu chuẩn quốc gia và quốc tế:

  • Phương pháp lấy mẫu: Mẫu nước rỉ rác được thu thập định kỳ theo mùa (mùa mưa từ tháng 5–10, mùa khô từ tháng 11–4 năm sau) trong suốt chu kỳ vận hành 2008–2017. Lấy mẫu đất lõi khoan nguyên dạng tại hai lỗ khoan chiến lược LK4 và LK5 xuyên qua các tầng đất sét pha nâu vàng, sét pha xám xanh và cát hạt nhỏ xám tro.
  • Thiết bị và phân tích phòng thí nghiệm: Toàn bộ mẫu được phân tích tại 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). Hàm lượng kim loại nặng ($As, Pb, Cd, Cr, Cu, Zn, Ni$) được phân tích bằng hệ thống quang phổ hấp thụ nguyên tử (AAS) và ICP-MS theo phương pháp chuẩn Standard Methods for the Examination of Water and Wastewater (APHA).
  • Tam giác đạc (Triangulation): Kết hợp đối chiếu chéo giữa dữ liệu hóa nghiệm thực địa, số liệu khí tượng thủy văn 14 năm (lượng mưa trung bình $1.612\ mm/năm$, độ ẩm trung bình $84,5%$) và mô phỏng số học.

Data và phân tích

Đặc tính mẫu đất hiện trường thể hiện rõ qua mặt cắt địa chất: Lớp 1 (đất bề mặt lẫn sét pha nâu hồng, dày $0,8\ m$), Lớp 2 (sét pha nâu vàng đáy hồ, dày $1,3\ m$, dung trọng $\rho = 2,0\ g/cm^3$), Lớp 3 (sét pha xám xanh hồng tro, dày $0,7\ m$), Lớp 4 (sét pha xám nâu, dày $0,6\ m$), Lớp 5 (cát hạt trung xám nâu, dày $4,0\ m$), Lớp 6 (cát hạt nhỏ xám tro, dày $3-4\ m$).

Phương pháp phần tử hữu hạn (PTHH) được lập trình trên phần mềm MATLAB. Phương trình vi phân đạo hàm riêng chi phối quá trình lan truyền 3D có dạng: $$R \frac{\partial C}{\partial t} = \frac{\partial}{\partial x}\left(D_x \frac{\partial C}{\partial x}\right) + \frac{\partial}{\partial y}\left(D_y \frac{\partial C}{\partial y}\right) + \frac{\partial}{\partial z}\left(D_z \frac{\partial C}{\partial z}\right) - \left(v_x \frac{\partial C}{\partial x} + v_y \frac{\partial C}{\partial y} + v_z \frac{\partial C}{\partial z}\right)$$

Miền tính toán được chia lưới phần tử tứ diện (tetrahedral elements) hoặc lăng trụ hữu hạn. Phân tích độ nhạy (robustness check) được thực hiện với các kịch bản hệ số hấp phụ làm chậm $R = 1$ (không tương tác hấp phụ) và $R = 2$ (có hấp phụ ion), cùng các kịch bản suy giảm lưu lượng nước rác sau khi phủ đỉnh bãi.

Phát hiện đột phá và implications

Những phát hiện then chốt

Dữ liệu phân tích thực nghiệm và mô phỏng số mang lại 5 phát hiện cốt lõi:

  1. Quy luật chuyển hóa sinh hóa của nước rỉ rác: Nước rác tại BCL Kiêu Kỵ trải qua sự chuyển pha rõ rệt. Trong giai đoạn 2008–2010 (pha axit), tỷ lệ $BOD_5/COD$ đạt mức $0,59 - 0,65$, nồng độ $COD$ vượt $20.000\ mg/L$. Giai đoạn sau 2012–2017 (pha metan hóa), $pH$ tăng lên $7,8 - 8,4$, tỷ lệ $BOD_5/COD$ giảm sâu xuống $0,1 - 0,2$, phản ánh sự chiếm ưu thế của các hợp chất hữu cơ mạch vòng khó phân hủy (axit humic và fulvic).
  2. Động thái hàm lượng kim loại nặng trong nước rỉ rác: Nồng độ các kim loại nặng hòa tan trong nước rỉ rác có xu hướng giảm theo thời gian do hiện tượng kiềm hóa bãi rác thúc đẩy quá trình kết tủa. Tuy nhiên, các anion kim loại như Asen ($As$) và ion Cadimi ($Cd$) vẫn duy trì nồng độ đáng kể trong pha lỏng ($As$ dao động từ $0,05 - 0,25\ mg/L$, $Cd$ từ $0,005 - 0,126\ mg/L$), tiềm ẩn nguy cơ thẩm thấu trường diễn.
  3. Mức độ tích lũy kim loại nặng trong nền đất: Kết quả phân tích mẫu đất tại hai lỗ khoan LK4 và LK5 cho thấy nồng độ tổng số của $As, Cr, Pb$ tại các lớp đất sát đáy ô chôn lấp đã vượt ngưỡng quy chuẩn kỹ thuật quốc gia về chất lượng đất (QCVN 03-MT:2015/BTNMT), chứng minh sự suy giảm cục bộ khả năng ngăn chặn của lớp màng chống thấm sau chu kỳ dài tiếp xúc với tải trọng rác và hóa chất.
  4. Quy luật lan truyền không gian 3D của Asen ($As$): Kết quả chạy mô hình PTHH trên MATLAB chỉ ra rằng tại tháng thứ 10 sau khi xảy ra sự cố rò rỉ, nồng độ $As$ lan truyền theo phương thẳng đứng ($Z$) qua lớp sét pha $6\ m$ với tốc độ chậm khi $R=2$, nhưng khi tiếp cận tầng cát hạt nhỏ ($qh_1$) có hệ số thấm cao, gradient nồng độ mở rộng nhanh chóng theo phương ngang ($X, Y$) hướng về phía sông Cầu Bây (cách bãi chôn lấp $120\ m$).
  5. Hiệu ứng kìm giữ của lớp đất á sét tự nhiên: So sánh giữa kịch bản $R=1$ và $R=2$ cho thấy dung tích hấp phụ trao đổi cation của tầng đất sét pha Kiêu Kỵ đóng vai trò như một màng lọc địa chất tự nhiên, làm giảm đỉnh nồng độ lan truyền của $Cd$ và $Pb$ tới $45-60%$ so với tính toán thuần túy đối lưu thủy lực.
                 NỒNG ĐỘ TƯƠNG ĐỐI CỦA ASEN THEO ĐỘ SÂU (Z) VÀ THỜI GIAN
               
      Nồng độ C/C0
            0           2           4           6 (Độ sâu Z - mét)
                     [Lớp Sét Pha]          [Tầng Cát Chứa Nước]

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh rằng mô hình vận chuyển ô nhiễm trong đất nhiệt đới bắt buộc phải tích hợp đồng thời trường gradient nồng độ $\nabla C$ và hệ số hấp phụ phi tuyến tính, bác bỏ các mô hình 1D đơn giản hóa vốn thường đánh giá sai lệch biên lan truyền ô nhiễm.
  • Về mặt phương pháp luận: Cung cấp bộ thuật toán PTHH tùy biến trên MATLAB có khả năng mở rộng để mô phỏng sự lan truyền của các hợp chất hữu cơ khó phân hủy (XOCs, Dioxin, Furan) tại các bãi rác công nghiệp.
  • Về mặt thực tiễn và chính sách:
    • Khuyến nghị Bộ Xây dựng và Bộ Tài nguyên & Môi trường bổ sung yêu cầu bắt buộc tính toán lan truyền 3D trong đánh giá tác động môi trường (ĐTM) khi đóng cửa các bãi chôn lấp quy mô cấp tỉnh.
    • Thiết lập vành đai cách ly vệ sinh tối thiểu $150 - 200\ m$ xung quanh BCL Kiêu Kỵ dọc theo hành lang thoát nước sông Cầu Bây nhằm ngăn ngừa nguy cơ nhiễm độc chuỗi thức ăn nông nghiệp.

Limitations và Future Research

Luận án thẳng thắn chỉ rõ 4 giới hạn nghiên cứu:

  1. Giới hạn số lượng lỗ khoan sâu: Do điều kiện kinh phí và tính chất phức tạp của bãi rác đã đóng cửa, số lượng lỗ khoan lõi sâu đại diện tập trung tại hai vị trí trọng điểm LK4 và LK5, chưa bao phủ toàn bộ lưới địa chất 5,4 ha.
  2. Giả định đồng nhất hóa cục bộ lớp đất: Mô hình số 3D giả định hệ số thấm và hệ số khuếch tán đồng nhất theo từng lớp địa tầng chính, chưa mô phỏng hết các vết nứt vi mô (macropores) phát sinh do hiện tượng lún không đều của khối rác.
  3. Chưa xét đến tương tác đa kim loại cạnh tranh hấp phụ: Nghiên cứu mô phỏng đơn lẻ từng kim loại ($As, Cd$) với hệ số $R$ không đổi, chưa tích hợp mô hình nhiệt động hóa học cạnh tranh vị trí hấp phụ phức tạp giữa $Fe, Mn, Pb, Zn$.
  4. Giới hạn chuỗi dữ liệu thời gian thực sau đóng bãi: Số liệu quan trắc tập trung đến năm 2017, cần tiếp tục cập nhật dữ liệu tự động liên tục trong giai đoạn ổn định hóa kéo dài (post-closure care).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào:

  • Tích hợp mô hình dòng chảy đa pha (khí sinh học $CH_4/CO_2$ - nước rỉ rác - hạt keo đất) trong khối chôn lấp.
  • Ứng dụng kỹ thuật đồng vị bền để truy xuất chính xác nguồn gốc kim loại nặng hòa tan trong tầng chứa nước ngầm Pleistocen.
  • Đánh giá hiệu quả xử lý tăng cường tại chỗ bằng công nghệ tường phản ứng thẩm thấu sinh học (Permeable Reactive Barrier - PRB) kết hợp thực vật phục hồi (Phytoremediation).

Tác động và ảnh hưởng

Nghiên cứu mang lại những giá trị tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật: Dự kiến đóng góp nguồn trích dẫn quan trọng cho các công trình nghiên cứu về địa kỹ thuật môi trường, quản lý chất thải rắn tại khu vực Đông Nam Á; cung cấp bộ dữ liệu thực chứng chuẩn hóa về thành phần nước rỉ rác bãi rác nhiệt đới.
  • Chuyển giao công nghệ ngành: Giúp các đơn vị tư vấn thiết kế hạ tầng đô thị (ngành cấp thoát nước, môi trường đô thị) ứng dụng công thức tính toán lượng nước rỉ rác mở rộng ($LLC$) để định cỡ chuẩn xác các trạm xử lý nước rác, tránh lãng phí vốn đầu tư hoặc quá tải cục bộ.
  • Ảnh hưởng chính sách quốc gia: Cung cấp cơ sở khoa học định lượng phục vụ công tác rà soát, sửa đổi các quy chuẩn thiết kế bãi chôn lấp (TCXDVN 261:2001 và TCVN 6696:2009), hỗ trợ các nhà quản lý xây dựng kịch bản đóng bãi an toàn cho hơn 900 bãi chôn lấp hiện hữu trên cả nước.
  • Lợi ích xã hội: Bảo vệ trực tiếp nguồn tài nguyên nước dưới đất - nguồn cấp nước sinh hoạt và tưới tiêu nông nghiệp trọng yếu cho hàng trăm nghìn hộ dân huyện Gia Lâm và lưu vực sông Cầu Bây, giảm thiểu nguy cơ phơi nhiễm kim loại nặng gây bệnh mãn tính.

Đối tượng hưởng lợi

  1. Nghiên cứu sinh và Giảng viên Chuyên ngành: Kế thừa phương pháp luận mô hình hóa số trị 3D bằng phương pháp PTHH trên MATLAB và hệ phương trình vi phân lan truyền đối lưu - phân tán để phát triển các hướng nghiên cứu chuyên sâu về ô nhiễm chất thải nguy hại.
  2. Kỹ sư Thiết kế và Doanh nghiệp Môi trường Đô thị (URENCO): Sử dụng trực tiếp công thức cân bằng nước mở rộng ($LLC$) và thông số địa kỹ thuật thực tế để thiết kế tối ưu hệ thống màng phủ đáy HDPE, tầng thu gom nước rác và hệ thống giếng quan trắc khí tượng - thủy văn.
  3. Cơ quan Hoạch định Chính sách (Bộ Xây dựng, Bộ Tài nguyên & Môi trường, Sở TN&MT Hà Nội): Sở hữu công cụ mô phỏng dự báo định lượng phục vụ công tác quy hoạch sử dụng đất sau đóng bãi, thẩm định báo cáo ĐTM và ban hành các ngưỡng an toàn khoảng cách môi trường.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc chuẩn hóa phương trình toán học mô tả quá trình lan truyền chất ô nhiễm trong không gian 3D dựa trên việc biểu diễn giải tích trường vectơ gradient nồng độ $\nabla C(x, y, z, t)$ tích hợp với phương pháp phần tử hữu hạn. Nghiên cứu đã mở rộng Lý thuyết Đối lưu - Phân tán Thủy động lực học cổ điển của Bear (1972) bằng cách bổ sung hệ số cản hấp phụ phi tuyến $R$ và giải thuật toán số thích ứng với cấu trúc địa tầng bất đẳng hướng phức tạp (xen kẹp sét pha và cát mịn) dưới điều kiện tải lượng ô nhiễm phân rã theo tuổi bãi.

2. Đột phá về phương pháp luận của luận án khi so sánh với hai nghiên cứu quốc tế tiền nhiệm?

So với mô hình định lượng thủy văn HELP một chiều của US EPA (Schroeder et al., 1994) vốn bỏ qua sự biến thiên cấu trúc vật liệu phủ khi nâng chiều cao rác, và mô hình giải tích 2D của Christensen et al. (2001) chỉ xét tầng đồng nhất, phương pháp luận của luận án vượt trội ở việc:

  • Tích hợp công thức cân bằng ẩm mở rộng $LLC$ tính đến tác động của dung dịch vi sinh bổ sung ($W_{EM}$) và đầm nén cơ học ($W_D$).
  • Xây dựng mạng lưới phần tử hữu hạn 3D giải đồng thời trường thủy lực và trường nồng độ kim loại nặng ($As, Cd$) trong điều kiện biên biến thiên theo thời gian thực 10 năm.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ?

Phát hiện bất ngờ nhất là sự đối lập sâu sắc về nồng độ giữa nước rỉ rác và môi trường đất đáy bãi: trong khi nồng độ kim loại nặng trong nước rỉ rác suy giảm mạnh ở giai đoạn sau năm 2012 do $pH$ tăng cao (pha metan hóa gây kết tủa), thì hàm lượng tích lũy $As, Cr, Pb$ trong các tầng đất sét pha tại lỗ khoan LK4 và LK5 lại vượt ngưỡng QCVN 03-MT:2015/BTNMT. Điều này chứng minh rằng đáy bãi chôn lấp đóng vai trò như một "bể tích tụ địa hóa", nơi kim loại nặng bị giữ lại trong bùn sét nhưng sẵn sàng tái hòa tan và phát tán tức thời nếu điều kiện oxy hóa khử ($Eh$) hoặc độ chua ($pH$) của nước ngầm thay đổi.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình tái lập hoàn chỉnh bao gồm:

  • Toàn bộ tham số đầu vào địa kỹ thuật: độ ẩm tự nhiên, dung tích tự nhiên ($2,0\ g/cm^3$), mô đun biến dạng, hệ số thấm $k$ của 6 lớp đất.
  • Thuật toán rời rạc hóa phương trình vi phân đạo hàm riêng, điều kiện biên Dirichlet/Neumann và các đoạn mã hàm MATLAB thực thi phương pháp PTHH.
  • Quy trình lấy mẫu đất nguyên dạng và phân tích kim loại nặng chuẩn hóa theo APHA và TCVN.

5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?

Chương trình 10 năm hướng tới ba trọng tâm:

  • Giai đoạn 1 (Năm 1-3): Tích hợp công nghệ cảm biến IoT quan trắc tự động liên tục thế oxy hóa khử ($Eh$), $pH$ và nồng độ ion kim loại tại các giếng ngầm xung quanh BCL Kiêu Kỵ.
  • Giai đoạn 2 (Năm 4-6): Phát triển mô hình mô phỏng dòng chảy ngầm kết hợp phản ứng sinh - địa - hóa đa cấu tử (Reactive Transport Modeling - PHREEQC/TOUGHREACT).
  • Giai đoạn 3 (Năm 7-10): Thử nghiệm và nhân rộng công nghệ xử lý ô nhiễm tầng ngậm nước bằng rào cản phản ứng thẩm thấu nano-sắt hóa trị 0 ($nZVI$) và phục hồi sinh thái cảnh quan bãi rác sau đóng cửa.

Kết luận

  1. Luận án đã giải quyết căn cơ bài toán cấp bách về quản lý môi trường bãi chôn lấp chất thải rắn thông qua việc phát triển công thức cân bằng nước mở rộng $LLC = I_C + W_{GC} + W_D + W_{EM}$, phản ánh sát thực tế vận hành tại Việt Nam.
  2. Thiết lập thành công mô hình toán học và giải thuật phần tử hữu hạn 3D trên nền tảng MATLAB, mô phỏng sinh động và định lượng trường lan truyền nồng độ kim loại nặng $\nabla C(x, y, z, t)$ trong lòng đất.
  3. Xác định rõ quy luật suy thoái chất lượng nước rỉ rác theo thời gian: tỷ lệ $BOD_5/COD$ giảm từ $>0,5$ xuống $<0,2$, nồng độ chất hữu cơ giảm mạnh nhưng các kim loại nguy hại ($As, Cd, Cr, Pb$) tích lũy cục bộ ở mức cao trong lớp đất nền.
  4. Cung cấp luận cứ khoa học thực chứng tại bãi chôn lấp Kiêu Kỵ (Gia Lâm, Hà Nội), chứng minh vai trò cản hấp phụ ($R=2$) của tầng sét pha $6\ m$ và cảnh báo nguy cơ lan truyền nhanh khi chất ô nhiễm xâm nhập tầng cát ngậm nước Holocen/Pleistocen.
  5. Mở ra 3 hướng nghiên cứu ứng dụng đột phá: mô hình hóa dòng thấm đa pha trong khối rác, thiết kế tường chắn địa kỹ thuật composite tối ưu, và quy hoạch vành đai cách ly an toàn cho toàn bộ hệ thống bãi chôn lấp đóng cửa tại Việt Nam.
  6. Kết quả nghiên cứu là tài liệu chuẩn mực về lý thuyết và thực tiễn, đóng góp trực tiếp vào mục tiêu kiểm soát ô nhiễm đất và nước dưới đất, kiến tạo giải pháp kỹ thuật bền vững cho công tác quản lý chất thải rắn quốc gia.