Luận văn thạc sĩ về lún mặt đất do khai thác nước ngầm tại đồng bằng sông Cửu Long

Luận văn thạc sĩ nghiên cứu kỹ thuật xây dựng nghiên cứu sự lún mặt đất do khai thác nước ngầm bằng mô hình toán áp dụng vào, đánh giá hiện trạng, phân tích vấn đề, đề xuất biện

• Tác giả: Nguyễn Quốc Thái
• Người hướng dẫn: PGS. Huỳnh Thanh Sơn
• Trường học: Đại Học Bách Khoa - ĐHQG TP.HCM
• Chuyên ngành: Kỹ Thuật Xây Dựng Công Trình Thủy
• Đề tài: Nghiên Cứu Sự Lún Mặt Đất Do Khai Thác Nước Ngầm Bằng Mô Hình Toán. Áp Dụng Vào Một Khu Vực Của Đồng Bằng Sông Cửu Long
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2022
• Địa điểm: TP. HCM

Mục lục chi tiết

• LỜI CẢM ƠN
• TÓM TẮT LUẬN VĂN
• ABSTRACT
• LỜI CAM ĐOAN
• CÁC TỪ VIẾT TẮT
• CÁC KÝ TỰ VÀ THUẬT NGỮ
• DANH MỤC HÌNH VẼ
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.2. Mục đích nghiên cứu
• 1.3. Đối tượng và phạm vi nghiên cứu
• 1.4. Nội dung nghiên cứu
• 1.5. Phương pháp nghiên cứu
• 2. CHƯƠNG 2: CƠ SỞ LÝ THUYẾT CỦA MÔ HÌNH TOÁN GIẢI TÍCH VÀ TOÁN SỐ
• 2.1. Thủy lực giếng khai thác nước ngầm [5], [6]
• 2.1.1. Trường hợp giếng đơn hoàn chỉnh
• 2.1.2. Trường hợp giếng đơn không hoàn chỉnh
• 2.1.3. Trường hợp nhiều giếng khai thác hoạt động đồng thời
• 2.2. Mô hình toán giải tích
• 2.2.1. Lý thuyết cố kết một chiều của lớp sét
• 2.2.2. Công thức tính toán lún do khai thác nước dưới đất
• 2.3. Mô hình toán số
• 2.3.1. Giới thiệu một số mô hình toán số tính lún mặt đất
• 2.3.2. Lựa chọn mô hình toán số
• 2.3.3. Cơ sở lý thuyết mô hình GEO SIGMA/W [13], [14]
• 2.3.4. Cơ sở lý thuyết mô hình GEO SEEP/W [15]
• 3. CHƯƠNG 3: ỨNG DỤNG SIGMA/W CHO MỘT SỐ BÀI TOÁN CƠ BẢN CỦA GIẾNG ĐƠN
• 3.1. Thiết lập mô hình tính toán
• 3.1.1. Xây dựng mô hình
• 3.1.2. Các trường hợp tính toán
• 3.1.3. Điều kiện ban đầu và điều kiện biên
• 3.1.4. Thông số địa chất (đặc tính vật liệu)
• 3.2. Các kết quả tính toán
• 3.2.1. Bài toán 1: Giếng đơn hoàn chỉnh – tầng chứa nước có áp
• 3.2.2. Bài toán 2: Giếng đơn hoàn chỉnh – tầng chứa nước không áp
• 3.2.3. Bài toán 3: Giếng đơn không hoàn chỉnh – tầng chứa nước có áp
• 3.2.4. Bài toán 4: Giếng đơn không hoàn chỉnh – tầng chứa nước không áp
• 3.3. Phân tích độ nhạy các thông số
• 3.3.1. Phương pháp đánh giá độ nhạy thông số [17]
• 3.3.2. Kết quả đánh giá độ nhạy thông số
• 3.3.3. Nhận xét về độ nhạy thông số
• 4. CHƯƠNG 4: ỨNG DỤNG SIGMA/W TÍNH TOÁN LÚN KHI KHAI THÁC NƯỚC NGẦM TẠI HUYỆN CHÂU THÀNH, TỈNH TRÀ VINH
• 4.1. Đặc điểm tự nhiên vùng nghiên cứu
• 4.1.1. Vị trí địa lý
• 4.1.2. Hiện trạng khai thác nước dưới đất
• 4.1.3. Đặc điểm địa chất thủy văn
• 4.1.4. Đặc điểm địa chất công trình
• 4.2. Đánh giá mức độ sụt lún mặt đất do khai thác nước dưới đất bằng mô hình toán số SIGMA/W
• 4.2.1. Xác định khu vực nghiên cứu và phạm vi nghiên cứu
• 4.2.2. Thiết lập mô hình toán số
• 4.2.3. Các trường hợp tính toán
• 4.2.4. Độ lún cố kết theo thời gian tại từng thời điểm tính toán
• 4.2.5. Kết quả tính toán độ lún mặt đất do khai thác nước ngầm
• 4.3. Đánh giá sơ bộ mức độ sụt lún mặt đất do khai thác nước dưới đất bằng mô hình giải tích
• 4.4. Phân tích kết quả tính toán
• 4.4.1. Độ lún trung bình và tốc độ lún trung bình phạm vi nghiên cứu
• 4.4.2. So sánh kết quả tính toán giữa mô hình toán số và mô hình giải tích
• 4.4.3. Dự báo độ lún do khai thác nước ngầm phạm vi nghiên cứu
• 5. CHƯƠNG 5: KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO
• PHỤ LỤC 1: MỘT SỐ CHI TIẾT TÍNH TOÁN ỨNG DỤNG SIGMA CHO BÀI TOÁN CƠ BẢN CỦA GIẾNG ĐƠN
• PHỤ LỤC 2: MỘT SỐ CHI TIẾT TÍNH TOÁN PHÂN TÍCH ĐỘ NHẠY CÁC THÔNG SỐ BÀI TOÁN CƠ BẢN CỦA GIẾNG ĐƠN
• PHỤ LỤC 3: MỘT SỐ CHI TIẾT TÍNH TOÁN BẰNG MÔ HÌNH GIẢI TÍCH
• PHỤ LỤC 4: MỘT SỐ HÌNH ẢNH KẾT QUẢ TÍNH TOÁN BẰNG MÔ HÌNH SIGMA/W
• PHỤ LỤC 5: MỘT SỐ HÌNH ẢNH KẾT QUẢ TÍNH TOÁN LÚN TẠI PHẠM VI NGHIÊN CỨU (H. CHÂU THÀNH, TỈNH TRÀ VINH)
• PHỤ LỤC 6: CHI TIẾT TÍNH TOÁN LÚN TẠI PHẠM VI NGHIÊN CỨU (H. TRÀ VINH) BẰNG MÔ HÌNH GIẢI TÍCH

