Tổng quan nghiên cứu

Nguồn nước dưới đất tại bán đảo Cà Mau đang chịu áp lực khai thác nặng nề với 137.988 giếng khoan hoạt động, tạo ra tổng lưu lượng khai thác lên đến 373.332 m³/ngày. Đáng chú ý, giếng khoan nhỏ lẻ quy mô hộ gia đình chiếm tới 99,7% tổng số công trình với lưu lượng rút nước đạt 275.180 m³/ngày. Thực trạng khai thác ồ ạt và thiếu kiểm soát này đã làm mực nước dưới đất sụt giảm nghiêm trọng với tốc độ từ 0,80 m/năm đến 1,0 m/năm tại các trạm quan trắc quốc gia. Trong bối cảnh nguồn nước mặt bị nhiễm mặn triệt để vào mùa khô trên toàn bộ diện tích 5.241 km² của tỉnh, các tầng chứa nước ngầm trở thành nguồn cấp nước sinh hoạt và sản xuất duy nhất nhưng lại đối mặt với nguy cơ suy thoái và xâm nhập mặn nghiêm trọng.

Vấn đề nghiên cứu trọng tâm là sự dịch chuyển không gian của ranh giới mặn - nhạt dưới tác động của áp lực khai thác nhân tạo kéo dài. Luận văn đặt mục tiêu xác định ranh giới mặn theo chỉ tiêu tổng chất rắn hòa tan (TDS), đánh giá định lượng mức độ biến động diện tích nhiễm mặn tại tầng chứa nước Pleistocen giữa - trên (qp2-3) giữa hai mốc thời gian: giai đoạn lịch sử (2004-2009) và thời điểm khảo sát hiện trạng (2017).

Nghiên cứu được triển khai trên phạm vi toàn bộ 8 huyện và 1 thành phố thuộc tỉnh Cà Mau, giới hạn trong tọa độ địa lý từ 8°30' đến 9°10' vĩ độ Bắc và 104°80' đến 105°05' kinh độ Đông. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp bộ dữ liệu không gian chuẩn hóa và bản đồ chuyên đề hỗ trợ các nhà quản lý môi trường hoạch định chính sách cấp phép khai thác, bảo vệ an ninh nguồn nước ngầm cho hơn 1,2 triệu người dân địa phươ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 vận dụng lý thuyết thủy địa hóa động lực công trình kết hợp với mô hình lan truyền chất tan trong môi trường xốp để giải thích cơ chế xâm nhập mặn. Khi áp lực khai thác vượt quá lượng bổ cập tự nhiên, phễu hạ thấp mực nước cục bộ hình thành làm thay đổi gradient thủy lực, thúc đẩy sự chuyển dịch bề mặt phân cách giữa nước ngọt và nước mặn. Đồng thời, nghiên cứu tích hợp lý thuyết địa thống kê không gian (Spatial Geostatistics) nhằm mô hình hóa sự biến thiên liên tục của các chỉ tiêu hóa lý trong không gian địa chất phức tạp.

Các khái niệm then chốt được chuẩn hóa bao gồm:

  • Tầng chứa nước Pleistocen giữa - trên (qp2-3): Phân vị địa chất thủy văn chủ đạo tại Cà Mau, phân bố ở chiều sâu mái trung bình 89,04 m và chiều sâu đáy trung bình 103,67 m, có bề dày chứa nước hữu dụng đạt 13,3 m.
  • Tổng chất rắn hòa tan (TDS): Chỉ số biểu thị tổng lượng khoáng chất, muối và kim loại hòa tan trong một đơn vị thể tích nước (tính bằng g/l hoặc mg/l).
  • Độ dẫn điện (EC): Thước đo gián tiếp khả năng dẫn dòng điện của nước, tỷ lệ thuận với nồng độ ion hòa tan.
  • Ranh giới mặn - nhạt: Đường phân định dựa trên ngưỡng giới hạn TDS = 1,0 g/l theo Quy chuẩn kỹ thuật quốc gia về chất lượng nước ăn uống (Thông tư số 04/2009/TT-BYT).

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

