Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu cùng hiện tượng nước biển dâng đang gây ra những áp lực chưa từng có đối với tài nguyên nước dưới đất tại các dải đồng bằng ven biển miền Trung Việt Nam. Tại tỉnh Hà Tĩnh, quá trình xâm nhập mặn đã làm suy thoái nghiêm trọng hơn 114 km2 diện tích đất và các tầng chứa nước ngọt ven bờ. Nhu cầu khai thác nước ngầm phục vụ sinh hoạt và phát triển kinh tế xã hội của hơn 1,2 triệu người dân địa phương ngày càng gia tăng, dẫn đến tình trạng khai thác tự phát quá mức, làm hạ thấp mực nước tĩnh và dịch chuyển ranh giới mặn sâu vào đất liền.

Nghiên cứu này được thực hiện nhằm mục tiêu làm sáng tỏ các nhân tố tự nhiên và nhân sinh tác động đến quá trình xâm nhập mặn, xác định cơ chế dịch chuyển của khối nước mặn trong các tầng trầm tích Đệ tứ, đồng thời xây dựng mô hình mô phỏng dự báo và đề xuất các giải pháp kỹ thuật, quản lý nhằm bảo vệ nguồn tài nguyên nước ngầm. Phạm vi không gian của nghiên cứu bao phủ diện tích khoảng 1.900 km2 trải dài qua 8 đơn vị hành chính cấp huyện từ Nghi Xuân, Đức Thọ, Can Lộc, Thạch Hà, Lộc Hà, Thành phố Hà Tĩnh, Cẩm Xuyên đến Kỳ Anh, dọc theo chiều dài 137 km bờ biển và 4 cửa sông lớn gồm Cửa Hội, Cửa Sót, Cửa Nhượng và Cửa Khẩu. Về mặt ý nghĩa thực tiễn, công trình cung cấp cơ sở khoa học định lượng với sai số hiệu chỉnh mô hình kiểm soát ở mức dưới 15%, phục vụ trực tiếp cho công tác quy hoạch cấp nước bền vững đến năm 2030.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu dựa trên nền tảng lý thuyết thủy lực nước dưới đất và các quy luật vận chuyển chất tan trong môi trường rỗng. Khung lý thuyết cốt lõi bao gồm:

  1. Nguyên lý cân bằng thủy tĩnh Ghyben - Herzberg: Thiết lập mối quan hệ độ sâu ranh giới mặn - nhạt dựa trên sự chênh lệch khối lượng riêng giữa nước biển (khoảng 1,025 g/cm3) và nước ngọt lục địa (1,000 g/cm3), theo công thức tính toán độ sâu đáy thấu kính nước ngọt gấp khoảng 40 lần độ cao cột nước ngầm trên mực nước biển.
  2. Phương trình vận chuyển chất tan ba chiều: Tích hợp định luật Fick về khuếch tán phân tử và quá trình phân tán thủy động lực học, biểu diễn sự biến thiên nồng độ ion tổng chất rắn hòa tan (TDS) theo thời gian dưới tác động của dòng thấm đối lưu.
  3. Định luật Archie: Xác định mối tương quan giữa điện trở suất của đất đá chứa nước bão hòa và độ dẫn điện của dung dịch lấp đầy lỗ rỗng, tạo cơ sở chuyển đổi dữ liệu địa vật lý thành chỉ số tổng khoáng hóa TDS.
  4. Cơ chế dòng thấm phụ thuộc mật độ: Phân tích sự bất ổn định đối lưu khi khối nước mặn có tỷ trọng lớn nằm kề cận hoặc phủ lên khối nước nhạt, được kiểm soát thông qua các chỉ số không thứ nguyên như số Peclet (trong dải chuyển tiếp từ 0,4 đến 6,0) và số Rayleigh tới hạn (ngưỡng 40).

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