Tóm tắt

I. Giới thiệu

Nghiên cứu về lún mặt đất do khai thác nước ngầm tại đồng bằng sông Cửu Long là một vấn đề cấp thiết, đặc biệt trong bối cảnh đô thị hóa và biến đổi khí hậu. Tình trạng sụt lún ảnh hưởng đến cơ sở hạ tầng và môi trường sống của người dân. Mục tiêu của nghiên cứu này là thiết lập mối quan hệ giữa mực nước ngầm và độ lún mặt đất, từ đó đưa ra các giải pháp quản lý tài nguyên nước hiệu quả. Theo nghiên cứu, sự biến động của mực nước ngầm không chỉ dẫn đến lún mặt đất mà còn ảnh hưởng đến các vấn đề môi trường khác như biến đổi khí hậu và tác động môi trường.

1.1 Tầm quan trọng của nghiên cứu

Nghiên cứu này không chỉ có ý nghĩa lý thuyết mà còn mang lại giá trị thực tiễn trong việc quản lý tài nguyên nước. Việc áp dụng mô hình toán học như SIGMA/W giúp dự đoán chính xác mức độ lún mặt đất, từ đó hỗ trợ cho các quyết định quy hoạch và xây dựng. Theo một nghiên cứu gần đây, lún mặt đất đã được ghi nhận tăng lên đáng kể tại nhiều khu vực thuộc đồng bằng sông Cửu Long, ảnh hưởng trực tiếp đến an toàn công trình và sinh hoạt của người dân.

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

Phương pháp nghiên cứu bao gồm việc sử dụng mô hình giải tích và mô hình toán số để phân tích lún mặt đất do khai thác nước ngầm. Mô hình SIGMA/W được áp dụng để tính toán sự thay đổi của mực nước ngầm và độ lún mặt đất. Nghiên cứu này cũng sử dụng dữ liệu địa chất và thủy văn để đảm bảo tính chính xác trong các tính toán. Các thông số như độ dẫn thủy lực, áp lực nước lỗ rỗng, và đặc tính của các lớp đất được xem xét kỹ lưỡng. Kết quả từ mô hình sẽ được so sánh với dữ liệu thực tế để đánh giá độ chính xác của mô hình.

2.1 Mô hình toán học

Mô hình toán học được thiết lập dựa trên các phương trình cơ bản của cơ học đất và thủy văn. Mô hình SIGMA/W cho phép tính toán các yếu tố ảnh hưởng đến lún mặt đất như mực nước ngầm và các yếu tố tác động khác. Nghiên cứu sẽ phân tích độ nhạy của các thông số đầu vào đối với kết quả tính toán, từ đó xác định được các yếu tố quan trọng nhất ảnh hưởng đến lún mặt đất.

III. Kết quả nghiên cứu

Kết quả nghiên cứu cho thấy có mối tương quan rõ ràng giữa sự giảm mực nước ngầm và lún mặt đất. Các khu vực có khai thác nước ngầm mạnh mẽ như huyện Châu Thành, tỉnh Trà Vinh đã ghi nhận mức độ lún cao hơn so với các khu vực khác. Phân tích cho thấy rằng việc giảm mực nước ngầm từ 1 đến 2 mét có thể dẫn đến lún mặt đất từ 10 đến 15 cm. Điều này chỉ ra rằng việc quản lý tài nguyên nước cần được thực hiện một cách nghiêm túc để giảm thiểu tác động tiêu cực đến môi trường và an toàn công trình.

