Luận án tiến sĩ nghiên cứu cơ chế dịch chuyển asen từ tầng chứa nước holocene vào tầng chứa nước pleistocene lấy ví dụ vùng thạch thất đan phượng hà nội

Luận án tiến sĩ nghiên cứu cơ chế dịch chuyển asen từ tầng chứa nước Holocene vào tầng chứa nước Pleistocene tại Thạch Thất, Đan Phượng, Hà Nội.

• Tác giả: Trần Vũ Long
• Người hướng dẫn: TS. Phạm Quý Nhân
• Trường học: Trường Đại Học Mỏ - Địa Chất
• Chuyên ngành: Kỹ thuật địa chất
• Đề tài: Nghiên Cứu Cơ Chế Dịch Chuyển Asen Trong Nước Dưới Đất Tại Thạch Thất - Đan Phượng, Hà Nội
• Loại tài liệu: luận án tiến sĩ địa chất
• Năm xuất bản: 2019
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• DANH MỤC CÁC KÝ HIỆU, CÁC TỪ VIẾT TẮT
• DANH MỤC CÁC BẢNG
• DANH MỤC CÁC HÌNH
• MỞ ĐẦU
• 1. CHƯƠNG 1: TỔNG QUAN VỀ NGHIÊN CỨU ASEN TRONG NƯỚC DƯỚI ĐẤT
• 1.1. Các công trình nghiên cứu trên thế giới
• 1.1.1. Nhóm đánh giá hiện trạng và xác định cơ chế giải phóng As vào NDĐ
• 1.1.2. Nhóm nghiên cứu con đường dịch chuyển As và khả năng tổn thương của TCN bên dưới trước ô nhiễm As
• 1.2. Nghiên cứu As ở Việt Nam và khu vực nghiên cứu
• 1.2.1. Nghiên cứu As ở Việt Nam
• 1.2.2. Các nghiên cứu As tại khu vực nghiên cứu
• 1.3. Các cơ chế dịch chuyển của As trong NDĐ
• 1.3.1. Khái niệm chung nghiên cứu cơ chế dịch chuyển As từ TCN qh vào TCN qp
• 1.3.2. Cơ chế thuỷ động lực
• 1.3.3. Cơ chế thuỷ địa hoá
• 2. CHƯƠNG 2: ĐẶC ĐIỂM PHÂN BỐ VÀ CÁC NHÂN TỐ ẢNH HƯỞNG ĐẾN SỰ TỒN TẠI VÀ DỊCH CHUYỂN ASEN TRONG NƯỚC DƯỚI ĐẤT KHU VỰC NGHIÊN CỨU
• 2.1. Khái quát về khu vực nghiên cứu
• 2.2. Đặc điểm phân bố As trong khu vực nghiên cứu
• 2.3. Các nhân tố ảnh hưởng đến sự tồn tại và dịch chuyển của As trong NDĐ
• 2.3.1. Thành phần khoáng vật - thạch học của trầm tích
• 2.3.2. Tuổi trầm tích và tuổi của NDĐ
• 2.3.3. Đặc điểm ĐCTV
• 2.3.4. Thuỷ địa hoá học
• 2.3.5. Vi sinh vật
• 2.3.6. Khai thác NDĐ
• 2.3.7. Thổ nhưỡng và sử dụng đất
• 3. CHƯƠNG 3: CƠ CHẾ DỊCH CHUYỂN ASEN TRONG NƯỚC DƯỚI ĐẤT
• 3.1. Vấn đề và phương pháp nghiên cứu
• 3.1.1. Vấn đề nghiên cứu
• 3.1.2. Phương pháp nghiên cứu
• 3.2. Cơ sở lý thuyết mô hình số dịch chuyển As hoà tan trong NDĐ
• 3.2.1. Lý thuyết về mô hình số dòng chảy MODFLOW
• 3.2.2. Lý thuyết mô hình số dịch chuyển As hoà tan trong NDĐ MT3D - USGS
• 3.3. Mô hình số dịch chuyển As trong NDĐ khu vực nghiên cứu
• 3.3.1. Mô hình số dòng chảy
• 3.3.2. Mô hình dịch chuyển As trong NDĐ khu vực nghiên cứu
• 3.4. Cơ chế dịch chuyển As từ TCN qh vào TCN qp
• 3.4.1. Mô hình dịch chuyển As tại Thạch Thất - Đan Phượng
• 3.4.2. So sánh với kết quả nghiên cứu mô hình dịch chuyển As tại Nam Dư
• 4. CHƯƠNG 4: CÁC GIẢI PHÁP GIẢM THIỂU Ô NHIỄM ASEN ĐỐI VỚI TẦNG CHỨA NƯỚC PLEISTOCEN
• 4.1. Nguyên tắc chung của các giải pháp giảm thiểu ô nhiễm As
• 4.2. Giải pháp khắc phục và hạn chế đối với khu vực nghiên cứu Thạch Thất - Đan Phượng
• 4.3. Công nghệ xử lý As trong nước
• 4.3.1. Các dạng công nghệ xử lý As trong nước
• 4.3.2. Công nghệ xử lý As ở Việt Nam và trên thế giới
• KẾT LUẬN VÀ KIẾN NGHỊ
• DANH MỤC CÁC CÔNG TRÌNH KHOA HỌC CỦA TÁC GIẢ LIÊN QUAN ĐẾN LUẬN ÁN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Nghiên Cứu Asen Trong Nước Ngầm Thách Thức Giải Pháp

Asen (As, arsenic, thạch tín) là một nguyên tố vi lượng, độc hại khi vượt quá ngưỡng cho phép. Ô nhiễm As trong nước dưới đất (NDĐ) là vấn đề nhức nhối, đặc biệt tại Đồng bằng Bắc Bộ (ĐBBB), nơi NDĐ từ các tầng chứa nước (TCN) Đệ Tứ là nguồn cung cấp chính cho sinh hoạt và sản xuất. Mặc dù không mới, vấn đề này vẫn thu hút sự quan tâm của giới khoa học trong và ngoài nước. Các nghiên cứu tập trung vào cơ chế hình thành, giải phóng As và đề xuất giải pháp giảm thiểu hiệu quả. Đối tượng nghiên cứu hàng đầu là As trong trầm tích và NDĐ. Các nghiên cứu trước đây chỉ ra rằng hàm lượng As trong TCN Holocen (qh) cao hơn nhiều so với TCN Pleistocen (qp).

1.1. Hiện Trạng Ô Nhiễm Asen Tại Đồng Bằng Sông Hồng

Nhiều nghiên cứu đã chỉ ra rằng As chủ yếu được giải phóng vào NDĐ trong môi trường khử mạnh. Cơ chế chính là khử hòa tan sắt oxi hydroxit có hấp phụ As trên trầm tích, do hoạt động phân giải vật chất hữu cơ trong TCN qh. Tuy nhiên, quá trình dịch chuyển NDĐ chứa As từ TCN qh vào TCN qp, và ảnh hưởng của các yếu tố tự nhiên, nhân tạo, đặc biệt là khai thác NDĐ quy mô công nghiệp, vẫn chưa được nghiên cứu chi tiết. Nguy cơ gia tăng hàm lượng As trong NDĐ của TCN qp vượt quá tiêu chuẩn là hiện hữu, đe dọa an ninh nguồn nước.

1.2. Tầm Quan Trọng Của Nghiên Cứu Cơ Chế Dịch Chuyển Asen

Việc nghiên cứu cơ chế dịch chuyển As từ TCN qh vào TCN qp, dự báo biến đổi theo thời gian và đề xuất giải pháp giảm thiểu nguy cơ là vô cùng cấp thiết. Nghiên cứu này sẽ đánh giá ảnh hưởng của các yếu tố tự nhiên và nhân tạo đến sự tồn tại và dịch chuyển As trong NDĐ. Đồng thời, nó sẽ nghiên cứu cơ chế và quá trình dịch chuyển As từ TCN qh vào TCN qp trong điều kiện tự nhiên, từ đó đề xuất giải pháp giảm thiểu và bảo vệ TCN qp.

