Tổng quan về luận án

Bán đảo Cà Mau (BĐCM) có diện tích tự nhiên $16.940\text{ km}^2$, bao gồm địa giới hành chính của 6 tỉnh, thành phố: Cần Thơ, Hậu Giang, Sóc Trăng, Bạc Liêu, Cà Mau và một phần tỉnh Kiên Giang. Với mật độ dân số trung bình đạt $377\text{ người/km}^2$ (quy mô dân số hơn $6,49\text{ triệu người}$), khu vực giữ vị trí hạt nhân trong chiến lược phát triển kinh tế nông - thủy sản của Đồng bằng sông Cửu Long (ĐBSCL). Tuy nhiên, mạng lưới kênh rạch tự nhiên thường xuyên bị xâm nhập mặn nghiêm trọng trong mùa khô kết hợp với tình trạng ô nhiễm nguồn nước mặt cục bộ, khiến nguồn nước dưới đất (NDĐ) trở thành nguồn tài nguyên duy nhất bảo đảm an ninh cấp nước sinh hoạt và sản xuất. Tình trạng khai thác nước ngầm quy mô lớn, phân tán và thiếu kiểm soát diễn ra đồng thời với hiện tượng biến đổi khí hậu (BĐKH) và nước biển dâng (NBD) đang đẩy hệ thống thủy văn ngầm vào nguy cơ cạn kiệt động lực và mặn hóa nghiêm trọng.

Khoảng trống nghiên cứu (research gap) then chốt trước đây nằm ở sự thiếu hụt các mô hình tích hợp đa cấp độ có khả năng định lượng tách biệt và đồng thời: tác động đơn lẻ của hoạt động khai thác nhân sinh so với tác động cộng hưởng của BĐKH tới từng tầng chứa nước (TCN) cụ thể; cũng như sự vắng bóng của một khung chỉ số bền vững cấp tiểu vùng để chuyển hóa các kết quả mô phỏng thủy địa động lực học thành công cụ quản trị thực tiễn. Nhằm giải quyết triệt để vấn đề này, luận án xác lập hệ thống câu hỏi và giả thuyết nghiên cứu:

  • Câu hỏi nghiên cứu 1 (RQ1): Mức độ suy giảm mực nước áp lực và lượng tích trữ tĩnh/động của 6 TCN chính biến động như thế nào dưới tác động đơn lẻ của hoạt động khai thác?
  • Câu hỏi nghiên cứu 2 (RQ2): Khi cộng hưởng với kịch bản BĐKH phát thải cao A2, lượng bổ cập khí tượng và cấu trúc biên mặn - nhạt của hệ thống TCN biến thiên ra sao theo thang thời gian 15 năm từ 2015 đến 2090?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ bền vững của tài nguyên NDĐ tại 43 huyện/thành phố thuộc BĐCM biến động như thế nào khi định lượng qua chỉ số cạn kiệt và chỉ số tổn thương xâm nhập mặn?
  • Giả thuyết nghiên cứu 1 (H1): Khai thác nước ngầm là tác nhân trực tiếp gây sụt giảm mực nước áp lực và suy thoái lượng tích trữ trong ngắn hạn, nhưng BĐKH đóng vai trò nhân tố kiểm soát chính lượng bổ cập ngấm bề mặt và làm dịch chuyển ranh giới mặn - nhạt dài hạn.
  • Giả thuyết nghiên cứu 2 (H2): Xu thế suy giảm tính bền vững của tài nguyên NDĐ không đồng nhất giữa các đơn vị hành chính mà chịu sự chi phối phân hóa mạnh mẽ bởi cấu trúc địa chất thủy văn (ĐCTV) và mật độ khai thác cục bộ.

