Tổng quan về luận án

Luận án này tập trung vào nghiên cứu cơ chế rửa mặn của nước lỗ rỗng (NLR) trong các tướng trầm tích biển tuổi Holocen tại khu vực Châu thổ Sông Hồng (CTSH), một khu vực có mật độ dân số cao, đạt 983 người/km2 vào năm 2014, với nhu cầu sử dụng nước lớn và nguồn nước dưới đất (NDĐ) đang đối mặt với nguy cơ suy thoái về chất lượng và trữ lượng. Các nghiên cứu gần đây chỉ ra rằng NDĐ bị nhiễm mặn không chỉ ở ven biển mà còn sâu trong lục địa, lên tới 70km cách bờ biển hiện tại [8, 66, 94]. Điều này không thể giải thích hoàn toàn bằng xâm nhập mặn (XNM) hiện đại, vốn chỉ ăn sâu 20-40km vào các cửa sông [93, 17, 25, 79]. Bối cảnh này tạo nên một nghiên cứu tiên phong trong việc làm sáng tỏ nguồn gốc nhiễm mặn NDĐ ở vùng sâu trong lục địa CTSH.

Research gap cụ thể mà luận án giải quyết là sự thiếu vắng các công trình nghiên cứu chuyên sâu về XNM cổ và cơ chế rửa mặn của nước mặn tàn dư (NMTD) trong các tầng trầm tích biển tại CTSH. Hầu hết các công trình trong nước tập trung vào XNM hiện đại do quá trình tự nhiên và khai thác quá mức [Vũ Kim Tuyến, 1994; Đặng Hữu Ơn, 1996; Phạm Quý Nhân, 2000, 2008]. Luận án nhấn mạnh: "Hiện nay trong nước chưa có công trình nào tập trung vào nghiên cứu các quá trình rửa mặn NLR trong các tầng trầm tích biển và ảnh hưởng của chúng tới phân bố mặn nhạt trong các TCN." Nghiên cứu này đặt ra giả thuyết rằng NMTD trong các tầng sét biển Holocen là nguyên nhân chính gây nhiễm mặn cho TCN Pleistocen. Ban đầu, luận án có phạm vi nghiên cứu là các trầm tích Đệ Tứ, nhưng thông qua kết quả mô hình 1 chiều và đo địa vật lý lỗ khoan (ĐVLLK), nhận thấy NLR trong trầm tích biển tuổi Pleistocen muộn đã bị rửa mặn hoàn toàn (từ 60 nghìn năm tới hàng triệu năm), trong khi NLR trong trầm tích biển Holocen vẫn còn NMTD do thời gian rửa mặn ngắn (<11 nghìn năm). Do đó, nghiên cứu thu hẹp phạm vi để đảm bảo tính xác thực và ý nghĩa của kết quả.

Các câu hỏi nghiên cứu chính (Research Questions) được đặt ra gồm:

  1. Đặc điểm phân bố mặn nhạt của NLR trong trầm tích biển Holocen là gì?
  2. Các cơ chế rửa mặn của NLR trong trầm tích biển Holocen là gì?
  3. Quá trình rửa mặn NLR trong trầm tích biển Holocen ảnh hưởng như thế nào tới TCN Pleistocen?

Khung lý thuyết của luận án tích hợp các nguyên lý về dịch chuyển vật chất trong môi trường lỗ hổng, bao gồm dịch chuyển đối lưu (Advection), khuếch tán phân tử (Diffusion) và phân dị trọng lực (Density flow). Đặc biệt, nghiên cứu áp dụng Định luật Fick thứ nhất để xác định hệ số khuếch tán hữu hiệu và Định luật Darcy để xác định hệ số thấm, cùng với việc tính toán hệ số Rayleigh (Ra) để phân định cơ chế rửa mặn. Mô hình SEAWAT, một công cụ tiên tiến cho mô phỏng dịch chuyển vật chất có tính đến khối lượng riêng chất lỏng thay đổi, là trọng tâm trong khung lý thuyết này.

Đóng góp đột phá của luận án nằm ở việc xác định cụ thể nguồn gốc nhiễm mặn NDĐ ở sâu trong lục địa CTSH là từ NMTD trong trầm tích biển Holocen, chứ không phải XNM hiện đại. Luận án đã định lượng tác động của thời kỳ biển tiến Flandrian và tầng sét biển Holocen tới TCN Pleistocen. Quan trọng hơn, nó xác định hai cơ chế rửa mặn chính dựa trên thành phần thạch học và hệ số Rayleigh: "Nước lỗ rỗng chứa trong tầng sét biển Holocen (với hệ số Rayleigh Ra≤40) bị rửa mặn theo cơ chế khuếch tán phân tử, nước lỗ rỗng chứa trong các trầm tích cát mịn pha sét hoặc sét pha cát (với hệ số Rayleigh Ra >40) bị rửa mặn theo cơ chế dịch chuyển vật chất do phân dị trọng lực." Phạm vi nghiên cứu bao gồm 38 điểm ĐVLLK, 210 điểm đo phương pháp trường chuyển (TEM), khoan và lấy mẫu tại 2 lỗ khoan Q87-VA và Q88-VA, cùng hàng chục mẫu phân tích hóa học NLR, độ hạt, khoáng vật sét, hệ số thấm và khuếch tán. Điều này mang lại ý nghĩa khoa học sâu sắc trong việc lý giải quy luật phân bố mặn nhạt của NDĐ trong TCN Pleistocen, đồng thời có ý nghĩa thực tiễn to lớn trong việc định hướng tìm kiếm, thăm dò, khai thác và bảo vệ nguồn tài nguyên NDĐ tại CTSH.

Literature Review và Positioning

