Tổng quan nghiên cứu

Khai thác nước dưới đất không kiểm soát đang đẩy tài nguyên nước tại vùng cực Nam bán đảo Cà Mau vào tình trạng báo động nghiêm trọng. Theo số liệu điều tra địa chất thủy văn, toàn tỉnh Cà Mau có 137.988 giếng khoan khai thác với tổng lưu lượng vượt 373.000 m³/ngày, trong đó giếng khoan hộ gia đình quy mô nhỏ chiếm tới 99,7% (137.590 giếng). Nguồn nước mặt hầu như bị nhiễm mặn trong suốt mùa khô do tác động chế độ triều, khiến nguồn nước dưới đất trở thành nguồn cung cấp duy nhất cho sinh hoạt và sản xuất.

Vấn đề nghiên cứu trọng tâm là sự sụt giảm mực nước áp lực liên tục với tốc độ hạ thấp từ 0,80 m/năm đến 1,0 m/năm, gây biến đổi cân bằng thủy lực và dịch chuyển ranh giới mặn vào các tầng chứa nước ngọt. Mục tiêu nghiên cứu cụ thể của đề tài là ứng dụng phương pháp địa thống kê Kriging và hệ thống thông tin địa lý GIS để xác định hiện trạng phân bố tổng hàm lượng chất rắn hòa tan (TDS), mô hình hóa và đánh giá mức độ dịch chuyển ranh giới mặn của tầng chứa nước Pleistocene giữa - trên ($qp_{2-3}$) qua hai mốc thời gian: giai đoạn 2004–2009 và năm 2017.

Phạm vi nghiên cứu bao quát toàn bộ không gian hành chính tỉnh Cà Mau với diện tích tự nhiên trên 5.294 km², tập trung vào các phức hệ tầng chứa nước Đệ Tứ và Neogen. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi xác lập cơ sở dữ liệu định lượng, phục vụ quy hoạch khai thác bền vững và bảo vệ vùng nước nhạt có diện tích phân bố khoảng 418,04 km² trước nguy cơ suy thoái và nhiễm mặn lan rộng.

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 được xây dựng dựa trên sự kết hợp chặt chẽ giữa lý thuyết động lực học dòng thấm - thủy địa chất ven biển và lý thuyết địa thống kê không gian (Geostatistics). Lý thuyết thủy lực nước dưới đất giải thích quy luật tương tác giữa pha nước ngọt và pha nước mặn theo mặt phân giới ranh mặn. Khi lưu lượng khai thác vượt quá trữ lượng bổ cập tự nhiên, cột áp thủy tĩnh suy giảm tạo ra nón hạ thấp mực nước, kích hoạt sự dịch chuyển của khối nước mặn từ biển và các tầng lân cận vào các thấu kính nước ngọt.

Các khái niệm chuyên ngành then chốt bao gồm: Tầng chứa nước lỗ hổng (Porous Aquifer), Tổng độ khoáng hóa/Tổng chất rắn hòa tan (TDS), Xâm nhập mặn nước ngầm (Groundwater Salinization), Ranh giới mặn - nhạt (xác định tại ngưỡng TDS = 1 g/l), và Hàm bán phương sai (Variogram). Mô hình hàm toán học Variogram cùng kỹ thuật chuyển đổi LogTDS được áp dụng để chuẩn hóa chuỗi dữ liệu biến động không gian có độ phân tán lớn.

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

Nguồn dữ liệu nghiên cứu kế thừa toàn diện từ mạng lưới trạm quan trắc quốc gia động thái nước dưới đất (tiêu biểu như các trạm Q177020, Q188020, Q199020, Q188030) và các lỗ khoan thăm dò địa chất thủy văn giai đoạn 2004–2009, kết hợp cùng đợt quan trắc, lấy mẫu phân tích thực tế trên địa bàn tỉnh vào năm 2017. Cỡ mẫu nghiên cứu bao gồm hàng chục vị trí quan trắc phân bố đại diện theo mạng lưới ô vuông không gian tại tầng chứa nước Pleistocene giữa - trên ($qp_{2-3}$) và Pleistocene dưới ($qp_1$).

Phương pháp chọn mẫu áp dụng theo tiêu chuẩn ngành địa chất thủy văn, ưu tiên các khu vực chịu áp lực khai thác cao như thành phố Cà Mau, huyện Năm Căn, huyện Cái Nước và vùng giáp ranh Bạc Liêu. Phương pháp phân tích dữ liệu chủ đạo là nội suy Kriging không gian tích hợp công nghệ GIS. Lý do lựa chọn Kriging vì đây là phương pháp ước lượng tuyến tính không chệch tối ưu (BLUE), có khả năng tính toán trọng số không gian thông qua hàm tương quan biến ngẫu nhiên, giúp tái tạo bề mặt đẳng nồng độ TDS với độ chính xác và độ tin cậy vượt trội so với các phương pháp nội suy hình học thông thường.

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 địa thống kê và bản đồ chuyên đề đã đưa ra các phát hiện cốt lõi về hiện trạng tài nguyên nước dưới đất Cà Mau:

