Tổng quan về luận án

Nghiên cứu địa chất thủy văn tại các vùng đồng bằng châu thổ ven biển luôn đối mặt với bài toán phức tạp về sự tương tác giữa các khối nước mặn và nước nhạt. 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 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à một công trình tiên phong tiếp cận đa chiều về động thái thấu kính nước nhạt trong không gian trầm tích Kainozoi ven biển Đông Nam đồng bằng Bắc Bộ.

Vùng nghiên cứu trải rộng trên địa bàn các huyện Hải Hậu, Nghĩa Hưng, Giao Thủy, Xuân Trường, Nam Trực, Trực Ninh thuộc tỉnh Nam Định và một phần tỉnh Ninh Bình (Kim Sơn, Yên Khánh, Hoa Lư). Tại đây tồn tại một thấu kính nước dưới đất nhạt quý giá phân bố ở độ sâu từ 60–70 m đến trên 120 m trong tầng chứa nước Pleistocen ($qp$). Trước bối cảnh áp lực khai thác gia tăng phục vụ dân sinh và nguy cơ từ biến đổi khí hậu, nước biển dâng, thấu kính này đang đứng trước nguy cơ nhiễm mặn nghiêm trọng.

Khoảng trống học thuật cốt lõi (research gap) mà luận án tập trung giải quyết là sự thiếu hụt các mô hình định lượng hóa cơ chế xâm nhập mặn phi truyền thống. Trước đây, xâm nhập mặn chủ yếu được diễn giải qua tương tác thủy lực ngang đơn thuần hoặc nêm mặn tĩnh Ghyben-Herzberg từ biển hiện đại. Luận án đã chứng minh hiện tượng xâm nhập mặn thẳng đứng do quá trình khuếch tán phân tử và phân dị trọng lực bắt nguồn từ các lớp trầm tích dính bám thấm nước yếu nguồn gốc biển (aquitards) chứa nước lỗ rỗng mặn cổ (connate water) nằm phủ bên trên.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  • RQ1: Quá trình hình thành và nguồn bổ cập động lực của thấu kính nước nhạt Pleistocen diễn ra như thế nào trong lịch sử tiến hóa địa chất Kainozoi?
    • Giả thuyết H1: Nước nhạt trong tầng Pleistocen được tích tụ từ các thời kỳ cổ khí hậu và liên tục nhận nguồn bổ cập thủy lực từ các thành tạo địa chất cổ hơn (Neogen, Triat) phía dưới theo hướng Tây và Tây Bắc.
  • RQ2: Cấu trúc không gian và ranh giới mặn – nhạt phân bố như thế nào dưới góc độ địa vật lý và thủy địa hóa?
    • Giả thuyết H2: Có sự tương quan hàm số chặt chẽ giữa điện trở suất, độ dẫn điện với hàm lượng ion $Cl^-$ và tổng chất rắn hòa tan (TDS), cho phép giải đoán chính xác ranh giới mặn – nhạt 3D.
  • RQ3: Động lực nào chi phối chính tốc độ và quy mô xâm nhập mặn vào thấu kính?
    • Giả thuyết H3: Sự chênh lệch cột áp do khai thác kết hợp cùng quá trình khuếch tán phân tử và phân dị trọng lực từ lớp sét biển Holocen là cơ chế kép dẫn đến sự thu hẹp thể tích thấu kính nước nhạt.

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 lý thuyết lan truyền chất tan trong môi trường lỗ rỗng (Advection-Dispersion-Diffusion Equation), thủy động lực học dòng chảy phụ thuộc tỷ trọng (Density-dependent flow), địa hóa đồng vị môi trường (Isotope Hydrology) và lý thuyết trầm tích luận Đệ tứ. Quy mô nghiên cứu bao quát dữ liệu quan trắc mực nước và chất lượng giai đoạn 1994–2014, kết hợp khối lượng khảo sát thực nghiệm chuyên sâu: 61 điểm đo trường chuyển (TEM), 22 lỗ khoan nông, 16 lỗ khoan địa vật lý carota, 03 cụm chùm lỗ khoan hút nước thí nghiệm, 59 mẫu hóa nước dưới đất, 27 mẫu ép nước lỗ rỗng, 87 mẫu đồng vị bền ($\delta^{18}\text{O}, \delta^{2}\text{H}, \delta^{13}\text{C}$), 32 mẫu đồng vị phóng xạ $^{14}\text{C}$ và 8 mẫu đồng vị khí trơ/Triti.