3.1 Phân tích độ lún

Đánh giá độ lún cho thấy rằng các yếu tố như độ dày lớp đất và điều kiện địa chất có ảnh hưởng lớn đến mức độ lún mặt đất. Các mô hình tính toán cho thấy sự tương tác giữa áp lực nước và độ lún là rất quan trọng trong việc dự đoán các biến động trong tương lai. Nghiên cứu cũng chỉ ra rằng việc duy trì mực nước ngầm ổn định là cần thiết để ngăn chặn tình trạng lún mặt đất gia tăng.

IV. Kết luận và kiến nghị

Nghiên cứu đã chỉ ra rằng lún mặt đất do khai thác nước ngầm là một vấn đề nghiêm trọng tại đồng bằng sông Cửu Long. Để giảm thiểu tình trạng này, cần có các biện pháp quản lý tài nguyên nước hiệu quả, bao gồm việc kiểm soát khai thác nước ngầm và áp dụng các công nghệ mới trong quản lý nước. Ngoài ra, cần tăng cường công tác tuyên truyền nhằm nâng cao nhận thức của cộng đồng về tác động của khai thác nước ngầm đến môi trường.

4.1 Kiến nghị

Đề xuất xây dựng các chính sách quản lý tài nguyên nước bền vững, khuyến khích sử dụng các nguồn nước thay thế và áp dụng công nghệ tiết kiệm nước. Cần thiết lập các khu vực bảo vệ nguồn nước ngầm và tổ chức các chương trình giám sát thường xuyên nhằm theo dõi tình trạng mực nước ngầm và lún mặt đất. Các nghiên cứu tiếp theo nên tập trung vào việc phát triển các mô hình dự báo chính xác hơn về lún mặt đất trong tương lai.

Nội dung chính

Tổng quan nghiên cứu

Khai thác tài nguyên nước dưới đất quá mức đang là nguyên nhân trực tiếp thúc đẩy tốc độ sụt lún bề mặt tại các đô thị và vùng đồng bằng châu thổ. Tại Thành phố Hồ Chí Minh, với khoảng 100.000 giếng khoan khai thác tổng lưu lượng 710.000 m³/ngày, nhiều khu vực ghi nhận mức lún nghiêm trọng như phường An Lạc (quận Bình Tân) lún tích lũy tới 81,4 cm trong 12 năm với tốc độ cực đại 6,78 cm/năm. Tại Đồng bằng sông Cửu Long, nếu tình trạng khai thác không được kiểm soát, các mô hình dự báo chỉ ra độ lún trung bình toàn vùng có thể đạt từ 68 cm đến 100 cm vào năm 2100. Hiện tượng này làm gia tăng nguy cơ ngập lụt, nứt vỡ đê kè, phá hủy kết cấu hạ tầng kỹ thuật và đẩy nhanh quá trình xâm nhập mặn trong bối cảnh biến đổi khí hậu cùng nước biển dâng.

Mục tiêu cốt lõi của nghiên cứu là xác lập mối quan hệ định lượng giữa việc hạ thấp mực nước ngầm và biến dạng lún mặt đất thông qua phương pháp mô hình toán giải tích kết hợp mô hình phần tử hữu hạn số. Nghiên cứu tập trung phân tích cơ chế nén chặt của các tầng địa chất dưới tác động của các chế độ bơm hút giếng đơn và giếng nhóm. Phạm vi thực nghiệm áp dụng số liệu thực tế tại huyện Châu Thành, tỉnh Trà Vinh với chuỗi quan trắc địa chất thủy văn giai đoạn 2000–2016, đồng thời mô phỏng dự báo các kịch bản biến dạng nền đất đến năm 2100. Kết quả nghiên cứu cung cấp cơ sở khoa học định lượng giúp các nhà quản lý quy hoạch tài nguyên nước, thiết kế công trình thủy công và kiểm soát biến dạng nền móng công trình ven biển.

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 trên nền tảng cơ học đất cổ điển và thủy lực ngầm hiện đại, tập trung vào ba lý thuyết trụ cột:

• Thuyết cố kết thấm một chiều của Terzaghi (1925): Giải thích quá trình nén chặt của tầng đất dính hạt mịn khi nước lỗ rỗng thoát dần ra ngoài dưới tác dụng của gradient thủy lực. Phương trình vi phân cố kết một chiều mô tả sự tiêu tán áp lực nước lỗ rỗng thặng dư theo chiều sâu và thời gian thông qua hệ số cố kết đứng $C_v$ và hệ số nén thể tích $m_v$.
• Nguyên lý ứng suất hữu hiệu ($\sigma' = \sigma - u$): Khi áp lực nước lỗ rỗng ($u$) suy giảm do mực nước ngầm bị hạ thấp, tổng ứng suất ($\sigma$) do trọng lượng bản thân đất không đổi sẽ làm gia tăng trực tiếp ứng suất có hiệu ($\sigma'$), nén chặt khung hạt khoáng vật và gây sụt lún bề mặt.
• Thủy lực giếng ngầm theo định luật Darcy và giả thiết Dupuit: Thiết lập phương trình dòng thấm ổn định và không ổn định hướng tâm về phía giếng khoan cho cả hai trường hợp tầng chứa nước có áp và không áp.