Nguồn dữ liệu của luận văn được tổng hợp từ hai giai đoạn quan trắc. Bộ dữ liệu lịch sử giai đoạn 2004-2009 gồm 172 mẫu nước phân tích chỉ tiêu TDS kế thừa từ các đề án địa chất thủy văn trước đây. Bộ dữ liệu hiện trạng năm 2017 bao gồm 233 mẫu, trong đó 203 mẫu được lấy trực tiếp tại thực địa từ các giếng khoan khai thác dân sinh có độ sâu từ 50 m trở lên và 30 mẫu bổ sung từ dự án quy hoạch tài nguyên nước chuyên ngành. Phương pháp chọn mẫu là chọn mẫu có chủ đích kết hợp lưới không gian, ưu tiên mật độ dày tại các khu vực chuyển tiếp ranh giới mặn để tối ưu hóa độ chính xác mô hình.

Chỉ tiêu TDS được xác định tại phòng thí nghiệm bằng phương pháp sấy khô bay hơi mẫu nước đã qua lọc chuẩn ở nhiệt độ cố định 180 ± 2°C đến khối lượng không đổi, kết hợp đo đạc nhanh độ dẫn điện EC ngoài hiện trường bằng thiết bị cầm tay chuyên dụng.

Về phương pháp xử lý dữ liệu, tác giả ứng dụng thuật toán nội suy Ordinary Kriging trên phần mềm ArcGIS 10.3 và MapInfo 15.0. Do tập dữ liệu TDS ban đầu có phân phối lệch trái, dữ liệu đã được chuẩn hóa qua hàm Log10(TDS) để đạt phân phối chuẩn với giá trị trung bình là -0,0608 và phương sai là 0,0497. Mô hình bán phương sai thực nghiệm (Semi-variogram) dạng hàm cầu (Spherical model) được lựa chọn với góc phương vị tối ưu 45° và giá trị Nugget xấp xỉ 0. Việc lựa chọn Kriging thay vì phương pháp nghịch đảo khoảng cách (IDW) là do Kriging có khả năng đánh giá cấu trúc không gian và hạn chế tối đa sai số ước lượng tại các điểm chưa có lỗ khoan quan trắc. Toàn bộ quy trình khảo sát, đo đạc và thành lập bản đồ được thực hiện xuyên suốt trong giai đoạn nghiên cứu hoàn thành vào năm 2018.

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

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

Phân tích không gian dữ liệu thủy địa hóa đã chỉ ra ba phát hiện quan trọng:

  1. Diện tích nhiễm mặn giai đoạn lịch sử (2004-2009): Kết quả nội suy 172 mẫu cho thấy khoảng 40% số mẫu nước có nồng độ TDS vượt ngưỡng quy chuẩn 1,0 g/l. Vùng nước mặn và lợ chiếm tổng diện tích khoảng 2.016 km², tương đương 37,8% diện tích tự nhiên của tỉnh Cà Mau. Nồng độ TDS cực đại ghi nhận được đạt tới 4,19 g/l tại khu vực phía Bắc huyện U Minh. Vùng mặn tập trung chủ yếu ở phía Nam thuộc các huyện Năm Căn, Ngọc Hiển và phân bố rải rác tại các huyện Thới Bình, Trần Văn Thời.

  2. Sự mở rộng và biến động ranh giới mặn hiện trạng (2017): Dựa trên 233 mẫu phân tích hiện trạng, ranh giới mặn 1,0 g/l đã dịch chuyển cục bộ sâu vào các tầng chứa nước ngọt. Phía Đông và Đông Bắc thành phố Cà Mau, cùng các dải chuyển tiếp tại huyện Cái Nước và huyện Trần Văn Thời, ghi nhận sự xâm lấn rõ rệt của pha nước lợ (TDS từ 1,0 đến 3,0 g/l). Tỷ lệ diện tích vùng nước ngọt thích hợp cho ăn uống sinh hoạt bị thu hẹp đáng kể so với tổng thể diện tích nghiên cứu.

  3. Mối tương quan giữa sụt giảm mực nước áp lực và xâm nhập mặn: Mực nước tĩnh tại các trạm quan trắc sụt giảm mạnh theo thời gian khai thác. Điển hình tại công trình quan trắc Q177020 (thành phố Cà Mau), mực nước giảm tổng cộng 6,50 m từ năm 1995 đến năm 2008. Tương tự, tại công trình Q188030 khai thác tầng Pleistocen dưới, mực nước sụt giảm 13,78 m trong giai đoạn 1991-2008. Sự suy giảm liên tục này tạo ra độ dốc thủy lực kéo nước mặn từ các thấu kính lân cận thâm nhập vào đới nước ngọt.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự biến động bất lợi của ranh giới mặn bắt nguồn từ áp lực khai thác tập trung quá mức tại các đô thị và sự phân tán dày đặc của 137.590 giếng khoan nông thôn. Mật độ giếng khoan quá dày nhưng không được thi công trám xi măng cách ly đúng kỹ thuật đã vô tình biến các ống vách giếng thành các "cửa sổ địa chất thủy văn", tạo điều kiện cho nước mặn từ tầng Holocen phía trên ngấm thẳng xuống tầng Pleistocen giữa - trên.