Literature Review và Positioning

Nghiên cứu xâm nhập mặn nước dưới đất trên thế giới đã trải qua nhiều giai đoạn phát triển với các trường phái tiếp cận khác nhau:

Trường phái địa chất cấu trúc và lịch sử tiến hóa trầm tích biển tiến – biển thoái được định hình qua các công trình kinh điển của J. De Vries (1981) tại Hà Lan, phân tích mối quan hệ giữa cấu trúc địa chất ven biển và sự phân bố thể nước nhạt. Tiếp nối hướng tiếp cận này, H. Groen (2000) cùng các cộng sự tại Đại học Vrije đã thực nghiệm mô hình máng thấm hai lớp, chỉ ra 4 cơ chế xâm nhập mặn trong thời kỳ biển tiến, khẳng định vai trò của quá trình khuếch tán từ trầm tích hạt mịn. Vincent E. Post (2004) tại Hà Lan cũng đã làm rõ mối quan hệ giữa tiến hóa địa chất Holocen với tuổi và nguồn gốc nước dưới đất thông qua đồng vị phóng xạ.

Trường phái cơ chế dịch chuyển vật chất phụ thuộc tỷ trọng và khuếch tán phân tử ghi nhận các đóng góp của Paschke và Hoopes (1984), Schincariol và Schwartz (1990), Koch và Zhang (1992). Đáng chú ý, George D. Wardlaw và David L. Valentine (2005) tại vùng Salton Sea (California, Mỹ) đã chứng minh sự gia tăng độ mặn theo chiều sâu trong lớp sét đáy hồ tuân theo định luật khuếch tán Fick, với thông lượng khuếch tán muối tính toán đạt $0{,}422 - 0{,}613\text{ g/cm}^2/\text{năm}$. Allen (2006) tại đồng bằng sông Fraser (Canada) xác lập mô hình nước mặn khuếch tán từ cửa sông vào tầng chứa nước và tiếp tục khuếch tán xuống tầng thấm yếu bên dưới.

Trường phái địa hóa đồng vị ứng dụng trong xâm nhập mặn cổ chứng minh rằng nguồn mặn không nhất thiết bắt nguồn từ biển hiện đại. J. Meesters (2000) sử dụng đồng vị $^{37}\text{Cl}$ tại Suriname chứng minh quá trình khuếch tán mặn từ lớp sét biển Holocen bên trên và trầm tích Creta bên dưới vào tầng chứa nước có nguồn gốc khí tượng. Tương tự, Serigne Faye và Piotr Maloszewski (2005) sử dụng $\delta^{18}\text{O}$ và $\delta^{2}\text{H}$ tại Senegal để truy xuất nguồn gốc nước ngầm.

Về phương pháp địa vật lý hiện đại, Peter Bauer Gottwein (2009) tại đồng bằng Okavango (Botswana) và Sung Ho Song (2007) tại Byunsan (Hàn Quốc) đã kết hợp phương pháp trường chuyển điện từ (TEM), carota lỗ khoan với mô hình số SEAWAT để thiết lập tương quan giữa điện trở suất và độ tổng khoáng hóa TDS.

Tồn tại hai luồng quan điểm tranh luận trong giới học thuật: một luồng quan điểm truyền thống xem xâm nhập mặn chủ yếu là quá trình thủy lực ngoại sinh do nêm mặn biển hiện đại xâm nhập theo phương ngang; luồng quan điểm hiện đại cho rằng tại các đồng bằng châu thổ, nguồn mặn nội tại bị bẫy trong các tầng sét trầm tích biển cổ đóng vai trò nguy hiểm và dai dẳng hơn thông qua khuếch tán thẳng đứng.

Luận án định vị nghiên cứu tại điểm giao thoa giữa trầm tích học Đệ tứ bồn trũng Sông Hồng, địa vật lý trường chuyển TEM độ phân giải cao, kỹ thuật ép nước lỗ rỗng độc đáo và mô hình hóa dòng chảy phụ thuộc tỷ trọng. So với công trình của Meesters (2000) tại Suriname và Allen (2006) tại Canada, công trình của Hoàng Văn Hoan giải quyết trọn vẹn bài toán từ nguồn gốc cổ địa lý, cơ chế hóa lý vi mô trong lỗ rỗng sét đến dự báo động lực học 3D dưới áp lực khai thác quy mô lớn.