II. Phương Pháp Nghiên Cứu Dịch Chuyển Asen Mô Hình Thực Nghiệm

Nghiên cứu này sử dụng kết hợp các phương pháp tiếp cận trực tiếp, kế thừa các kết quả nghiên cứu trước đây và áp dụng các phương pháp điều tra, đánh giá hiện đại. Cách tiếp cận trực tiếp bao gồm khảo sát thực địa để nắm bắt đặc điểm địa hình, địa mạo, địa chất, địa chất thủy văn (ĐCTV), hiện trạng ô nhiễm và khai thác nước. Việc kế thừa các kết quả nghiên cứu trước đây giúp giảm khối lượng công việc và định hướng đánh giá nguồn gốc, phân bố As và các con đường dịch chuyển As trong NDĐ. Các phương pháp nghiên cứu hiện đại được áp dụng một cách hợp lý để đạt hiệu quả cao.

2.1. Thu Thập Dữ Liệu Địa Chất Thủy Văn Khu Vực Nghiên Cứu

Các phương pháp nghiên cứu bao gồm thu thập tài liệu, điều tra khảo sát sơ bộ và thành lập các điểm nghiên cứu phân bố As trên tuyến nghiên cứu. Thu thập và quan trắc các thông số về thủy văn, ĐCTV: động thái NDĐ, lượng mưa, bốc hơi, mực nước sông trong khu vực nghiên cứu. Tiến hành các công tác đo Địa vật lý và giải đoán tài liệu: Phương pháp đo điện đa cực, phương pháp trường chuyển (TEM) và phương pháp địa vật lý lỗ khoan.

2.2. Phân Tích Mẫu Trầm Tích Và Nước Ngầm Xác Định Tuổi Thành Phần

Thành lập các mặt cắt địa chất, ĐCTV. Khoan lấy mẫu nguyên dạng theo chiều sâu và phân tích mẫu trầm tích, xác định tuổi trầm tích. Xác định các thông số ĐCTV, lấy mẫu NDĐ phân tích đặc điểm thủy địa hóa và hàm lượng As. Quan trắc lâu dài động thái NDĐ tại khu vực nghiên cứu để đánh giá biến đổi hàm lượng As theo thời gian. Lấy mẫu đồng vị phóng xạ, đồng vị bền tại các lỗ khoan xác định nguồn gốc và tuổi của NDĐ.

2.3. Xây Dựng Mô Hình Dòng Chảy Và Dịch Chuyển Asen

Xây dựng các mô hình dòng chảy, dịch chuyển As trong khu vực nghiên cứu và dự báo diễn biến hàm lượng As trong NDĐ theo thời gian với điều kiện khai thác trong khu vực nghiên cứu. Đánh giá sự khác nhau giữa vận tốc di chuyển chậm và nhanh của NDĐ đối với quá trình dịch chuyển của As. Luận án sẽ đi sâu nghiên cứu cơ chế dịch chuyển của As trong NDĐ từ TCN qh tới TCN qp, ảnh hưởng của khả năng hấp phụ và giải hấp phụ của TCN qp đối với quá trình dịch chuyển này.

III. Đặc Điểm Phân Bố Asen Ảnh Hưởng Từ Yếu Tố Tự Nhiên Nhân Tạo

Khu vực nghiên cứu là Thạch Thất - Đan Phượng, Hà Nội, nơi có hàm lượng As trong NDĐ, đặc biệt trong TCN qh, rất cao so với các khu vực khác của ĐBBB. Về vị trí, khu vực nghiên cứu phân bố từ vùng rìa vào đến khu trung tâm cấu trúc bồn trũng Hà Nội. Trầm tích TCN qh tại đây khá trẻ và rất giàu hàm lượng vật chất hữu cơ hoạt động, chỉ khoảng vài trăm năm. Các TCN của khu vực nghiên cứu có mối quan hệ thủy lực chặt chẽ với nhau và với sông Hồng.

3.1. Vị Trí Địa Lý Và Đặc Điểm Địa Chất Khu Vực Thạch Thất Đan Phượng

Đồng thời trong khu vực nghiên cứu hầu như không có các công trình khai thác NDĐ tập trung với quy mô công nghiệp nên động thái của NDĐ là tự nhiên. Một tuyến mặt cắt đặc trưng được chọn để tập trung nghiên cứu. Tuyến mặt cắt này kéo dài qua các điểm nghiên cứu chính từ vùng rìa đồng bằng Phú Kim (Thạch Thất), Phụng Thượng (Phúc Thọ), Vân Cốc (Phúc Thọ) tới khu vực sát sông Hồng tại Trung Châu (Đan Phượng).

3.2. Tuyến Nghiên Cứu Đặc Trưng Phân Tích Hàm Lượng Asen Theo Chiều Sâu

Tuyến nghiên cứu chính này được lựa chọn dựa trên cơ sở đặc điểm hàm lượng As phân bố trong NDĐ theo chiều sâu, đặc điểm ĐCTV, đặc điểm và tuổi của trầm tích, đặc điểm thủy địa hóa. Tuyến nghiên cứu vuông góc với đường dòng tổng quát trong khu vực nghiên cứu, chảy từ vùng rìa tới sông Hồng.

IV. Cơ Chế Dịch Chuyển Asen Từ Tầng Holocen Vào Tầng Pleistocene

Các nghiên cứu trước đây đều chỉ ra hàm lượng As trong TCN Holocen (qh) cao hơn trong TCN Pleistocen (qp) rất nhiều và As chủ yếu được giải phóng vào NDĐ trong môi trường khử mạnh theo cơ chế khử hoà tan sắt oxi hydroxit có hấp phụ As trên trầm tích do hoạt động phân giải vật chất hữu cơ xảy ra trong TCN qh. Tuy nhiên quá trình dịch chuyển của NDĐ có hàm lượng As cao tại TCN qh vào TCN qp có thể gây ra biến đổi như thế nào đến hàm lượng As của TCN qp dưới ảnh hưởng của các nhân tố tự nhiên và nhân tạo, của các hoạt động khai thác NDĐ với quy mô công nghiệp và đặc biệt là cơ chế khống chế chính quá trình dịch chuyển này vẫn chưa được nghiên cứu một cách chi tiết.

4.1. Ảnh Hưởng Của Hoạt Động Khai Thác Nước Ngầm Đến Dịch Chuyển Asen

Nguy cơ làm gia tăng hàm lượng As trong NDĐ của TCN qp vượt quá tiêu chuẩn cho phép là hiện hữu và đe dọa tới an ninh nguồn nước cho các hoạt động kinh tế - xã hội. Chính vì lý do đó, việc nghiên cứu cơ chế dịch chuyển của As từ TCN qh vào TCN qp, dự báo biến đổi theo thời gian và đề xuất các giải pháp giảm thiểu nguy cơ đối với an ninh nguồn nước là một vấn đề rất cấp thiết.

4.2. Vai Trò Của Hệ Số Trễ Retardation Trong Quá Trình Dịch Chuyển Asen

Luận án sẽ đi sâu nghiên cứu cơ chế dịch chuyển của As trong NDĐ từ TCN qh tới TCN qp, ảnh hưởng của khả năng hấp phụ và giải hấp phụ của TCN qp đối với quá trình dịch chuyển này. Đánh giá sự khác nhau giữa vận tốc di chuyển chậm và nhanh của NDĐ đối với quá trình dịch chuyển của As.

V. Giải Pháp Giảm Thiểu Ô Nhiễm Asen Công Nghệ Xử Lý Quản Lý