Bên cạnh đó, các hoạt động nuôi tôm thâm canh và chế biến thủy sản (chiếm 94% giá trị sản xuất công nghiệp toàn tỉnh) đòi hỏi khối lượng nước ngầm rất lớn, làm gia tăng tốc độ suy kiệt nguồn nước. Đối chiếu với các nghiên cứu tương tự tại vùng đồng bằng ven biển châu thổ sông Mekong, diễn biến xâm nhập mặn tại Cà Mau có tốc độ phức tạp hơn do cấu trúc địa tầng bị phân cắt và chịu tác động kép của chế độ bán nhật triều Biển Đông lẫn nhật triều Biển Tây.

Để phục vụ công tác quản lý trực quan, các dữ liệu biến động này được thể hiện hiệu quả qua hệ thống bản đồ số hóa chuyên đề: bản đồ phân bố đẳng trị TDS, bản đồ chồng lớp không gian đa thời gian (Overlay Map) hiển thị sự dịch chuyển ranh giới mặn, biểu đồ tần suất Histogram kiểm định phân phối Log(TDS), và các bảng tổng hợp diện tích biến động mặn - nhạt theo từng đơn vị hành chính cấp huyện.

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

Dựa trên kết quả đánh giá thực nghiệm, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm kiểm soát biến động ranh giới mặn và bảo vệ tầng chứa nước:

  • Ban hành quy chế kiểm soát và trám lấp giếng khoan hư hỏng: Cơ quan quản lý cần tiến hành tổng kiểm kê và thực hiện trám lấp kỹ thuật bắt buộc 100% các giếng khoan khai thác không còn sử dụng hoặc bị hỏng theo đúng quy định kỹ thuật của Bộ Tài nguyên và Môi trường. Mục tiêu là triệt tiêu hoàn toàn các phễu thấm xuyên tầng, giảm 20% nguy cơ lan truyền mặn nhân tạo trước năm 2025. Chủ thể thực hiện: Sở Tài nguyên và Môi trường tỉnh Cà Mau phối hợp cùng Ủy ban nhân dân các huyện, thị xã.

  • Hiện đại hóa mạng lưới quan trắc và giới hạn ngưỡng khai thác: Lắp đặt bổ sung 15 đến 20 trạm quan trắc tự động liên tục đo đạc mực nước và độ dẫn điện EC tại tầng Pleistocen giữa - trên. Áp dụng hạn ngạch khai thác, khống chế tổng lưu lượng khai thác toàn tỉnh dưới mức 300.000 m³/ngày trong giai đoạn 2020-2030. Chủ thể thực hiện: Trung tâm Quy hoạch và Điều tra tài nguyên nước quốc gia kết hợp cơ quan tài nguyên địa phương.

  • Chuyển đổi mô hình cấp nước từ phân tán sang tập trung: Đầu tư mở rộng công suất cho 145 trạm cấp nước sinh hoạt nông thôn tập trung và 43 giếng khoan công nghiệp đô thị, từng bước thay thế hệ thống 137.590 giếng khoan tự phát nhỏ lẻ. Mục tiêu là giảm tỷ lệ khai thác nước ngầm hộ gia đình từ 99,7% xuống dưới 65% vào năm 2030. Chủ thể thực hiện: Trung tâm Nước sạch và Vệ sinh môi trường nông thôn cùng Công ty Cấp thoát nước Cà Mau.

  • Triển khai giải pháp trữ nước mặt và bổ cập nhân tạo: Xây dựng các cụm hồ trữ nước mưa và nước mặt ngọt hóa quy mô từ 500.000 m³ đến 1.000.000 m³ tại các vùng sinh thái rừng tràm U Minh. Áp dụng thử nghiệm công nghệ bổ cập nhân tạo (Managed Aquifer Recharge - MAR) cho các tầng ngậm nước sâu trong khung thời gian 2025-2035. Chủ thể thực hiện: Sở Nông nghiệp và Phát triển nông thôn phối hợp với các viện nghiên cứu chuyên ngành.

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

