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 đóng vai trò sống còn trong việc đảm bảo an ninh nguồn nước trước các tác động dị thường của biến đổi khí hậu (BĐKH) và nước biển dâng (NBD). Luận án tiến sĩ chuyên ngành Địa chất thủy văn với đề tài "Nghiên cứu và dự báo ảnh hưởng của Biến đổi khí hậu và nước biển dâng đến nước dưới đất tỉnh Thái Bình" được thực hiện tại Trường Đại học Mỏ - Địa chất, đặt nền móng tiên phong trong việc lượng hóa tương tác đa chiều giữa khí tượng, hải văn, thủy thạch học và nhân sinh đối với tài nguyên nước dưới đất (NDĐ). Bối cảnh địa lý đặc thù chỉ ra: "Thái Bình được bao bọc bởi hệ thống sông, biển khép kín với đường bờ biển chạy dài trên 50 km... Các tầng chứa nước có đặc điểm thuỷ địa hoá rất phức tạp. Nước mặn và nước nhạt phân bố xen nhau không có quy luật".

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu hụt các mô hình mô phỏng 3D mật độ biến đổi (variable-density flow) có khả năng phân tách rành mạch cơ chế lan truyền mặn giữa tầng chứa nước không áp tầng nông và tầng áp lực tầng sâu. Các nghiên cứu trước đây tại Việt Nam (Lại Đức Hùng, 1996; Đặng Đình Phúc, 1997; Nguyễn Như Trung, 2007) chủ yếu tiếp cận theo phương pháp giải tích hoặc mô hình 2D đẳng mật độ đơn lẻ, chưa tích hợp đồng thời kịch bản dâng cao mực nước biển với hệ thống công trình kiểm soát lũ và đê biển.

Để lấp đầy khoảng trống tri thức này, luận án tập trung giải quyết ba câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Đặc trưng cấu trúc địa chất thủy văn và nguồn gốc cổ địa lý chi phối quy luật phân bố mặn - nhạt trong các tầng chứa nước Holocen ($qh$) và Pleistocen ($qp$) như thế nào?
  2. RQ2: Mối quan hệ thủy lực và biên truyền mặn giữa NDĐ với nước mưa, nước sông - biển diễn ra theo quy luật tuyến tính hay phi tuyến trong phạm vi không gian bao xa tính từ bờ biển?
  3. RQ3: Dưới các kịch bản BĐKH&NBD (A2, B2) đến năm 2100, diện tích và trữ lượng nước nhạt tại các tầng chứa nước biến đổi ra sao trong hai trường hợp: giữ nguyên đê biển hiện tại và nâng cấp toàn diện hệ thống đê biển?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • H1: Tầng chứa nước Holocen ($qh$) chịu tác động trực tiếp và tức thời từ bổ cập khí tượng và xâm nhập mặn biển - sông theo tương quan tuyến tính trong đới đệm ven bờ ($< 3,0\text{ km}$).
  • H2: Tầng chứa nước Pleistocen ($qp$) được cách ly tương đối bởi các lớp sét $Q_1^{2}vp_2$, sự thu hẹp diện tích nước nhạt không do NBD trực tiếp mà chịu kiểm soát bởi áp lực khai thác nhân tạo và hiện tượng truyền áp thủy lực ven duyên.