Đó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 thủy địa hóa trầm tích và cơ học chất lưu trong môi trường lỗ hổng qua các luận điểm then chốt:

Thứ nhất, chứng minh sự tồn tại của cơ chế "thấu mặn ngược" do khuếch tán phân tử và phân dị trọng lực từ tầng thấm nước yếu nguồn gốc biển (hệ tầng Hải Hưng $m\text{Q}_{IV}^{1-2}$) xuống tầng chứa nước Pleistocen ($qp$). Quá trình này diễn ra độc lập với dòng chảy đối lưu thủy lực và trở thành nguồn cung cấp muối liên tục khi có sự xáo trộn áp lực.

Thứ hai, xác lập bản chất cổ khí hậu của thấu kính nước nhạt Nam Định. Dựa trên giá trị đồng vị bền $\delta^{18}\text{O}$ dao động từ $-8‰$ đến $-10‰$ kết hợp với hoạt độ phóng xạ $^{14}\text{C}$ xác định tuổi nước từ vài nghìn năm đến trên 20.000 năm trước hiện tại (BP), luận án bác bỏ nhận định thấu kính nước nhạt được hình thành từ nguồn bổ cập nước mưa tại chỗ cận đại, khẳng định đây là nguồn nước cổ được bảo tồn qua các chu kỳ biển thoái - biển tiến Kainozoi và được nuôi dưỡng liên tục từ các tầng chứa nước Neogen ($N$) và Triat ($T$) phía dưới.

Thứ ba, hoàn thiện mô hình tương quan thực nghiệm địa vật lý - thủy địa hóa. Thiết lập các phương trình tương quan tuyến tính và phi tuyến có hệ số xác định $R^2 > 0{,}85$ giữa các thông số điện trở suất ($\rho$), độ dẫn điện tầng đá ($\sigma$) với hàm lượng ion chỉ thị chloride ($\text{Cl}^-$) và tổng chất rắn hòa tan ($\text{TDS}$).

[Khí quyển / Biển cổ Holocen]
          │
          ▼ (Lắng đọng trầm tích hạt mịn mQIV)
[Lớp sét biển Hải Hưng (Aquitard)] ──(TDS nước lỗ rỗng: 10 - 25 g/l)
          │
          ├─► Khuếch tán phân tử (Fickian Diffusion)
          ├─► Phân dị trọng lực (Density Sinking)
          ▼
[Thấu kính nước nhạt Pleistocen (qp)] ◄─── Bổ cập áp lực sâu từ tầng Neogen (N) & Triat (T)
          ▲                                 (δ18O: -8‰ đến -10‰; Tuổi 14C > 20.000 BP)
          │
   (Hạ thấp mực nước do giếng khai thác)

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

Khung phân tích của luận án tích hợp bốn cấu phần khoa học:

  1. Thuyết dòng chảy và chuyển khối phụ thuộc tỷ trọng (Density-Dependent Groundwater Flow): Sử dụng phương trình bảo toàn khối lượng chất tan kết hợp biến thiên mật độ chất lưu $\rho = \rho_0 (1 + \epsilon C)$.
  2. Định luật khuếch tán Fick trong môi trường xốp: $$J = -D_e \frac{\partial C}{\partial z}$$ với hệ số khuếch tán hiệu dụng $D_e$ được chuẩn hóa qua độ lỗ rỗng và độ uốn khúc của trầm tích sét.
  3. Định luật Archie và chuyển đổi địa vật lý điện trở: Xác lập mối liên hệ giữa điện trở suất dung dịch lỗ rỗng ($R_w$), điện trở suất tầng đá ($R_0$) và hệ số thành hệ ($F$).
  4. Địa tầng phân loại thứ tự (Sequence Stratigraphy): Phân chia các chu kỳ trầm tích đầm hồ, vũng vịnh, biển nông kỷ Đệ tứ vùng châu thổ Sông Hồng.