Các khái niệm cơ bản xuyên suốt mô hình gồm:

• Độ dẫn suất thủy lực ($K$): Khả năng truyền dẫn nước qua môi trường rỗng, tính bằng m/s hoặc m/ngày.
• Áp lực nước lỗ rỗng ($u$): Áp suất thủy tĩnh của chất lưu trong các khe rỗng của đất đá.
• Chỉ số nén ($C_c$) và nén lại ($C_r$): Thông số đặc trưng cho độ biến dạng thể tích của mẫu đất sét chịu tải.
• Độ lún cố kết sơ cấp ($S_c$) và đàn hồi ($S_e$): Các thành phần biến dạng chính của nền đất nhiều tầng.

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

Nghiên cứu tích hợp phương pháp giải tích lý thuyết và phương pháp số phần tử hữu hạn (FEM 2D) thông qua bộ công cụ GeoStudio:

• Quy mô và phương pháp chọn mẫu: Thu thập và chuẩn hóa dữ liệu địa chất công trình từ 7 cụm giếng khoan đại diện phân bố trên địa bàn huyện Châu Thành, tỉnh Trà Vinh cùng dữ liệu đo áp liên tục tại trạm quan trắc quốc gia Q404020 ở tầng chứa nước Pleistocen dưới ($qp_3$). Phương pháp chọn mẫu dựa trên mật độ khai thác tập trung và sự đồng nhất của các tầng trầm tích sét - cát đan xen điển hình ở độ sâu từ 40 m đến hơn 120 m.
• Quy trình mô hình hóa: Sử dụng module SEEP/W để mô phỏng trường dòng thấm không ổn định và xác định phân bố hạ thấp mực nước ngầm $S(r,t)$ quanh giếng bơm. Kết quả áp lực nước lỗ rỗng từ SEEP/W được ghép nối liên hoàn (uncoupled analysis) sang module SIGMA/W nhằm phân tích trạng thái ứng suất - biến dạng và xác định trường lún cố kết $\Delta B$ theo thời gian.
• Lý do lựa chọn phương pháp: Mô hình số phần tử hữu hạn 2D cho phép xử lý chính xác các điều kiện biên thủy lực phức tạp, tính phi tuyến của vật liệu đất yếu nhiều lớp và tương tác địa tầng thực tế mà các công thức giải tích đơn giản không thể phản ánh đầy đủ.

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

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

• Mối tương quan giữa lưu lượng bơm và độ sụt lún: Khi lưu lượng khai thác ($Q$) tăng từ 25 m³/h lên 205 m³/h ở bài toán giếng đơn hoàn chỉnh, độ hạ thấp mực nước trong giếng ($S_0$) tăng từ 0,36 m lên 2,93 m (tăng hơn 713%), kéo theo độ lún cực đại tại tâm giếng tăng tuyến tính từ 0,42 cm lên 3,45 cm.
• Mức độ nhạy cảm của các thông số địa kỹ thuật: Phân tích độ nhạy cho thấy lưu lượng bơm ($Q$) và hệ số dẫn suất thủy lực ($K$) là hai yếu tố chi phối mạnh nhất đến biến dạng lún, chiếm hơn 65% mức độ biến thiên tổng thể của độ lún mặt đất. Bán kính giếng ($r_0$) và khoảng cách ống lọc có độ nhạy thấp hơn, chỉ làm thay đổi biên độ lún dưới 8% khi kích thước ống thay đổi từ 0,125 m đến 0,5 m.
• Tốc độ lún thực tế tại khu vực Châu Thành (Trà Vinh): Trong giai đoạn quan trắc 2000–2016, tổng độ hạ thấp mực nước trung bình tại tầng $qp_3$ đạt khoảng 6,8 m. Mô hình tính toán cho thấy độ lún bề mặt tích lũy trong 16 năm dao động từ 19,2 cm đến 38,6 cm tùy vị trí, tương ứng với tốc độ lún trung bình từ 1,2 cm/năm đến 2,41 cm/năm.
• Mức độ tương quan giữa mô hình số và giải tích: Sai số tính toán độ lún giữa phần mềm SIGMA/W và công thức giải tích Terzaghi - Girinski dao động trong phạm vi từ 5,2% đến 11,8%, xác nhận độ tin cậy cao của phương pháp mô phỏng số.

Thảo luận kết quả

Cơ chế gây lún chủ yếu xuất phát từ quá trình nén ép cố kết của lớp sét dẻo mềm dày khoảng 40 m nằm phía trên tầng chứa nước cát hạt mịn dày 80 m. Khi nước ngầm bị hút ra liên tục, áp lực nước lỗ rỗng suy giảm tạo ra gradient áp lực thẳng đứng, khiến nước từ lớp bùn sét thoát dần vào tầng cát có tính thấm cao hơn. Quá trình này chuyển phần lớn tải trọng thủy tĩnh sang tải trọng hạt khoáng, gây biến dạng thể tích không phục hồi.