Nghiên cứu này sẽ đề xuất các giải pháp giảm thiểu và bảo vệ TCN qp. Các giải pháp này bao gồm công nghệ xử lý As trong nước và các biện pháp quản lý nguồn nước hiệu quả. Việc áp dụng các giải pháp này sẽ giúp giảm thiểu nguy cơ ô nhiễm As trong NDĐ và bảo vệ an ninh nguồn nước cho các hoạt động kinh tế - xã hội.

5.1. Các Dạng Công Nghệ Xử Lý Asen Trong Nước Hiện Nay

Các dạng công nghệ xử lý As trong nước bao gồm công nghệ hấp phụ, công nghệ oxy hóa, công nghệ trao đổi ion và công nghệ màng lọc. Mỗi công nghệ có ưu và nhược điểm riêng, phù hợp với các điều kiện khác nhau. Việc lựa chọn công nghệ phù hợp cần dựa trên các yếu tố như hàm lượng As trong nước, chi phí đầu tư và vận hành, và hiệu quả xử lý.

5.2. Ứng Dụng Công Nghệ Xử Lý Asen Tại Việt Nam Và Trên Thế Giới

Công nghệ xử lý As đã được ứng dụng rộng rãi tại Việt Nam và trên thế giới. Tại Việt Nam, nhiều dự án xử lý As trong nước đã được triển khai, đặc biệt tại các khu vực có hàm lượng As cao. Các dự án này đã góp phần cải thiện chất lượng nước và bảo vệ sức khỏe người dân.

VI. Kết Luận Kiến Nghị Hướng Nghiên Cứu Tiếp Theo Về Asen

Nghiên cứu này đã đánh giá ảnh hưởng của các yếu tố tự nhiên và nhân tạo đến sự tồn tại và dịch chuyển As trong NDĐ. Đồng thời, nó đã nghiên cứu cơ chế và quá trình dịch chuyển As từ TCN qh vào TCN qp trong điều kiện tự nhiên, từ đó đề xuất giải pháp giảm thiểu và bảo vệ TCN qp. Tuy nhiên, vẫn còn nhiều vấn đề cần được nghiên cứu sâu hơn.

6.1. Đánh Giá Rủi Ro Ô Nhiễm Asen Trong Tương Lai

Cần có các nghiên cứu đánh giá rủi ro ô nhiễm As trong tương lai, đặc biệt trong bối cảnh biến đổi khí hậu và gia tăng khai thác NDĐ. Các nghiên cứu này sẽ giúp đưa ra các biện pháp phòng ngừa và ứng phó kịp thời.

6.2. Nghiên Cứu Ảnh Hưởng Của Vi Sinh Vật Đến Quá Trình Dịch Chuyển Asen

Cần có các nghiên cứu về ảnh hưởng của vi sinh vật đến quá trình dịch chuyển As. Vi sinh vật có thể đóng vai trò quan trọng trong việc giải phóng và hấp phụ As, ảnh hưởng đến quá trình dịch chuyển của As trong NDĐ.

Nội dung chính

Tổng quan về luận án

Ô nhiễm asen (As) tự nhiên trong môi trường nước dưới đất (NDĐ) tại các bồn trũng châu thổ Đệ Tứ là một trong những thách thức địa chất môi trường và sức khỏe cộng đồng nghiêm trọng nhất trên phạm vi toàn cầu. Tại Đồng bằng Bắc Bộ (ĐBBB), nguồn nước dưới đất khai thác từ các tầng chứa nước lỗ hổng đóng vai trò huyết mạch đối với nhu cầu sinh hoạt và phát triển kinh tế - xã hội. Luận án tiến sĩ kỹ thuật địa chất của tác giả Trần Vũ Long với đề tài "Nghiên cứu cơ chế dịch chuyển Asen từ tầng chứa nước Holocen vào tầng chứa nước Pleistocen, lấy ví dụ vùng Thạch Thất - Đan Phượng, Hà Nội" (Mã số: 9520501; Người hướng dẫn: PGS.TS Phạm Quý Nhân và PGS.TS Flemming Larsen) là công trình khoa học tiên phong giải quyết bài toán cơ chế lan truyền chất ô nhiễm trong điều kiện thủy động lực tự nhiên chưa bị xáo trộn bởi khai thác công nghiệp quy mô lớn.

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trước công trình này là: Mặc dù các nghiên cứu trước đây (như Berg et al., 2001, 2008; Postma et al., 2007, 2012; Đỗ Trọng Sự, 1993; Đỗ Văn Bình, 2007) đã xác định rõ hiện trạng ô nhiễm và cơ chế giải phóng As vào tầng chứa nước Holocen ($q_h$) thông qua phản ứng khử hòa tan oxit hydroxit sắt, nhưng quá trình dịch chuyển thẳng đứng của As từ tầng $q_h$ xuống tầng chứa nước Pleistocen ($q_p$) – đối tượng khai thác nước chính – dưới sự kiểm soát đồng thời của trường thủy động lực học tự nhiên và các phản ứng thủy địa hóa hấp phụ/giải hấp phụ vẫn chưa được định lượng và mô hình hóa chi tiết. Các câu hỏi và giả thuyết nghiên cứu được xác lập:

• Câu hỏi nghiên cứu 1 (RQ1): Các nhân tố tự nhiên và cấu trúc địa chất thủy văn chi phối sự tồn tại và phân bố nồng độ As trong hai tầng chứa nước $q_h$ và $q_p$ tại khu vực Thạch Thất - Đan Phượng như thế nào?
• Câu hỏi nghiên cứu 2 (RQ2): Cơ chế nào đóng vai trò khống chế chính đối với sự dịch chuyển của As qua ranh giới giữa hai tầng chứa nước: đối lưu - phân tán thủy động lực hay tương tác thủy địa hóa bề mặt?
• Giả thuyết 1 (H1): Sự dịch chuyển của As hòa tan từ $q_h$ vào $q_p$ tuân theo cơ chế thủy động lực, trong đó vận tốc dòng thấm chậm trong điều kiện tự nhiên cho phép các phản ứng hấp phụ trên pha rắn của trầm tích $q_p$ đạt trạng thái cân bằng triệt để, tạo ra hiệu ứng trễ thủy địa hóa rất cao ($R \gg 1$), giúp làm chậm đáng kể quá trình xâm nhập ô nhiễm.
• Giả thuyết 2 (H2): Tuổi trầm tích và mức độ thoái hóa của vật chất hữu cơ hoạt động là động lực sinh địa hóa cốt lõi quyết định mức độ tích tụ và giải phóng As(III) vào pha lỏng.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa Thuyết dòng chảy trong môi trường rỗng Darcy, Thuyết vận chuyển chất tan đối lưu - phân tán (Advection-Dispersion Theory), Động học phản ứng oxy hóa khử trầm tích, và Thuyết đẳng nhiệt hấp phụ Langmuir/Freundlich. Về phạm vi, nghiên cứu triển khai trên tuyến mặt cắt ĐCTV dài hơn 20 km kéo dài từ vùng rìa đồng bằng tại Phú Kim (Thạch Thất), qua Phụng Thượng, Vân Cốc (Phúc Thọ) đến đới sát sông Hồng tại Trung Châu (Đan Phượng), bao phủ 64 lỗ khoan quan trắc, 121 mẫu phân tích hóa học nước, 269 mẫu đồng vị bền/phóng xạ/khí trơ và 26 mẫu định tuổi trầm tích OSL trong giai đoạn 2009–2013 thuộc dự án quốc tế VietAS pha II (DANIDA tài trợ).

Literature Review và Positioning

Nghiên cứu ô nhiễm và vận chuyển Asen trên thế giới phân hóa thành hai trường phái chính:

1. Trường phái địa hóa học trầm tích và sinh địa hóa: Tập trung vào cơ chế giải phóng As vào nước lỗ rỗng. Đại diện tiêu biểu như Michael Berg et al. (2007) tại đồng bằng sông Bassac - Mekong (Campuchia), Shamsudduha & Uddin (2007) tại đồng bằng Bengal, và Benjamin D. Kocar et al. (2008) tại Campuchia đều khẳng định quá trình vi sinh vật phân giải hợp chất cacbon hữu cơ hoạt động (DOC) trong điều kiện yếm khí là động lực khử oxit/hydroxit sắt (Fe(III) $\rightarrow$ Fe(II)), kéo theo sự giải phóng As(III) hấp phụ vào pha lỏng. Ngược lại, một số nghiên cứu tại Nam Á như Ashraf et al. (2008) lại đề xuất quá trình phong hóa hóa học khoáng vật biotit là cơ chế nền tảng ban đầu trước khi diễn ra các pha khử tiếp theo.
2. Trường phái thủy động lực học và tính tổn thương tầng chứa nước: Tập trung vào con đường vận chuyển As dưới tác động nhân tạo. Holly A. Michael & Clifford I. Voss (2008, 2009) khi mô phỏng bồn trũng Bengal (Ấn Độ - Bangladesh) đã chỉ ra rằng khai thác công nghiệp quy mô lớn phá vỡ thế cân bằng tự nhiên, kéo nước ô nhiễm từ tầng nông xuống tầng sâu thông qua tính bất đồng nhất của trầm tích với hệ số bất đẳng hướng thấm ($K_h/K_v$) lên đến $10^4$. Khan et al. (2016) tại Dhaka nhấn mạnh các dòng chảy ưu thế (preferential flow pathways) có thể làm ô nhiễm tầng sâu (>150m) chỉ trong vòng một thập kỷ.

Tại Việt Nam, các nghiên cứu của Elisabeth Eiche (2008), Đặng Mai et al. (2012), Postma et al. (2007, 2012, 2016) tại Vạn Phúc và Đan Phượng đã xác lập cơ sở sinh địa hóa về sự tích tụ As trong trầm tích Holocen trẻ (<1000 năm). Van Geen et al. (2013) tính toán hệ số trễ của As tại bãi giếng Vạn Phúc đạt $R = 16 - 20$ dưới tác động khai thác. Nghiên cứu thực nghiệm của Trần Vũ Long & Phạm Quý Nhân (2019) tại bãi giếng Nam Dư (Hà Nội) cho thấy dưới động thái khai thác công nghiệp cực mạnh gần sông Hồng, vận tốc thấm thẳng đứng đạt tới ~19 m/năm, khiến hệ số trễ thực tế bị suy giảm mạnh xuống chỉ còn $R = 2 - 20$.

Luận án của Trần Vũ Long định vị chính xác khoảng trống khoa học bằng việc chọn vùng Thạch Thất - Đan Phượng: nơi có nồng độ As trong $q_h$ cao bậc nhất ĐBBB nhưng hoàn toàn duy trì chế độ thủy động lực tự nhiên (không có bãi giếng công nghiệp công suất lớn). Luận án đặt tiền đề so sánh đối sánh trực tiếp với trường hợp khai thác áp lực cao tại Nam Dư (Việt Nam) và các nghiên cứu của Stollenwerk et al. (2007), Radloff et al. (2011) tại bồn trũng Bengal (Bangladesh), qua đó làm sáng tỏ vai trò cốt lõi của vận tốc thấm đối với trạng thái cân bằng hấp phụ và độ lớn của hệ số trễ $R$.

Đó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 Thuyết vận chuyển chất tan phản ứng (Reactive Solute Transport Theory) của Bear (1979) và Appelo & Postma (2005) khi áp dụng cho á kim Asen trong môi trường xáo trộn tự nhiên liên tầng:

• Xác lập mối quan hệ phi tuyến giữa vận tốc thấm thực tế ($v$) và hiệu quả trễ hấp phụ ($R$): Khác với các giả định truyền thống xem $R$ là một đại lượng hằng số địa hóa học tĩnh, luận án chứng minh $R$ là một biến số động phụ thuộc chặt chẽ vào thời gian lưu (residence time) của dòng thấm. Khi vận tốc dòng chảy chậm trong điều kiện tự nhiên ($v$ nhỏ), thời gian tiếp xúc giữa dung dịch chứa As(III) và pha khoáng vật (goethite, biotit, khoáng sét) đủ dài để phản ứng hấp phụ bề mặt đạt trạng thái cân bằng nhiệt động học hoàn toàn, đẩy hệ số trễ lên mức cực đại $R = 69 - 162$.
• Bổ sung luận điểm về sự phân dị hóa trị và tính thuận nghịch hấp phụ: Xác định As tồn tại tuyệt đối ở dạng $As(III)$ ($H_3AsO_3^0$) chiếm 90–100% trong môi trường khử mạnh ($Eh < 0, pH = 6 - 8$). Luận án chứng minh rằng mặc dù $As(V)$ hấp phụ mạnh hơn $As(III)$ gấp 6–8 lần trên hydroxit sắt, nhưng phản ứng hấp phụ của $As(III)$ trên trầm tích $q_p$ là phản ứng thuận nghịch được kiểm soát bởi dung lượng trao đổi bề mặt và nồng độ các ion cạnh tranh ($PO_4^{3-}, HCO_3^-$).

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

Khung phân tích của luận án tích hợp đa ngành gồm 4 trụ cột lý thuyết: Thủy văn địa chất cấu trúc, Địa hóa đồng vị môi trường, Thí nghiệm cân bằng pha rắn - lỏng và Mô hình hóa số 3 chiều dòng chảy - dịch chuyển chất.

Điều kiện biên và giới hạn quy chuẩn: Khung mô hình áp dụng cho hệ thống hai tầng chứa nước cát - sỏi bở rời có quan hệ thủy lực trực tiếp hoặc gián tiếp qua thấu kính sét - bột thấm nước yếu; môi trường địa hóa khử sắt - metan; bỏ qua quá trình oxy hóa $As(III) \rightarrow As(V)$ do môi trường khử sâu được duy trì ổn định.

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

Thiết kế nghiên cứu

Luận án áp dụng quan điểm triết học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism) thông qua phương pháp tiếp cận đa mức (Multi-level empirical design): từ quy mô phân tử (cơ chế trao đổi bề mặt ion), quy mô lõi khoan (thí nghiệm độ hạt, chiết mẫu, định tuổi OSL), quy mô điểm quan sát (slugtest, quan trắc mực nước tự động) đến quy mô khu vực toàn bồn trũng (mô hình số 3D MODFLOW/MT3D-USGS).

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu thực nghiệm tuân thủ các tiêu chuẩn quốc tế nghiêm ngặt thuộc dự án VietAS:

• Khảo sát địa vật lý: Đo 21 điểm trường chuyển (TEM) và 4 tuyến mặt cắt ảnh điện 2D với 320 điểm đo sâu điện nhằm xác định chính xác cấu trúc ranh giới địa tầng và phát hiện các "cửa sổ địa chất thủy văn" nơi lớp sét ngăn cách giữa $q_h$ và $q_p$ bị vắng mặt.
• Khoan và quan trắc động thái: Thi công 64 lỗ khoan địa chất - ĐCTV; đo địa vật lý trong 10 lỗ khoan; thực hiện 12 thí nghiệm bơm hút nước và slugtest để xác định hệ số thấm ($K$). Quan trắc liên tục mực nước tại 26 lỗ khoan với hơn 300 lượt đo trong giai đoạn 2009–2012.
• Phân tích hóa học và đồng vị:
  • Lấy 121 mẫu nước phân tích chỉ tiêu lý hóa tại hiện trường ($pH, Eh, DO, EC$) và trong phòng thí nghiệm (tổng $As, As(III), Fe(II), NH_4^+, CH_4, DOC, DIC$, các cation và anion chính).
  • Phân tích 226 mẫu đồng vị bền ($\delta^{18}\text{O}, \delta\text{D}$), 35 mẫu đồng vị phóng xạ ($^{14}\text{C}$), và 08 mẫu đồng vị khí trơ ($^3\text{He}/^4\text{He}$) tại Viện Nghiên cứu Địa chất Đan Mạch và Greenland (GEUS).
  • Thu thập 26 mẫu trầm tích nguyên dạng để định tuổi quang tích huỳnh quang (OSL).