Khung lý thuyết của luận án tích hợp phương trình dòng chảy mật độ biến đổi Darcy mở rộng, định luật bảo toàn khối lượng và phương trình đối lưu - phân tán (Advection-Dispersion Equation - ADE), kết hợp lý thuyết chu kỳ trầm tích biển tiến - biển thoái Đệ Tứ. Phạm vi nghiên cứu bao quát toàn bộ tỉnh Thái Bình ($1.570\text{ km}^2$), tập trung phân tích 2 tầng chứa nước chính với chuỗi số liệu quan trắc khí tượng - thủy văn từ năm 1960 đến 2014, kết hợp chiến dịch thực nghiệm đo đạc tự động chu kỳ cao và mô phỏng dự báo động lực học đến năm 2100.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự chuyển dịch mô hình nghiên cứu xâm nhập mặn từ các giải pháp giải tích Ghyben-Herzberg cổ điển sang các mô hình số phân tán mật độ đa chiều. Sherif and Singh (1999) khi ứng dụng mô hình FEFLOW tại đồng bằng sông Nile (Ai Cập) và vịnh Bengal (Ấn Độ) đã chứng minh tính dễ bị tổn thương dị biệt: mực nước biển dâng 50 cm khiến biên mặn tại sông Nile lấn sâu 9 km, trong khi tại Bengal chỉ là 0,4 km do độ dốc thủy lực nội địa lớn. Nghiên cứu của Ghosh Bobba (2002) tại đồng bằng Godavari (Ấn Độ) bằng mô hình SUTRA 2D nhấn mạnh tương tác kép giữa NBD 50 cm và chế độ khai thác nước ngầm quá mức làm dịch chuyển liên tục ranh giới mặn - nhạt trong 20 năm. Tại châu Âu, Rasmussen et al. (2013) tại đảo Falster (Đức) và Van Baaren & De Louw tại Hà Lan chỉ ra rằng hệ thống kênh thoát nước nhân tạo và đê bao triều có thể làm suy giảm mực nước ngầm từ 1 đến 3 m, lấn át cả tác động trực tiếp của NBD trong phạm vi 10 km bờ biển. Dausman and Langevin (2005) tại hạt Broward, Florida đã chuẩn hóa việc tích hợp MODFLOW và SEAWAT để dự báo sự di chuyển ion $\text{Cl}^-$ khi nước biển dâng trên 48 cm trong thế kỷ 21.

Tại Việt Nam, các công trình của Phạm Quý Nhân, Đỗ Trọng Sự (2000) tại Nam Định, Nguyễn Sơn (2007, 2010) tại lưu vực sông Đáy và ven biển Quảng Bình, Trịnh Hoài Thu (2014) tại đồng bằng sông Hồng, và Nguyễn Văn Lâm (2015) tại dải ven biển Bắc Bộ đã từng bước ứng dụng mô hình số dòng chảy. Tuy nhiên, một cuộc tranh luận học thuật kéo dài xoay quanh hai luồng quan điểm đối lập:

  • Luồng quan điểm 1 (Biến động ngoại sinh chiếm ưu thế): Cho rằng sự suy thoái trữ lượng và mặn hóa NDĐ đồng bằng ven biển Bắc Bộ chịu điều phối áp đảo bởi nước biển dâng và sự thay đổi biên bổ cập khí hậu (Bùi Trần Vượng, 2013; Nguyễn Thị Hạ, 2016).
  • Luồng quan điểm 2 (Động lực nội sinh và nhân sinh chiếm ưu thế): Khẳng định quá trình hạ thấp mực nước và thu hẹp thấu kính nước nhạt chủ yếu do khai thác tập trung cục bộ và cấu trúc trầm tích cổ kiểm soát (Trịnh Hoài Thu, 2014; Nguyễn Như Trung, 2007).

Luận án này định vị chính xác tại giao điểm của hai luồng tranh luận, cung cấp lời giải xác đáng thông qua việc phân tách cấu trúc địa chất hai tầng riêng biệt: khẳng định tính đúng đắn của Luồng 1 đối với tầng chứa nước không áp $qh$ nông, đồng thời chứng minh Luồng 2 giữ vai trò chi phối đối với tầng chứa nước có áp $qp$ sâu. Đây là bước tiến tri thức quan trọng vượt lên trên các nghiên cứu đơn tầng trước đó.

