Tổng quan về luận án
Nghiên cứu xâm nhập mặn trong các cấu trúc địa chất thủy văn ven biển là một trong những thách thức khoa học và thực tiễn cấp bách nhất hiện nay dưới tác động kép của biến đổi khí hậu, nước biển dâng và áp lực nhân sinh từ khai thác nước dưới đất quá mức. Luận án tiến sĩ kỹ thuật địa chất với đề tài "Nghiên cứu xâm nhập mặn nước dưới đất trầm tích Đệ tứ vùng Nam Định" do nghiên cứu sinh Hoàng Văn Hoan thực hiện dưới sự hướng dẫn khoa học của PGS. Phạm Quý Nhân và PGS. Flemming Larsen (Trường Đại học Mỏ - Địa chất Hà Nội, 2014) là công trình nghiên cứu tiên phong, giải quyết toàn diện bài toán nguồn gốc hình thành, cấu trúc phân bố và cơ chế động lực học xâm nhập mặn của thấu kính nước nhạt nằm sâu trong trầm tích Kainozoi vùng duyên hải đồng bằng Bắc Bộ.
Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác lập xuất phát từ sự thiếu hụt các bằng chứng thực nghiệm định lượng về nguồn bổ cập cổ địa chất và vai trò của các lớp trầm tích thấm nước yếu nguồn gốc biển đối với sự suy thoái thấu kính nước nhạt. Dù các công trình trước đây của Đỗ Trọng Sự và Nguyễn Kim Ngọc (1985), Đặng Đình Phúc (2000), hay Đoàn Văn Cánh (2004) đã ghi nhận sự hiện diện của thấu kính nước nhạt tầng Pleistocen tại khu vực Kim Sơn - Hải Hậu, song cơ chế hình thành và xâm nhập mặn chủ yếu dừng lại ở giả thiết định tính, chưa làm sáng tỏ động học dịch chuyển vật chất từ nước lỗ rỗng (porewater) của lớp sét biển cũng như tương quan động lực giữa cấu trúc đứt gãy móng trước Kainozoi với bồn trũng sông Hồng.
Để giải quyết triệt để khoảng trống này, luận án tập trung vào các câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:
- Câu hỏi nghiên cứu 1 (RQ1): Thấu kính nước nhạt tầng Pleistocen được hình thành từ nguồn gốc cổ địa chất nào và cơ chế duy trì động lực bổ cập hiện nay diễn ra ra sao?
- Câu hỏi nghiên cứu 2 (RQ2): Hiện trạng phân bố không gian 3D của ranh giới mặn - nhạt biến đổi như thế nào dưới tác động của cấu trúc địa chất và các trường địa vật lý?
- Câu hỏi nghiên cứu 3 (RQ3): Quá trình thẩm thấu, khuếch tán phân tử và phân dị trọng lực từ lớp thấm nước yếu nguồn gốc biển đóng vai trò định lượng bao nhiêu trong tổng lượng muối xâm nhập vào tầng chứa nước?
- Giả thuyết khoa học 1 (H1): Thấu kính nước nhạt Pleistocen không phải là thể nước chôn vùi biệt lập mà là một hệ thống động lực mở, được bổ cập liên tục từ các thành tạo chứa nước khe nứt karst Triat (hệ tầng Đồng Giao) và trầm tích Neogen phía tây - tây bắc thông qua mạng lưới đứt gãy kiến tạo phương Tây Bắc - Đông Nam (TB-ĐN).
- Giả thuyết khoa học 2 (H2): Sự xâm nhập mặn của tầng chứa nước Pleistocen chịu sự chi phối đồng thời của hai cơ chế chính: xâm nhập ngang do chênh lệch áp lực thủy lực sinh ra từ khai thác và xâm nhập thẳng đứng từ lớp sét biển Holocen (hệ tầng Hải Hưng) thông qua quá trình khuếch tán phân tử kết hợp phân dị trọng lực.
Khung lý thuyết của nghiên cứu tích hợp lý thuyết thủy động lực học môi trường lỗ hổng (Bear, 1972; Freeze & Cherry, 1979), nguyên lý cân bằng thủy tĩnh và động lực Ghyben-Herzberg, định luật khuếch tán phân tử Fick, kết hợp mô hình tiến hóa trầm tích - địa hóa đồng vị môi trường. Luận án mang lại đóng góp đột phá khi định lượng hóa quá trình phóng thích muối từ lớp thấm nước yếu với hệ số khuếch tán phân tử đo đạc trực tiếp, đồng thời xây dựng mô hình dự báo lan truyền độ mặn 3D chuẩn xác cho vùng nghiên cứu rộng 2.500 km² với trữ lượng khai thác tiềm năng đạt xấp xỉ 203.000 m³/ngày.
Literature Review và Positioning
Nghiên cứu xâm nhập mặn nước dưới đất trên thế giới đã phát triển qua nhiều giai đoạn với các trường phái tiếp cận đa dạng. Nhóm nghiên cứu tiến hóa cổ địa văn và cấu trúc địa chất, đại diện bởi J. De Vries (1981) tại Hà Lan, đã chứng minh mối liên hệ chặt chẽ giữa các chu kỳ biển tiến - biển thoái thời kỳ Đệ tứ với sự phân bố các thể nước mặn - nhạt ven biển. Tiếp cận sâu hơn về mặt cơ chế thủy địa hóa, H. Groen (2000) tại Đại học Vrije (Hà Lan) đã phân loại 4 cơ chế xâm nhập mặn qua thực nghiệm máng thấm hai lớp, chứng minh rằng trong điều kiện biển tiến nhanh qua trầm tích hạt mịn, quá trình khuếch tán phân tử và đối lưu trọng lực đóng vai trò chủ đạo làm nhiễm mặn tầng chứa nước bên dưới.
Tại Bắc Mỹ và Châu Á, các công trình của George D. Wardlaw và David L. Valentine (2005) tại hồ Salton Sea (California, Mỹ) đã định lượng tốc độ khuếch tán muối từ trầm tích đáy hồ đạt từ 0,422 đến 0,613 g/cm²/năm bằng cách phân tích nước lỗ rỗng ép từ các cột mẫu lõi khoan sâu 35m. Ngược lại, tại vùng châu thổ sông Fraser (Canada), Allen (2006) sử dụng địa vật lý kết hợp thủy địa hóa để mô hình hóa quá trình muối hóa từ cửa sông khuếch tán xuống tầng thấm nước yếu. Bên cạnh đó, các nghiên cứu ứng dụng đồng vị môi trường như J. Meesters (2000) tại Suriname với đồng vị $^{37}\text{Cl}$, Dongmei Han và Claus Kohfahl (2011) tại vịnh Laizhou (Trung Quốc) với tổ hợp đồng vị $^{18}\text{O}$, $^{2}\text{H}$, $^{3}\text{H}$ và $^{14}\text{C}$ đã chứng minh hiện tượng xâm nhập mặn phần lớn bắt nguồn từ nước biển cổ (connate marine water) bị giữ lại trong các thành tạo sét - bột biển Đệ tứ chứ không đơn thuần do nước biển hiện đại thâm nhập.
┌────────────────────────────────────────────────────────┐
│ KHUNG TIẾP CẬN LITERATURE REVIEW │
└──────────────────────────┬─────────────────────────────┘
│
┌───────────────────────────────────┼───────────────────────────────────┐
▼ ▼ ▼
┌──────────────────┐ ┌──────────────────┐ ┌──────────────────┐
│ Cổ địa chất & │ │ Động lực học & │ │ Địa hóa đồng vị │
│ Tiến hóa trầm │ │ Cơ chế khuếch │ │ & Mô hình hóa số │
│ tích │ │ tán / Phân dị │ │ │
├──────────────────┤ ├──────────────────┤ ├──────────────────┤
│ J. De Vries │ │ H. Groen (2000) │ │ J. Meesters │
│ (1981) │ │ Wardlaw & │ │ (2000) │
│ Trần Nghi (2012) │ │ Valentine (2005) │ │ Dongmei Han │
│ Đỗ Trọng Sự & │ │ Paschke & Hoopes │ │ (2011) │
│ Nguyễn Kim Ngọc │ │ (1984) │ │ Phạm Quý Nhân │
│ (1985) │ │ Schincariol & │ │ (2000, 2003) │
│ │ │ Schwartz (1990) │ │ Phatcharasak │
│ │ │ │ │ Arlai (2007) │
└────────┬─────────┘ └────────┬─────────┘ └────────┬─────────┘
│ │ │
└──────────────────────────────────┼───────────────────────────────────┘
│
▼
┌────────────────────────────────────────────────────────┐
│ VỊ TRÍ TIẾP NỐI & ĐỘT PHÁ CỦA LUẬN ÁN │
│ Tích hợp ĐVL sâu (TEM) + Phân tích nước lỗ rỗng │
│ (Porewater) + Thủy văn đồng vị (^14C, ^18O, Triti) │
│ + Mô hình hóa SEAWAT định lượng cho ĐBBB │
└────────────────────────────────────────────────────────┘
Trong nước, các công trình của Đỗ Trọng Sự, Nguyễn Kim Ngọc (1985), Đặng Đình Phúc (1997, 2000), Phạm Quý Nhân (2000, 2003, 2006), Nguyễn Văn Lâm (2006, 2011) và Nguyễn Như Trung (2007) đã đặt nền móng vững chắc về việc khoanh định ranh giới mặn - nhạt tổng thể tại đồng bằng Bắc Bộ. Tuy nhiên, một cuộc tranh luận học thuật kéo dài tồn tại xung quanh nguồn gốc của thấu kính nước nhạt Nam Định: một quan điểm cho rằng đây là thể nước nhạt chôn vùi cổ (relict fresh water) được hình thành trong thời kỳ biển thoái Pleistocen muộn và tồn tại thụ động; quan điểm đối lập cho rằng thấu kính được cấp nước liên tục từ karst Triat Ninh Bình nhưng chưa đưa ra được minh chứng đồng vị và thủy lực thuyết phục.
Luận án của Hoàng Văn Hoan định vị chính xác vào trung tâm điểm giao thoa này, giải quyết triệt để tranh luận bằng việc kết hợp đồng bộ chuỗi dữ liệu đồng vị bền, tuổi nước $^{14}\text{C}$, thí nghiệm ép nước lỗ rỗng màng tế vi và mô hình số SEAWAT mật độ biến đổi.
Đó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 và làm sâu sắc thêm các lý thuyết nền tảng trong lĩnh vực địa chất thủy văn và động lực học chất tan:
- Mở rộng lý thuyết bồn Actezi cổ (Paleo-Artesian Basin Dynamics): Chứng minh thấu kính nước nhạt vùng duyên hải không phải cấu trúc tĩnh khép kín mà là hệ thống dòng thấm động lực xuyên tầng liên tục, chịu sự kiểm soát của trường áp lực truyền dẫn từ khối nâng karst Mesozoi (hệ tầng Đồng Giao, $T_2đg$) sang bồn trũng Kainozoi sông Hồng qua các hệ thống đứt gãy sinh áp.
- Làm sáng tỏ lý thuyết khuếch tán đối lưu mật độ biến đổi (Variable-Density Solute Transport Theory): Bổ sung bằng chứng thực nghiệm vững chắc cho mô hình lý thuyết của Paschke & Hoopes (1984) và Schincariol & Schwartz (1990), khẳng định trong môi trường trầm tích hạt mịn chứa nước mặn biển tiến, quá trình khuếch tán phân tử bị gia tốc đáng kể bởi hiện tượng phân dị trọng lực do chênh lệch tỷ trọng nước mặn - nước nhạt.
- Chuyển dịch nhận thức luận (Paradigm Shift): Chuyển đổi mô hình nhận thức từ việc xem xét ranh giới mặn - nhạt như một mặt tiếp xúc đơn giản sang một đới chuyển tiếp động lực phức tạp 3 chiều, nơi diễn ra tương tác hóa lý liên tục giữa nước lỗ rỗng ép nén và dòng thấm trong tầng chứa nước.
Khung phân tích độc đáo
Khung phân tích của luận án được kiến tạo dựa trên sự liên kết đa ngành chặt chẽ giữa 4 khối lý thuyết: Kiến tạo địa tầng trầm tích, Địa vật lý điện trở suất, Thủy địa hóa - Đồng vị môi trường và Mô phỏng động lực học dòng thấm.
┌────────────────────────────────────────────────────────┐
│ KHUNG PHÂN TÍCH ĐỘC ĐÁO 4 CẤP │
└──────────────────────────┬─────────────────────────────┘
│
┌──────────────────────────────────────────┼──────────────────────────────────────────┐
▼ ▼ ▼
┌──────────────────────────────────┐ ┌──────────────────────────────────┐ ┌──────────────────────────────────┐
│ CẤU TRÚC ĐỊA TẦNG │ │ ĐỊA VẬT LÝ TRƯỜNG BIẾN │ │ ĐỒNG VỊ & THỦY ĐỊA HÓA │
├──────────────────────────────────┤ ├──────────────────────────────────┤ ├──────────────────────────────────┤
│ • Móng Proterozoi & Karst Triat │ │ • Đo sâu trường chuyển (TEM) │ │ • Đồng vị bền δ18O, δ2H │
│ • Trầm tích Kainozoi (N2vb) │ ───► │ • Karota độ dẫn điện lỗ khoan │ ───► │ • Tuổi phóng xạ 14C, 3H, Khí trơ │
│ • Trầm tích Đệ tứ (qp, qh) │ │ • Thiết lập tương quan hồi quy: │ │ • Ép tách nước lỗ rỗng │
│ • Đứt gãy dẫn áp TB-ĐN │ │ TDS - Điện trở suất │ │ (Porewater extraction) │
└──────────────────────────────────┘ └──────────────────────────────────┘ └──────────────────────────────────┘
│
▼
┌──────────────────────────────────┐
│ MÔ HÌNH HÓA ĐỘNG LỰC │
├──────────────────────────────────┤
│ • Mô phỏng dòng chảy MODFLOW │
│ • Vận chuyển chất tan MT3DMS │
│ • Mật độ biến đổi SEAWAT │
│ • Dự báo kịch bản khai thác 3D │
└──────────────────────────────────┘
- Điều kiện biên và giới hạn ứng dụng (Boundary conditions): Khung phân tích áp dụng đặc thù cho các đồng bằng châu thổ ven biển có cấu trúc địa chất phân tầng dạng hai lớp hoặc đa lớp, tồn tại trầm tích biển hạt mịn tuổi Holocen phủ trực tiếp lên tầng chứa nước áp lực Pleistocen, tiếp giáp với các khối nâng chứa nước karst biên rìa.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu được thiết kế theo triết học thực chứng phản biện (post-positivism / critical realism), kết hợp định lượng thực nghiệm đa cấp độ (multi-level empirical design):
- Cấp độ khu vực (Regional Scale): Phân tích cấu trúc kiến tạo, giải đoán địa chấn 2D thềm lục địa bể sông Hồng và mạng lưới đo sâu điện từ trường chuyển (TEM).
- Cấp độ địa phương (Local Scale): Mạng lưới giếng khoan quan trắc ĐCTV quốc gia và dự án quốc tế (VietAS, IGPVN).
- Cấp độ vi mô / vật liệu (Micro/Core Scale): Thí nghiệm ép tách nước lỗ rỗng từ mẫu lõi khoan nguyên dạng trong phòng thí nghiệm chuyên sâu.
┌────────────────────────────────────────────────────────┐
│ SƠ ĐỒ THIẾT KẾ PHƯƠNG PHÁP NGHIÊN CỨU │
└──────────────────────────┬─────────────────────────────┘
│
┌──────────────────────────────────────────┼──────────────────────────────────────────┐
▼ ▼ ▼
┌──────────────────┐ ┌──────────────────┐ ┌──────────────────┐
│ ĐỊA VẬT LÝ DIỆN │ │ KHOAN ĐCTV CHÙM │ │ ĐỊA HÓA ĐỒNG VỊ │
├──────────────────┤ ├──────────────────┤ ├──────────────────┤
│ • 61 điểm TEM │ │ • 02 chùm khoan │ │ • 59 mẫu hóa NDĐ │
│ • 22 giếng nông │ │ VietAS_ND01-02 │ │ • 27 mẫu pore- │
│ • 16 giếng Karota│ │ • 03 chùm hút │ │ water ép màng │
│ điện trở suất │ │ nước thí nghiệm│ │ • 87 mẫu δ18O/δ2H│
│ • Tuyến chấn 2D │ │ • Xác định T, k, │ │ • 32 mẫu 14C DIC │
│ thềm lục địa │ │ mực nước tĩnh │ │ • 08 mẫu Triti/He│
└────────┬─────────┘ └────────┬─────────┘ └────────┬─────────┘
│ │ │
└─────────────────────────────────────────┼──────────────────────────────────────────┘
│
▼
┌────────────────────────────────────────────────────────┐
│ TÍCH HỢP MÔ HÌNH HÓA SỐ GMS │
│ MODFLOW + SEAWAT 3D (Chuỗi quan trắc 1994-2014) │
└────────────────────────────────────────────────────────┘
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý dữ liệu được tiến hành với độ tin cậy khoa học cao:
- Phương pháp trường chuyển (TEM): Đo đạc tại 61 điểm khảo sát phân bố đồng đều trên diện tích nghiên cứu, kết hợp giải đoán bằng phần mềm chuyên dụng để bóc tách điện trở suất theo các lát cắt độ sâu từ 25–30m (hệ tầng Hải Hưng) đến 55–60m (hệ tầng Vĩnh Phúc).
- Phương pháp địa vật lý lỗ khoan: Đo Karota độ dẫn điện và điện trở tại 16 lỗ khoan sâu, thiết lập tương quan hồi quy giữa tổng hàm lượng chất rắn hòa tan ($\text{TDS}$) với điện trở suất đất đá ($\rho$).
- Kỹ thuật ép nước lỗ rỗng (Porewater extraction): Khoan lấy mẫu đất nguyên dạng bằng ống mẫu chuyên dụng tại 02 lỗ khoan nghiên cứu sâu VietAS_ND01 và VietAS_ND02 thuộc lớp sét dẻo hệ tầng Hải Hưng ($mQ_{2}^{1-2}hh$), sử dụng thiết bị ép cơ học áp lực cao của Trung tâm CETASD kết hợp màng lọc vi mô để tách chiết 27 mẫu nước lỗ rỗng nguyên chất dọc theo chiều sâu cột địa tầng.
- Phân tích đồng vị môi trường: Phân tích thành phần đồng vị bền ($\delta^{18}\text{O}$, $\delta^{2}\text{H}$) trên 87 mẫu nước; xác định tuổi tuyệt đối qua hoạt độ phóng xạ $^{14}\text{C}$ trong cacbon vô cơ hòa tan ($\text{DIC}$) trên 32 mẫu tại Viện Khoa học và Kỹ thuật Hạt nhân; phân tích 8 mẫu đồng vị khí trơ và Triti ($^{3}\text{H}$) để xác định thời gian lưu giữ và nhiệt độ cổ khí hậu.
Data và phân tích
Dữ liệu quan trắc mực nước và chất lượng nước liên tục trong 20 năm (1994–2014) từ mạng lưới quan trắc quốc gia được chuẩn hóa. Các mô hình dòng chảy và lan truyền chất tan 3 chiều được xây dựng bằng hệ thống GMS (Groundwater Modeling System) tích hợp module MODFLOW và SEAWAT. Các bài toán hiệu chỉnh tham số (calibration) và kiểm chứng (validation) mực nước tính toán so với mực nước thực tế tại các cụm công trình quan trắc chuẩn (Q108a, Q109a, Q110a) đạt hệ số tương quan rất cao ($R^2 > 0,92$), đảm bảo độ tin cậy tuyệt đối cho các kịch bản mô phỏng.
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Nguồn gốc và cơ chế bổ cập cổ - hiện đại của thấu kính nước nhạt:
Kết quả phân tích đồng vị bền $\delta^{18}\text{O}$ dao động từ -8,5‰ đến -7,2‰ kết hợp tuổi nước dưới đất xác định qua hoạt độ $^{14}\text{C}$ đạt từ 3.000 đến trên 18.000 năm BP chứng minh rằng nước nhạt trong tầng Pleistocen ($qp$) được tích tụ trong giai đoạn băng hà Pleistocen muộn và liên tục nhận nguồn bổ cập xuyên tầng từ các thành tạo đá vôi karst Triat ($T_2đg$) và cát kết Neogen ($N_2vb$) phía tây và tây bắc Ninh Bình. Dòng áp lực từ khối nâng karst di chuyển theo phương Tây Bắc - Đông Nam bổ cập trực tiếp vào bồn trũng trũng sụt, tạo ra động lực đẩy khối nước mặn nguyên thủy ra biển.
-
Quy luật phân bố nồng độ muối dị thường trong lớp thấm nước yếu:
Số liệu phân tích hóa học nước lỗ rỗng tại hai lỗ khoan VietAS_ND01 và VietAS_ND02 cho thấy một phát hiện then chốt: hàm lượng $\text{Cl}^-$ và độ dẫn điện ($\text{EC}$) trong lớp sét biển Holocen đạt cực đại ở phần giữa của lớp sét (hàm lượng $\text{Cl}^-$ lên tới $4.200 - 5.800 \text{ mg/l}$) và giảm dần về phía tiếp xúc với tầng chứa nước nhạt Pleistocen bên dưới. Đồ thị phân bố tạo nên đường cong khuếch tán điển hình theo định luật Fick, minh chứng nước mặn bị giữ lại trong trầm tích biển đang liên tục khuếch tán thẳng đứng xuống tầng nước nhạt.
-
Xác lập hệ phương trình tương quan thực nghiệm có độ tin cậy cao:
Luận án đã thành lập mối tương quan hồi quy hàm mũ và tuyến tính chặt chẽ giữa hàm lượng ion $\text{Cl}^-$, $\text{TDS}$ với độ dẫn điện của tầng đất đá ($\sigma_b$) và điện trở suất ($\rho$), cho phép giải đoán tức thời chất lượng nước ngầm từ tài liệu đo địa vật lý giếng khoan và trường chuyển:
$$\text{TDS} = a \cdot \rho^{-b} \quad (R^2 = 0,89 - 0,94)$$
-
Định lượng động lực xâm nhập mặn do khai thác nhân tạo:
Dữ liệu lịch sử khai thác từ 6.000 giếng khoan UNICEF (khoan từ năm 1982 đến 1995) cùng hàng loạt bãi giếng công nghiệp tập trung đã gây suy giảm mực nước áp lực tầng Pleistocen nghiêm trọng (tốc độ hạ thấp mực nước từ 0,2–0,4 m/năm giai đoạn 1994–2014). Chênh lệch cột áp đã làm đảo chiều đường dòng thủy lực, khiến biên mặn phía bắc và đông bắc dịch chuyển sâu vào trung tâm thấu kính nước nhạt với vận tốc trung bình từ $15 - 35 \text{ m/năm}$.
ĐỘ SÂU (m)
0 ┌────────────────────────────────────────────────────────┐
│ TẦNG CHỨA NƯỚC HOLOCEN NÔNG (Nhiễm mặn/nhạt cục bộ) │
15 ├────────────────────────────────────────────────────────┤
│ LỚP THẤM NƯỚC YẾU HỆ TẦNG HẢI HƯNG (mQ2^1-2 hh) │
│ Nồng độ Cl- cực đại: 4.200 - 5.800 mg/l (Nước biển cổ│
│ quá trình Flandrian còn lưu giữ trong porewater) │
│ │
│ ↓↓ Dòng khuếch tán phân tử + Phân dị trọng lực ↓↓ │
45 ├────────────────────────────────────────────────────────┤
│ TẦNG CÁCH NƯỚC VĨNH PHÚC (Loang lổ laterit phong hóa) │
60 ├────────────────────────────────────────────────────────┤
│ THẤU KÍNH NƯỚC NHẠT TẦNG PLEISTOCEN (a,am Q1^2-3 hn) │
│ (Mực nước hạ thấp do >6.000 giếng khai thác liên tục) │
│ │
│ ◄─── Xâm nhập ngang từ biên mặn Đông Bắc (15-35 m/năm) │
│ ══════════════════════════════════════════════════════ │
│ ▲▲ Bổ cập xuyên tầng từ Karst Triat ($T_2đg$) & Neogen │
120 └────────────────────────────────────────────────────────┘
Implications đa chiều
- Về mặt học thuật: Cung cấp phương pháp luận chuẩn hóa tích hợp giữa địa vật lý trường chuyển (TEM), địa hóa nước lỗ rỗng ép nén và đồng vị thủy văn để nghiên cứu các cấu trúc thấu kính nước nhạt ven biển trên toàn thế giới.
- Về mặt quản lý và quy hoạch: Cảnh báo nguy cơ xóa sổ hoàn toàn thấu kính nước nhạt Nam Định trong vòng 30-50 năm tới nếu duy trì công suất khai thác tự phát hiện nay; đề xuất chuyển dịch vùng khai thác tập trung về phía tây nam (vùng đón đầu nguồn bổ cập từ karst Ninh Bình) và thiết lập hành lang bảo vệ nghiêm ngặt ranh giới mặn - nhạt.
- Về mặt kỹ thuật công trình: Đề xuất các giải pháp kỹ thuật cụ thể như xây dựng mạng lưới giếng bổ cập nhân tạo (artificial recharge wells) bằng nguồn nước mặt mùa lũ đã qua xử lý sơ bộ nhằm duy trì áp lực đẩy mặn; kiểm soát độ sâu đặt ống lọc giếng khai thác để tránh hút phải nước mặn khuếch tán từ nóc tầng chứa nước.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 hạn chế khoa học do điều kiện kỹ thuật và nguồn lực tại thời điểm nghiên cứu:
- Mật độ mẫu nước lỗ rỗng: Thí nghiệm ép tách nước lỗ rỗng chỉ mới thực hiện chi tiết tại 02 vị trí lỗ khoan nghiên cứu sâu (VietAS_ND01 và VietAS_ND02), chưa đủ độ bao phủ không gian để lập bản đồ biến thiên hệ số khuếch tán 3D trên toàn bộ diện tích 2.500 km².
- Mô phỏng đứt gãy trong mô hình số: Các đứt gãy kiến tạo TB-ĐN mới được mô hình hóa dưới dạng các đới dẫn thủy lực tương đương (hydraulic conductance zones) mà chưa tích hợp đầy đủ tính dị hướng phi tuyến và đặc tính cơ nhiệt thủy lực (THMC coupling) của các khe nứt đá móng sâu.
- Tác động biến đổi khí hậu dài hạn: Kịch bản mô hình hóa chưa tích hợp đầy đủ biến số nước biển dâng kết hợp với hiện tượng sụt lún mặt đất cục bộ do nén chặt tầng trầm tích khi khai thác nước ngầm quá mức.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Mở rộng kỹ thuật đo phân tích nước lỗ rỗng bằng phương pháp ly tâm cao tốc và phân tích đồng vị khí hiếm ($^3\text{He}/^4\text{He}$) để hiệu chỉnh chính xác hơn thời gian lưu của nước dưới đất.
- Thiết lập mô hình số tích hợp dòng chảy bề mặt - nước dưới đất - tương tác biển (Coupled Surface-Subsurface & Ocean Model) để đánh giá ảnh hưởng tổng hợp của biến đổi khí hậu toàn cầu.
- Thử nghiệm các mô hình thí nghiệm hiện trường (pilot project) về nạp nhân tạo tầng chứa nước (Aquifer Storage and Recovery - ASR) tại các khu vực ranh giới mặn - nhạt.
Tác động và ảnh hưởng
Công trình luận án tạo ra chuỗi giá trị tác động sâu rộng:
- Giá trị học thuật quốc tế: Đóng góp dữ liệu thực nghiệm quý giá về đồng bằng sông Hồng vào kho tàng nghiên cứu châu thổ thế giới; kết quả đã được công bố trên các tạp chí chuyên ngành quốc tế uy tín thuộc dự án hợp tác Đan Mạch (VietAS - DANIDA) và Cộng hòa Liên bang Đức (IGPVN).
- Tác động kinh tế - xã hội: Cung cấp luận cứ khoa học chuẩn xác để bảo vệ nguồn cấp nước sinh hoạt duy nhất cho hơn 1,5 triệu dân cư tại các huyện duyên hải Nam Định và Ninh Bình, ngăn chặn nguy cơ nhiễm mặn hạ tầng nông nghiệp và nuôi trồng thủy sản.
- Định hình chính sách công: Là tài liệu nền tảng phục vụ Quy hoạch phân bổ tài nguyên nước tỉnh Nam Định đến năm 2020, tầm nhìn 2030; hỗ trợ Cục Quản lý Tài nguyên nước (Bộ Tài nguyên và Môi trường) ban hành các quy định khoanh định vùng cấm, vùng hạn chế khai thác nước dưới đất.
Đối tượng hưởng lợi
┌────────────────────────────────────────────────────────┐
│ ĐỐI TƯỢNG HƯỞNG LỢI │
└──────────────────────────┬─────────────────────────────┘
│
┌────────────────────────────────────────────┼────────────────────────────────────────────┐
▼ ▼ ▼
┌──────────────────────────────────┐ ┌──────────────────────────────────┐ ┌──────────────────────────────────┐
│ GIỚI HỌC THUẬT & │ │ CƠ QUAN QUẢN LÝ NHÀ │ │ DOANH NGHIỆP & CỘNG │
│ NGHIÊN CỨU SINH │ │ NƯỚC & ĐỊA PHƯƠNG │ │ ĐỒNG DÂN CƯ │
├──────────────────────────────────┤ ├──────────────────────────────────┤ ├──────────────────────────────────┤
│ • Tiếp nhận phương pháp luận │ │ • Sở TN&MT Nam Định, Ninh Bình │ │ • Các công ty cấp nước sạch: │
│ tích hợp ĐVL + Porewater + │ │ • Cục Quản lý Tài nguyên Nước │ │ Tối ưu hóa vị trí bãi giếng │
│ Đồng vị 14C/δ18O │ │ • Xây dựng quy hoạch khai thác │ │ • 1,5 triệu dân vùng ven biển: │
│ • Cơ sở dữ liệu chuẩn hóa │ │ bền vững, thiết lập vùng cấm/ │ │ Đảm bảo an ninh nguồn nước ăn │
│ cho các nghiên cứu châu thổ │ │ hạn chế khai thác ngầm │ │ uống, sinh hoạt lâu dài │
└──────────────────────────────────┘ └──────────────────────────────────┘ └──────────────────────────────────┘
- Nghiên cứu sinh và Giới học thuật Địa chất thủy văn: Thừa hưởng phương pháp tiếp cận đa kỹ thuật chuẩn mực và hệ số tương quan thực nghiệm đã được kiểm chứng.
- Nhà hoạch định chính sách và Quản lý tài nguyên: Sở Tài nguyên và Môi trường tỉnh Nam Định, Ninh Bình có căn cứ khoa học rõ ràng để cấp phép khai thác hợp lý.
- Doanh nghiệp cấp nước và Cộng đồng dân cư: Hưởng lợi trực tiếp từ việc duy trì ổn định chất lượng nguồn nước ngầm phục vụ sinh hoạt, nâng cao sức khỏe cộng đồng và phát triển kinh tế bền vữ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 lý thuyết độc đáo nhất là chứng minh và lượng hóa cơ chế bổ cập xuyên tầng từ cấu trúc karst Mesozoi ($T_2đg$) sang bồn trũng Kainozoi qua hệ thống đứt gãy dẫn áp TB-ĐN, bác bỏ hoàn toàn giả thuyết "thấu kính nước nhạt chôn vùi tĩnh". Luận án đã mở rộng Lý thuyết bồn Actezi cổ (Paleo-Artesian Hydrodynamics) và Lý thuyết dòng chảy tỷ trọng biến đổi (Variable-Density Flow Theory) bằng bằng chứng thực nghiệm đồng vị phóng xạ $^{14}\text{C}$ và đồng vị bền $\delta^{18}\text{O}$.
2. Điểm mới vượt trội về mặt phương pháp luận so với các nghiên cứu trước đây là gì?
So với các nghiên cứu chỉ dùng phương pháp đo sâu điện đối xứng đơn thuần (như Sung Ho Song, 2007; Edet, 1993) hoặc mô hình hóa thuần lý thuyết (Đặng Đình Phúc, 1997), luận án đã thiết lập quy trình tích hợp 4 bước đột phá: (1) Phương pháp trường chuyển (TEM) 61 điểm đo độ phân giải cao; (2) Khoan lấy mẫu nguyên dạng và ép tách nước lỗ rỗng (porewater extraction) màng tế vi; (3) Phân tích tổ hợp đa đồng vị ($^{18}\text{O}$, $^{2}\text{H}$, $^{14}\text{C}$, $^{3}\text{H}$, khí trơ); (4) Xây dựng mô hình số SEAWAT mật độ biến đổi được chuẩn hóa bằng chuỗi quan trắc thực tế 20 năm.
3. Phát hiện nào gây bất ngờ nhất (counter-intuitive) và có số liệu thực nghiệm nào chứng minh?
Phát hiện bất ngờ nhất là hiện tượng đỉnh nồng độ mặn (salinity peak) không nằm ở đáy tầng chứa nước mặn nông mà nằm tập trung cục bộ ở phần giữa của lớp sét dẻo hệ tầng Hải Hưng ($mQ_{2}^{1-2}hh$) với hàm lượng ion $\text{Cl}^-$ đạt tới $4.200 - 5.800 \text{ mg/l}$ tại lỗ khoan VietAS_ND01 và VietAS_ND02. Điều này chứng minh rằng lớp "cách nước" thực chất đóng vai trò là "nguồn phát tán muối liên tục" xuống tầng chứa nước nhạt thông qua cơ chế khuếch tán phân tử và phân dị trọng lực theo thời gian địa chất.
MẶT CẮT THỰC NGHIỆM BIẾN THIÊN NỒNG ĐỘ MUỐI THEO ĐỘ SÂU
(Dữ liệu ép nước lỗ rỗng tại lỗ khoan VietAS_ND01)
Độ sâu (m) Hàm lượng Cl- (mg/l)
0 |======== 800 mg/l (Tầng nông)
10 |================== 1.800 mg/l
20 |================================================== 5.800 mg/l (ĐỈNH DỊ THƯỜNG)
30 |====================================== 3.900 mg/l
40 |======================== 2.200 mg/l (Khuếch tán Fick đi xuống)
50 |============= 1.100 mg/l
60 |==== 250 mg/l (Thấu kính nước nhạt Pleistocen)
4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?
Có. Luận án cung cấp chi tiết toàn bộ thông số thiết kế mạng lưới TEM (kích thước khung phát, thời gian trễ, dải tần số), quy trình ép tách cơ học nước lỗ rỗng trong phòng thí nghiệm đạt chuẩn quốc tế của dự án VietAS/CETASD, quy chuẩn lấy và bảo quản mẫu đồng vị phóng xạ $^{14}\text{DIC}$, cũng như toàn bộ ma trận ô lưới, điều kiện biên thủy lực và bộ thông số địa chất thủy văn ($K, T, S, \alpha_L, \alpha_T, D_m$) sử dụng trong mô hình số GMS/SEAWAT.
5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
Chương trình nghiên cứu 10 năm tập trung vào 3 trụ cột: (1) Ứng dụng công nghệ đồng vị vi lượng và đồng vị bền đa nguyên tố để lập bản đồ chi tiết thời gian lưu giữ nước ngầm toàn đồng bằng Bắc Bộ; (2) Thiết lập trạm quan trắc tự động kết hợp truyền tin viễn thám ranh giới mặn - nhạt thời gian thực; (3) Nghiên cứu ứng dụng công nghệ bơm ép nạp nhân tạo tầng chứa nước (ASR) quy mô công nghiệp để chủ động phục hồi trữ lượng thấu kính nước nhạt.
Kết luận
Công trình luận án tiến sĩ của Hoàng Văn Hoan là một đóng góp học thuật xuất sắc, toàn diện và có giá trị thực tiễn to lớn cho ngành địa chất thủy văn Việt Nam với 6 kết luận cốt lõi:
- Xác lập bản chất động lực học của thấu kính nước nhạt Nam Định: được tích tụ trong Pleistocen muộn và nhận bổ cập liên tục từ các thành tạo karst Triat ($T_2đg$) và Neogen phía tây - tây bắc qua hệ thống đứt gãy kiến tạo TB-ĐN.
- Khoanh định chính xác cấu trúc hình học 3D của thấu kính nước nhạt và ranh giới mặn - nhạt với diện tích phân bố trong tầng Pleistocen ở độ sâu 60–120m, trữ lượng tiềm năng khoảng 203.000 m³/ngày.
- Phát hiện và chứng minh thực nghiệm cơ chế xâm nhập mặn thẳng đứng do khuếch tán phân tử và phân dị trọng lực từ nước lỗ rỗng của lớp sét biển Holocen ($mQ_{2}^{1-2}hh$) xuống tầng nước nhạt.
- Lượng hóa quy luật dịch chuyển biên mặn theo phương ngang (15–35 m/năm) dưới áp lực khai thác nước ngầm qua chuỗi quan trắc 20 năm và mô hình số SEAWAT mật độ biến đổi.
- Thiết lập bộ phương trình tương quan hồi quy chuẩn xác giữa các thông số địa vật lý (trường chuyển, Karota) với các chỉ tiêu hóa học ($\text{TDS}, \text{Cl}^-$), tạo công cụ giải đoán gián tiếp hiệu quả cao.
- Đề xuất hệ thống giải pháp công trình và phi công trình có tính khả thi cao, phục vụ trực tiếp cho công tác quy hoạch, quản lý và khai thác bền vững tài nguyên nước dưới đất vùng duyên hải Bắc Bộ trước bối cảnh biến đổi khí hậu toàn cầu.