Nghiên cứu áp dụng quy trình tích hợp liên ngành giữa khảo sát thực địa, thăm dò địa vật lý và mô hình hóa toán học số trị:

  • Thu thập và chuẩn hóa dữ liệu: Kế thừa và tổng hợp số liệu từ 34 lỗ khoan thăm dò địa chất thủy văn với tổng chiều sâu 2.685,5 m, 1.039 mẫu nước phân tích thành phần hóa học toàn diện, cùng chuỗi số liệu quan trắc khí tượng thủy văn tại 107 trạm đo trong khu vực nghiên cứu.
  • Thăm dò địa vật lý chuyên sâu: Áp dụng phương pháp đo sâu điện đối xứng dòng một chiều (VES) với khoảng cách phân bố cực dòng AB/2 mở rộng từ 100 m đến 1.000 m. Việc lựa chọn phương pháp này cho phép quét liên tục điện trở suất theo chiều sâu mặt cắt từ 10 m đến 100 m mà không phá hủy cấu trúc địa tầng, giúp xác định chính xác ranh giới mặn nhạt tại các vùng thiếu mạng lưới giếng khoan quan trắc.
  • Phân tích hệ thống thông tin địa lý (GIS) và mô hình số 3D: Dữ liệu địa chất được số hóa trên nền GIS kết hợp phần mềm mô phỏng dòng thấm Visual MODFLOW và module lan truyền chất tan SEAWAT. Lý do lựa chọn SEAWAT vì đây là công cụ tối ưu mô phỏng đồng thời dòng chảy và sự dịch chuyển muối có xét đến biến thiên mật độ chất lỏng. Quá trình hiệu chỉnh bài toán ngược ổn định và không ổn định đạt chỉ số sai số tuyệt đối trung bình (MAE) và sai số bình quân phương (RMS) dưới 12%, đảm bảo độ tin cậy khoa học cao cho các kịch bản dự báo.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích dữ liệu thực địa và kết quả chạy mô hình số trị đã làm sáng tỏ 3 phát hiện khoa học quan trọng:

  1. Phân bố không gian của thấu kính nước ngọt: Tầng cát biển - gió Holocen muộn (mvQ23) phân bố dọc theo bờ biển tạo thành các cồn cát cao từ 5 m đến 15 m (cá biệt đạt 22 m tại Kỳ Phương), có bề dày trung bình 11,0 m và độ rỗng hữu hiệu từ 35% đến 40%. Đây là bể chứa nước nhạt tự do có chất lượng rất tốt với chỉ số TDS dưới 0,5 g/l, đóng vai trò tấm chắn thủy lực tự nhiên làm chậm quá trình xâm nhập mặn từ biển vào đất liền.
  2. Thực trạng suy thoái tầng chứa nước Pleistocen (qp): Do các hoạt động khai thác công nghiệp và dân sinh tập trung, ranh giới mặn nhạt (ngưỡng TDS = 1,0 g/l) tại tầng qp đã tiến sâu vào lục địa từ 1,5 km đến 3,8 km ở các khu vực hạ lưu sông La, Cửa Sót và Cẩm Xuyên. Diện tích nước ngầm bị nhiễm mặn trong tầng qp chiếm khoảng 28% tổng diện tích phân bố tầng chứa nước toàn tỉnh.
  3. Dự báo động thái lan truyền mặn đến năm 2020 và 2030: Nếu tốc độ khai thác nước ngầm gia tăng bình quân từ 5% đến 8% mỗi năm mà không có giải pháp bổ cập, phễu hạ thấp mực nước sẽ mở rộng với tốc độ hạ thấp mực nước từ 0,25 m đến 0,45 m/năm. Ranh giới mặn 1,0 g/l sẽ lấn sâu thêm từ 350 m đến 900 m về phía Tây, đe dọa trực tiếp đến hơn 45 công trình cấp nước tập trung nông thôn.

Thảo luận kết quả

Cơ chế xâm nhập mặn tại vùng ven biển Hà Tĩnh là kết quả tổng hợp của 3 động lực chính: xâm nhập ngang từ biển do chênh lệch áp lực thủy tĩnh khi mực nước ngầm hạ thấp, thẩm thấu thẳng đứng từ mạng lưới sông ngòi bị nhiễm mặn trong mùa kiệt qua các cửa sổ địa chất thủy văn, và khuếch tán phân tử từ các lớp sét biển tuổi Holocen chứa muối cổ sang các vỉa cát chứa nước ngọt liền kề.

Khi so sánh với các nghiên cứu tương tự tại dải ven biển Nam Định hay bán đảo Cà Mau, tốc độ dịch chuyển ranh giới mặn tại Hà Tĩnh có tính cục bộ cao hơn do sự phân cắt địa hình của các đứt gãy kiến tạo á vĩ tuyến dài trên 30 km và sự phân hóa mạnh của hệ số thấm (dao động từ 2,5 m/ngày ở trầm tích sông biển amQ13 đến 38,0 m/ngày ở trầm tích cát biển mQ21-2). Trên các đồ thị tương quan và bản đồ mặt cắt mô phỏng 3D, sự phân bố độ mặn thể hiện rõ dạng nêm mặn đáy cong hướng về phía lục địa, minh chứng rằng các dòng đối lưu tỷ trọng đóng vai trò chi phối cấu trúc phân tầng thủy địa hóa trong khu vực.

Đề xuất và khuyến nghị