Khung lý thuyết của công trình được định hình trên sự giao thoa giữa lý thuyết cân bằng ẩm khí quyển - đất - thảm phủ (WetSpass framework), lý thuyết dòng thấm thủy địa động lực học trong môi trường tầng chứa nước xốp bão hòa đa tầng (Darcy-Continuity Theory) và khung đánh giá chỉ số bền vững tài nguyên nước ngầm của UNESCO/IAEA/IAH. Luận án mang lại đóng góp đột phá khi định lượng chính xác sự suy giảm tích trữ của toàn bồn trũng qua 6 lát cắt thời gian dự báo ($2015, 2030, 2045, 2060, 2075, 2090$), cung cấp luận cứ khoa học chuẩn xác để hoạch định lại toàn bộ mạng lưới khai thác nước ngầm trên phạm vi $16.940\text{ km}^2$.

Literature Review và Positioning

Nghiên cứu tương tác giữa BĐKH, khai thác nhân sinh và động thái nước dưới đất là chủ đề trọng điểm của thủy địa chất học đương đại. Luận án đã tổng hợp có hệ thống các trường phái nghiên cứu quốc tế tiêu biểu:

  • Trường phái mô hình hóa lượng bổ cập liên kết dòng chảy ngầm: Riasat Ali et al. (2012) khi nghiên cứu vùng Tây Nam nước Úc đã áp dụng mô hình liên kết VFM (Vertical Flux Model) và PRAMS (Perth Regional Aquifer Modelling System), kết luận rằng các TCN có áp chịu ảnh hưởng rất chậm từ BĐKH do độ trễ thấm xuyên qua các lớp sét cách nước. Ngược lại, Woldeamlak et al. (2007) ứng dụng WetSpass kết hợp MODFLOW tại lưu vực Grote-Nete (Bỉ) chứng minh kịch bản khô hạn làm sụt giảm mực nước ngầm trung bình $0,5\text{ m}$ (cực đại $3,1\text{ m}$), đe dọa trực tiếp đến thảm thực vật mặt đất. Scibek & Allen (2006) tại Grand Forks (Canada) sử dụng mô hình giáng cấp thống kê SDSM kết hợp mô hình HELP và MODFLOW, chỉ ra rằng tương tác thủy lực giữa sông và TCN không áp chi phối mạnh hơn lượng bổ cập ngấm trực tiếp.
  • Trường phái động lực học xâm nhập mặn ven bờ: Sherif & Singh (1999) áp dụng định luật thủy tĩnh Ghyben-Herzberg tại đồng bằng sông Nile (Ai Cập) và vịnh Bengal (Ấn Độ), chứng minh khi mực nước biển dâng $50\text{ cm}$, ranh giới mặn tại đồng bằng sông Nile tiến sâu vào nội địa tới $9\text{ km}$, trong khi tại vịnh Bengal chỉ tăng thêm $0,4\text{ km}$ do độ dốc thủy lực và cấu trúc địa chất khác biệt. Ghosh Bobba (2002) sử dụng mô hình SUTRA mô phỏng 20 năm tại đồng bằng Godavari (Ấn Độ), khẳng định duy trì cường độ khai thác hiện hữu sẽ đẩy nêm mặn tiến sâu không thể đảo ngược. Sun Woo Chang et al. (2011) sử dụng MODFLOW và SEAWAT để phân tích, chỉ ra trong trạng thái dòng chảy không ổn định dài hạn, cả tầng có áp lẫn không áp đều bị xâm nhập mặn nghiêm trọng dưới tác động của NBD. Antoine Armandine Les Landes et al. (2013) phân tách định lượng tại vùng đất ngập nước Tây Bắc nước Pháp, cho thấy BĐKH làm suy giảm 5,3–13,6% diện tích đất ngập nước, trong khi khai thác tăng 100% chỉ làm giảm tối đa 3,7%.