Đó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à hoàn thiện khung lý thuyết Địa chất thủy văn động lực học và Thủy địa hóa học trầm tích Đệ Tứ thông qua 4 đóng góp bản lề:

  1. Mở rộng lý thuyết lan truyền ranh giới mặn - nhạt mật độ biến đổi: Chứng minh hiện tượng xâm nhập mặn không diễn ra đồng nhất mà bị chia cắt không gian bởi các "cửa sổ địa chất thủy văn" (hydrogeological windows) tại các lỗ khoan $LK2B$ (Việt Thuận - Vũ Thư), $LK22$ (Quỳnh Côi - Quỳnh Phụ) và $LK711$ (Đông Hoàng - Đông Hưng), nơi lớp thấm nước yếu $Q_1^{2}vp_2$ bị vát mỏng hoàn toàn, tạo dòng thấm xuyên trực tiếp nối liền hai tầng $qh$ và $qp$.
  2. Xác lập nguyên lý phân hóa tác nhân chi phối theo tầng sinh quyển: Làm sáng tỏ bản chất nguồn gốc nước nhạt tầng $qp$ phía Bắc tỉnh có tuổi trầm tích trẻ hơn ($20.000\text{ năm BP}$) so với phía Nam ($40.000 - 100.000\text{ năm BP}$). Năng lượng dòng chảy cổ Tây Bắc - Đông Nam kết hợp đứt gãy kiến tạo Vĩnh Ninh đã đẩy lùi nước biển cổ và đưa khối nước nhạt từ Neogen lên, tạo nên dị thường thủy địa hóa độc nhất vô nhị.
  3. Định lượng hóa cơ chế truyền áp thủy lực (Hydraulic Pressure Transmission): Khẳng định dao động triều không đưa trực tiếp ion muối vào tầng $qp$ ở khoảng cách sâu, mà tác động qua sóng áp lực truyền trong môi trường lỗ hổng bão hòa trong phạm vi $1,5 - 2,0\text{ km}$ tính từ đường bờ.

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

Khung phân tích của luận án tích hợp liên ngành 3 trụ cột lý thuyết:

  • Lý thuyết dòng chảy biến đổi mật độ Darcy-SEAWAT: Mô hình hóa toán học phương trình trạng thái chất lưu phụ thuộc nồng độ muối: $$\rho = \rho_f + E \cdot C$$ Trong đó $\rho_f$ là khối lượng riêng nước nhạt ($1000\text{ kg/m}^3$), $C$ là nồng độ $\text{TDS}$, $E \approx 0,7143$.
  • Lý thuyết trầm tích học biến cố Đệ Tứ: Giải thích tính phân dị tướng thạch học aluvi, vũng vịnh, châu thổ ảnh hưởng đến hệ số dẫn nước $Km$ và hệ số nhả nước trọng lực $\mu$.
  • Khung phân tích tương quan tuyến tính biên động lực: Thiết lập phương trình hồi quy thực nghiệm đánh giá mức độ tương tác giữa biên triều Hòn Dấu, mực nước sông Hồng, sông Trà Lý với mực nước các lỗ khoan quan trắc.

Điều kiện biên giới hạn (boundary conditions) được kiểm soát chặt chẽ: Biên trên là biên bổ cập lượng mưa biến đổi theo chu kỳ thực nghiệm; biên đáy là tầng sét Neogen xem như cách nước tương đối; biên biển là biên cột nước không đổi ($CHD$) biến thiên theo kịch bản dâng cao mực nước biển và triều ký; biên khai thác là các giếng khoan công nghiệp và sinh hoạt với lưu lượng trích xuất xác định.

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 hiện thực phê phán (Critical Realism) và phương pháp luận thực chứng (Positivism) nghiêm ngặt trong khoa học Trái đất. Thiết kế nghiên cứu đa cấp độ (multi-level design) kết hợp đồng thời 5 bậc tiếp cận: tiếp cận lý thuyết cơ bản, tiếp cận kế thừa lịch sử địa chất, tiếp cận khảo sát thực địa hiện trường, tiếp cận phân tích thí nghiệm chuyên sâu và tiếp cận mô hình hóa số trị 3 chiều.

Quy trình nghiên cứu rigorous