Điều kiện biên áp dụng khung phân tích: Vùng đồng bằng ven biển ngập triều, nơi có tầng chắn sét biển chứa mặn cổ nằm trực tiếp trên tầng chứa nước áp lực hạt thô bở rời, chịu tác động của áp lực khai thác nước ngầm phân tán và tập trung.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) với phương pháp tiếp cận diễn dịch thực nghiệm kết hợp tam giác đạc phương pháp (methodological triangulation):

  • Thiết kế đa cấp độ (Multi-level design): Khảo sát địa mạo - cấu trúc khu vực $\rightarrow$ Khảo sát địa vật lý trường chuyển bề mặt $\rightarrow$ Đo đạc carota lỗ khoan sâu $\rightarrow$ Thí nghiệm ép nước lỗ rỗng quy mô micromet $\rightarrow$ Phân tích hóa học và đồng vị phân tử $\rightarrow$ Mô phỏng số động lực học 3D.
  • Mẫu nghiên cứu: Chọn mẫu đại diện theo các tuyến mặt cắt vuông góc và song song với đường bờ biển, xuyên qua các tâm hạ thấp mực nước tại Hải Hậu, Nghĩa Hưng, Giao Thủy.
                  ┌──────────────────────────────────────────────┐
                  │          THIẾT KẾ NGHIÊN CỨU TÍCH HỢP        │
                  └──────────────────────┬───────────────────────┘
                                         │
        ┌────────────────────────────────┼────────────────────────────────┐
        ▼                                ▼                                ▼
┌───────────────────────┐ ┌──────────────────────────────┐ ┌──────────────────────────────┐
│  ĐỊA VẬT LÝ DIỆN RỘNG │ │ THÍ NGHIỆM ĐẤT ĐÁ & NƯỚC LỖ  │ │  ĐỊA HÓA ĐỒNG VỊ & THỦY LỰC  │
│ 61 điểm TEM           │ │ RỖNG                         │ │ 87 mẫu δ18O, δ2H, δ13C       │
│ 16 lỗ khoan Carota    │ │ Ép nước lõi khoan VietAS     │ │ 32 mẫu 14C, 8 mẫu khí trơ    │
│ 22 lỗ khoan nông      │ │ Phân tích Cl-, TDS, EC       │ │ 3 chùm hút nước thí nghiệm   │
└───────────┬───────────┘ └──────────────┬───────────────┘ └──────────────┬───────────────┘
            │                            │                                │
            └────────────────────────────┼────────────────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │       MÔ HÌNH HÓA SỐ (GMS / SEAWAT)          │
                  │ Mô phỏng lan truyền mặn 3D (1994 - 2014)     │
                  │ Dự báo các kịch bản khai thác tương lai      │
                  └──────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn quốc tế:

  • Phương pháp trường chuyển điện từ (TEM): Đo tại 61 điểm bằng thiết bị kỹ thuật số chuyên dụng của Cục Địa chất Đan Mạch (GEUS), xử lý nghịch đảo 1D/2D bằng thuật toán Occam, loại trừ triệt để nhiễu công nghiệp để thu được mặt cắt điện trở suất từ độ sâu 0 đến trên 150 m.
  • Carota địa vật lý giếng khoan: Thực hiện trên 16 lỗ khoan sâu, đo độ dẫn điện tầng đá ($\sigma_a$), điện trở suất biểu kiến, thế tự nhiên ($SP$) và bức xạ gamma tự nhiên ($GR$).
  • Khoan lấy mẫu nguyên dạng và ép nước lỗ rỗng (Porewater extraction): Tiến hành tại 2 lỗ khoan nghiên cứu trọng điểm VietAS_ND01 và VietAS_ND02. Mẫu sét nguyên dạng được lấy bằng ống mẫu chuyên dụng, bảo quản kín và ép trích xuất dung dịch nước lỗ rỗng tại phòng thí nghiệm của Trung tâm CETASD (Đại học Khoa học Tự nhiên Hà Nội) theo quy trình chuẩn hóa của Đan Mạch.
  • Phân tích đồng vị: Đồng vị bền ($\delta^{18}\text{O}, \delta^{2}\text{H}$) đo bằng máy quang phổ khối lượng laser; đồng vị phóng xạ $^{14}\text{C}$ phân tích qua phương pháp đo hoạt độ phóng xạ của tổng cacbon vô cơ hòa tan (DIC) tại Viện Năng lượng Nguyên tử Quốc tế (IAEA); đồng vị khí trơ phân tích tại các phòng thí nghiệm chuyên biệt châu Âu.
  • Thí nghiệm ĐCTV: Tiến hành bơm hút nước chùm tại 03 cụm lỗ khoan thí nghiệm để xác định chính xác hệ số dẫn nước ($T$), hệ số thấm ($K$) và hệ số nhả nước ($S$) của tầng chứa nước Pleistocen.