Công trình tiêu biểu Khu vực nghiên cứu Phương pháp & Mô hình Kết quả chính
Riasat Ali et al. (2012) Tây Nam nước Úc VFM liên kết PRAMS TCN có áp ít nhạy cảm với BĐKH ngắn hạn; tầng nông suy giảm mạnh ở vùng khô hạn.
Sherif & Singh (1999) ĐB Sông Nile & Vịnh Bengal Giải tích Ghyben-Herzberg NBD $50\text{ cm}$ làm mặn xâm nhập $9\text{ km}$ tại Nile nhưng chỉ $0,4\text{ km}$ tại Bengal.
Woldeamlak et al. (2007) Lưu vực Grote-Nete (Bỉ) WetSpass liên kết MODFLOW BĐKH làm giảm mực nước trung bình $0,5\text{ m}$ (cực đại $3,1\text{ m}$), hủy hoại đồng cỏ.
Đoàn Văn Cánh et al. (2011) ĐB Bắc Bộ & ĐB Nam Bộ Phân tích ĐCTV & Thấm xuyên Thiết lập 5 tiêu chí và 9 chỉ số; xác định trữ lượng khai thác phụ thuộc thấm xuyên.
Luận án (Đào Hồng Hải, 2018) Bán đảo Cà Mau WetSpass - GMS MODFLOW/SEAWAT Định lượng suy giảm 6 TCN; dự báo cạn kiệt và tổn thương mặn cho 43 đơn vị hành chính.

Tại Việt Nam, các công trình của Đoàn Văn Cánh et al. (2011) và Bùi Học et al. (2010) đã tiên phong thiết lập hệ tiêu chí phân vùng khai thác bền vững và chứng minh nguồn hình thành trữ lượng khai thác chủ yếu bắt nguồn từ dòng thấm xuyên giữa các TCN. Tuy nhiên, các nghiên cứu này chủ yếu dừng lại ở việc minh chứng thủy lực trên các mặt cắt khái lược, chưa xây dựng được mô hình 3D dòng chảy kết hợp vận chuyển chất mô phỏng chi tiết 7 tầng chứa nước trong điều kiện biên biển kép (Biển Đông và Biển Tây) dưới tác động đồng thời của kịch bản phát thải cao SRES A2. Luận án định vị chính xác khoảng trống này để phát triển một nghiên cứu toàn diện, tích hợp công cụ mô phỏng toán - thủy văn với bộ chỉ số bền vững quốc tế của UNESCO.

Đó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 dòng thấm thủy địa động lực học đa tầng trong môi trường địa chất trầm tích Đệ tứ và Neogen bị chia cắt phức tạp. Bằng việc lượng hóa sự truyền áp lực và dịch chuyển ranh mặn dưới điều kiện khai thác tập trung, công trình đã chứng minh:

  1. Mệnh đề lý thuyết 1 (Proposition 1): Tốc độ suy giảm áp lực thủy tĩnh trong các TCN có áp sâu ($qp_1, n_2^2, n_2^1$) tỷ lệ thuận trực tiếp với tổng lưu lượng khai thác tập trung tại các phễu hạ thấp đô thị, nhưng quá trình mặn hóa lại diễn ra theo cơ chế kép: dòng thấm xuyên thẳng đứng qua các "cửa sổ địa chất" trầm tích biển $mQ$ bị bào mòn và sự dịch chuyển ngang của nêm mặn cổ từ đới ven biển.
  2. Mệnh đề lý thuyết 2 (Proposition 2): Hiện tượng BĐKH đóng vai trò nhân tố tái phân bố lại lượng bổ cập ngấm bề mặt theo không gian, làm giảm đáng kể lượng nước bổ cập tự nhiên trong mùa khô, gián tiếp triệt tiêu gradien áp lực tự nhiên hướng ra biển và gia tốc tốc độ xâm nhập mặn vào các thấu kính nước ngọt.

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