Chiến dịch thu thập và xử lý số liệu đạt độ tin cậy và chuẩn mực học thuật cao:

  • Mạng lưới mẫu thực địa quy mô lớn: Khảo sát chi tiết tại hơn $200$ điểm trên toàn tỉnh Thái Bình, đo đạc trực tiếp nồng độ $\text{TDS}$ bằng thiết bị quan trắc tự động và lấy $60$ mẫu phân tích thành phần hóa học toàn diện trong phòng thí nghiệm tại Liên đoàn Quy hoạch và Điều tra Tài nguyên nước miền Bắc.
  • Thí nghiệm thấm hiện trường và thạch học: Triển khai $20$ điểm thí nghiệm đổ nước hố đào trên bề mặt kết hợp phân tích thành phần hạt thạch học của $35$ mẫu đất đá để xây dựng bản đồ phân vùng khả năng thấm cho tầng $qh$.
  • Quan trắc động lực học tự động liên tục: Thiết lập $20$ trạm quan trắc tự động đo mực nước và độ mặn với tần suất $1\text{ giờ/lần}$ liên tục suốt $30\text{ ngày}$ (từ ngày 11/01/2015 đến 15/02/2015) đồng bộ theo chu kỳ tuần trăng và biên độ triều tại các huyện ven biển Thái Thụy, Tiền Hải; kết hợp phân tích chuỗi số liệu tại các lỗ khoan sâu nội địa $Q155, Q156, Q158, Q159$.
  • Phân tích đồng vị phóng xạ (Isotope Hydrology): Sử dụng đồng vị phóng xạ Triti ($^3\text{H}$) đạt giá trị $11,5 \pm 1,0\text{ TU}$, xác nhận tuyệt đối nguồn gốc nước tầng $qh$ là nước trẻ khí tượng - thủy văn hiện đại.

Data và phân tích

Các thông số thủy địa chất đầu vào được hiệu chỉnh chính xác:

  • Tầng chứa nước $qh_2$: Hệ số thấm $K = 1,49\text{ m/ngày}$, hệ số nhả nước trọng lực $\mu = 0,1238$, mô đun dòng ngầm $M_d = 8,6927\text{ l/s}\cdot\text{km}^2$.
  • Tầng chứa nước $qh_1$: Chiều dày trung bình $14,2\text{ m}$, $K = 0,61 - 1,29\text{ m/ngày}$, $\mu = 0,1188$, độ khoáng hóa cực đại đạt $M = 27,42\text{ g/l}$ (loại hình $\text{Cl-Na-Mg}$).
  • Tầng chứa nước $qp$: Chiều dày trung bình $57\text{ m}$ (dao động $15 - 127\text{ m}$), $K = 22\text{ m/ngày}$, $\mu = 0,182$, hệ số nhả nước đàn hồi $\mu^* = 0,0069$.
  • Tầng chứa nước $n$: $K = 3,05 - 16,23\text{ m/ngày}$, hệ số dẫn nước $Km = 105 - 432\text{ m}^2\text{/ngày}$, $\text{TDS} = 0,38 - 20,18\text{ g/l}$.

Công cụ phần mềm chủ đạo bao gồm hệ thống GMS (Groundwater Modeling System), Visual MODFLOW kết hợp module động học SEAWAT chuyên dụng cho dòng chảy biến đổi mật độ và MT3DMS cho truyền khối chất tan.

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

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

Trích dẫn luận điểm bảo vệ cốt lõi của tác giả mang lại các chỉ số thực nghiệm đột phá: "Tầng chứa nước Holocen chịu ảnh hưởng trực tiếp bởi nước mưa và nước sông, biển trong phạm vi 1,5 đến 3,0 km so với đường bờ với quan hệ tuyến tính. Lượng mưa tăng, lượng bổ cập cho tầng chứa nước qh trung bình ước tính khoảng 0,0003 m/ng đã góp phần rửa mặn nước dưới đất với diện tích mặn thu hẹp khoảng 180 km2 so với năm 1996."