Data và phân tích

Số liệu được tích hợp và xử lý trên nền tảng phần mềm chuyên dụng AquaVEO GMS (Groundwater Modeling System):

• Mô hình dòng chảy (MODFLOW): Không gian mô phỏng được phân chia thành lưới 3 chiều với các mắt lưới được làm dày tại bãi giếng; mô hình được hiệu chỉnh dựa trên tương quan mực nước quan trắc và mực nước tính toán tại các lỗ khoan chuẩn ($Q56, VC_qp$).
• Mô hình lan truyền chất tan (MT3D-USGS): Tích hợp phương trình đối lưu - phân tán kết hợp hấp phụ phi tuyến và tuyến tính: $$\frac{\partial c}{\partial t} = \nabla \cdot (D \nabla c) - v \nabla c - \frac{\rho_b}{\theta} \frac{\partial q}{\partial t}$$ Trong đó $\rho_b$ là khối lượng thể tích khô của đất đá, $\theta$ là độ lỗ rỗng hiệu dụng, và $q$ tuân theo phương trình hấp phụ đẳng nhiệt Langmuir: $$q = \frac{s_{tot} c}{K_L + c}$$ Số Peclet ($Pe = \frac{v_x d}{D^*}$) được tính toán để phân định miền ưu thế giữa khuếch tán phân tử và phân tán cơ học.

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

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

"As tồn tại trong NDĐ dưới dạng As hoá trị III với hàm lượng rất cao thay đổi lớn tùy theo TCN và theo vị trí; từ 77µg/L (Phụng Thượng) đến 450µg/L (Đan Phượng) trong TCN qh và từ 55µg/L (Phụng Thượng) đến 170µg/L (Đan Phượng) trong TCN qp." (Trích luận án, tr. 8)

1. Dạng tồn tại và nồng độ bất thường: Asen trong nước dưới đất vùng nghiên cứu tồn tại chủ yếu ở dạng khử $As(III)$ chiếm trên 90% (nhiều mẫu đạt 98–100% tổng As hòa tan). Nồng độ As tổng số trong tầng $q_h$ đạt đỉnh tại Đan Phượng (lên tới 450–550 $\mu\text{g/L}$), vượt tiêu chuẩn WHO (10 $\mu\text{g/L}$) tới 45–55 lần. Ở tầng $q_p$, nồng độ As thấp hơn đáng kể nhưng vẫn dao động từ 55 đến 170 $\mu\text{g/L}$.
2. Mối quan hệ nhân quả nghịch đảo giữa tuổi trầm tích và nồng độ As: Dữ liệu định tuổi OSL và đồng vị cho thấy sự tương quan rõ nét: trầm tích $q_h$ tại Đan Phượng sát sông Hồng rất trẻ (~460 năm), chứa hàm lượng cacbon hữu cơ hoạt động cao, tạo môi trường khử cực mạnh (tương quan đồng biến chặt chẽ giữa $As(III)$ với $Fe(II), NH_4^+, CH_4$). Ngược lại, tại vùng rìa Phú Kim, trầm tích cổ hơn rất nhiều (~5.900 năm), vật chất hữu cơ đã thoái hóa gần như hoàn toàn, khiến nồng độ $As(III)$ trong NDĐ giảm rõ rệt.

"As(V) có khả năng hấp phụ mạnh hơn rất nhiều so với As(III) từ 6 - 8 lần tuy nhiên lại không có tính thuận nghịch... số liệu phân tích trên các mẫu NDĐ được lấy từ các lỗ khoan tại các điểm nghiên cứu cho thấy hàm lượng As(III) chiếm >90% (có nhiều mẫu đến 98 - 100%) của tổng hàm lượng As hòa tan trong NDĐ." (Trích luận án, tr. 22, 35-36)

3. Xác lập giá trị hệ số trễ cao trong điều kiện tự nhiên:

"Sự dịch chuyển của As từ TCN qh vào TCN qp chủ yếu tuân theo cơ chế thuỷ động lực. Vận tốc di chuyển của NDĐ chậm, phản ứng hấp phụ trên trầm tích có điều kiện thuận lợi đạt tới trạng thái cân bằng và triệt để thì hàm lượng As gia tăng trong TCN qp theo thời gian là rất thấp tương ứng với hệ số trễ lớn R = 69 - 162." (Trích luận án, tr. 8)

Kết quả mô hình số MT3D-USGS chứng minh rằng tại Thạch Thất - Đan Phượng, do vận tốc thấm tự nhiên nhỏ, hệ số trễ đạt mức $R = 69 - 162$. Điều này giải thích tại sao dù tầng $q_h$ có nồng độ As cực cao và tiếp xúc thủy lực trực tiếp với tầng $q_p$ qua các cửa sổ ĐCTV, tầng $q_p$ vẫn duy trì được nồng độ As tương đối thấp trong suốt hàng trăm năm qua.

4. Sự tương phản cơ chế giữa điều kiện tự nhiên và khai thác công nghiệp: Khi so sánh với bãi giếng Nam Dư ($R = 2 - 20$), sự khác biệt lên đến một bậc độ lớn về hệ số trễ khẳng định rằng chế độ thủy động lực bị cưỡng bức bởi bơm hút lưu lượng lớn sẽ làm giảm mạnh thời gian tiếp xúc của dung dịch với bề mặt hạt, phá vỡ trạng thái cân bằng Langmuir và làm suy giảm khả năng tự bảo vệ của tầng $q_p$.

Chỉ số so sánh	Vùng Thạch Thất - Đan Phượng (Luận án)	Bãi giếng Nam Dư (Trần Vũ Long & Phạm Quý Nhân, 2019)	Bồn trũng Bengal, Bangladesh (Radloff et al., 2011; Van Geen et al., 2013)
Chế độ Thủy động lực	Dòng chảy tự nhiên, chưa bị xáo trộn	Khai thác công nghiệp quy mô lớn sát sông	Khai thác phân tán kết hợp thủy lợi nông nghiệp
Vận tốc thấm thẳng đứng	Rất nhỏ ($v \approx 0.05 - 0.2\text{ m/năm}$)	Cực lớn ($v \approx 19\text{ m/năm}$)	Trung bình đến cao
Hệ số trễ hấp phụ ($R$)	$R = 69 - 162$ (Cân bằng triệt để)	$R = 2 - 20$ (Mất cân bằng động)	$R = 16 - 35$
Nồng độ As(III) trong $q_h$	77 – 450 $\mu\text{g/L}$ (cực đại 550 $\mu\text{g/L}$)	100 – 400 $\mu\text{g/L}$	100 – 700 $\mu\text{g/L}$
Khả năng tự bảo vệ của $q_p$	Rất cao trong điều kiện tự nhiên	Bị phá vỡ, As xâm nhập nhanh	Bị suy thoái theo thời gian

Implications đa chiều

• Về lý thuyết: Khẳng định tính tất yếu của việc ghép nối động học hấp phụ phi tuyến với phương trình thủy động lực học 3D trong việc dự báo lan truyền ô nhiễm á kim trong môi trường địa chất Đệ Tứ.
• Về phương pháp luận: Cung cấp quy trình tích hợp chuẩn xác giữa Địa vật lý (TEM/Điện 2D), Địa hóa đồng vị ($^{14}\text{C}, ^3\text{He}/^4\text{He}, \delta^{18}\text{O}, \delta\text{D}$), Định tuổi trầm tích (OSL) và Mô hình số (GMS/MODFLOW/MT3D-USGS) có thể chuyển giao cho các bồn trũng châu thổ khác (như Châu thổ sông Cửu Long, đồng bằng sông Hằng).
• Về quản lý và chính sách: Đưa ra cảnh báo khoa học xác đáng: Tuyệt đối không quy hoạch các bãi giếng khai thác tập trung công suất lớn tại các vùng có "cửa sổ địa chất thủy văn" tiếp xúc trực tiếp giữa $q_h$ và $q_p$ (như Đan Phượng), vì sẽ làm giảm hệ số trễ $R$, kéo sụt As từ tầng nông xuống phá hủy nguồn cấp nước an toàn của tầng sâu.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn tồn tại một số giới hạn khoa học:

1. Giả định đẳng nhiệt hấp phụ Langmuir/Freundlich đơn thành phần: Mô hình dịch chuyển chưa tích hợp hoàn toàn mô hình phức hóa bề mặt đa cấu tử (Surface Complexation Modeling - SCM) để mô phỏng tương tác cạnh tranh vị trí hấp phụ phức tạp giữa $As(III)$ với các anion cạnh tranh mạnh như $PO_4^{3-}, SiO_4^{4-}, HCO_3^-$ và ion $Fe(II)$ hòa tan.
2. Độ phân giải không gian của thấu kính sét: Mặc dù đã đo 4 tuyến địa điện 2D và 21 điểm TEM, tính dị hướng vi mô của các thấu kính bùn sét phân bố rải rác trong bồn trũng bồi tích vẫn chưa được phản ánh tuyệt đối trong mô hình số 3D.
3. Thời gian quan trắc thực địa: Chuỗi số liệu quan trắc động thái mực nước và chất lượng nước tập trung trong giai đoạn 3 năm (2009–2012), cần được tiếp tục duy trì quan trắc dài hạn đa thập kỷ để kiểm chứng độ tin cậy của các kịch bản dự báo 50–100 năm.

Hướng nghiên cứu tiếp theo bao gồm:

• Ứng dụng mô hình thủy địa hóa tương tác dòng chảy nhiệt động lực học 3D (như PHT3D hoặc TOUGHREACT) kết hợp mô hình phức hóa bề mặt (SCM).
• Nghiên cứu cơ chế giải hấp phụ As khi có sự xâm nhập của nguồn nước bổ cập nhân tạo giàu oxy hoặc nitrat nhằm phục hồi tầng chứa nước.
• Mở rộng nghiên cứu biến đổi vi sinh vật học phân tử (metagenomics) liên quan đến các chủng vi khuẩn khử sắt/khử asen trong trầm tích $q_h$.

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

• Học thuật: Đóng góp dữ liệu thực nghiệm chuẩn xác cho mạng lưới nghiên cứu asen quốc tế; là tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh, học viên cao học ngành Địa chất thủy văn, Kỹ thuật địa chất và Khoa học môi trường.
• Chính sách và quy hoạch tài nguyên nước: Cung cấp cơ sở khoa học trực tiếp cho Bộ Tài nguyên và Môi trường, UBND Thành phố Hà Nội và Liên đoàn Quy hoạch và Điều tra tài nguyên nước miền Bắc trong việc phân vùng cấm/hạn chế khai thác nước dưới đất, bảo vệ tầng chứa nước Pleistocen theo Luật Tài nguyên nước.
• Xã hội: Giúp các nhà máy nước địa phương và người dân nông thôn hiểu rõ nguy cơ ô nhiễm, định hướng chuyển đổi nguồn cấp nước hoặc áp dụng các công nghệ xử lý nước giếng khoan thích hợp, giảm thiểu bệnh tật do phơi nhiễm asen mãn tính (như tổn thương da, ung thư).

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

• Nghiên cứu sinh & Nhà khoa học Địa chất: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa địa chất thủy văn hiện trường, địa hóa đồng vị và mô hình số 3D.
• Cơ quan quản lý Nhà nước (Bộ TN&MT, Sở TN&MT): Có cơ sở định lượng để ban hành quy chuẩn khai thác an toàn và khoanh định hành lang bảo vệ nguồn nước ngầm.
• Doanh nghiệp cấp nước & R&D công nghệ môi trường: Nắm bắt chính xác nồng độ, hóa trị As(III) và các ion đi kèm ($Fe, NH_4^+$) để thiết kế dây chuyền xử lý tối ưu (kết hợp oxy hóa - keo tụ - lọc cát).

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?

Trả lời: Đóng góp lý thuyết độc đáo nhất là việc chứng minh và định lượng hóa sự phụ thuộc của hệ số trễ hấp phụ ($R$) vào vận tốc dòng thấm thủy động lực trong môi trường tự nhiên liên tầng. Luận án đã mở rộng Thuyết vận chuyển chất tan phản ứng (Reactive Transport Theory) khi chứng minh rằng trong điều kiện dòng thấm tự nhiên vận tốc chậm ($v \approx 0.05 - 0.2\text{ m/năm}$), phản ứng hấp phụ của $As(III)$ trên pha rắn trầm tích đạt trạng thái cân bằng nhiệt động học hoàn toàn, tạo ra giá trị $R = 69 - 162$. Điều này giải thích cơ chế "tự bảo vệ tự nhiên" của tầng chứa nước Pleistocen ($q_p$) trước nguồn ô nhiễm khổng lồ từ tầng Holocen ($q_h$).

2. Điểm cải tiến phương pháp luận của luận án so với các nghiên cứu trước đây là gì?

Trả lời: So với các nghiên cứu của Postma et al. (2010, 2012) vốn chỉ dùng mô hình địa hóa 1D hoặc Larsen et al. (2008) chỉ dừng lại ở mô hình dòng chảy đơn thuần, luận án đã:

1. Thiết lập mô hình số không gian 3 chiều hoàn chỉnh (MODFLOW kết hợp MT3D-USGS) trên phần mềm AquaVEO GMS mô phỏng đồng thời dòng thấm 3D và sự lan truyền Asen có xét đến hấp phụ phi tuyến.
2. Tích hợp chuỗi phương pháp đa dạng chưa từng có tại khu vực: 21 điểm TEM, 320 điểm ảnh điện 2D, 64 lỗ khoan, 26 mẫu OSL, 269 mẫu đồng vị ($\delta^{18}\text{O}, \delta\text{D}, ^{14}\text{C}, ^3\text{He}/^4\text{He}$) để hiệu chỉnh mô hình số.

3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm là gì?

Trả lời: Phát hiện bất ngờ nhất là nồng độ As trong nước dưới đất tại điểm nghiên cứu Đan Phượng đạt mức cực cao (lên tới 450–550 $\mu\text{g/L}$) nhưng gần như 100% tồn tại dưới dạng $As(III)$, đi kèm nồng độ cực cao của $Fe(II), NH_4^+, CH_4$. Mặc dù tầng $q_h$ và tầng $q_p$ tại đây tiếp xúc trực tiếp qua cửa sổ ĐCTV không có tầng sét ngăn cách, nồng độ As trong tầng $q_p$ vẫn thấp hơn 3–4 lần so với $q_h$ nhờ hệ số trễ $R$ cực lớn ($69 - 162$), chứng minh pha khoáng vật của tầng $q_p$ sở hữu dung lượng hấp phụ Asen rất lớn chưa bị bão hòa.

4. Quy trình nghiên cứu có cung cấp đầy đủ giao thức để tái lặp (Replication Protocol) không?

Trả lời: Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình thực nghiệm: từ phương pháp bố trí mạng lưới giếng quan trắc, kỹ thuật lấy mẫu nước yếm khí bảo toàn hóa trị $As(III)/As(V)$, phương pháp phân tích quang phổ khối plasma (ICP-MS), quy trình đo đạc địa vật lý TEM/2D điện trở, giao thức chuẩn bị mẫu định tuổi OSL, cho đến bộ thông số thủy lực đầu vào ($K_x, K_y, K_z, S_s, \alpha_L, \alpha_T, K_d$) sử dụng trong tệp dữ liệu MODFLOW/MT3D-USGS, cho phép các nhà nghiên cứu độc lập tái lập hoàn toàn mô hình.