Khung phân tích của luận án là sự tích hợp chặt chẽ giữa 3 nền tảng lý thuyết và công cụ mô hình hóa:

                                   (SEAWAT/MT3DMS: Vận chuyển Muối)
                             [Bộ chỉ số Bền vững UNESCO/IAEA/IAH]
  1. Lý thuyết cân bằng nước bề mặt (WetSpass): Xác định lượng ngấm bổ cập theo không gian thực dựa trên sự phân bố lượng mưa, nhiệt độ, bốc thoát hơi tiềm năng, loại đất và thảm phủ.
  2. Lý thuyết dòng thấm không ổn định 3D (MODFLOW): Thiết lập phương trình vi phân đạo hàm riêng mô tả dòng chảy qua môi trường rỗng dị hướng, đa lớp thấm và bán thấm.
  3. Lý thuyết đánh giá chỉ số bền vững (UNESCO/IAEA/IAH Framework): Lựa chọn và chuẩn hóa 2 chỉ số cốt lõi:
    • Chỉ số cạn kiệt tài nguyên NDĐ ($I_d$): Tỷ lệ giữa biến thiên lưu lượng khai thác, lượng tích trữ và tốc độ sụt giảm mực nước áp lực.
    • Chỉ số tổn thương do xâm nhập mặn ($I_s$): Tỷ lệ phần trăm diện tích bị nhiễm mặn ($\text{TDS} > 1\text{ g/L}$) so với tổng diện tích phân bố tự nhiên của từng TCN trên phạm vi từng đơn vị hành chính.
    • Điều kiện biên (Boundary conditions): Biên mực nước xác định trước (Constant Head) mô phỏng chế độ triều Biển Đông - Biển Tây; biên dòng chảy tổng hợp (General Head Boundary) mô phỏng sông Hậu; biên không dòng chảy (No-flow boundary) tại ranh giới kết thúc địa tầng hoặc đứt gãy kiến tạo.

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

Thiết kế nghiên cứu

Nghiên cứu quán triệt thế giới quan thực chứng (positivism) với phương pháp tiếp cận định lượng diễn dịch (quantitative deductive approach). Quy trình tích hợp đa phương pháp bao gồm: chuẩn hóa số liệu ĐCTV thực địa $\rightarrow$ chi tiết hóa kịch bản BĐKH $\rightarrow$ mô phỏng bổ cập ẩm bề mặt $\rightarrow$ xây dựng mô hình dòng chảy ngầm 3 chiều $\rightarrow$ mô phỏng dịch chuyển khối chất hòa tan $\rightarrow$ tính toán chỉ số bền vững không gian.

Thiết kế nghiên cứu đa tầng bao phủ toàn bộ 7 phân vị địa chất thủy văn theo chiều sâu:

  1. Tầng chứa nước Holocen ($qh$)
  2. Tầng chứa nước Pleistocen trên ($qp_3$)
  3. Tầng chứa nước Pleistocen giữa - trên ($qp_{2-3}$)
  4. Tầng chứa nước Pleistocen dưới ($qp_1$)
  5. Tầng chứa nước Pliocen giữa ($n_2^2$)
  6. Tầng chứa nước Pliocen dưới ($n_2^1$)
  7. Tầng chứa nước Miocen trên ($n_1^3$)

Quy trình nghiên cứu rigorous

Cơ sở dữ liệu đầu vào phục vụ mô hình hóa được kiểm định nghiêm ngặt với quy mô lớn:

  • Số liệu địa hình & địa mạo: $7.779\text{ điểm}$ cao độ số hóa bề mặt đất.
  • Số liệu địa chất - ĐCTV: $268\text{ lỗ khoan}$ địa chất sâu; $1.737\text{ kết quả}$ phân tích thành phần hóa học mẫu nước; $3.318\text{ điểm}$ đo sâu điện trường (VES); $268\text{ biểu đồ}$ carota địa vật lý giếng khoan; $5\text{ tuyến}$ mặt cắt ĐCTV chuẩn hóa khu vực.
  • Số liệu thủy văn - khí tượng: Quan trắc liên tục tại $5\text{ trạm}$ khí tượng bề mặt, $35\text{ trạm}$ đo mưa nhân dân, $13\text{ trạm}$ thủy văn và hệ thống trạm hải văn ven biển.
  • Số liệu khai thác: $572\text{ lỗ khoan}$ khai thác công nghiệp quy mô lớn ($>200\text{ m}^3/\text{ngày}$) và $4.516\text{ giếng khoan}$ hộ gia đình phân tán.
  • Số liệu quan trắc động thái: Chuỗi dữ liệu quan trắc mực nước tự động tại hàng chục trạm quốc gia (như các giếng Q17701T, Q40102T, Q403020, Q401030, Q19904T, Q40104Z, Q017050).

Data và phân tích