Các phát hiện cụ thể được tổng hợp với bằng chứng định lượng chính xác:

  1. Cơ chế rửa mặn tự nhiên giai đoạn 1996 - 2014: Khẳng định xu thế gia tăng lượng mưa ($0,9 - 1,0%/\text{năm}$ từ năm 1995) đóng vai trò rửa mặn tích cực cho tầng $qh$, thu hẹp $180\text{ km}^2$ diện tích nhiễm mặn so với bản đồ phân bố của Lại Đức Hùng (1996).
  2. Phân hóa tác động của BĐKH&NBD theo kịch bản phát thải A2 đến năm 2100: Trích dẫn số liệu dự báo từ luận án: "Tùy theo đặc trưng của tầng chứa nước, hiện trạng sử dụng đất và đê điều, diện tích vùng nước nhạt TCN qh tính đến năm 2100 theo kịch bản phát thải A2 bị thu hẹp khoảng 109,7 km2 trong trường hợp đê biển hiện tại và 42,9 km2 khi nâng cấp hệ thống đê biển..."
  3. Tính độc lập tương đối của tầng Pleistocen trước NBD: Tầng $qp$ hầu như không chịu ảnh hưởng trực tiếp từ nước biển dâng mặt bờ; diện tích nước nhạt bị thu hẹp khoảng $25\text{ km}^2$ tính đến năm 2100 chủ yếu do áp lực khai thác nước ngầm tập trung tại khu vực phía Bắc của tỉnh.
  4. Hiệu ứng chắn triều của hệ thống đê biển: Lượng hóa lần đầu tiên vai trò của $71,4\text{ km}$ đê biển (tuyến số 5, 6, 7, 8 có cao trình $4,6 - 5,5\text{ m}$), chứng minh việc nâng cấp đê biển giúp giảm thiểu tới $60,9%$ diện tích nước nhạt bị mất mát ($66,8\text{ km}^2$ được bảo vệ) ở tầng $qh$.
Tầng chứa nước Tác nhân chi phối chính Phạm vi ảnh hưởng ven bờ Biến động diện tích nước nhạt đến 2100 (Kịch bản A2) Cơ chế dịch chuyển ranh giới mặn - nhạt
Holocen ($qh$) Bổ cập khí tượng + Mực nước biển dâng $1,5 - 3,0\text{ km}$ (tuyến tính) • Thu hẹp $109,7\text{ km}^2$ (Đê hiện tại)
• Thu hẹp $42,9\text{ km}^2$ (Đê nâng cấp)		 Xâm nhập trực tiếp qua mặt cắt ngang ven duyên và dòng mặt sông nhiễm triều
Pleistocen ($qp$) Khai thác nhân tạo + Truyền áp thủy triều $1,5 - 2,0\text{ km}$ (truyền áp) Thu hẹp $25,0\text{ km}^2$ (Tập trung tại phía Bắc tỉnh) Dịch chuyển do phễu hạ thấp mực nước và thấm xuyên qua cửa sổ ĐCTV ($2,03\text{ m/ng}$)