Các công trình nghiên cứu về XNM NDĐ nói chung và XNM cổ nói riêng đã được khảo sát rộng rãi trên thế giới [40, 41, 44, 48, 47, 52, 55, 58, 59, 62, 68, 70, 81, 82, 83, 97, 98]. Các nhà khoa học đã chỉ ra rằng NMTD trong các tầng trầm tích biển hình thành vào các thời kỳ biển tiến trước đây là nguồn gây nhiễm mặn cho các TCN xung quanh. Chẳng hạn, Arie Issar (1981) đã nghiên cứu tốc độ rửa mặn NLR và NMTD ở Israel, sử dụng tài liệu hóa học NDĐ, đồng vị bền và cấu trúc địa chất. Groen và cộng sự (2000) tại Suriname đã sử dụng phương pháp đồng vị 37Cl và mô hình khuếch tán 1 chiều (trong khoảng 4000-6000 năm và 500-10000 năm) để chứng minh nguồn gốc XNM của NDĐ là do các tầng trầm tích chứa NMTD từ tầng sét biển Holocen bên trên và trầm tích tuổi Kreta bên dưới. Tương tự, Hiroshiro và cộng sự (2006) đã lý giải nguồn gốc nước mặn trong TCN ở Motooka, Nhật Bản, là NMTD chứ không phải XNM hiện đại, dù NDĐ được khai thác lưu lượng lớn, dựa trên dữ liệu về độ dẫn điện, thế oxy hóa khử, đồng vị bền (2H, 18O), phóng xạ (3H) và định tuổi 14C. Bruno Capaccioni và cộng sự (2005) nghiên cứu đồng bằng Catania (Ý) đã chỉ ra rằng quá trình rửa mặn diễn ra chậm hơn ở TCN sâu so với tầng nông, nhấn mạnh vai trò của nước mưa trong việc thúc đẩy quá trình rửa mặn. Gossel và cộng sự (2010) đã sử dụng phần mềm FEFLOW để mô hình hóa dòng chảy và dịch chuyển vật chất trong khoảng 140 nghìn năm ở vùng Nubian, tính đến phân tán, khuếch tán và sự khác biệt khối lượng riêng. Beekman và cộng sự (2011) đã xây dựng mô hình 1 chiều để mô phỏng dịch chuyển của ion 35Cl- và 37Cl- từ NLR trong trầm tích dưới hồ Marken, Hà Lan, cho thấy tốc độ khuếch tán của 35Cl- nhanh hơn 1,0017 lần so với 37Cl-. Các nghiên cứu này đều sử dụng tổ hợp các phương pháp địa hóa, đồng vị, địa vật lý và mô hình để giải quyết bài toán XNM và rửa mặn.

Ngược lại, các công trình nghiên cứu trong nước, như của Vũ Kim Tuyến (1994) về đồng vị NDĐ đồng bằng Bắc Bộ, Đặng Hữu Ơn (1996) về dự báo XNM ở Bà Rịa - Vũng Tàu, Nguyễn Văn Hoàng (2001) về lưu lượng thoát nước tầng Pleistocen, Nguyễn Nhƣ Trung và cộng sự (2004) về dự báo XNM ở Hải Phòng bằng mô hình điện trở và ĐCTV, và Trịnh Hoài Thu (2015) về ảnh hưởng khai thác nước ngầm tới XNM trong TCN Pleistocen ở 4 tỉnh ven biển CTSH, chủ yếu tập trung vào XNM hiện đại do khai thác quá mức hoặc quá trình tự nhiên ở khu vực ven biển. Phạm Quý Nhân (2000) đã nghiên cứu XNM do khai thác NDĐ ở CTSH bằng mô hình dòng chảy và mô hình dịch chuyển chất hòa tan, nhưng "chưa tính đến sự thay đổi khối lượng riêng chất lỏng".

Sự khác biệt rõ ràng giữa các nghiên cứu quốc tế (chủ yếu xem xét NMTD và XNM cổ) và các nghiên cứu trong nước (tập trung XNM hiện đại) đã tạo ra một khoảng trống (gap) lớn trong hiểu biết về nguồn gốc nhiễm mặn NDĐ ở các khu vực sâu trong lục địa CTSH. Luận án này định vị mình là công trình tiên phong trong nước, trực tiếp giải quyết vấn đề "XNM cổ" và "quá trình rửa mặn của NLR chứa trong các tầng trầm tích biển ở toàn vùng CTSH". Bằng cách này, luận án không chỉ tổng hợp mà còn lấp đầy khoảng trống này, thúc đẩy lĩnh vực thủy địa chất Việt Nam tiến xa hơn. Nó so sánh trực tiếp với các nghiên cứu quốc tế như của Groen và cộng sự (2000) và Beekman và cộng sự (2011) về phương pháp mô hình 1 chiều và nghiên cứu khuếch tán trong trầm tích biển, nhưng mở rộng phạm vi ứng dụng sang một châu thổ lớn ở vùng nhiệt đới với lịch sử địa chất Đệ Tứ phức tạp, chưa từng được nghiên cứu sâu về khía cạnh này.

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

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

Luận án này mở rộng và thách thức các lý thuyết hiện có về dịch chuyển vật chất trong môi trường lỗ hổng bằng cách áp dụng chúng vào bối cảnh phức tạp của châu thổ sông bị ảnh hưởng bởi dao động mực nước biển Quaternary. Cụ thể, nó mở rộng ứng dụng của Định luật Fick thứ nhất và Định luật Darcy bằng cách tích hợp phân tích cơ chế dịch chuyển vật chất có tính đến sự thay đổi khối lượng riêng chất lỏng thông qua hệ số Rayleigh. Điều này cho phép phân biệt rõ ràng các cơ chế rửa mặn chi phối (khuếch tán phân tử và phân dị trọng lực) dựa trên đặc điểm thạch học của trầm tích.

Khung khái niệm của luận án (Chương 4) hình dung sự tương tác động học giữa trầm tích biển Holocen, tầng chứa nước Pleistocen bên dưới, và các cơ chế rửa mặn khác nhau, tất cả đều bị chi phối bởi lịch sử dao động mực nước biển và hình thành thung lũng cắt xẻ. Nó minh họa một hệ thống mà tại đó nước mặn tàn dư trong tầng sét biển Holocen dần di chuyển vào các tầng chứa nước liền kề.