Hệ thống phần mềm chuyên dụng được vận hành gồm: WetSpass (tính toán lượng bổ cập trên lưới ô vuông GIS), GMS (Groundwater Modeling System) tích hợp MODFLOW (mô phỏng dòng chảy ngầm ba chiều trạng thái không ổn định giai đoạn lịch sử 2000–2010 và dự báo 2015–2090), SEAWAT/MT3DMS (mô phỏng dịch chuyển ranh giới mặn - nhạt với ngưỡng nồng độ khoáng hóa $\text{TDS} = 1\text{ g/L}$).

Thông số thấm của các lớp đất đá được hiệu chỉnh chuẩn xác qua tài liệu hút nước thí nghiệm tại $234\text{ lỗ khoan}$:

  • TCN $qp_3$: Hệ số thấm $k = 0,22 \div 65,82\text{ m/ngày}$ (trung bình $20,62\text{ m/ngày}$).
  • TCN $qp_{2-3}$: $k = 0,89 \div 55,07\text{ m/ngày}$ (trung bình $21,24\text{ m/ngày}$).
  • TCN $qp_1$: $k = 0,76 \div 53,28\text{ m/ngày}$ (trung bình $24,73\text{ m/ngày}$).

Mô hình được hiệu chỉnh (calibration) và kiểm định (validation) thông qua đối sánh chuỗi mực nước tính toán và mực nước thực đo tại các lỗ khoan quan trắc chuẩn quốc gia, kiểm soát sai số căn phương trung bình tuyệt đối (Root Mean Square Error - RMSE) và sai số cân bằng nước toàn bồn trũng luôn dưới $1%$.

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

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

[!IMPORTANT] Dữ liệu phân tích định lượng từ mô hình dòng chảy và vận chuyển chất đã làm sáng tỏ cơ chế tác động độc lập lẫn tương hỗ giữa khai thác nước ngầm và biến đổi khí hậu trên toàn bán đảo Cà Mau:

  1. Suy giảm mực nước áp lực và trữ lượng do khai thác thuần túy (giai đoạn 2000–2010): Văn bản luận án khẳng định:

    "Hoạt động khai thác nước dưới đất trong điều kiện không có biến đổi khí hậu đã làm mực nước áp lực của các tầng chứa nước $qp_3$; $qp_{2-3}$; $qp_1$; $n_2^2$; $n_2^1$; $n_1^3$ suy giảm lần lượt là: 0,33; 0,31; 1,0; 0,91; 0,52; 1,1 m/năm, và làm lượng tích trữ nước dưới đất trung bình hàng năm suy giảm lần lượt là: 3,13; 31,07; 7,01; 6,21; 1,1; và 0,46 triệu m³/năm." Tầng $qp_{2-3}$ gánh chịu mức suy giảm tích trữ nghiêm trọng nhất ($31,07\text{ triệu m}^3/\text{năm}$) do đây là đối tượng khai thác chủ lực của mạng lưới giếng khoan hộ gia đình và công nghiệp quy mô vừa.

  2. Quy luật suy giảm lượng bổ cập dưới kịch bản BĐKH phát thải cao A2: Khi cố định lưu lượng khai thác tại mốc nền 2015, lượng bổ cập ngấm tự nhiên vào bồn trũng suy giảm mạnh mẽ theo thời gian:

    "Tác động của hoạt động khai thác có tính đến sự cộng hưởng của BĐKH làm thay đổi lượng bổ cập theo các khoảng thời gian 2015; 2030; 2045; 2060; 2075; 2090 lần lượt là: 1.548.505; 1.549.563; 1.408.663; 1.281.480; 1.045.515; 936.591 m³/năm, với giả thiết là lượng khai thác giữ nguyên như năm 2015." Mức bổ cập năm 2090 giảm tới $39,5%$ so với năm 2015, chứng minh BĐKH triệt tiêu nghiêm trọng nguồn cấp nước ngập ẩm tái tạo.

  3. Cộng hưởng khai thác và BĐKH giai đoạn 2015–2090: Dưới tác động tích hợp, mực nước áp lực tiếp tục tụt giảm sâu ở các tầng: $qp_3$ ($0,114 \div 0,137\text{ m/năm}$), $qp_{2-3}$ ($0,194 \div 0,232\text{ m/năm}$), $qp_1$ ($0,061 \div 0,064\text{ m/năm}$), $n_2^2$ ($0,133 \div 0,495\text{ m/năm}$), $n_2^1$ ($0,018 \div 0,020\text{ m/năm}$), $n_1^3$ ($0,012 \div 0,248\text{ m/năm}$). Lượng tích trữ nước trung bình năm sụt giảm tương ứng: $0,34$; $2,50$; $0,69$; $0,66$; $0,12$ và $0,12\text{ triệu m}^3/\text{năm}$.