Dựa trên kết quả mô phỏng và đánh giá trữ lượng khai thác an toàn, luận văn đề xuất 4 nhóm giải pháp công nghệ và quản trị mang tính đồng bộ:

  1. Chuyển đổi công nghệ khai thác nước ngầm tầng nông: Thay thế dần các giếng khoan sâu cục bộ bằng hệ thống giếng tia thu nước hoặc hành lang thu nước nằm ngang đặt ở độ sâu từ 3 m đến 7 m trong các cồn cát biển gió mvQ23. Mục tiêu giảm thiểu tối đa sự hình thành phễu hạ thấp mực nước sâu, bảo tồn đáy thấu kính nước ngọt và nâng công suất thu nước an toàn thêm 20% đến 30% tại các trạm cấp nước sinh hoạt huyện Cẩm Xuyên và Kỳ Anh. Thời gian triển khai: Giai đoạn 2015 - 2022. Chủ thể thực hiện: Sở Xây dựng phối hợp với Trung tâm Nước sạch và Vệ sinh môi trường nông thôn Hà Tĩnh.
  2. Xây dựng công trình bổ cập nhân tạo nguồn nước: Thiết lập các hồ trữ nước mưa và hào thấm trên bề mặt các dải cát ven biển với diện tích hứng nước khoảng 50 ha đến 100 ha, kết hợp điều tiết xả nước từ các hồ chứa thượng nguồn (như hồ Kẻ Gỗ, hồ Sông Rác) để duy trì lưu lượng đẩy mặn tối thiểu 15 m3/s tại các cửa sông trong mùa khô. Thời gian triển khai: Giai đoạn 2016 - 2025. Chủ thể thực hiện: Công ty TNHH MTV Thủy lợi Nam Hà Tĩnh và Sở Nông nghiệp và Phát triển nông thôn.
  3. Thiết lập mạng lưới quan trắc tự động và phân vùng cấm khai thác: Lắp đặt 25 trạm quan trắc tự động liên tục đo mực nước, nhiệt độ và độ dẫn điện (EC/TDS) tại 2 tầng chứa nước qh và qp; cắm mốc hành lang bảo vệ nghiêm ngặt cách ranh giới mặn 1,0 g/l tối thiểu 500 m. Nghiêm cấm khoan khai thác mới với lưu lượng trên 50 m3/ngày trong phạm vi dải cách bờ biển dưới 1.000 m. Thời gian triển khai: Giai đoạn 2015 - 2020. Chủ thể thực hiện: Sở Tài nguyên và Môi trường tỉnh Hà Tĩnh.
  4. Khôi phục thảm thực vật và rừng phòng hộ chắn cát: Triển khai trồng mới và bảo vệ 2.500 ha rừng phi lao chắn cát ven biển và rừng ngập mặn tại các cửa sông Cửa Sót, Cửa Nhượng nhằm tăng cường lượng nước ngấm bổ cập tự nhiên thêm 15% và hạn chế quá trình bốc thoát hơi nước trong mùa khô hạn. Thời gian triển khai: Giai đoạn 2015 - 2030. Chủ thể thực hiện: Chi cục Kiểm lâm phối hợp với Ủy ban nhân dân 8 huyện, thị xã ven biển.

Đối tượng nên tham khảo luận văn

Nội dung luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên nước: Lãnh đạo và cán bộ chuyên môn tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn Hà Tĩnh sử dụng số liệu phân vùng mặn và trữ lượng khai thác an toàn để phê duyệt quy hoạch cấp phép thăm dò, khai thác khoáng sản nước dưới đất và xây dựng kịch bản ứng phó biến đổi khí hậu cấp tỉnh.
  2. Đơn vị tư vấn, thiết kế hạ tầng cấp thoát nước: Các kỹ sư công trình ứng dụng sơ đồ thiết kế giếng tia, hành lang thu nước ngang và thông số hệ số thấm thực nghiệm để tối ưu hóa mạng lưới giếng khai thác nước sạch công suất từ 500 m3/ngày đến 5.000 m3/ngày cho các khu dân cư nông thôn và đô thị ven biển.
  3. Viện nghiên cứu và giảng viên, học viên sau đại học: Các nhà khoa học thuộc chuyên ngành Địa chất thủy văn, Quản lý tài nguyên và Khoa học môi trường tham khảo quy trình tích hợp mô hình SEAWAT và phương pháp đo sâu điện đối xứng để phát triển các đề tài nghiên cứu chuyên sâu về thủy địa hóa tại các vùng đồng bằng ven biển nhiệt đới.
  4. Ban quản lý các khu công nghiệp và doanh nghiệp ven biển: Các chủ đầu tư hạ tầng tại Khu kinh tế Vũng Áng, Khu công nghiệp Gia Lách tham khảo bản đồ dự báo ranh giới nhiễm mặn để lựa chọn vị trí xây dựng công trình, thiết kế phương án xử lý nước ngầm phục vụ sản xuất công nghiệp và nuôi trồng thủy sản bền vững.