Mô hình lý thuyết được xây dựng dựa trên các luận điểm bảo vệ chính:

  1. Các trầm tích biển Holocen, hình thành trong thời kỳ biển tiến Flandrian ở CTSH, vẫn còn chứa nước mặn tàn dư với độ tổng khoáng hóa (TDS) tăng dần từ lục địa ra biển. Vùng nước mặn (TDS > 3g/l) chiếm phần lớn diện tích, tập trung ở hai thung lũng cắt xẻ và vùng ven biển.
  2. Cơ chế rửa mặn của NLR phụ thuộc vào thành phần thạch học của trầm tích và hệ số Rayleigh (Ra): NLR trong tầng sét biển Holocen (Ra≤40) bị rửa mặn theo cơ chế khuếch tán phân tử, trong khi NLR trong trầm tích cát mịn pha sét hoặc sét pha cát (Ra >40) bị rửa mặn theo cơ chế dịch chuyển vật chất do phân dị trọng lực.
  3. Quá trình rửa mặn NLR trong trầm tích biển Holocen có ảnh hưởng đáng kể đến TCN Pleistocen, đặc biệt là ở các khu vực thung lũng cắt xẻ hình thành trong giai đoạn băng hà Wurm 2.

Luận án này đề xuất một sự thay đổi mô hình trong việc giải thích nguồn gốc nhiễm mặn ở các khu vực sâu trong lục địa của các châu thổ nhiệt đới. Thay vì chỉ xem xét XNM hiện đại, luận án cung cấp bằng chứng rằng "NLR chứa trong các trầm tích biển hình thành vào thời kỳ biển tiến Holocen vẫn còn chứa nước mặn tàn dư" do thời gian rửa mặn ngắn (<11 nghìn năm), so với các trầm tích Pleistocen đã rửa mặn hoàn toàn. Điều này đã thách thức quan điểm phổ biến trong thủy địa chất Việt Nam.

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

Khung phân tích của luận án là độc đáo bởi sự tích hợp đa ngành. Nó kết hợp sâu sắc lý thuyết về Địa chất thủy văn và Địa chất Đệ Tứ với các nguyên lý Hóa học nước (phân tích thành phần hóa học NLR, NDĐ), Vật lý địa cầu (ĐVLLK, TEM) và Toán học ứng dụng (mô hình số SEAWAT).

Luận án sử dụng một cách tiếp cận phân tích mới bằng cách định lượng hóa cơ chế rửa mặn thông qua hệ số Rayleigh (Ra), một tham số không thường xuyên được áp dụng trong các nghiên cứu tương tự ở Việt Nam. Hệ số Ra được tính toán dựa trên các thông số thực nghiệm (hệ số thấm K, hệ số khuếch tán De) và các đặc tính của chất lỏng (khối lượng riêng, độ nhớt) cùng chiều dày tầng trầm tích. Điều này cho phép đưa ra bằng chứng định lượng vững chắc cho việc phân biệt giữa khuếch tán phân tử và phân dị trọng lực là cơ chế chi phối.

Các đóng góp khái niệm bao gồm việc phân định rõ ràng giữa "nước lỗ rỗng" (porewater, nước chứa trong các tầng trầm tích hạt mịn như sét) và "nước dưới đất" (groundwater, nước trong các tầng chứa nước Holocen và Pleistocen). Điều này quan trọng trong việc xây dựng mô hình và giải thích các quá trình dịch chuyển vật chất.

Các điều kiện biên được nêu rõ: nghiên cứu tập trung vào trầm tích biển Holocen và ảnh hưởng của chúng đến TCN Pleistocen, với vùng nghiên cứu cụ thể là CTSH, giới hạn bởi các thung lũng cắt xẻ và vùng ven biển. Luận án cũng thừa nhận các điều kiện biên về thời gian, chỉ tập trung vào quá trình rửa mặn trong Holocen (<11 nghìn năm) thay vì toàn bộ Quaternary.

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 triết lý nghiên cứu thực chứng hậu kỳ (Post-positivism) hoặc chủ nghĩa hiện thực phê phán (Critical Realism). Nó tìm cách khám phá các cơ chế và mối quan hệ nhân quả tiềm ẩn trong các quá trình thủy địa chất thông qua các phương pháp định lượng nghiêm ngặt, nhưng cũng thừa nhận tính phức tạp và đa chiều của thực tế, yêu cầu sự giải thích và tổng hợp.

Nghiên cứu sử dụng thiết kế phương pháp hỗn hợp (Mixed Methods) với sự kết hợp cụ thể của điều tra thực địa, thí nghiệm trong phòng và mô hình hóa số. Điều này được chứng minh bằng việc tổng hợp dữ liệu địa vật lý (ĐVLLK, TEM), hóa học (phân tích NLR và NDĐ), địa chất (phân tích độ hạt, khoáng vật sét) và các thí nghiệm thủy lực (hệ số thấm, khuếch tán) để cung cấp một bức tranh toàn diện và đa chiều về quá trình rửa mặn. Ví dụ, ĐVLLK và TEM cung cấp thông tin không gian về phân bố độ mặn, trong khi thí nghiệm và phân tích hóa học cung cấp dữ liệu về cơ chế dịch chuyển vật chất ở quy mô vi mô.

Thiết kế đa cấp (Multi-level design) được áp dụng khi nghiên cứu tương tác giữa nước lỗ rỗng (cấp vi mô trong trầm tích hạt mịn) và nước dưới đất (cấp vĩ mô trong tầng chứa nước) trong một khung thời gian địa chất rộng (Pleistocen-Holocen).