Data và phân tích

Dữ liệu được phân tích qua hệ thống công cụ định lượng tiên tiến:

  • Thiết lập tương quan hồi quy giữa độ dẫn điện của tầng với hàm lượng $\text{Cl}^-$ và $\text{TDS}$ trong nước lỗ rỗng và nước dưới đất, cho kết quả tương quan tuyến tính đồng biến chặt chẽ.
  • Mô hình hóa số: Ứng dụng hệ thống mô hình hóa nước dưới đất chuyên dụng GMS (Groundwater Modeling System) tích hợp module MODFLOW giải bài toán dòng chảy và SEAWAT giải phương trình truyền khối phụ thuộc tỷ trọng 3 chiều.
  • Hiệu chỉnh và kiểm định (Calibration & Validation): Mô hình được hiệu chỉnh dựa trên chuỗi quan trắc liên tục 20 năm (1994–2014) tại các trạm quan trắc quốc gia $Q92, Q108a, Q109a, Q110a, Q225b$. Chỉ số sai số căn phương trung bình (RMSE) đạt mức tối ưu dưới $0{,}05\text{ m}$ về mực nước và dưới $5%$ về nồng độ chất tan.

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

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

  1. Bản chất nguồn gốc cổ và nguồn nuôi dưỡng thấu kính nước nhạt: Kết quả phân tích đồng vị bền cho thấy $\delta^{18}\text{O}$ trong nước nhạt tầng Pleistocen dao động từ $-8{,}0‰$ đến $-10{,}2‰$, âm hơn đáng kể so với nước mưa hiện đại ($\delta^{18}\text{O}$ trung bình khoảng $-5{,}5‰$ đến $-6{,}5‰$) và nước sông Hồng. Tuổi phóng xạ $^{14}\text{C}$ xác định tuổi nước đạt từ $4.500$ đến $23.000\text{ năm BP}$. Điều này khẳng định nước nhạt được hình thành vào thời kỳ băng hà Pleistocen muộn khi mực nước biển thấp hơn hiện nay hơn $100\text{ m}$. Hướng dòng chảy khu vực chứng minh thấu kính được bổ cập áp lực liên tục từ phía Tây và Tây Bắc thông qua các cửa sổ địa chất thủy văn liên thông với tầng chứa nước Neogen và đá gốc Triat.

  2. Lớp thấm nước yếu nguồn gốc biển là "túi muối" lưu trữ mặn cổ: Kết quả ép nước lỗ rỗng tại các lỗ khoan VietAS_ND01 và VietAS_ND02 phát hiện tổng chất rắn hòa tan (TDS) trong nước lỗ rỗng của tầng sét Hải Hưng ($m\text{Q}_{IV}^{1-2}$) và Vĩnh Phúc ($m\text{Q}_I^{3vp}$) đạt từ $10{,}5\text{ g/l}$ đến $25{,}8\text{ g/l}$, hàm lượng $\text{Cl}^-$ đạt trên $12.000\text{ mg/l}$. Biểu đồ biến thiên theo chiều sâu cho thấy nồng độ muối tăng dần từ mặt đất xuống đáy tầng sét, tạo nên một gradient nồng độ thẳng đứng rất lớn hướng thẳng vào nóc tầng chứa nước Pleistocen.

  3. Xác lập cơ chế xâm nhập mặn thẳng đứng do khuếch tán và phân dị trọng lực: Thí nghiệm xác định hệ số khuếch tán hiệu dụng phân tử ($D_e$) của ion $\text{Cl}^-$ trong trầm tích sét biển đạt giá trị $1{,}8 \times 10^{-10}\text{ m}^2/\text{s}$ đến $3{,}2 \times 10^{-10}\text{ m}^2/\text{s}$. Luận án chứng minh rằng ngay cả khi chưa có sự khai thác nước, dòng khuếch tán phân tử kết hợp với xu thế chìm lắng của các dung dịch có tỷ trọng cao (hiệu ứng phân dị trọng lực) đã làm mặn hóa phần trên của tầng Pleistocen.

  4. Biến động dịch chuyển biên mặn dưới tác động khai thác tập trung: Tổng lưu lượng khai thác trong vùng tăng vọt với sự tham gia của các nhà máy nước tập trung và hàng chục nghìn giếng khoan nông hộ gia đình (bao gồm chương trình nước sạch UNICEF). Mực nước áp lực tầng Pleistocen tại giếng quan trắc $Q109a$ hạ thấp liên tục qua các năm, làm suy giảm cột áp thủy lực của thấu kính nước nhạt. Dưới tác động này, 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 với tốc độ $15 - 35\text{ m/năm}$, đồng thời gia tăng lưu lượng thấm xuyên kéo nước mặn từ tầng sét biển bên trên xuống tầng chứa nước.