Câu hỏi thường gặp

  1. Phương pháp đo sâu điện đối xứng (VES) đóng vai trò gì trong việc xác định ranh giới mặn nhạt tại vùng ven biển Hà Tĩnh? Phương pháp đo sâu điện đối xứng với khoảng cách cực AB/2 từ 100 m đến 1.000 m cho phép ghi nhận sự thay đổi điện trở suất biểu kiến của các tầng đất đá chứa nước. Dựa trên định luật Archie và tương quan thực nghiệm với nồng độ TDS từ 1.039 mẫu nước, phương pháp giúp khoanh định nhanh chóng ranh giới đới nước nhạt có điện trở suất cao và đới nước mặn có điện trở suất thấp mà không cần khoan thăm dò tốn kém.

  2. Tại sao mô hình số SEAWAT lại ưu việt hơn mô hình MODFLOW truyền thống trong nghiên cứu này? Phần mềm SEAWAT tích hợp chặt chẽ giữa MODFLOW và MT3DMS, có khả năng giải đồng thời phương trình dòng chảy và phương trình vận chuyển chất tan có xét đến sự biến thiên mật độ chất lỏng (khoảng 2,5% giữa nước biển 1,025 g/cm3 và nước ngọt 1,000 g/cm3). Điều này giúp mô phỏng chính xác hiện tượng nêm mặn và chuyển động đối lưu tự do vốn không thể tính toán bằng MODFLOW đơn thuần.

  3. Dải cồn cát ven biển tuổi Holocen muộn (mvQ23) có ý nghĩa như thế nào đối với nguồn nước ngầm Hà Tĩnh? Tầng trầm tích cát biển gió mvQ23 có bề dày trung bình 11,0 m và độ rỗng lớn từ 35% đến 40%, tạo nên các thấu kính nước ngọt có chất lượng cao nổi trên nền nước mặn. Đây là nguồn cung cấp nước sinh hoạt chính cho hàng trăm nghìn hộ dân ven biển và là vùng bổ cập nước mưa trọng yếu bảo vệ các tầng chứa nước sâu bên dưới.

  4. Vì sao công nghệ giếng tia và hành lang thu nước ngang lại được khuyến nghị thay thế giếng khoan đứng truyền thống? Giếng khoan đứng truyền thống thường hút nước tập trung tại một điểm sâu, dễ tạo phễu hạ thấp mực nước cục bộ khiến nước mặn từ phía dưới bốc lên theo dạng nón mặn. Ngược lại, giếng tia và hành lang ngang thu nước trên diện rộng ở độ sâu nông từ 3 m đến 7 m, duy trì áp lực thủy tĩnh ổn định và ngăn ngừa hiện tượng kéo nêm mặn vào công trình thu.

  5. Tốc độ xâm nhập mặn tại khu vực nghiên cứu được dự báo thay đổi ra sao trong giai đoạn 2020 đến 2030? Dưới kịch bản khai thác gia tăng từ 5% đến 8%/năm kết hợp tác động của nước biển dâng, mô hình dự báo ranh giới nồng độ muối TDS = 1,0 g/l sẽ tiếp tục tịnh tiến vào phía đất liền thêm từ 350 m đến 900 m tại các huyện Thạch Hà, Cẩm Xuyên và Kỳ Anh, làm thu hẹp khoảng 15% diện tích các tầng chứa nước nhạt còn lại.

Kết luận

  • Công trình đã làm sáng tỏ cấu trúc địa chất thủy văn phức tạp của 1.900 km2 dải ven biển Hà Tĩnh, xác định chính xác cơ chế lan truyền muối phụ thuộc mật độ trong trầm tích Đệ tứ.
  • Dữ liệu địa vật lý với 34 lỗ khoan sâu 2.685,5 m và mạng lưới đo sâu điện phân bố rộng đã giúp xây dựng thành công bản đồ phân vùng ranh giới mặn nhạt với độ tin cậy sai số RMS dưới 12%.
  • Mô hình số SEAWAT đã dự báo chi tiết động thái dịch chuyển nêm mặn đến năm 2030, chứng minh nguy cơ suy thoái nghiêm trọng nguồn nước nhạt nếu không kiểm soát việc bơm hút tự phát.
  • Hệ giải pháp công nghệ gồm giếng tia, hành lang thu nước ngang, bổ cập nhân tạo và quản lý vùng cấm khai thác 500 m cung cấp lời giải toàn diện bảo vệ an ninh nguồn nước vùng bờ.
  • Kế hoạch triển khai mạng lưới quan trắc tự động 25 trạm và phục hồi 2.500 ha rừng phòng hộ ven biển cần được chính quyền địa phương phê duyệt và đưa vào quy hoạch phát triển kinh tế xã hội giai đoạn tới. Các nhà quản lý và doanh nghiệp cần sớm khai thác cơ sở dữ liệu của luận văn để định hướng đầu tư hạ tầng cấp nước bền vững cho tỉnh Hà Tĩnh.