Nội dung và kết quả của luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng sau:

  • 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, Chi cục Bảo vệ Môi trường tỉnh Cà Mau và các tỉnh Đồng bằng sông Cửu Long có thể sử dụng trực tiếp bản đồ phân vùng TDS và ranh giới mặn để xây dựng quy hoạch phân vùng cấm, hạn chế khai thác nước dưới đất đến năm 2030.

  • Các viện nghiên cứu và nhà khoa học chuyên ngành: Các chuyên gia địa chất thủy văn, tài nguyên nước và biến đổi khí hậu có thể ứng dụng quy trình xử lý bán phương sai Semi-variogram và thuật toán Ordinary Kriging làm tài liệu tham khảo cho các nghiên cứu mô hình hóa chất lượng nước ngầm ven biển.

  • Doanh nghiệp cấp nước và các đơn vị sản xuất công nghiệp: Các công ty cấp thoát nước đô thị và nhà máy chế biến thủy sản có căn cứ khoa học chuẩn xác để lựa chọn vị trí đặt giếng khai thác công suất lớn, tránh các khu vực có nguy cơ nhiễm mặn cao và mực nước tĩnh đang tụt dốc.

  • Giảng viên, học viên cao học và sinh viên đại học: Phục vụ học tập, giảng dạy các học phần Quản lý tài nguyên nước, Hệ thống thông tin địa lý (GIS) trong môi trường, và Đánh giá tác động môi trường tại các cơ sở đào tạo đại học khối ngành khoa học tự nhiên.

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

Tầng chứa nước Pleistocen giữa - trên (qp2-3) tại Cà Mau có đặc trưng địa chất thủy văn ra sao?
Tầng qp2-3 phân bố rộng rãi khắp tỉnh Cà Mau ở chiều sâu mái trung bình 89,04 m và chiều sâu đáy 103,67 m, bề dày chứa nước thực tế biến đổi từ 2,0 m đến 31,0 m (trung bình 13,3 m). Thành phần trầm tích chủ yếu là cát hạt mịn đến trung xen kẹp bột sét. Đây là tầng ngậm nước quan trọng bậc nhất, cung cấp lưu lượng cho hàng trăm nghìn giếng khai thác nhưng đang có tốc độ hạ thấp mực nước khoảng 0,80 m/năm.

Nghiên cứu sử dụng tiêu chí và ngưỡng kỹ thuật nào để phân loại nước mặn và nước ngọt?
Nghiên cứu căn cứ vào chỉ tiêu tổng chất rắn hòa tan (TDS) theo Thông tư 04/2009/TT-BYT: nước ngọt đạt chuẩn ăn uống sinh hoạt có TDS ≤ 1,0 g/l, nước lợ có TDS trong khoảng 1,0 - 3,0 g/l, và nước mặn có TDS > 3,0 g/l. Ranh giới mặn - nhạt được xác định chính xác tại đường đẳng trị nồng độ TDS = 1,0 g/l.