Kích thước mẫu và tiêu chí lựa chọn:

  • ĐVLLK: 38 lỗ khoan từ mạng quan trắc Quốc gia, cung cấp thông tin về độ dẫn điện của tầng và độ phát xạ gamma tự nhiên.
  • TEM: 210 điểm đo dọc theo 4 tuyến nghiên cứu, tạo ra các mặt cắt điện trở suất để khoanh vùng mặn nhạt.
  • Khoan lấy mẫu nguyên dạng: 2 lỗ khoan (Q87-VA, Q88-VA) tại Hà Nam, độ sâu 55-57m. Vị trí này được chọn dựa trên sự tương phản rõ rệt về độ mặn của NLR và NDĐ giữa LK Q88 (nước lợ/mặn) và LK Q87 (nước nhạt), phục vụ việc nghiên cứu nguyên nhân sự khác biệt.
  • Mẫu NLR hóa học: 40 mẫu chiết từ trầm tích nguyên dạng (tần suất 2 mét khoan/1 mẫu).
  • Mẫu thành phần độ hạt: 16 mẫu.
  • Mẫu thành phần khoáng vật sét: 8 mẫu.
  • Mẫu thí nghiệm cột thấm: 6 mẫu.
  • Mẫu thí nghiệm khuếch tán: 6 mẫu.
  • Mẫu TPHH NDĐ tầng Pleistocen: 10 mẫu từ năm 2011.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu nguyên dạng được thực hiện bằng phương pháp khoan đập (máy khoan XJ 100A) và lấy mẫu liên tục bằng ống mẫu inox dài 3m, với tỉ lệ thu hồi mẫu cao (60-95%). Các mẫu được bịt kín cẩn thận. Việc chiết NLR được thực hiện bằng thiết bị tạo áp suất cao (2-8 atm tùy thuộc loại trầm tích), đảm bảo thu được NLR từ trầm tích hạt mịn.

Giao thức thu thập dữ liệu rất chi tiết:

  • Phân tích hóa học NLR: Độ dẫn điện (EC), pH, HCO3-, Fe2+ được đo ngay tại chỗ để đảm bảo độ chính xác. Các chỉ tiêu khác (ion chính, TDS, đồng vị bền 18O, 2H) được phân tích tại Cục Địa chất Đan Mạch (GEUS-Đan Mạch) vào năm 2012. Mẫu cation được axit hóa đến pH<2 và bảo quản lạnh.
  • Phân tích độ hạt: Dùng phương pháp rây cho cát và tỷ trọng cho bột sét, phân chia thành 6 nhóm cỡ hạt theo USGS.
  • Phân tích khoáng vật sét: Sử dụng phương pháp nhiễu xạ tia X (XRD) cho mẫu định hướng ngẫu nhiên và mẫu định hướng (<2µm). Các mẫu định hướng được xử lý qua nhiều điều kiện (bão hòa Mg/khô, bão hòa Mg/Glycerol, bão hòa K/khô, nung 300oC) để nhận diện illit, smectit, vermiculit, chlorit dựa trên vị trí các peak (ví dụ, smectit từ 17Å/18Å chuyển về 10Å sau khi nung K-300oC).
  • Thí nghiệm cột thấm: Xác định hệ số thấm K bằng thiết bị cột thấm dài 7cm, đường kính 3cm, duy trì áp suất thành để đảm bảo dòng chảy thẳng đứng, đo độ dẫn điện đầu vào/ra.
  • Thí nghiệm khuếch tán: Xác định hệ số khuếch tán phân tử hữu hiệu De bằng hệ thống 2 bình thủy tinh (600ml) nối qua mẫu trầm tích dày 1cm. Bình nhận chứa nước deion (EC~4-10S/cm) đặt trên, bình cho chứa nước muối NaCl (EC~51000S/cm) đặt dưới để loại trừ phân dị trọng lực. Độ dẫn điện bình nhận được quan trắc hàng ngày trong 50 ngày.

Phương pháp tam giác (Triangulation) được áp dụng rộng rãi:

  • Tam giác dữ liệu: Độ mặn của NLR được đánh giá từ TPHH, ĐVLLK (độ dẫn điện tầng), và TEM (điện trở suất).
  • Tam giác phương pháp: Kết hợp điều tra thực địa, thí nghiệm phòng và mô hình số.
  • Tam giác lý thuyết: Sử dụng các lý thuyết về advection, diffusion, density flow để giải thích các cơ chế rửa mặn.

Về tính hợp lệ (Validity) và độ tin cậy (Reliability), thí nghiệm khuếch tán ghi nhận "độ biến đổi độ dẫn điện tăng tuyến tính theo thời gian", cho phép áp dụng định luật Fick một cách tin cậy. Các phân tích tại GEUS-Đan Mạch, một cơ quan uy tín quốc tế, đảm bảo độ chính xác của dữ liệu. Độ tin cậy trong diễn giải địa chất được củng cố bằng việc tham chiếu các tài liệu định tuổi 14C của Tanabe và cộng sự (2006) và Doan Dinh Lam & Boyd WE (2003).

Data và phân tích

Đặc điểm mẫu: 2 lỗ khoan Q87-VA và Q88-VA được chọn do vị trí gần nhau (~2km) nhưng có sự tương phản rõ rệt về độ mặn của NLR và NDĐ. Q88-VA cho thấy NLR và NDĐ lợ/mặn, trong khi Q87-VA cho thấy nước nhạt. Dữ liệu địa chất Đệ Tứ cung cấp thông tin về các hệ tầng Lệ Chi (Q11lc), Hà Nội (Q12-3hn), Vĩnh Phúc (Q13vp), Hải Hưng (Q21-2hh), Thái Bình (Q23tb) và Đệ Tứ không phân chia (Q), là cơ sở để giải thích sự phân bố trầm tích và quá trình rửa mặn.

Các kỹ thuật phân tích tiên tiến bao gồm:

  • Tính toán hệ số Rayleigh (Ra): Dựa trên hệ số thấm, khuếch tán, độ rỗng, khối lượng riêng của nước mặn và nước nhạt, độ nhớt động học, và chiều dày tầng.
  • Mô hình hóa số: Sử dụng phần mềm SEAWAT, một mô hình dịch chuyển vật chất với khối lượng riêng chất lỏng thay đổi.
    • Mô hình 1 chiều: Mô phỏng tốc độ rửa mặn của NLR trong tầng sét biển theo thời gian, giúp xác định ảnh hưởng của thời gian rửa mặn. Kết quả mô hình 1 chiều chỉ ra rằng NLR trong trầm tích biển Pleistocen muộn đã rửa mặn hoàn toàn, trong khi NLR trong tầng sét biển Holocen vẫn chứa NMTD.
    • Mô hình 2 chiều: Mô phỏng lịch sử tiến hóa trầm tích Holocen trong mối quan hệ với dao động mực nước biển và quá trình rửa mặn NLR, được xây dựng cho tuyến nghiên cứu 3 từ Hà Nội đến Nam Định.
  • Kiểm tra độ bền vững (Robustness checks): Kết quả mô hình rửa mặn NLR thời kỳ hiện tại được đối sánh với các kết quả quan trắc thực địa để đánh giá độ chính xác của mô hình.