[CHỈ SỐ CẠN KIỆT]

[CHỈ SỐ TỔN THƯƠNG MẶN]
  1. Sự sụp đổ của các chỉ số bền vững cấp huyện/thành phố:
    • Về chỉ số cạn kiệt: Năm 2015 có $21/43\text{ đơn vị}$ ở mức bền vững, $10/43$ tương đối bền vững, $12/43$ không bền vững. Đến năm 2090, số đơn vị bền vững giảm mạnh xuống chỉ còn $7/43$, số đơn vị không bền vững tăng gấp đôi lên $24/43\text{ huyện/thành phố}$.
    • Về chỉ số tổn thương xâm nhập mặn: Năm 2015 có $11/43\text{ đơn vị}$ bền vững ($<50%$ diện tích tầng bị mặn hóa), $10/43$ bền vững trung bình, $22/43$ không bền vững. Đến năm 2090, chỉ còn $7/43\text{ đơn vị}$ bền vững, $11/43$ tương đối bền vững và có tới $25/43\text{ đơn vị}$ rơi vào tình trạng không bền vững.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định cơ chế lan truyền phễu hạ thấp áp lực trong bồn trũng trầm tích trẻ châu thổ không diễn ra đồng tâm đẳng hướng mà bị biến dạng bởi ranh giới mặn - nhạt nguyên sinh và thứ sinh.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình liên kết động WetSpass $\rightarrow$ GMS-MODFLOW $\rightarrow$ SEAWAT $\rightarrow$ UNESCO Sustainability Indicators, tạo mẫu hình tính toán có thể chuyển giao cho các bồn trũng ven biển chịu áp lực kép tại Đông Nam Á.
  • Về ứng dụng thực tiễn và chính sách:
    • Giải pháp phi công trình: Ban hành hạn ngạch khai thác nước dưới đất nghiêm ngặt tại $24\text{ huyện/thành phố}$ nằm trong vùng cạn kiệt nghiêm trọng; chuyển đổi mô hình cấp nước nông thôn phân tán sang hệ thống cấp nước mặt tập trung liên vùng.
    • Giải pháp công trình: Xây dựng hành lang bổ cập nhân tạo nguồn nước (Managed Aquifer Recharge - MAR) trong mùa lũ tại các vùng cồn cát ven biển; vận hành đồng bộ hệ thống kiểm soát cống ngăn mặn ven biển (như cống Cái Lớn, Cái Bé, cống Xẻo Rô) để duy trì áp lực nước mặt, giảm thiểu tốc độ xâm nhập mặn vào các tầng chứa nước nông.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn học thuật:

  • Giả định kịch bản khai thác: Mô hình dự báo 2015–2090 giả định lưu lượng khai thác được giữ cố định ở mốc năm 2015 để cô lập tác động của BĐKH, chưa tích hợp các kịch bản bùng nổ dân số và biến động nhu cầu công nghiệp đột biến.
  • Dữ liệu phân tán tầng sâu: Số liệu địa chất và mẫu nước tại tầng Miocen trên ($n_1^3$) và Miocen giữa - trên ($n_1^{2-3}$) còn tương đối thưa thớt do chiều sâu phân bố lớn ($>300\text{ m}$), dẫn đến mức độ chi tiết hóa ở các tầng sâu này thấp hơn các tầng Pleistocen nông.
  • Quan trắc đồng vị thủy văn: Nghiên cứu chưa áp dụng phương pháp phân tích đồng vị bền ($^{18}\text{O}, ^2\text{H}$) và đồng vị phóng xạ ($^{14}\text{C}, ^3\text{H}$) trên diện rộng để kiểm chứng độc lập tuổi tuyệt đối và nguồn gốc bổ cập cổ thủy văn của các thấu kính nước ngọt.