5. Luận án phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?

Trả lời: Chương trình nghiên cứu dài hạn được định hướng gồm 3 trọng tâm:

1. Xây dựng mạng lưới quan trắc tự động liên tục (telemetry monitoring) về áp lực nước và nồng độ Asen tại các cửa sổ ĐCTV dọc sông Hồng.
2. Phát triển mô hình thủy địa hóa nhiệt động lực học 3D đa cấu tử (SCM-PHT3D) đánh giá tác động cạnh tranh hấp phụ của ion $PO_4^{3-}$ từ phân bón nông nghiệp ngấm xuống.
3. Thử nghiệm công nghệ nạp nước nhân tạo bổ cập (Managed Aquifer Recharge - MAR) có kiểm soát để tạo màng chắn địa hóa ngăn chặn Asen dịch chuyển vào các bãi giếng cấp nước đô thị.

Kết luận

1. Luận án đã xác định toàn diện hiện trạng ô nhiễm Asen trong nước dưới đất tại vùng Thạch Thất - Đan Phượng: nồng độ As trong tầng $q_h$ dao động từ 77 đến 450 $\mu\text{g/L}$ và tầng $q_p$ từ 55 đến 170 $\mu\text{g/L}$, tồn tại chủ yếu ở dạng $As(III)$ (>90%) trong môi trường khử mạnh.
2. Làm sáng tỏ mối quan hệ giữa tuổi trầm tích và ô nhiễm As: Trầm tích Holocen trẻ (~460 năm) giàu chất hữu cơ hoạt tính là nguồn sinh địa hóa giải phóng As mạnh nhất; trầm tích càng cổ (vùng rìa Phú Kim ~5.900 năm), tiềm năng giải phóng As càng suy giảm.
3. Chứng minh cơ chế dịch chuyển Asen từ $q_h$ vào $q_p$ tuân theo cơ chế thủy động lực đối lưu - phân tán nhưng bị kiểm soát mạnh mẽ bởi quá trình hấp phụ trên pha rắn trầm tích.
4. Định lượng thành công hệ số trễ $R = 69 - 162$ trong điều kiện thủy động lực tự nhiên, khẳng định vai trò của vận tốc thấm chậm trong việc duy trì trạng thái cân bằng hấp phụ triệt để, bảo vệ tầng chứa nước Pleistocen.
5. Cảnh báo rủi ro suy giảm hệ số trễ ($R \rightarrow 2 - 20$) nếu khai thác nước ngầm quy mô công nghiệp tại các vùng cửa sổ ĐCTV sát sông Hồng, cung cấp luận cứ khoa học then chốt cho quy hoạch và bảo vệ an ninh tài nguyên nước dưới đất vùng Thủ đô Hà Nội.

Trích đoạn nội dung tài liệu

CHƯƠNG 1 - TỔNG QUAN VỀ NGHIÊN CỨU ASEN TRONG NƯỚC DƯỚI ĐẤT Việc gia tăng hàm lượng và dịch chuyển của As trong NDĐ đang là vấn đề được quan tâm liên quan tới chất lượng nước và sức khoẻ người dân tại một số đồng bằng châu thổ lớn trên thế giới như Tây Bengal (Ấn Độ), Bangladesh, Đài Loan, một số vùng ở Achentina, Canađa, Mỹ và ĐBBB - Việt Nam. Rất nhiều nghiên cứu đã được thực hiện tại các khu vực này [47], [54], [60], [67], [77], [94], [95]. Tuy nhiên các nghiên cứu này mới chỉ tập trung vào các vấn đề như phân bố của As trong các TCN, các quá trình địa hoá và thuỷ địa hoá liên quan đến giải phóng As vào NDĐ và mô hình số thuỷ địa hoá mô phỏng quá trình giải phóng As từ trầm tích vào NDĐ. Các công trình nghiên cứu trên thế giới Trên cơ sở tổng hợp và phân tích các công trình khoa học nghiên cứu As đã công bố cho thấy, các công trình nghiên cứu này đều sử dụng đa dạng các phương pháp nghiên cứu khác nhau để giải quyết vấn đề hàm lượng cao và sự dịch chuyển của As trong NDĐ.

Đồng thời đề xuất các giải pháp ngăn chặn, giảm thiểu nguy cơ, hạn chế ô nhiễm As đối với các TCN là đối tượng khai thác chính. Như vậy, có thể phân ra nhóm các kết quả công bố như sau: 1. Nhóm đánh giá hiện trạng và xác định cơ chế giải phóng As vào NDĐ Việc gia tăng hàm lượng As trong NDĐ là vấn đề nghiêm trọng đối với chất lượng nước tại một số đồng bằng châu thổ trên thế giới như Tây Bengal (Ấn Độ), Băngladesh, Đài Loan, một số vùng ở Achentina, Canađa, Mỹ. Kết quả nghiên cứu cho thấy hàm lượng As trong NDĐ nhiều nơi đã vượt quá giới hạn cho phép của WHO (10µg/L) như ở Manikaganj, Harirampar, Faridpur, Gopalganj (Bangladesh) có 14/19 mẫu vượt quá tiêu chuẩn cho phép của Bangladesh (50µg/L), riêng vùng Harirampar cả 4 mẫu đều trên 100µg/L [47], [90].

Hàm lượng cao của As có thể tìm thấy lên tới 1. Còn ở phía tây nam Đài Loan hàm lượng As trung bình từ 147 - 671µg/L và người dân sử dụng nước ở đây đã bị bệnh đen chân 11 (blackfoot) [83], [93]. Hai khu vực của Achentina là San Antonio delos Codres và Taco Poro, mỗi nơi có hàm lượng As khoảng 200µg/L [86]. Ô nhiễm As trong NDĐ ở phía đông sông Hoogky, một nhánh của sông Hằng phía tây Bengal đã được báo cáo từ đầu năm 1978.

Nhóm bệnh nhân đầu tiên được phát hiện vào 7/1983 [67]. Kể từ đó phạm vi ảnh hưởng và số bệnh nhân mới ngày càng tăng. Khu vực ảnh hưởng rộng 3.400km2, xấp xỉ 30 triệu dân, số người sử dụng nước nhiễm độc As lên tới 1 triệu người, trong đó hơn 200.000 người đã được xác nhận là có triệu chứng nhiễm độc As. Đây là vụ nhiễm độc As lớn nhất trong lịch sử và ô nhiễm As trong NDĐ đã trở thành một vấn đề quốc gia [67].

Một cuộc khảo sát của Michael Berg và cộng sự (2007) [34] thực hiện ở tỉnh Kandal nằm ở ồng bằng châu thổ giữa hai con sông Bassac và Mekong, Campuchia. Kết quả thu được cho thấy: NDĐ ở khu vực có hàm lượng As nằm trong khoảng từ 1 - 1. Thêm vào đó, hàm lượng As cao tương ứng với điều kiện khử trong TCN qh, ở đó hàm lượng trung bình của NH4+ (6,0mg/L), DOC (3,9mg/L), Fe(II) (2,8mg/L), Mn(II) (0,62mg/L) đều cao. Tác giả cũng tìm thấy mối tương quan thuận giữa As và NH4+, DOC được xem như là một bằng chứng cho thấy NDĐ có môi trường khử và diễn ra quá trình khử hòa tan các khoáng vật oxit có chứa As.