Độ sâu (m)  Mặt cắt thạch học          TDS nước lỗ rỗng (g/l)      Cơ chế chi phối
 0 ┬─────── ┌──────────────────────┐   0,5 - 1,5 g/l              [Thấm ngấm tự nhiên]
   │        │ Sét bột Holocen nông │
20 ┼─────── ├──────────────────────┤   10,5 - 18,2 g/l            [Bẫy mặn biển cổ]
   │        │ Sét biển Hải Hưng    │                              │
40 ┼─────── │ (mQIV1-2)            │   20,1 - 25,8 g/l (Max)      ▼ Khuếch tán Fick &
   │        ├──────────────────────┤                              ▼ Phân dị trọng lực
60 ┼─────── │ Sét Vĩnh Phúc (Q13vp)│   Giảm dần xuống 3,0 g/l     │
   │        ├──────────────────────┤                              ▼
80 ┼─────── │ Cát cuội Pleistocen  │   < 1,0 g/l (Thấu kính nhạt) ◄── Bổ cập áp lực sâu
   │        │ (qp - Hà Nội)        │                              (Tuổi 14C > 20.000 BP)
   ┴─────── └──────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào kho tàng địa chất thủy văn quốc tế dẫn liệu thực nghiệm về cơ chế mặn hóa thẳng đứng từ trầm tích biển cổ tại các đồng bằng châu thổ sông lớn châu Á.
  • Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa phương pháp địa vật lý bề mặt (TEM), carota lỗ khoan, kỹ thuật ép nước lỗ rỗng và mô hình số SEAWAT, tạo cẩm nang phương pháp có thể chuyển giao cho các nghiên cứu tại đồng bằng sông Cửu Long và ven biển miền Trung.
  • Về thực tiễn và chính sách:
    • Khoanh định ranh giới không gian 3D của thấu kính nước nhạt (ngưỡng $\text{TDS} < 1{,}0\text{ g/l}$) với độ chính xác cao.
    • Đề xuất giới hạn mực nước hạ thấp cho phép ($H_{cp}$) cho từng tiểu vùng khai thác tại Nam Định.
    • Khuyến nghị quy hoạch chuyển đổi vị trí các bãi giếng khai thác nước ngầm tập trung dịch chuyển về phía Tây - Tây Nam (khu vực có nguồn bổ cập dồi dào từ tầng Neogen), tránh xa đới nén mặn phía Đông Bắc và dải ven biển.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn ghi nhận một số giới hạn nghiên cứu:

  • Số lượng lỗ khoan ép nước lỗ rỗng chuyên sâu còn hạn chế (02 vị trí), chưa phủ kín toàn bộ diện tích biến tính thạch học phức tạp của hệ tầng Hải Hưng trên toàn tỉnh Nam Định.
  • Số lượng mẫu phân tích đồng vị khí trơ ($n=8$) còn khiêm tốn do chi phí kỹ thuật cao, chưa cho phép xây dựng mô hình nhiệt độ cổ khí hậu liên tục theo thang thời gian chi tiết.
  • Mô hình mô phỏng SEAWAT chưa tích hợp đầy đủ hiện tượng lún sụt đất (land subsidence) do suy giảm áp lực nước lỗ rỗng trong các lớp sét dẻo mịn.

Các định hướng nghiên cứu tiếp nối trong tương lai:

  1. Mở rộng ứng dụng đồng vị phi truyền thống như Chlorine-36 ($^{36}\text{Cl}$), Strontium ($^{87}\text{Sr}/^{86}\text{Sr}$) và Boron ($\delta^{11}\text{B}$) để phân tách tuyệt đối các nguồn mặn hỗn hợp (nước biển hiện đại, nước biển cổ ngưng đọng, và hòa tan khoáng vật).
  2. Xây dựng mô hình ghép đôi (coupled hydro-mechanical model) mô phỏng đồng thời dòng chảy tỷ trọng biến đổi, lan truyền chất tan và biến dạng nén lún cơ học của các tầng trầm tích Đệ tứ.
  3. Nghiên cứu thực nghiệm các giải pháp công trình kiểm soát xâm nhập mặn: thiết lập đập chắn thủy lực (hydraulic barrier) bằng giếng phun nhân tạo, và bổ cập nhân tạo tầng chứa nước (Managed Aquifer Recharge - MAR) trong mùa mưa lũ.
  4. Mở rộng nghiên cứu sang tầng chứa nước Neogen sâu hơn nhằm đánh giá tổng thể tiềm năng tài nguyên nước dưới đất toàn bồn trũng Sông Hồng.

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