Chương trình nghiên cứu tiếp nối (Future Research Agenda):

  1. Ứng dụng phương pháp đồng vị môi trường để định lượng chính xác tỷ lệ pha trộn giữa nước bổ cập hiện đại và nước ngầm nguyên sinh trong các TCN Pliocen và Miocen.
  2. Tích hợp mô hình sụt lún mặt đất (Land Subsidence Modeling) liên kết với trường suy giảm áp lực nước lỗ rỗng để cảnh báo nguy cơ ngập lụt đô thị.
  3. Nghiên cứu thực nghiệm các giải pháp bổ cập nhân tạo tầng chứa nước (MAR) sử dụng nguồn lũ sông Hậu.
  4. Ứng dụng trí tuệ nhân tạo và mạng nơ-ron tích hợp vật lý (Physics-Informed Neural Networks - PINNs) để nâng cao tốc độ dự báo động thái nước dưới đất theo thời gian thực.

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

Luận án là công trình nghiên cứu toàn diện, đóng góp cơ sở dữ liệu quan trọng cho chuyên ngành Kỹ thuật Địa chất và Thủy địa chất học tại Việt Nam. Về mặt học thuật, các kết quả mô phỏng và phương pháp liên kết mô hình cung cấp nguồn tư liệu chuẩn mực phục vụ đào tạo sau đại học và các nghiên cứu chuyên sâu về bồn trũng Nam Bộ.

Về phương diện chính sách và kinh tế - xã hội, các bản đồ phân vùng bền vững tại 43 huyện/thành phố là luận cứ trực tiếp giúp Bộ Tài nguyên và Môi trường cùng UBND các tỉnh Cần Thơ, Cà Mau, Bạc Liêu, Sóc Trăng, Hậu Giang, Kiên Giang phê duyệt quy hoạch tài nguyên nước đến năm 2030, tầm nhìn 2050; đồng thời hiện thực hóa các mục tiêu thích ứng sinh kế bền vững theo Nghị quyết 120/NQ-CP của Chính phủ về phát triển bền vững ĐBSCL thích ứng với biến đổi khí hậu.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận liên kết đa mô hình hoàn chỉnh (WetSpass - MODFLOW - SEAWAT - UNESCO Indicators) cùng bộ dữ liệu địa chất thủy văn chi tiết 7 tầng chứa nước.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, Sở TN&MT các tỉnh ĐBSCL): Sử dụng trực tiếp danh mục phân loại $24/43\text{ huyện/thành phố}$ nguy cơ cạn kiệt cao để khoanh định vùng hạn chế và cấp phép khai thác nước dưới đất chuẩn xác.
  • Đơn vị Cấp thoát nước & R&D Doanh nghiệp: Điều chỉnh thiết kế chiều sâu giếng khoan công nghiệp, chuyển dịch khai thác từ các tầng bị cạn kiệt mạnh ($qp_{2-3}, qp_1$) sang các phương án khai thác hoặc tích trữ an toàn hơn.
  • Cộng đồng địa phương: Hưởng lợi gián tiếp thông qua việc bảo vệ nguồn nước ngầm sinh hoạt, ngăn chặn nguy cơ mặn hóa và sụt lún nền đất cục bộ.

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ì?
    Trả lời: Luận án đã mở rộng lý thuyết thủy địa động lực học trong môi trường địa chất trầm tích Đệ tứ bồn trũng Nam Bộ bằng cách chứng minh hiện tượng suy giảm mực nước áp lực cục bộ do khai thác tập trung làm biến đổi gradien thủy lực tự nhiên, từ đó kích hoạt dòng thấm xuyên qua các lớp cách nước bán thấm $mQ$ và gia tốc quá trình xâm nhập mặn ngang từ đới duyên hải vào các thấu kính nước ngọt.

  2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế trước đây?
    Trả lời: So với các nghiên cứu đơn lẻ như Riasat Ali (2012) (chỉ mô phỏng bổ cập và dòng chảy mà bỏ qua ranh mặn) hay Sherif & Singh (1999) (chỉ áp dụng giải tích thủy tĩnh Ghyben-Herzberg 2D), luận án đã tích hợp chuỗi công cụ 3D hoàn chỉnh: WetSpass $\rightarrow$ GMS-MODFLOW $\rightarrow$ SEAWAT/MT3DMS, đồng thời lượng hóa kết quả đầu ra thành bộ chỉ số bền vững không gian theo chuẩn mực UNESCO/IAEA/IAH cho 43 đơn vị hành chính.

  3. Phát hiện bất ngờ và có giá trị thực nghiệm cao nhất từ tập dữ liệu?
    Trả lời: Tầng chứa nước Pleistocen giữa - trên ($qp_{2-3}$) tuy có diện tích phân bố rộng lớn và hệ số thấm cao ($k_{\text{tb}} = 21,24\text{ m/ngày}$) nhưng lại là tầng có lượng tích trữ bị suy giảm nghiêm trọng nhất trong giai đoạn lịch sử ($31,07\text{ triệu m}^3/\text{năm}$), do chịu áp lực khai thác chiếm hơn $50%$ tổng lưu lượng toàn vùng ($628.000\text{ m}^3/\text{ngày}$).

  4. Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) không?
    Trả lời: Hoàn toàn đầy đủ. Luận án chi tiết hóa toàn bộ các bước thiết lập lưới mô hình, cao độ đáy - mái 7 TCN, ma trận hệ số thấm, điều kiện biên (Constant Head, General Head, No-flow), danh mục tọa độ và lưu lượng của $572\text{ giếng}$ công nghiệp, $4.516\text{ giếng}$ hộ gia đình, cùng quy trình tính toán chỉ số cạn kiệt và tổn thương mặn.

  5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?
    Trả lời: Trọng tâm 10 năm tới tập trung vào: (1) Tích hợp mô hình biến dạng đàn - dẻo mô phỏng sụt lún mặt đất do hạ thấp mực nước ngầm; (2) Nghiên cứu công nghệ bổ cập nhân tạo tầng chứa nước (MAR) quy mô công nghiệp; (3) Ứng dụng kỹ thuật đồng vị thủy văn xác định tuổi nước và kiểm chuẩn ranh mặn; (4) Xây dựng hệ thống giám sát và cảnh báo sớm tài nguyên nước ngầm tự động thời gian thực trên nền tảng IoT và GIS.