Tại sao phương pháp nội suy Kriging lại ưu việt hơn trong việc thành lập bản đồ phân bố mặn?
Thuật toán Ordinary Kriging kết hợp phân tích mô hình hàm cầu (Spherical semi-variogram model) ở góc phương vị 45° cho phép ước lượng tối ưu không gian dựa trên cấu trúc tương quan địa thống kê. Nhờ chuyển đổi dữ liệu Log10(TDS) với phương sai 0,0497, mô hình Kriging xử lý 233 mẫu hiện trạng đạt sai số nội suy nhỏ nhất, phản ánh trung thực gradient nồng độ muối.

Đâu là nguyên nhân chính thúc đẩy sự dịch chuyển ranh giới mặn tại Cà Mau?
Nguyên nhân chủ yếu là việc khai thác nước dưới đất quá mức đạt 373.332 m³/ngày từ 137.988 giếng khoan, vượt quá khả năng bổ cập tự nhiên. Mực nước áp lực bị hạ thấp liên tục từ 6,5 m đến 13,78 m qua các giai đoạn, tạo phễu hạ thấp mực nước và kích hoạt sự xâm nhập của nước mặn vào các tầng chứa nước ngọt.

Làm cách nào để ngăn chặn hiện tượng nhiễm mặn thông tầng do các giếng khoan dân sinh?
Giải pháp cấp thiết là thực hiện quy trình trám lấp kỹ thuật bằng vữa sét hoặc xi măng chuyên dụng cho các giếng khoan hư hỏng nhằm khóa kín các đường rò rỉ xuyên tầng. Đồng thời, địa phương cần đầu tư hệ thống cấp nước tập trung để từng bước thay thế 137.590 giếng khoan nhỏ lẻ, giảm thiểu việc khoan ngầm tự phát qua lớp cách nước.

Kết luận

  • Luận văn đã phân tích toàn diện bộ dữ liệu gồm 172 mẫu nước lịch sử (2004-2009) và 233 mẫu nước hiện trạng (2017) để định lượng biến động mặn tại tầng chứa nước Pleistocen giữa - trên (qp2-3).
  • Xác định diện tích vùng nước mặn và lợ (TDS > 1,0 g/l) chiếm tới 2.016 km², tương đương 37,8% tổng diện tích tự nhiên của tỉnh Cà Mau trong giai đoạn trước, với nồng độ mặn đỉnh điểm đạt 4,19 g/l.
  • Chứng minh áp lực khai thác từ 137.988 giếng khoan với lưu lượng 373.332 m³/ngày là tác nhân chính gây sụt giảm mực nước ngầm với tốc độ 0,80 - 1,0 m/năm, làm thay đổi gradient thủy lực của tầng chứa nước.
  • Ứng dụng thành công phương pháp Ordinary Kriging kết hợp công nghệ GIS để xây dựng hệ thống bản đồ số hóa ranh giới mặn có độ tin cậy và tính trực quan cao.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật, công nghệ và quản lý đồng bộ nhằm bảo vệ và khai thác bền vững nguồn tài nguyên nước dưới đất tỉnh Cà Mau.

Đóng góp khoa học lớn nhất của luận văn là thiết lập phương pháp luận chuẩn xác trong việc tích hợp phân tích địa thống kê với công nghệ GIS để đánh giá rủi ro xâm nhập mặn tầng sâu tại vùng ven biển. Trong giai đoạn 2025-2030, các bước nghiên cứu tiếp theo cần tập trung xây dựng mô hình mô phỏng 3D dòng thấm và lan truyền chất ô nhiễm để dự báo chính xác các kịch bản biến đổi khí hậu. Hãy chủ động liên hệ các đơn vị quản lý chuyên trách để tiếp cận trọn bộ tài liệu nghiên cứu và áp dụng bản đồ phân vùng vào công tác quy hoạch tài nguyên nước thực tế.