Dữ liệu mô phỏng được trực quan hóa thông qua bảng tổng hợp ứng suất - độ lún và đồ thị tương quan đa biến $Q - K - \Delta B$. Đồng thời, phương pháp nội suy không gian nghịch đảo khoảng cách (IDW) đã xây dựng thành công 5 bản đồ đẳng trị độ lún bề mặt qua các mốc năm 2001, 2004, 2008, 2012 và 2016. Khi đối chiếu với các công bố khoa học quốc tế sử dụng công nghệ viễn thám InSAR tại Đồng bằng sông Cửu Long (tốc độ lún trung bình 1,5–3,0 cm/năm), kết quả mô hình SIGMA/W cho thấy sự trùng khớp chặt chẽ về biên độ biến dạng. Điều này khẳng định nếu duy trì kịch bản khai thác tự do với mức tăng trưởng 2% mỗi năm, độ lún tích lũy đến năm 2100 có thể vượt ngưỡng 85 cm, biến nhiều khu vực trũng thấp thành vùng ngập triều thường trực.

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

• Áp dụng hạn ngạch và khoanh vùng hạn chế khai thác nước dưới đất: Sở Tài nguyên và Môi trường tỉnh Trà Vinh cùng các tỉnh lân cận cần ban hành quy chế kiểm soát lưu lượng, cắt giảm 30% đến 40% công suất khai thác tại các khu công nghiệp và vùng đô thị tập trung trước năm 2028, nhằm giữ tốc độ lún nền dưới ngưỡng 1,0 cm/năm.
• Hiện đại hóa mạng lưới trạm quan trắc tự động liên tầng: Cục Quản lý Tài nguyên Nước phối hợp cùng chính quyền địa phương triển khai lắp đặt 15–20 cụm cảm biến đo áp lực nước lỗ rỗng và mốc chuẩn độ lún sâu tự động (extensometer) đa tầng tại các tầng chứa nước $qp_3$ và $qp_{2-3}$ trong giai đoạn 2025–2027.
• Thực hiện giải pháp bổ cập nhân tạo nguồn nước ngầm (Managed Aquifer Recharge - MAR): Sở Nông nghiệp và Phát triển Nông thôn chủ trì xây dựng các công trình lắng lọc và tái nạp nước mặt mùa lũ vào các tầng chứa nước với lưu lượng thử nghiệm 10.000–15.000 m³/ngày, hoàn thành kiểm định kỹ thuật vào năm 2029 để phục hồi áp lực thủy tĩnh.
• Tích hợp bản đồ cảnh báo rủi ro lún vào quy hoạch xây dựng hạ tầng: Các viện quy hoạch và doanh nghiệp tư vấn xây dựng bắt buộc áp dụng chỉ số lún dự báo từ mô hình SIGMA/W vào tiêu chuẩn tính toán móng cọc, đê bao ngăn mặn và hệ thống thoát nước đô thị với tầm nhìn quy hoạch đến năm 2050, triển khai ngay từ năm 2026.

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

• Cơ quan quản lý nhà nước về tài nguyên và môi trường: Cán bộ chuyên trách tại các Sở Tài nguyên và Môi trường khu vực Đồng bằng sông Cửu Long có thể sử dụng kết quả dự báo và phương pháp phân tích độ nhạy để xây dựng chính sách cấp phép, xác định vùng hạn chế khai thác nước dưới đất sát với thực tế địa chất.
• Kỹ sư tư vấn thiết kế và địa kỹ thuật công trình: Các chuyên gia thiết kế hạ tầng thủy lợi, đê kè, cầu đường và nhà cao tầng có thể áp dụng quy trình mô phỏng liên hoàn SEEP/W - SIGMA/W để dự báo chính xác độ lún phụ do suy giảm mực nước ngầm, từ đó tối ưu hóa chiều sâu móng cọc và cao độ đỉnh công trình.
• Giảng viên, nghiên cứu sinh và học viên cao học: Luận văn là tài liệu tham khảo chuyên sâu thuộc nhóm ngành Kỹ thuật Xây dựng Công trình Thủy, Địa kỹ thuật Xây dựng và Quản lý Tài nguyên Nước, cung cấp chi tiết phương pháp giải bài toán cố kết thấm 2D và cách thiết lập mô hình số.
• Doanh nghiệp cấp nước và quản lý khu công nghiệp: Đơn vị vận hành các nhà máy nước ngầm có thêm công cụ định lượng để tính toán lưu lượng bơm an toàn, bố trí khoảng cách giữa các giếng khoan tối ưu nhằm hạn chế tối đa nguy cơ sụt lún công trình lân cận.

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

• Cơ chế vật lý nào khiến việc hạ thấp mực nước ngầm gây lún mặt đất?
  Khi nước ngầm bị hút cạn, áp lực nước lỗ rỗng bên trong các lỗ rỗng của đất giảm xuống. Theo định luật ứng suất hữu hiệu Terzaghi, tổng tải trọng bản thân đất không đổi khiến toàn bộ phần áp lực nước mất đi chuyển thành ứng suất hữu hiệu nén trực tiếp lên khung hạt đất sét, làm giảm thể tích lỗ rỗng và gây sụt lún bề mặt với độ lún có thể vượt 30 cm sau nhiều năm.