Kết luận

  1. Thiết lập thành công mô hình cấu trúc ĐCTV 3 chiều đồng bộ gồm 7 tầng chứa nước chính ($qh, qp_3, qp_{2-3}, qp_1, n_2^2, n_2^1, n_1^3$) trên diện tích $16.940\text{ km}^2$ bán đảo Cà Mau.
  2. Định lượng chính xác tốc độ hạ thấp mực nước áp lực ($0,31 \div 1,1\text{ m/năm}$) và tốc độ suy giảm trữ lượng tích trữ ($0,46 \div 31,07\text{ triệu m}^3/\text{năm}$) do hoạt động khai thác giai đoạn 2000–2010.
  3. Dự báo xu thế suy thoái lượng bổ cập ngấm bề mặt dưới kịch bản BĐKH A2, sụt giảm từ $1,55\text{ triệu m}^3/\text{năm}$ (năm 2015) xuống còn $0,94\text{ triệu m}^3/\text{năm}$ (năm 2090).
  4. Phân loại và dự báo mức độ suy thoái tính bền vững tại 43 huyện/thành phố đến năm 2090: số đơn vị không bền vững về cạn kiệt tăng từ 12 lên 24; số đơn vị tổn thương mặn cao tăng lên 25 đơn vị.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, chuyển dịch mô hình khai thác nước ngầm thích ứng linh hoạt với BĐKH và nước biển dâng.
  6. Mở ra 3 hướng nghiên cứu mới về cơ chế thấm xuyên biến dạng cấu trúc trầm tích, bổ cập nhân tạo nguồn nước (MAR) và tích hợp mô hình sụt lún địa tầng đô thị tại đồng bằng châu thổ.