Implications đa chiều

  • Về mặt học thuật: Khẳng định tính tất yếu của việc mô hình hóa phân tầng trong nghiên cứu thủy văn đồng bằng châu thổ, bác bỏ giả định xem toàn bộ hệ thống chứa nước Đệ Tứ là một thể đồng nhất.
  • Về phương pháp luận: Cung cấp quy trình tích hợp mô hình SEAWAT 3D với kịch bản đê điều công trình, có thể chuyển giao áp dụng cho toàn bộ các tỉnh duyên hải Bắc Bộ và đồng bằng sông Cửu Long.
  • Về ứng dụng thực tiễn: Xác lập bản đồ phân vùng khai thác an toàn, chỉ định chuyển dịch các bãi giếng công nghiệp tầng $qp$ xa khỏi hành lang đệm $2\text{ km}$ ven biển và các đới nhạy cảm phía Bắc.
  • Về chính sách công: Cung cấp luận cứ khoa học chuẩn xác để UBND tỉnh Thái Bình và Bộ Tài nguyên và Môi trường phê duyệt dự án nâng cấp $71,4\text{ km}$ đê biển và ban hành quy chế bảo vệ hành lang nước ngầm.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 3 giới hạn khoa học nội tại:

  1. Độ phân giải mắt lưới không gian mô hình: Do giới hạn tính toán, kích thước ô lưới chưa đủ mịn để mô phỏng chi tiết các thấu kính nước nhạt quy mô hộ gia đình ($< 50\text{ m}$) trong tầng $qh$.
  2. Tính bất định của kịch bản kinh tế - xã hội: Lưu lượng khai thác ngầm trong tương lai được giả định không đổi hoặc tăng tuyến tính cục bộ, chưa tích hợp mô hình quy hoạch đô thị hóa phi tuyến.
  3. Độ dài chuỗi quan trắc liên tục cao tần: Chuỗi đo tự động 1 giờ/lần kéo dài 1 tháng (mùa khô 2015), tuy bắt trọn chu kỳ tuần trăng nhưng cần được kiểm chứng thêm bằng chuỗi số liệu đa năm liên tục qua cả mùa mưa lũ cực đoan.

Chương trình nghiên cứu tiếp nối trong 5-10 năm tới bao gồm:

  • Tích hợp mô hình biến dạng đàn - dẻo mô phỏng hiện tượng lún sụt mặt đất (Land Subsidence) do khai thác nước ngầm quá mức kết hợp SEAWAT.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Physics-Informed Neural Networks - PINNs) để tự động hiệu chỉnh thông số thấm $K$ phi đồng nhất.
  • Mở rộng phân tích đồng vị bền ($\delta^{18}\text{O}, \delta^2\text{H}$) và đồng vị $\text{Carbon-14}$ để định tuổi chính xác tuyệt đối các pha bổ cập cổ của tầng Neogen.

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

Công trình tạo lập sức ảnh hưởng sâu rộng trên 4 phương diện:

  • Giá trị học thuật quốc tế: Đóng góp bộ dữ liệu thực địa hiếm có về tương tác nước mặt - nước ngầm tại châu thổ sông Hồng, dự kiến gia tăng trích dẫn trong các nghiên cứu tổng quan về đồng bằng Đông Nam Á.
  • Chuyển đổi thực tiễn quản trị tài nguyên: Định hình lại quy hoạch cấp nước nông thôn tỉnh Thái Bình, chuyển từ khai thác phân tán tự phát sang mạng lưới cấp nước tập trung từ nước mặt sông Hồng kết hợp giếng khai thác $qp$ an toàn.
  • Chính sách thích ứng BĐKH cấp quốc gia: Đóng góp cơ sở dữ liệu nền cho Kịch bản Biến đổi khí hậu và Nước biển dâng cập nhật của Bộ Tài nguyên và Môi trường.
  • Lợi ích kinh tế - xã hội lượng hóa: Giúp tiết kiệm hàng trăm tỷ đồng ngân sách thông qua việc tối ưu hóa vị trí tuyến đê biển và ngăn chặn nguy cơ nhiễm mặn không thể phục hồi tại các tầng chứa nước ngọt sinh hoạt của $1,8\text{ triệu}$ dân.

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

  • Nghiên cứu sinh & Nhà khoa học Địa chất thủy văn: Tiếp cận phương pháp luận xây dựng mô hình SEAWAT 3D kết nối nhiều tầng chứa nước và bộ thông số địa chất thủy văn chuẩn hóa ($K, \mu, Md$).
  • Viện nghiên cứu & Trường Đại học: Kế thừa khung lý thuyết trầm tích - thủy địa hóa Đệ Tứ đồng bằng châu thổ sông Hồng phục vụ đào tạo sau đại học.
  • Doanh nghiệp R&D & Cấp nước đô thị: Định vị chính xác tọa độ tầng nước nhạt $qp$ và $n$ để thiết kế giếng khoan công nghiệp hiệu suất cao, tránh rủi ro nhiễm mặn phễu hút.
  • Cơ quan quản lý Nhà nước (Sở TN&MT, Cục Quản lý Tài nguyên nước): Sở hữu công cụ số trực quan để cấp phép khai thác, phân vùng cấm/hạn chế khai thác NDĐ theo luật định.

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ì? Đó là việc khám phá và chứng minh quy luật phân dị tác nhân mặn hóa: Tầng nông ($qh$) chịu kiểm soát bởi biến số khí hậu ngoại sinh và ranh giới động ven duyên; tầng sâu ($qp$) chịu kiểm soát bởi cấu trúc cổ địa lý, hoạt động nhân sinh và truyền áp thủy triều, qua đó bác bỏ mô hình đơn tầng truyền thống.