• Mô hình toán số SIGMA/W có ưu điểm gì vượt trội so với công thức giải tích?
  SIGMA/W sử dụng phương pháp phần tử hữu hạn 2D cho phép mô phỏng địa tầng nhiều lớp phi tuyến tính, tích hợp điều kiện biên thấm thay đổi linh hoạt từ SEEP/W và phản ánh đúng quá trình cố kết theo thời gian. Sai số giữa mô hình số và thực nghiệm chỉ dao động quanh mức 5% đến 12%, giải quyết được những giới hạn mà công thức giải tích đồng nhất không làm được.

• Thông số nào ảnh hưởng mạnh nhất đến tốc độ lún mặt đất quanh giếng khoan?
  Kết quả phân tích độ nhạy khẳng định lưu lượng bơm ($Q$) và độ dẫn suất thủy lực ($K$) là hai thông số quyết định nhất. Tăng lưu lượng bơm lên 4 lần có thể làm tăng độ lún cục bộ lên hơn 5 lần, trong khi thay đổi bán kính ống vách giếng từ 0,125 m lên 0,5 m chỉ làm thay đổi độ lún dưới 8%.

• Tầng chứa nước nào chịu tác động khai thác nặng nề nhất tại Châu Thành, Trà Vinh?
  Tầng chứa nước Pleistocen dưới ($qp_3$) ở độ sâu trung bình từ 80 m đến 120 m là tầng bị khai thác tập trung nhất phục vụ sinh hoạt và tưới tiêu. Mực nước đo áp tại tầng này đã hạ thấp trung bình hơn 0,42 m/năm trong giai đoạn 2000–2016, đóng góp hơn 70% vào tổng độ lún toàn khu vực.

• Hậu quả dài hạn đối với Đồng bằng sông Cửu Long nếu không kiểm soát việc khai thác nước ngầm là gì?
  Theo các kịch bản mô phỏng đến năm 2100, tốc độ lún trung bình toàn vùng có thể đạt 1,0–2,5 cm/năm, dẫn tới tổng mức lún tích lũy từ 68 cm đến 100 cm. Điều này khiến hơn 40% diện tích đồng bằng có nguy cơ chìm dưới mực nước biển triều cường, gây vỡ đê bao và nhiễm mặn nguồn nước mặt nghiêm trọng.

Kết luận

• Luận văn đã xây dựng thành công quy trình mô phỏng liên hoàn giữa mô hình dòng thấm SEEP/W và mô hình biến dạng SIGMA/W để giải quyết bài toán lún do khai thác nước ngầm.
• Phân tích độ nhạy xác định lưu lượng khai thác ($Q$) và hệ số thấm ($K$) là hai yếu tố chi phối then chốt đến biên độ sụt lún nền đất.
• Dữ liệu thực nghiệm tại huyện Châu Thành (Trà Vinh) giai đoạn 2000–2016 ghi nhận độ lún tích lũy lên tới 38,6 cm với tốc độ trung bình 1,2–2,41 cm/năm.
• Mô hình dự báo cảnh báo nguy cơ sụt lún ĐBSCL có thể đạt ngưỡng 68–100 cm vào năm 2100 nếu tiếp tục duy trì mức khai thác nước ngầm hiện tại.
• Các giải pháp giảm thiểu cần tập trung cắt giảm 30%–40% sản lượng khai thác giếng khoan, lắp đặt hệ thống quan trắc áp lực lỗ rỗng tự động và triển khai bổ cập nhân tạo nguồn nước trước năm 2030.

Công trình là đóng góp học thuật và thực tiễn quan trọng trong công tác quản lý tài nguyên nước và địa kỹ thuật xây dựng vùng đồng bằng ven biển. Các cơ quan quản lý, viện nghiên cứu và doanh nghiệp hạ tầng cần nhanh chóng ứng dụng các kết quả mô hình số này vào việc hoạch định chính sách cấp phép khai thác nước ngầm và tiêu chuẩn thiết kế công trình nhằm đảm bảo sự phát triển bền vững cho Đồng bằng sông Cửu Long.

Trích đoạn nội dung tài liệu

Đặt vấn đề Nước là nguồn tài nguyên không thể thiếu, cần thiết cho mọi hoạt động sống. Đối với con người, nước được dùng trong các hoạt động nông nghiệp, công nghiệp, môi trường, giải trí và hơn hết là duy trì sự sống. Nước bao phủ 71% diện tích Trái Đất, trong đó có 97% là nước mặn, còn lại là nước ngọt. Nước ngầm là một dạng nước dưới đất, là nước ngọt được chứa trong các lỗ rỗng của đất hoặc đá.