Luận án tạo ra những tác động khoa học và xã hội rõ nét:

  • Tác động học thuật: Đặt nền móng cho chuỗi nghiên cứu chuyên sâu về thủy văn đồng vị và địa vật lý điện từ tại Việt Nam; tạo nguồn trích dẫn quan trọng cho các công bố quốc tế về xâm nhập mặn vùng châu thổ sông Hồng trên các tạp chí chuyên ngành như Journal of Hydrology, Hydrogeology Journal.
  • Tác động công nghiệp và cấp nước: Cung cấp thông số địa chất thủy văn chuẩn xác ($T, K, S, D_e$) cho các công ty tư vấn khảo sát thiết kế cấp thoát nước, giúp tối ưu hóa chiều sâu đặt ống lọc giếng khai thác, tránh đới nước lợ chuyển tiếp.
  • Tác động chính sách: Cung cấp cơ sở khoa học cốt lõi phục vụ Sở Tài nguyên và Môi trường tỉnh Nam Định và Bộ Tài nguyên và Môi trường trong việc xây dựng Quy hoạch phân bổ và bảo vệ tài nguyên nước dưới đất, xác định vùng hạn chế khai thác nước dưới đất theo Nghị định của Chính phủ.
  • Lợi ích xã hội: Bảo vệ nguồn nước ngọt duy nhất của hơn 1 triệu cư dân nông thôn vùng ven biển Nam Định - Ninh Bình, ngăn chặn nguy cơ mặn hóa không thể phục hồi của nguồn tài nguyên thiên nhiên chiến lược.

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

  • Nghiên cứu sinh và học viên cao học ngành Địa chất Thủy văn, Địa kỹ thuật: Tiếp cận phương pháp luận mẫu mực về kết hợp đa công cụ địa vật lý - đồng vị - mô hình hóa số.
  • Các nhà khoa học và chuyên gia tài nguyên nước: Kế thừa tập dữ liệu quan trắc và thực nghiệm quý giá về đồng vị, thành phần hóa học nước lỗ rỗng và các thông số thủy động lực học trầm tích Đệ tứ.
  • Kỹ sư và doanh nghiệp khai thác nước ngầm: Ứng dụng bản đồ ranh giới mặn - nhạt và các phương trình tương quan điện trở suất - TDS để tối ưu hóa công tác khoan thăm dò và khai thác an toàn.
  • Cơ quan quản lý nhà nước (Bộ TN&MT, Cục Quản lý Tài nguyên nước, UBND các tỉnh ven biển): Sở hữu công cụ dự báo và luận cứ khoa học tin cậy để ban hành chính sách quản lý, cấp phép khai thác và bảo vệ tầng chứa nước trước thách thức biến đổi khí hậu.

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 định lượng hóa 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ừ tầng sét biển chứa nước mặn cổ xuống tầng chứa nước ngầm áp lực. Nghiên cứu đã mở rộng lý thuyết dịch chuyển chất tan truyền thống (Advection-Dispersion Equation) bằng việc bổ sung vai trò quyết định của gradient nồng độ nội tại trong môi trường thấm nước yếu (aquitard-driven salinization) mà không cần dòng đối lưu thủy lực cưỡng bức.

2. Đột phá phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế trước đây?
Luận án kết hợp sáng tạo phương pháp trường chuyển điện từ (TEM) độ phân giải sâu với kỹ thuật ép nước lỗ rỗng (porewater extraction) từ mẫu sét nguyên dạng và phân tích đồng vị môi trường ($\delta^{18}\text{O}, \delta^{2}\text{H}, \delta^{13}\text{C}, ^{14}\text{C}$). So với nghiên cứu của Sung Ho Song (2007) chỉ dùng đo sâu điện thông thường hay Meesters (2000) chỉ tập trung vào mô hình hóa 1D, quy trình tích hợp của luận án giải quyết trọn vẹn bài toán từ cấu trúc không gian diện rộng đến vi cấu trúc lỗ rỗng và kiểm chứng bằng mô hình số 3D SEAWAT.