Thứ nhất, tốc độ hạ thấp mực nước áp lực diễn ra liên tục và nghiêm trọng. Tại trạm quan trắc Q177020 (TP. Cà Mau, tầng $qp_{2-3}$), mực nước hạ thấp tổng cộng 6,50 m trong giai đoạn 1995–2008, đạt tốc độ trung bình 0,80 m/năm. Tại trạm Q188030 (tầng $qp_1$), mực nước sụt giảm 13,78 m từ năm 1991 đến năm 2008, tương đương tốc độ suy giảm 1,0 m/năm.

Thứ hai, diện tích phân bố nước nhạt của tầng $qp_{2-3}$ đạt khoảng 418,04 km² (chiếm khoảng 3/4 diện tích khu vực có nước ngọt ngầm toàn tỉnh), tập trung chủ yếu ở khu vực trung tâm thành phố Cà Mau và các xã lân cận như An Xuyên, Tân Thành, Hòa Thành, Định Bình, cùng một phần huyện Thới Bình và U Minh.

Thứ ba, vùng nước lợ và mặn có sự mở rộng cục bộ rõ rệt. Phân vùng nước lợ phía Bắc chiếm diện tích 118 km², phía Nam chiếm 39,8 km² và phía Đông chiếm 2,96 km². Kết quả chồng xếp bản đồ GIS giữa hai giai đoạn 2004–2009 và 2017 cho thấy ranh giới mặn (ngưỡng TDS = 1 g/l) dịch chuyển lấn sâu vào vùng nước nhạt từ 200 m đến 800 m tại các khu vực tập trung đông dân cư, làm thu hẹp khoảng 8% diện tích trữ lượng nước ngọt khả khai.

Thảo luận kết quả

Nguyên nhân cốt lõi của sự dịch chuyển ranh giới mặn bắt nguồn từ cơ cấu khai thác thiếu quy hoạch. Với 137.590 giếng khoan nhỏ lẻ bơm hút 275.180 m³/ngày, mật độ giếng dày đặc đã vô tình tạo nên hàng loạt phễu hạ thấp mực nước cục bộ. Mạng lưới 8 sông chính và 3 kênh cấp I chiếm gần 3% diện tích tự nhiên chịu tác động thủy triều biển Đông (biên độ triều chênh lệch 0,5 m) và biển Tây, làm gia tăng quá trình thẩm thấu mặn qua các cửa sổ địa chất thủy văn.

Kết quả biến động được trực quan hóa rõ nét qua đồ thị biến thiên mực nước theo thời gian (Hydrograph) kết hợp biểu đồ phân bố tần suất LogTDS dạng chuẩn và bề mặt lưới Raster GIS 2D. So với các nghiên cứu trước đây vốn chỉ đánh giá chất lượng nước đơn lẻ hoặc dùng chỉ số WQI tổng quát, nghiên cứu này đã lượng hóa chính xác tọa độ không gian của ranh giới mặn, chứng minh rằng sự suy giảm áp lực tĩnh là động lực chính thúc đẩy mặn xâm nhập đa chiều.

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

Nhằm bảo vệ nguồn nước dưới đất và ngăn chặn đà xâm nhập mặn tại tỉnh Cà Mau, các giải pháp kỹ thuật và quản lý cần được triển khai đồng bộ:

  • Thực hiện quy hoạch và phân vùng khai thác nghiêm ngặt: Sở Tài nguyên và Môi trường tỉnh Cà Mau cần ban hành ngay danh mục vùng cấm, vùng hạn chế khai thác nước dưới đất; mục tiêu giảm 20% lưu lượng bơm hút tại các vệt giáp ranh mặn - nhạt của tầng $qp_{2-3}$ trong lộ trình 2024–2026.
  • Triển khai đề án trám lấp giếng khoan hư hỏng: Tiến hành kiểm tra và xử lý kỹ thuật trám lấp đối với hơn 5.000 giếng khoan không sử dụng hoặc bị nhiễm mặn mỗi năm, loại bỏ hoàn toàn các đường dẫn trực tiếp đưa nước mặn tầng nông xuống tầng sâu.
  • Nâng cấp mạng lưới quan trắc tự động: Lắp đặt hệ thống cảm biến đo mực nước và độ dẫn điện tự động tại 100% các lỗ khoan quan trắc trọng điểm, tích hợp cảnh báo sớm biến động TDS theo thời gian thực trước năm 2025 do Trung tâm Quan trắc Tài nguyên và Môi trường vận hành.
  • Ứng dụng công nghệ bổ cập nhân tạo: Nghiên cứu thí nghiệm bổ cập nước ngầm nhân tạo (MAR) bằng nguồn nước mưa thu gom trong mùa mưa (lượng mưa trung bình đạt từ 218,1 mm đến 385,9 mm/tháng), tạo đệm thủy lực nhân tạo đẩy lùi ranh giới mặn tại các vùng trũng TP. Cà Mau.

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

Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên và môi trường: Sở Tài nguyên và Môi trường Cà Mau, Chi cục Thủy lợi, Cục Quản lý Tài nguyên Nước ứng dụng kết quả để thiết lập hành lang pháp lý, cấp phép thăm dò và xây dựng bản đồ quy hoạch cấp nước an toàn.
  2. Các viện nghiên cứu, trường đại học chuyên ngành: Giảng viên, nhà khoa học thuộc lĩnh vực Địa chất thủy văn, Kỹ thuật Môi trường và Biến đổi khí hậu sử dụng phương pháp luận Kriging - GIS làm tài liệu giảng dạy và mở rộng mô hình hóa 3D dòng thấm.
  3. Doanh nghiệp cấp thoát nước và quy hoạch hạ tầng đô thị: Các công ty cấp nước địa phương định vị an toàn các bãi giếng công nghiệp công suất lớn, tránh nguy cơ nhiễm mặn công trình khai thác sau vài năm vận hành.
  4. Học viên cao học và nghiên cứu sinh: Cung cấp quy trình xử lý dữ liệu địa thống kê thực chiến, cách chuẩn hóa dữ liệu giếng khoan và kỹ năng phân tích bản đồ không gian chuyên sâu.

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