Nó cũng có thể là nước chứa trong các tầng ngậm nước bên dưới mực nước ngầm. Đôi khi người ta còn phân biệt nước ngầm nông, nước ngầm sâu và nước chôn vùi. Hiện nay, với sự tiến bộ của công nghệ khai thác (công nghệ khoan, bơm, …) cũng như nhu cầu dùng nước ngày càng tăng của con người, nước ngầm đang bị khai thác quá mức, vượt quá giới hạn nguồn tự nhiên. Điều này dẫn đến nhiều hệ lụy khác nhau, như: ô nhiễm nguồn nước, xâm nhập mặn, giảm về lượng và chất nguồn nước ngầm, … trong đó sự sụt, lún đất do khái thác nước ngầm quá mức là vấn đề cần quan tâm, nhất là ở các thành phố lớn, vùng châu thổ đông dân cư trên thế giới.

Hiện tượng sụt, lún đất do khai thác quá mức nước ngầm diễn ra tại nhiều nơi trên thế giới. Các sự cố về sụt, lún được ghi nhận tại New Orleans, Arizona, Nevada, Texas, California – Hoa Kỳ (Hình 1.2); Bangkok – Thái Lan; Suzhou – Trung Quốc; Osaka – Nhật Bản; Jakarta – Indonesia (Hình 1.3),… đều liên quan đến khai thác quá mức nước ngầm. Điều này kéo theo lũ lụt cục bộ, ngập úng đô thị, làm biến dạng hư hỏng các công trình xây dựng.1 Khai thác nước ngầm khu vực gần hồ Lucerne Lake, trong mùa khô tại San Bernardino County, Mojave Desert, California làm sụt nền đất, gây nứt lòng hồ. Nguồn: https://www.gov/ Trang 1 Hình 1.2 Sụt lún đất tại San Joaquin Valley, California.