Shamsudduha và Uddin (2007) [84] nghiên cứu ảnh hưởng của sự dịch chuyển đường bờ biển trong Đệ Tứ và điều kiện ĐCTV đến sự phân bố As trong các TCN ở đồng bằng Bengal. Quá trình biển tiến trong Holocen đã làm hình thành nên các tầng trầm tích giàu vật liệu hữu cơ với điều kiện môi trường khử là nguồn cung cấp chủ yếu As cho NDĐ. Cao độ tuyệt đối và dốc địa hình cũng có ảnh hưởng lớn tới đặc điểm phân bố As trong NDĐ bởi các kết quả cho thấy hàm lượng As cao được tìm thấy ở những vị trí có địa hình thấp và bằng phẳng. Một nghiên cứu của Benjamin D.

Kocar và nnk (2008) [55] được tiến hành ở vùng đồng bằng châu thổ sông Mê-Kong gần thủ đô Phnom Pênh, Căm-pu-chia với mục đích tìm hiểu các quá trình thủy văn và sinh địa hóa ảnh hưởng đến hàm lượng 12 hòa tan cũng như quá trình dịch chuyển của As trong TCN qh cho thấy: Hàm lượng As hòa tan thay đổi theo không gian, khoảng trên 1.300µg/L trong NDĐ đến 600µg/L ở nước lỗ rỗng chiết từ lớp sét sát bề mặt. NDĐ có hàm lượng As cao có môi trường khử với hàm lượng oxy hoà tan không đáng kể và hàm lượng Fe(II), NH4+, Carbon hữu cơ hòa tan rất cao. Ở môi trường gần bề mặt nhưng nằm dưới các vùng ngập nước lâu dài, các chỉ thị này càng dễ nhận thấy. Các tác giả cho rằng chính các Carbon hữu cơ không bền, lắng đọng cùng các sắt oxy hydroxit mang As, dưới điều kiện khử mạnh đã giải phóng As từ dạng rắn thành dạng hòa tan.

Tác giả cũng đưa ra giả thuyết As ở tầng gần bề mặt, nằm dưới các vùng đất ngập nước được giải phóng như vậy được coi là nguồn gây ô nhiễm As chính đối với các TCN phía dưới trong thời gian dài (100 đến 1.000 năm) ở vùng Nam và Đông Nam Á. Phân bố As trong NDĐ tại Bangladesh [54] Ashraf và nnk (2008) [81] đã nghiên cứu về sự giải phóng As từ khoáng vật biotit vào TCN qh ở Băngladesh. Kết quả chỉ ra rằng, ở độ sâu 18 - 29m trong TCN qh, nơi các lỗ khoan nông được người dân sử dụng để khai thác nước phục vụ sinh hoạt và tưới tiêu thì As tồn tại chủ yếu trong pha biotit. Mặc dù, quá trình khử hòa tan các sắt oxy hydroxit được cho là nguyên nhân chính gây nên hàm lượng As cao 13 trong NDĐ, tuy nhiên các tác giả đã kết luận rằng quá trình phong hóa hóa học khoáng vật biotit là cơ chế nền tảng ban đầu và các điều kiện khử tiếp theo đóng góp vào việc giải phóng As làm tăng cao hàm lượng As trong NDĐ.

Yilong Zhang và nnk (2013) [98] nghiên cứu mối quan hệ giữa hàm lượng As trong NDĐ và gradient thủy lực dọc theo hướng dòng chảy NDĐ TCN nằm nông ở đồng bằng Hetao, Bắc Trung Quốc. Điều kiện thủy động lực được chia theo 4 nhóm từ Bắc tới Nam. Phương trình hàm mũ được xây dựng thể hiện cho mối tương quan giữa hàm lượng As và gradient thủy lực ở vùng được bổ cập từ nước mưa. Các kết quả cho thấy gradient thủy lực thấp trong điều kiện tự nhiên là điều kiện quan trọng làm tăng hàm lượng As trong NDĐ.

Qi Guo và nnk (2014) [45] nghiên cứu về sự khác biệt giữa điều kiện thủy địa hóa giữa NDĐ có hàm lượng As cao và NDĐ có hàm lượng As thấp trong TCN qh đồng bằng Yinchuan, Trung Quốc. Mục đích của công trình này là đánh giá mối liên quan giữa phân bố As trong NDĐ và các quá trình địa hóa làm di chuyển As ở 3 khu vực: Khu vực phân bố của nón phóng vật, khu vực canh tác khô hạn, và khu vực trồng lúa. Kết quả cho thấy ở khu vực phân bố nón phóng vật, hàm lượng As trong NDĐ thấp. Tuy nhiên, ở 2 khu vực còn lại hàm lượng As lại cao và liên quan chặt chẽ với điều kiện môi trường khử.

Ở khu vực canh tác lúa, nước tưới tiêu đóng vai trò là nguồn bổ cập cho NDĐ dẫn đến thế oxy hóa khử của NDĐ ở vùng này cao hơn so với vùng canh tác khô hạn. Các đồng vị bền 18O và 2H cho thấy bốc hơi do nước tưới đóng vai trò không đáng kể trong việc tăng cao hàm lượng As. Mối tương quan giữa hàm lượng As và sắt trong NDĐ cho thấy As có nguồn gốc chủ yếu từ quá trình khử hòa tan các sắt oxy hydroxit. Bên cạnh đó, phốt pho hòa tan có thể đóng vai trò là chất cạnh tranh với As.

Nhận xét chung: Như vậy các nghiên cứu trên chủ yếu tập trung vào nghiên cứu phân bố của ô nhiễm As trong NDĐ của các TCN, làm sáng tỏ cơ chế giải phóng As vào NDĐ của TCN bên trên. Các kết quả nghiên cứu đều chỉ ra hàm lượng As cao chủ yếu được phát hiện tại TCN trầm tích aluvi trẻ bên trên và cơ chế giải phóng 14 chính là khử hoà tan sắt oxy hydroxit có hấp phụ As vào NDĐ dưới ảnh hưởng của quá trình phân giải vật chất hữu cơ hoạt động trầm tích đồng thời với thành tạo TCN. Nhóm nghiên cứu con đường dịch chuyển As và khả năng tổn thương của TCN bên dưới trước ô nhiễm As Holly A. Michael và Cliffor I.

Voss (2008) [61] đã nghiên cứu đánh giá về tính bền vững của NDĐ trong TCN nằm sâu cho mục đích cung cấp nước với hàm lượng As thấp tại đồng bằng Bengal. Trong nghiên cứu này, TCN nằm sâu hơn 150m đã được đánh giá chi tiết bằng việc sử dụng mô hình số cho toàn đồng bằng Bengal. Nghiên cứu cũng chỉ ra việc khai thác nước tại TCN bên dưới với mối quan hệ thuỷ lực với TCN nằm bên trên sẽ kéo theo quá trình di chuyển As hoà tan trong NDĐ từ TCN bên trên và phá huỷ TCN này. Kết quả của mô hình số cho thấy quá trình khai thác NDĐ quy mô nhỏ ở hộ gia đình ở TCN bên dưới không ảnh hưởng nhiều đến hệ thống dòng chảy ngầm tại đây.

Đồng thời với quá trình khai thác tại TCN bên trên phục vụ cho nông nghiệp cũng góp phần bảo vệ TCN bên dưới khỏi các tác động của quá trình dịch chuyển As. Các phân tích từ mô hình số còn chỉ ra việc khai thác được kiểm soát chặt chẽ có thể giúp cung cấp nguồn nước an toàn cho >90% nhu cầu trong hơn 1. Michael và Cliffor I. Voss (2009) [62] đã nghiên cứu đánh giá về các đặc tính dòng chảy ngầm trên quy mô khu vực tại đồng bằng Bengal, Ấn Độ và Bangladesh.

Nghiên cứu này nhằm đánh giá các chiến lược quản lý trong dài hạn đối với cung cấp nước sinh hoạt an toàn không có As cho đồng bằng Bengal. Tại đồng bằng này, các trầm tích có mức độ bất đồng nhất cao với các tập thấu kính thấm nước yếu cục bộ chia tách các TCN.