3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là hàm lượng muối trong nước lỗ rỗng của tầng sét Hải Hưng ở độ sâu 25–40 m cao gấp nhiều lần độ mặn trung bình của nước dưới đất trong vùng (TDS đạt cực đại trên $25\text{ g/l}$), cùng với giá trị đồng vị $\delta^{18}\text{O}$ của tầng Pleistocen âm sâu ($-8‰$ đến $-10‰$), chứng minh rằng thấu kính nước nhạt hoàn toàn không phải nước mưa hiện đại ngấm xuống mà là nguồn nước cổ tích tụ từ kỷ băng hà được bảo tồn nguyên vẹn qua hàng vạn năm.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (replication protocol) không?
Luận án cung cấp chi tiết toàn bộ quy trình thực nghiệm: từ sơ đồ bố trí mạng lưới đo TEM (khoảng cách điểm đo, cấu hình vòng dây), quy trình carota giếng khoan, phương pháp bảo quản và kỹ thuật ép mẫu sét trích xuất nước lỗ rỗng, đến các bước thiết lập ô lưới, gán điều kiện biên và hiệu chỉnh thông số trong phần mềm GMS/SEAWAT.

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?
Chương trình nghiên cứu dài hạn hướng tới: (1) Ứng dụng kỹ thuật đồng vị bền của clo ($^{37}\text{Cl}$) và đồng vị phóng xạ Clor-36 ($^{36}\text{Cl}$) để định tuổi nước mặn cổ chính xác tuyệt đối; (2) Tích hợp mô hình biến dạng cơ học đất và sụt lún mặt đất với mô hình dòng chảy mặn 3D; (3) Triển khai thử nghiệm công nghệ bổ cập nhân tạo (MAR) tại các cấu trúc địa chất thuận lợi phía Tây Bắc Nam Định.

Kết luận

Luận án tiến sĩ của Hoàng Văn Hoan là một công trình khoa học xuất sắc, giải quyết thấu đáo và toàn diện bài toán xâm nhập mặn nước dưới đất tại vùng ven biển Nam Định. Những đóng góp cốt lõi bao gồm:

  1. Xác lập cơ sở khoa học về lịch sử hình thành thấu kính nước nhạt Pleistocen Nam Định: nguồn gốc cổ khí hậu Pleistocene muộn, được bảo tồn và bổ cập liên tục từ các thành tạo Neogen và Triat phía Tây - Tây Bắc.
  2. Định lượng hóa thành công 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ừ tầng sét biển Holocen (hệ tầng Hải Hưng), mở ra bước chuyển paradigm từ tư duy "xâm nhập mặn ngang từ biển hiện đại" sang "xâm nhập mặn nội tại từ bẫy trầm tích mặn cổ".
  3. Thiết lập hệ thống phương trình tương quan thực nghiệm có độ tin cậy cao giữa các thông số địa vật lý (điện trở suất, độ dẫn điện) với các chỉ số thủy địa hóa ($\text{TDS}, \text{Cl}^-$).
  4. Khoanh định chính xác ranh giới mặn – nhạt không gian 3 chiều của thấu kính nước nhạt, đánh giá định lượng tốc độ dịch chuyển biên mặn dưới tác động khai thác thông qua mô hình số SEAWAT.
  5. Đề xuất hệ thống giải pháp kỹ thuật, quy hoạch phân vùng khai thác và ngưỡng hạ thấp mực nước an toàn phục vụ bảo vệ và khai thác bền vững tài nguyên nước dưới đất tỉnh Nam Định.
  6. Mở ra 3 nhánh nghiên cứu học thuật mới: thủy địa hóa nước lỗ rỗng trầm tích hạt mịn, mô hình hóa dòng chảy tỷ trọng biến đổi kết hợp sụt lún đất, và quản trị bổ cập tầng chứa nước ven biển (MAR).

Công trình khẳng định vị thế khoa học ngang tầm với các nghiên cứu thủy địa hóa quốc tế tại Hà Lan, Bắc Mỹ và vùng Nam Á, để lại di sản thực tiễn quan trọng cho công tác quản lý tài nguyên nước và thích ứng biến đổi khí hậu tại các vùng đồng bằng châu thổ Việt Nam.