Tại sao nghiên cứu lại tập trung chi tiết vào tầng chứa nước Pleistocene giữa - trên ($qp_{2-3}$)? Tầng $qp_{2-3}$ có trữ lượng khai thác lớn nhất và phân bố rộng khắp tỉnh Cà Mau với bề dày trung bình 13,68 m. Đây là tầng chứa nước chịu áp lực khai thác dân sinh cao nhất với mật độ giếng dày đặc, cung cấp chuỗi số liệu đầy đủ và có độ tin cậy thống kê cao nhất để mô hình hóa ranh giới mặn.

Tốc độ suy giảm mực nước ngầm tại Cà Mau hiện nay ở mức bao nhiêu? Chuỗi số liệu quan trắc dài hạn chỉ ra tốc độ hạ thấp mực nước ngầm trung bình dao động từ 0,80 m/năm (tại trạm Q177020, tầng $qp_{2-3}$) đến 1,0 m/năm (tại trạm Q188030, tầng $qp_1$). Mức sụt giảm tích lũy qua các năm đã vượt quá 6,5 m đến 13,78 m, gây sụt giảm áp lực nước nghiêm trọng.

Phương pháp Kriging mang lại ưu thế kỹ thuật gì nổi bật trong phân tích ranh giới mặn? Kriging tính toán hàm bán phương sai không gian, cho phép nội suy nồng độ TDS không thiên lệch và tối thiểu hóa phương sai sai số ước lượng. Phương pháp này tối ưu hơn hẳn nội suy khoảng cách nghịch đảo (IDW), phản ánh trung thực quy luật phân tán của chất tan trong môi trường trầm tích không đồng nhất.

Nguyên nhân chính gây nên biến động dịch chuyển ranh giới mặn tại địa phương là gì? Nguyên nhân chủ yếu do hoạt động khai thác quá mức từ 137.590 giếng khoan nhỏ lẻ (chiếm 99,7% tổng số giếng, lưu lượng 275.180 m³/ngày) phục vụ nông nghiệp và nuôi trồng thủy sản. Việc suy giảm cột áp tĩnh đã phá vỡ thế cân bằng tự nhiên, khiến nước mặn xâm nhập ngang từ biển và thẩm thấu dọc qua các tầng chứa nước.

Giải pháp nào cấp thiết nhất để bảo vệ vùng nước nhạt 418,04 km² của Cà Mau? Giải pháp cấp bách hàng đầu là đóng cửa và trám lấp các giếng khoan hư hỏng, chuyển đổi mô hình cấp nước sang hệ thống nhà máy nước mặt tập trung, kết hợp quy định hạn ngạch khai thác nghiêm ngặt tại các dải tiếp giáp ranh giới mặn nhằm phục hồi áp lực thủy tĩnh cho tầng chứa nước.

Kết luận

  • Luận văn đã lượng hóa toàn diện đặc điểm địa chất thủy văn và cơ cấu khai thác nước dưới đất tại Cà Mau với 137.988 giếng khoan, xác định tốc độ tụt mực nước ngầm đạt từ 0,80 m/năm đến 1,0 m/năm.
  • Phân định ranh giới và xác lập diện tích vùng nước nhạt tầng Pleistocene giữa - trên ($qp_{2-3}$) đạt 418,04 km², cùng các vùng nước lợ phía Bắc (118 km²) và phía Nam (39,8 km²).
  • Chứng minh tính hiệu quả của phương pháp nội suy Kriging tích hợp GIS trong việc lập bản đồ phân bố TDS và theo dõi biến động dịch chuyển ranh mặn 1 g/l qua các giai đoạn 2004–2009 và 2017.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật, chính sách quản lý và bổ cập nhân tạo khả thi nhằm bảo vệ an ninh nguồn nước tỉnh Cà Mau trong bối cảnh biến đổi khí hậu.
  • Nghiên cứu là nguồn tài liệu chuẩn mực về phương pháp luận và số liệu thực địa, kêu gọi các cấp chính quyền triển khai ngay quy hoạch tài nguyên nước dưới đất đến năm 2030 để ngăn chặn đà suy thoái nguồn nước ngọt quý giá.