Thông số thống kê như effect sizes và confidence intervals không được trình bày cụ thể trong phần tóm tắt này, nhưng chúng là thành phần không thể thiếu trong báo cáo đầy đủ của các phân tích định lượng để đánh giá ý nghĩa thống kê của các phát hiện.

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

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

Luận án đã đạt được một số phát hiện then chốt với bằng chứng cụ thể từ dữ liệu:

  1. Sự tồn tại của Nước Mặn Tàn Dư Holocen: Các trầm tích biển Holocen trong CTSH vẫn còn chứa NMTD đáng kể, với độ tổng khoáng hóa (TDS) > 3g/l, phân bố rộng rãi ở hai thung lũng cắt xẻ (dọc theo mặt cắt nghiên cứu 3 và 4) và các khu vực ven biển. "Vùng nƣớc mặn (TDS > 3g/l) chiếm phần lớn diện tích, phân bố tập trung ở khu vực 2 thung lũng cắt xẻ và ở vùng ven biển CTSH." Ngược lại, nước nhạt (TDS < 1g/l) chỉ tập trung thành khoảnh nhỏ ở phía nam Hà Nội. Điều này được chứng minh qua kết quả phân tích TPHH của 40 mẫu NLR và dữ liệu điện trở suất từ 38 lỗ khoan ĐVLLK và 210 điểm TEM.
  2. Phân định Cơ chế Rửa Mặn: Cơ chế rửa mặn của NLR được xác định dựa trên thành phần thạch học và hệ số Rayleigh (Ra). "Nước lỗ rỗng chứa trong tầng sét biển Holocen (với hệ số Rayleigh Ra≤40) bị rửa mặn theo cơ chế khuếch tán phân tử, nước lỗ rỗng chứa trong các trầm tích cát mịn pha sét hoặc sét pha cát (với hệ số Rayleigh Ra >40) bị rửa mặn theo cơ chế dịch chuyển vật chất do phân dị trọng lực." Phát hiện này được hỗ trợ bởi các thí nghiệm xác định hệ số thấm (K) và hệ số khuếch tán (De) từ 6 mẫu trầm tích đại diện.
  3. Tác động của Biển Tiến Flandrian đến TCN Pleistocen: Quá trình biển tiến Flandrian và sự tồn tại của tầng sét biển Holocen đã ảnh hưởng đáng kể đến phân bố độ mặn của TCN Pleistocen bên dưới. "Bản đồ phân bố mặn nhạt NDĐ tầng Pleistocen (Hình 5.1) cho thấy vị trí của các khu vực phân bố NDĐ mặn/lợ nằm tương đối trùng khớp với vị trí các thung lũng cắt xẻ." Điều này gợi ý một cơ chế XNM cổ từ các tầng biển Holocen vào TCN Pleistocen, đặc biệt là ở những khu vực có sự bào mòn trầm tích Pleistocen muộn do quá trình hình thành thung lũng cắt xẻ trong thời kỳ băng hà Wurm 2.
  4. Sự khác biệt về thời gian rửa mặn: Kết quả mô hình 1 chiều mô phỏng quá trình rửa mặn theo thời gian cho thấy NLR trong các trầm tích biển tuổi Pleistocen muộn đã bị rửa mặn hoàn toàn (trải qua thời gian rửa mặn dài hơn 60 nghìn năm), trong khi NLR trong tầng sét biển Holocen vẫn còn chứa NMTD do thời gian rửa mặn ngắn hơn (<11 nghìn năm).

Kết quả này so sánh với các nghiên cứu trước đây ở Việt Nam, vốn thường tập trung vào XNM hiện đại. Phát hiện này là phản trực giác (counter-intuitive) đối với một vùng châu thổ nhiệt đới ẩm, nơi lượng mưa dồi dào có thể được kỳ vọng sẽ rửa mặn nhanh chóng. Tuy nhiên, bằng chứng từ dữ liệu (ví dụ, biểu đồ Piper của NLR tại LK Q88-VA và Q87-VA) và mô hình hóa chỉ ra rằng các yếu tố địa chất lịch sử và thành phần thạch học đã trì hoãn đáng kể quá trình rửa mặn ở một số khu vực.

Implications đa chiều

  • Advances lý thuyết: Luận án đóng góp vào lý thuyết về dịch chuyển vật chất trong môi trường lỗ hổng bằng cách định lượng hóa các cơ chế khuếch tán và phân dị trọng lực trong điều kiện môi trường tự nhiên phức tạp của một châu thổ. Nó mở rộng khung lý thuyết về XNM cổ, làm rõ hơn mối quan hệ giữa biến đổi mực nước biển, trầm tích học và thủy địa chất.
  • Đổi mới phương pháp: Tổ hợp các phương pháp nghiên cứu (ĐVLLK, TEM, thí nghiệm thực địa, thí nghiệm phòng thí nghiệm chi tiết tại GEUS-Đan Mạch, mô hình SEAWAT 1D và 2D) tạo thành một khuôn khổ có thể áp dụng cho các hệ thống tầng chứa nước ven biển và châu thổ khác trên thế giới, đặc biệt là những nơi có lịch sử địa chất Quaternary tương tự.
  • Ứng dụng thực tiễn: Các phát hiện cung cấp thông tin quan trọng cho việc định hướng thăm dò, khai thác NDĐ. Việc xác định rõ các khu vực chứa NMTD Holocen và cơ chế rửa mặn giúp tối ưu hóa vị trí khoan giếng nước sạch và tránh các khu vực có nguy cơ nhiễm mặn cao. Ví dụ, khuyến nghị tập trung khai thác ở những khu vực không chịu ảnh hưởng của hai thung lũng cắt xẻ.
  • Khuyến nghị chính sách: Luận án cung cấp bằng chứng khoa học vững chắc để các nhà hoạch định chính sách có thể xây dựng các quy định quản lý tài nguyên nước hiệu quả hơn, bao gồm việc lập bản đồ phân vùng rủi ro nhiễm mặn chi tiết và phát triển các chiến lược bảo vệ tầng chứa nước Pleistocen, vốn là nguồn nước quan trọng cho dân cư với mật độ cao tại CTSH.
  • Điều kiện tổng quát hóa: Các kết quả nghiên cứu có thể được tổng quát hóa cho các châu thổ sông khác trên thế giới, đặc biệt là các châu thổ ở vùng nhiệt đới và cận nhiệt đới, nơi trải qua các dao động mực nước biển lớn trong kỷ Đệ Tứ và có sự hiện diện của các trầm tích biển hạt mịn Holocen.