(Nguồn: https://phys.org/news/) Hình 1.3 Gần 40% thành phố Jakarta, Indonesia chìm dưới mực nước biển, hậu quả của việc khai thác nước ngầm, gây ngập lụt nghiêm trọng (The New York Times, 2017). Nguồn: https://storymaps.com/ Tại TP. Hồ Chí Minh, hiện nay có khoảng 100.000 giếng khoan khai thác nước ngầm có đường kính và độ sâu khác nhau, phân bố không đều trên các khu vực với tổng lượng khai thác nước ngầm là 710.000 m3/ngày do các hộ khai thác đơn lẻ, các đơn vị, doanh nghiệp khai thác 225.000 m3/ngày, còn lại là Tổng công ty Cấp nước Sài Gòn (SAWACO) [1]. Khối lượng khai thác này đã gần tiệm cận với trữ lượng khai thác an toàn và có nguy cơ thiếu an toàn.

Hiện tượng lún xảy ra nghiêm trọng, cụ thể: tại phường An Lạc, quận Bình Tân (TPHCM) trong 12 năm lún tới 81,4 cm. Đây cũng là nơi có tốc độ lún lớn nhất, lên tới 6,78 cm/năm [2], tại Quận 6 là từ 10÷15 cm, Quận 1 là từ 15÷20 cm, thậm chí tại Quận 2 và 7 là hơn 20 cm tại thời điểm ghi nhận một số vị trí trong năm 2003 và 2004 [3]. Trang 2 Đối với khu vực Đồng bằng sông Cửu Long (ĐBSCL), hoạt động khai thác nước ngầm diễn ra phức tạp gây sụt lún đất nghiêm trọng, đặc biệt ở những vùng hoạt động sản xuất, sinh hoạt phụ thuộc nhiều vào tài nguyên nước dưới đất. Theo nghiên cứu về tình trạng lún do khai thác nước ngầm vừa được thực hiện, nếu tình trạng khai thác nước ngầm vẫn tiếp tục tăng cao như hiện nay, cho đến năm 2100, ĐBSCL có thể lún trung bình 68 cm đến 100 cm phụ thuộc vào mức khai thác nước ngầm [4].

Sự biến động của mực nước ngầm dẫn đến sụt lún nền và công trình xây dựng đã và đang tiếp diễn tại các thành phố lớn. Vì vậy, nghiên cứu thiết lập mối tương quan giữa hạ thấp cao độ mực nước ngầm và sụt lún nền sẽ giúp đưa ra những qui hoạch và phương thức quản lý hợp lý, nhất là trong bối cảnh của sự đô thị hóa như hiện nay. Cùng với sự phát triển của công nghệ, các mô hình toán số dùng để giải quyết bài toán sụt lún đất đã ra đời, một trong những mô hình đó là bộ phần mềm GeoStudio, mà cụ thể là phần mềm con SIGMA/W, do công ty GEOSLOPE International Ltd (Canada) phát triển.2 Mục đích nghiên cứu Nghiên cứu sự lún mặt đất khi khai thác nước ngầm bằng mô hình toán (giải tích và toán số). Áp dụng vào một khu vực cụ thể thuộc ĐBSCL.3 Đối tượng và phạm vi nghiên cứu - Đối tượng nghiên cứu: Nghiên cứu độ sụt lún mặt đất khi khai thác nước ngầm bằng mô hình toán (giải tích và toán số).

- Phạm vi nghiên cứu: Vùng Đồng bằng sông Cửu Long.4 Nội dung nghiên cứu - Cơ sở lý thuyết về nước dưới đất, giếng khai thác nước ngầm và biến dạng lún. - Cơ sở lý thuyết của phần mềm SIGMA/W. - Nghiên cứu đề xuất một biểu thức giải tích cho phép ước tính độ lún mặt đất do khai thác nước ngầm. - Ứng dụng phần mềm SIGMA/W vào nghiên cứu một số bài toán cơ bản khi khai thác nước ngầm và lún mặt đất.

- Ứng dụng phần mềm SIGMA/W vào nghiên cứu độ lún mặt đất do khai thác nước ngầm tại một khu vực cụ thể thuộc ĐBSCL.5 Phương pháp nghiên cứu Các phương pháp nghiên cứu chính được sử dụng trong đề tài này bao gồm: - Phương pháp kế thừa: kế thừa kết quả nghiên cứu của các đề tài, dự án trên thế giới cũng như tại Việt Nam liên quan đến dòng thấm nước dưới đất cũng như vấn đề sụt lún đất do khai thác nước ngầm. - Phương pháp thu thập, tổng hợp thống kê tài liệu nghiên cứu trước đây. - Phương pháp mô hình toán giải tích và toán số. Trang 3 CHƯƠNG 2 CƠ SỞ LÝ THUYẾT CỦA MÔ HÌNH TOÁN GIẢI TÍCH VÀ TOÁN SỐ 2.1 Thủy lực giếng khai thác nước ngầm [5], [6] 2.1 Trường hợp giếng đơn hoàn chỉnh 2.1 Tầng thấm nước có áp Xét giếng hoàn chỉnh bán kính r0 được khoan tới tận tầng không thấm của một tầng thấm nước có áp đồng chất, đẳng hướng có chiều dày B và độ dẫn suất thuỷ lực K (Hình 2.

Q - lưu lượng bơm từ giếng r - bán kính tính từ tâm giếng đến vị trí đo áp R1 - bán kính ảnh hưởng r0 - bán kính giếng H0 - cột nước đo áp trong giếng H - cột nước đo áp tại vị trí cách tâm giếng một khoảng cách r H1 - cột nước đo áp ban đầu B - bề dày tầng thấm nước S0 - độ hạ thấp mực nước trong giếng S - độ hạ thấp mực nước tại vị trí cách tâm giếng một khoảng cách r Hình 2.1 Sơ đồ tính trường hợp giếng đơn hoàn chỉnh – tầng thấm nước có áp. Theo [5], sử dụng các giả thiết của Dupuit, phương trình tổng quát xác định lưu lượng bơm Q (hệ toạ độ cực, lỗ khoan là gốc) tại một khoảng cách r bất kỳ như sau: Trang 4  2 T  H1  H  2 TS Q    R1 R  ln ln 1  r r (2-1)  H1  H  S  Q R  ln 1 2 T r Tại vị trí thành giếng  r  r0 , H  H 0 , S  S0  , ta có:   H1  H 0  S0  2 T ln r Q R1   2 T  H1  H 0  2 TS0 0    (2-2) Q  R R  ln 1 ln 1 r0 r0 Trong đó: Q (m³/s) là lưu lượng bơm từ giếng; r (m) là bán kính tính từ tâm giếng đến vị trí đo áp; H (m) là cột nước đo áp tại vị trí cách tâm giếng một khoảng cách r; H1 (m) là cột nước đo áp ban đầu; H0 (m) là cột nước đo áp trong giếng; T  K  B (m²/s) là độ dẫn truyền; K (m/s) là độ dẫn suất thủy lực; B (m) là bề dày tầng thấm nước; R1 (m) là bán kính ảnh hưởng; S0 (m) là độ hạ thấp mực nước trong giếng; r0 (m) là bán kính giếng; Theo [7], một vài biểu thức kinh nghiệm xác định bán kính ảnh hưởng R1:  Công thức Siechardt: R1  10S0 K  Công thức Kusakin: R1  2 S 0 H1 K Với K (m/ngày) là độ dẫn suất thủy lực.2 Tầng thấm nước không áp Xét giếng hoàn chỉnh bán kính ro được khoan tới tận tầng không thấm của một TTN không áp đồng chất, đẳng hướng có chiều dày B và ĐDSTL K (Hình 2. Trang 5 Q - lưu lượng bơm từ giếng r - bán kính tính từ tâm giếng đến vị trí đo áp R1 - bán kính ảnh hưởng d - khoảng cách tính từ MĐTN đến đường mực nước ban đầu H0 - cột nước đo áp trong giếng H - cột nước đo áp tại vị trí cách tâm giếng một khoảng cách r H1 - cột nước đo áp ban đầu B - bề dày tầng thấm nước; r0 - bán kính giếng S0 - độ hạ thấp mực nước trong giếng S - độ hạ thấp mực nước tại vị trí cách tâm giếng một khoảng cách r Hình 2.2 Sơ đồ tính trường hợp giếng đơn hoàn chỉnh – tầng thấm nước không áp.