2. Điểm mới phương pháp luận so với các nghiên cứu tiền nhiệm? So với nghiên cứu giải tích của Đặng Đình Phúc (1997) và mô hình 2D của Nguyễn Sơn (2007), luận án đã thiết lập mô hình 3D biến đổi mật độ SEAWAT tích hợp dữ liệu đồng vị Triti ($11,5\text{ TU}$), $200$ điểm khảo sát và $20$ điểm đo triều tự động chu kỳ cao, kết hợp định lượng kịch bản có và không có công trình đê biển $71,4\text{ km}$.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất? Diện tích nhiễm mặn tầng $qh$ năm 2014 thu hẹp $180\text{ km}^2$ so với năm 1996 nhờ lượng mưa gia tăng trung bình đạt $0,0003\text{ m/ngày}$ giúp rửa mặn tự nhiên, chứng minh BĐKH không chỉ mang lại tác động tiêu cực mà còn có hiệu ứng tích cực cục bộ theo mùa.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) không? Có. Toàn bộ thông số thủy thạch học ($K = 1,49\text{ m/ng}$ cho $qh_2$, $K = 22\text{ m/ng}$ cho $qp$), tọa độ lỗ khoan kiểm chứng ($LK2B, LK22, Q155-159$), điều kiện biên SEAWAT và phương trình trạng thái đều được công khai chi tiết trong các phụ lục kỹ thuật.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Tập trung vào phát triển mô hình ghép bậc cao (Coupled Flow-Transport-Geomechanics) giải quyết đồng thời bài toán xâm nhập mặn và lún chìm mặt đất dưới áp lực kép của đô thị hóa và nước biển dâng cực hạn.

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc mục tiêu khoa học với 5 đóng góp cốt tử:

  1. Xác lập bản đồ phân vùng cấu trúc địa chất thủy văn chi tiết tỉnh Thái Bình, làm rõ vai trò của 3 tầng chứa nước ($qh, qp, n$) và các cửa sổ địa chất thủy văn kết nối.
  2. Định lượng hóa quy luật tương quan tuyến tính của nước sông - biển với tầng $qh$ trong đới $1,5 - 3,0\text{ km}$ và quy luật truyền áp thủy triều với tầng $qp$ trong đới $1,5 - 2,0\text{ km}$.
  3. Chứng minh quy luật rửa mặn tự nhiên $180\text{ km}^2$ giai đoạn 1996 - 2014 dưới ảnh hưởng của xu thế tăng lượng bổ cập khí tượng.
  4. Dự báo chính xác suy giảm diện tích nước nhạt đến năm 2100 theo kịch bản A2: tầng $qh$ mất $109,7\text{ km}^2$ (giảm xuống còn $42,9\text{ km}^2$ nếu nâng cấp đê), tầng $qp$ mất $25\text{ km}^2$ do khai thác.
  5. Đề xuất quy hoạch khai thác không gian hợp lý và giải pháp công trình đê điều chắn mặn bền vững, đóng góp trực tiếp vào an ninh nguồn nước vùng kinh tế trọng điểm Bắc Bộ.