Limitations và Future Research

Luận án trung thực thừa nhận một số hạn chế cụ thể. Thứ nhất, dù ban đầu có phạm vi nghiên cứu là trầm tích Đệ Tứ, luận án đã thu hẹp trọng tâm vào trầm tích biển Holocen do các phát hiện ban đầu cho thấy trầm tích Pleistocen đã được rửa mặn hoàn toàn trong thời gian dài (trên 60 nghìn năm), dẫn đến việc bỏ qua một số khía cạnh chi tiết của quá trình này. Thứ hai, các mô hình 1 chiều và 2 chiều, dù tiên tiến, vẫn phải dựa trên một số giả định về các thông số đầu vào và điều kiện biên, có thể ảnh hưởng đến mức độ chính xác của kết quả mô phỏng. Ví dụ, "Các thông số đầu vào mô hình 1 chiều khảo sát tốc độ rửa mặn của NLR" và "Các thông số tổng hợp của mô hình 2 chiều" dù được xác định kỹ lưỡng nhưng vẫn là những giá trị trung bình hoặc ước tính. Thứ ba, độ phân giải về không gian và thời gian của một số dữ liệu lịch sử địa chất có thể hạn chế khả năng tái tạo hoàn hảo các quá trình xảy ra trong hàng nghìn năm.

Các điều kiện biên của nghiên cứu được xác định rõ: giới hạn địa lý trong khu vực CTSH, tập trung vào trầm tích biển Holocen và tầng chứa nước Pleistocen. Điều này có nghĩa là các kết luận về cơ chế rửa mặn và ảnh hưởng của chúng có thể không trực tiếp áp dụng cho các môi trường địa chất hoặc khí hậu khác.

Dựa trên những hạn chế này, một chương trình nghiên cứu tương lai 4-5 hướng cụ thể được đề xuất:

  1. Nghiên cứu đồng vị chuyên sâu: Ngoài 18O và 2H, việc sử dụng các đồng vị phóng xạ khác (ví dụ 14C) và đồng vị của các nguyên tố hòa tan (ví dụ Sr, B) có thể cung cấp dấu vết chính xác hơn về nguồn gốc nước, quá trình trộn lẫn và thời gian lưu giữ của NLR.
  2. Mô hình hóa 3D động lực: Phát triển mô hình 3 chiều về dòng chảy và dịch chuyển vật chất có tính đến khối lượng riêng thay đổi (coupled flow-transport) sẽ cải thiện đáng kể độ chính xác trong việc dự báo phân bố mặn nhạt theo không gian và thời gian.
  3. Mở rộng phạm vi địa lý: Áp dụng khuôn khổ phương pháp luận tích hợp này cho các châu thổ khác ở Việt Nam (ví dụ: châu thổ sông Cửu Long) hoặc các châu thổ quốc tế với điều kiện địa chất và thủy địa chất tương tự để kiểm tra tính tổng quát của các phát hiện.
  4. Đánh giá tác động biến đổi khí hậu: Nghiên cứu ảnh hưởng của mực nước biển dâng và thay đổi lượng mưa do biến đổi khí hậu đến tốc độ rửa mặn NLR và sự ổn định của TCN trong tương lai.
  5. Nghiên cứu phản ứng hóa học-sinh học: Điều tra các phản ứng hóa học và sinh học diễn ra trong trầm tích hạt mịn và tác động của chúng đến thành phần NLR trong quá trình rửa mặn.

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

Luận án có tiềm năng tạo ra tác động và ảnh hưởng sâu rộng trên nhiều lĩnh vực:

  • Tác động học thuật: Nghiên cứu này có vị thế trở thành một công trình nền tảng cho lĩnh vực thủy địa chất cổ ở Việt Nam, đặc biệt là trong các nghiên cứu về hệ thống châu thổ. Với tính tiên phong trong việc làm rõ cơ chế rửa mặn và nguồn gốc XNM cổ, luận án dự kiến sẽ được trích dẫn rộng rãi bởi các nhà khoa học nghiên cứu về thủy địa chất ven biển, biến đổi mực nước biển Đệ Tứ và dịch chuyển vật chất trong môi trường lỗ hổng. Ước tính có khả năng đạt được 50-100 trích dẫn trong vòng 5-10 năm tới.
  • Chuyển đổi công nghiệp: Các phát hiện sẽ cung cấp dữ liệu và mô hình quan trọng cho ngành công nghiệp thăm dò và khai thác nước dưới đất. Việc hiểu rõ hơn về phân bố NMTD Holocen và cơ chế rửa mặn sẽ giúp các công ty tài nguyên nước tối ưu hóa chiến lược khoan và quản lý giếng, giảm rủi ro gặp phải nước mặn và tăng hiệu quả khai thác nước ngọt. Cụ thể, nó có thể ảnh hưởng đến các công ty cấp nước địa phương và các nhà phát triển bất động sản trong việc lập kế hoạch nguồn nước cho các khu đô thị mới.
  • Ảnh hưởng chính sách: Luận án cung cấp bằng chứng khoa học cần thiết cho các cơ quan chính phủ (ví dụ: Bộ Tài nguyên và Môi trường, Sở Tài nguyên và Môi trường các tỉnh) để xây dựng và thực thi các chính sách quản lý tài nguyên nước hiệu quả hơn. Các khuyến nghị về phân vùng bảo vệ tầng chứa nước, kiểm soát khai thác và chiến lược thích ứng với nhiễm mặn sẽ có đường lối triển khai cụ thể, giúp bảo vệ nguồn nước ngầm chiến lược cho một khu vực có mật độ dân số cao.
  • Lợi ích xã hội: Thông qua việc hỗ trợ quản lý NDĐ bền vững, luận án đóng góp vào an ninh nguồn nước ngọt cho hàng triệu người dân sinh sống tại CTSH. Điều này giúp giảm thiểu các vấn đề sức khỏe liên quan đến việc sử dụng nước nhiễm mặn và bảo vệ năng suất nông nghiệp, đặc biệt là trong bối cảnh biến đổi khí hậu. Lợi ích có thể được định lượng bằng việc bảo vệ một nguồn nước ngọt tiềm năng ước tính hàng tỷ m3 cho hàng chục triệu người.
  • Liên quan quốc tế: Các kết quả nghiên cứu về ảnh hưởng của dao động mực nước biển Đệ Tứ đến thủy địa chất deltaic và cơ chế rửa mặn có ý nghĩa toàn cầu. Chúng cung cấp các bài học quý giá cho các châu thổ khác trên thế giới, đặc biệt là ở Đông Nam Á, nơi đang đối mặt với những thách thức tương tự về XNM và khan hiếm nước ngọt.

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

  • Các nhà nghiên cứu tiến sĩ: Luận án cung cấp một khuôn mẫu toàn diện cho các nghiên cứu thủy địa chất phức tạp, đặc biệt là trong môi trường châu thổ. Nó chỉ ra các khoảng trống nghiên cứu cụ thể như cần phát triển mô hình 3D chi tiết hơn, tích hợp phân tích đồng vị chuyên sâu, và mở rộng ứng dụng ra các khu vực khác, tạo tiền đề cho các đề tài tiến sĩ tiếp theo.
  • Các học giả cấp cao: Nghiên cứu này làm phong phú thêm hiểu biết lý thuyết về dịch chuyển chất hòa tan trong môi trường lỗ rỗng, vai trò của các kiểm soát địa chất đối với hóa địa chất thủy văn, và tác động lâu dài của biến đổi mực nước biển đối với các tầng chứa nước. Nó cung cấp các trường hợp điển hình và dữ liệu mạnh mẽ để kiểm định và mở rộng các lý thuyết hiện có.
  • Bộ phận Nghiên cứu & Phát triển của ngành công nghiệp: Luận án cung cấp dữ liệu và mô hình cần thiết để tối ưu hóa việc khai thác và quản lý nước dưới đất, giúp các công ty cấp nước và khai thác tài nguyên giảm thiểu rủi ro nhiễm mặn và chi phí xử lý nước. Ví dụ, nó có thể giúp định lượng rủi ro XNM cho các dự án cơ sở hạ tầng ven biển mới.
  • Các nhà hoạch định chính sách: Nghiên cứu cung cấp cơ sở bằng chứng vững chắc cho việc xây dựng các chính sách quản lý và bảo vệ nguồn nước bền vững, đặc biệt là ở các khu vực ven biển và châu thổ đông dân cư. Các khuyến nghị chính sách về phân vùng, kiểm soát khai thác và các biện pháp giảm thiểu XNM sẽ trực tiếp hỗ trợ việc đạt được các mục tiêu an ninh nguồn nước quốc gia.

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

  1. Đóng góp lý thuyết độc đáo nhất: Đóng góp lý thuyết độc đáo nhất của luận án là việc định lượng hóa và phân biệt hai cơ chế rửa mặn chính – khuếch tán phân tử và phân dị trọng lực – trong các trầm tích biển Holocen của CTSH, dựa trên hệ số Rayleigh (Ra) và thành phần thạch học cụ thể. Luận án mở rộng lý thuyết về dịch chuyển chất hòa tan bằng cách tích hợp sâu sắc lịch sử địa chất Đệ Tứ và đặc điểm trầm tích học vào các quá trình thủy địa chất, đặc biệt làm rõ nguồn gốc NMTD Holocen như một yếu tố chính gây nhiễm mặn cho TCN Pleistocen. Điều này khác biệt với nhiều nghiên cứu trước đây chỉ xem xét khuếch tán hoặc đối lưu mà không định lượng mối quan hệ với đặc tính địa chất chi tiết.

  2. Đổi mới phương pháp luận: Đổi mới phương pháp luận nằm ở việc xây dựng một khuôn khổ tích hợp đa ngành. Luận án kết hợp đồng thời 38 điểm ĐVLLK (đo gamma tự nhiên và độ dẫn điện tầng), 210 điểm đo TEM (phương pháp trường chuyển), khoan lấy mẫu trầm tích nguyên dạng từ 2 lỗ khoan Q87-VA và Q88-VA, 40 mẫu phân tích hóa học NLR chi tiết, 16 mẫu phân tích độ hạt, 8 mẫu phân tích khoáng vật sét bằng XRD với các xử lý đặc biệt (Mg-saturated/glycerol, K-saturated/300oC) tại GEUS-Đan Mạch, cùng 6 thí nghiệm cột thấm và 6 thí nghiệm xác định hệ số khuếch tán. Dữ liệu này sau đó được tích hợp vào mô hình số SEAWAT 1 chiều và 2 chiều để mô phỏng dịch chuyển vật chất có tính đến khối lượng riêng chất lỏng thay đổi.

    • So sánh với các nghiên cứu trước đây: Phạm Quý Nhân (2000) đã xây dựng mô hình dịch chuyển chất hòa tan ở CTSH nhưng "chưa tính đến sự thay đổi khối lượng riêng chất lỏng". Luận án này cải tiến đáng kể bằng cách sử dụng SEAWAT để giải quyết vấn đề khối lượng riêng, vốn rất quan trọng trong XNM và rửa mặn. So với các nghiên cứu quốc tế như Groen và cộng sự (2000) ở Suriname hay Beekman và cộng sự (2011) ở Hà Lan, luận án này áp dụng một phạm vi rộng hơn các kỹ thuật địa vật lý (ĐVLLK, TEM) và phân tích phòng thí nghiệm phức tạp cho trầm tích, bên cạnh mô hình hóa, cung cấp cái nhìn toàn diện hơn về một hệ thống châu thổ lớn ở vùng nhiệt đới với sự phức tạp địa chất đặc thù.

  3. Phát hiện đáng ngạc nhiên nhất: Phát hiện đáng ngạc nhiên nhất là sự tồn tại dai dẳng của nước mặn tàn dư (NMTD) với độ tổng khoáng hóa (TDS) > 3g/l trong các trầm tích biển Holocen tương đối trẻ (hình thành dưới 11 nghìn năm trước) ở CTSH, đặc biệt là ở các khu vực sâu trong lục địa. Điều này mâu thuẫn với kỳ vọng ban đầu rằng trong một môi trường nhiệt đới ẩm với lượng mưa dồi dào, các trầm tích này lẽ ra đã được rửa mặn hoàn toàn. Tuy nhiên, kết quả từ mô hình 1 chiều và dữ liệu thực địa đã chứng minh rằng thời gian rửa mặn tương đối ngắn của Holocen (<11 nghìn năm) là không đủ để rửa sạch hoàn toàn NMTD, trong khi các trầm tích biển Pleistocen cũ hơn (>60 nghìn năm) đã được rửa mặn. Phát hiện này cung cấp một lời giải thích quan trọng cho nguồn gốc nhiễm mặn NDĐ ở các khu vực nội địa của CTSH, vốn không thể được giải thích bằng XNM hiện đại.

  4. Giao thức tái tạo được cung cấp?: Có, luận án cung cấp một giao thức tái tạo chi tiết. Các chương 1 và 4 mô tả cặn kẽ các phương pháp nghiên cứu, bao gồm các thông số cụ thể của thiết bị (máy khoan XJ 100A, thiết bị chiết nước áp suất cao 2-8 atm, thiết bị cột thấm dài 7cm đường kính 3cm, hệ thống 2 bình thí nghiệm khuếch tán 600ml), quy trình lấy mẫu (ống mẫu inox, tỉ lệ thu hồi 60-95%), giao thức phòng thí nghiệm (ví dụ: các điều kiện xử lý mẫu XRD cho khoáng vật sét, thời gian quan trắc 50 ngày cho thí nghiệm khuếch tán), và các thông số đầu vào cho mô hình số (SEAWAT). Mức độ chi tiết này cho phép các nhà nghiên cứu khác có thể tái tạo hoặc kiểm định lại các phương pháp và kết quả chính của luận án.

  5. Chương trình nghiên cứu 10 năm được phác thảo?: Có, luận án đã phác thảo một chương trình nghiên cứu tương lai cụ thể, bao gồm 4-5 hướng chính như: nghiên cứu đồng vị chuyên sâu (để xác định nguồn gốc và thời gian lưu giữ nước), mô hình hóa 3D động lực (để cải thiện dự báo không gian), mở rộng phạm vi địa lý (để kiểm tra tính tổng quát của phát hiện), đánh giá tác động biến đổi khí hậu (đến quá trình rửa mặn và an ninh nguồn nước), và nghiên cứu các phản ứng hóa học-sinh học trong quá trình rửa mặn. Các hướng này đủ để định hình một chương trình nghiên cứu 10 năm, cho phép phát triển và mở rộng các phát hiện hiện tại.

Kết luận

Luận án này đã tạo ra một bước tiến quan trọng trong việc thấu hiểu cơ chế rửa mặn của nước lỗ rỗng và ảnh hưởng của nó đến chất lượng nước dưới đất trong các châu thổ sông. Có thể tóm tắt 5 đóng góp cụ thể sau:

  1. Xác định nguồn gốc nhiễm mặn: Luận án đã khẳng định rằng nước mặn tàn dư chứa trong các trầm tích biển Holocen là một nguồn gây nhiễm mặn chính cho tầng chứa nước Pleistocen ở CTSH, đặc biệt ở các khu vực sâu trong lục địa, khác biệt với XNM hiện đại.
  2. Phân định cơ chế rửa mặn: Nghiên cứu đã định lượng và phân định hai cơ chế rửa mặn chính yếu – khuếch tán phân tử (khi Ra≤40) và dịch chuyển vật chất do phân dị trọng lực (khi Ra>40) – dựa trên hệ số Rayleigh và đặc điểm thạch học của trầm tích.
  3. Đánh giá tác động lịch sử địa chất: Luận án đã định lượng ảnh hưởng của thời kỳ biển tiến Flandrian và sự hiện diện của tầng sét biển Holocen đến phân bố mặn nhạt của tầng chứa nước Pleistocen, đặc biệt ở các thung lũng cắt xẻ.
  4. Phát triển phương pháp luận tích hợp: Luận án đã xây dựng và áp dụng thành công một tổ hợp phương pháp nghiên cứu tiên tiến, bao gồm địa vật lý, hóa địa chất, thí nghiệm phòng và mô hình số (SEAWAT) để giải quyết bài toán thủy địa chất phức tạp.
  5. Cung cấp cơ sở khoa học cho quản lý tài nguyên: Các kết quả nghiên cứu cung cấp cơ sở dữ liệu và mô hình khoa học vững chắc để định hướng tìm kiếm, thăm dò, khai thác và bảo vệ nguồn tài nguyên nước dưới đất một cách hiệu quả và bền vững tại CTSH.

Nghiên cứu này đã tạo ra một sự thúc đẩy mô hình (paradigm advancement) trong thủy địa chất Việt Nam, chuyển trọng tâm từ chỉ XNM hiện đại sang một cách tiếp cận toàn diện hơn, có tính đến XNM cổ và các quá trình rửa mặn lâu dài. Điều này mở ra ít nhất 3 luồng nghiên cứu mới: nghiên cứu đồng vị chuyên sâu để truy vết nguồn nước, mô hình hóa 3D phức hợp dòng chảy-dịch chuyển vật chất, và đánh giá tác động của biến đổi khí hậu lên các quá trình rửa mặn. Với tầm quan trọng của các châu thổ trên toàn cầu, đặc biệt là trong bối cảnh biến đổi khí hậu, các phát hiện của luận án có tính liên quan quốc tế cao. Di sản của nó là những hiểu biết sâu sắc và công cụ thực tiễn, có khả năng đo lường tác động thông qua việc tăng cường an ninh nguồn nước, bảo vệ nông nghiệp và thông tin chính sách cho hàng chục triệu người dân sinh sống tại các khu vực châu thổ trên thế giới.