Tổng quan nghiên cứu
Thành phố Cẩm Phả thuộc tỉnh Quảng Ninh là trung tâm kinh tế - công nghiệp trọng điểm với hoạt động khai thác than và tốc độ đô thị hóa ven biển diễn ra nhanh chóng. Tuy nhiên, giai đoạn từ năm 2010 đến 2016 đã ghi nhận sự gia tăng đột biến của hiện tượng sụt lún đất với diễn biến phức tạp. Thống kê thực tế ghi nhận số lượng hố sụt tăng mạnh qua các năm: năm 2010 xuất hiện 2 hố sụt, năm 2012 có 9 hố sụt, đỉnh điểm năm 2013 ghi nhận tới 19 hố sụt và năm 2014 xảy ra 10 hố sụt. Các hố sụt có đường kính miệng hố phổ biến từ 1,5 đến 4 m, độ sâu sụt lún dao động từ 4 đến 8,5 m, phân bố trên diện rộng gây hư hại hạ tầng kỹ thuật, làm nứt đổ nhà cửa và ảnh hưởng trực tiếp đến đời sống dân sinh.
Vấn đề nghiên cứu trọng tâm là làm rõ mối quan hệ tương tác giữa cấu trúc địa chất karst ngầm, đặc tính cơ lý tầng phủ bở rời và sự mất cân bằng động lực nước dưới đất do tác động nhân sinh. Mục tiêu cụ thể của luận văn tập trung vào việc xác định điều kiện, nguyên nhân, cơ chế hình thành tai biến sụt đất tại dải ven biển Cẩm Phả; từ đó ứng dụng mô hình chỉ tiêu tích hợp để phân vùng nguy cơ sụt đất trên bản đồ tỷ lệ 1/25.000.
Phạm vi nghiên cứu bao quát các khu vực trọng điểm thuộc phường Cẩm Đông, Cẩm Sơn, Cẩm Phú và Cẩm Tây. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp hệ thống cơ sở dữ liệu định lượng, giúp các nhà quản lý quy hoạch xây dựng công trình an toàn và chủ động giảm thiểu rủi ro tai biến địa chất trong phát triển đô thị bền vững.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Luận văn xây dựng trên nền tảng lý thuyết phát triển karst và tai biến xói ngầm ngầm hóa của các nhà khoa học trên thế giới như Waltham, Fookes và Williams. Karst chiếm khoảng 15% diện tích bề mặt Trái Đất; tại Việt Nam, diện tích đá carbonat đạt tới 50.000 km² (tương đương khoảng 20% tổng diện tích tự nhiên). Quá trình hình thành hố sụt karst ven biển đòi hỏi sự hội tụ của ba điều kiện cốt lõi:
- Sự tồn tại của khoảng không gian ngầm karst dạng hang hốc, đới dập vỡ kiến tạo trong đá móng carbonat.
- Sự hiện diện của lớp đất tầng phủ xung yếu có tính bất đồng nhất cao, dễ bị rửa trôi khi bão hòa nước.
- Tác nhân kích hoạt từ sự biến động mực nước dưới đất kết hợp tải trọng tĩnh và động trên bề mặt.
Mô hình chỉ tiêu tích hợp các yếu tố phát triển tai biến được ứng dụng làm khung tính toán xuyên suốt. Bản chất của mô hình là xác định hàm mục tiêu $K_y = f(a_1, a_2, ..., a_n)$ đại diện cho cường độ hoặc mật độ tai biến, sau đó chuẩn hóa các tham số địa chất công trình - địa chất thủy văn về cùng thứ nguyên (thang đo từ 0 đến 1). Chỉ số tích hợp $I_{\sum}$ tại từng ô lưới không gian được tính toán theo công thức:
$$I_{\sum} = \sum_{i=1}^{n} g_i R_{iH}$$
Trong đó $g_i$ là tỷ trọng ảnh hưởng của yếu tố thứ $i$ và $R_{iH}$ là giá trị lượng hóa đã chuẩn hóa của yếu tố đó.
Phương pháp nghiên cứu
Nghiên cứu sử dụng tổ hợp các phương pháp khảo sát và phân tích hiện đại nhằm đảm bảo độ tin cậy khoa học:
- Nguồn dữ liệu và cỡ mẫu khảo sát: Thu thập và xử lý toàn diện số liệu từ 76 lỗ khoan địa chất công trình, hồ sơ đo đạc thực tế của 22 vị trí sụt đất điển hình trên địa bàn thành phố Cẩm Phả. Dữ liệu động thái nước ngầm được quan trắc tự động tại 3 trạm đo chuyên dụng với tần suất 15 phút/lần liên tục trong 6 tháng (từ tháng 7/2014 đến tháng 1/2015).
- Phương pháp chọn mẫu: Chọn mẫu đất thí nghiệm cơ lý theo phương pháp mẫu nguyên dạng tại các độ sâu từ 2 m đến 18 m; khảo sát thực địa đối chứng toàn diện các khu vực trọng điểm Tây Nam Đá Mài, Thủy Sơn và Quảng Hồng vào tháng 7/2017.
- Phương pháp phân tích: Ứng dụng hệ thống thông tin địa lý (GIS) để số hóa bản đồ chuyên đề, thành lập mô hình trường biến đổi (MHTBĐ) dạng raster cho từng lớp thông tin. Phân tích thống kê đa biến được lựa chọn để thiết lập ma trận tương quan cặp đôi ($r_{ij}$) và giải hệ phương trình xác định hệ số tiêu chuẩn hóa ($\beta$). Việc chọn mô hình chỉ tiêu tích hợp giúp triệt tiêu tính định tính chủ quan, đảm bảo hệ số tương quan nhiều chiều đạt $R > 0,7$, đáp ứng đầy đủ tiêu chuẩn kiểm định thống kê địa kỹ thuật.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Kết quả khoan thăm dò và phân tích tính chất cơ lý đất đá đã xác lập cấu trúc địa tầng khu vực nghiên cứu gồm 5 phân lớp rõ rệt:
- Lớp 1 (Đất lấp bề mặt): Nguồn gốc đất đá thải mỏ và phế thải xây dựng, thành phần hỗn tạp, chiều dày dao động từ 2,5 m đến 7,1 m (trung bình 3,5 m). Đất có trạng thái rời rạc, chỉ số xuyên tiêu chuẩn SPT đạt $N_{30} = 6$, cường độ kháng nén tiêu chuẩn $R_0 = 1,2\text{ kG/cm}^2$, mô đun biến dạng $E_0 = 60\text{ kG/cm}^2$. Tầng này chứa nước dưới đất nông cách mặt đất chỉ từ 0,5 m đến 2 m.
- Lớp 2 (Cát lẫn sạn sỏi và bụi sét): Phân bố cục bộ ở độ sâu 3-4 m, chiều dày trung bình 1,0 m, góc ma sát trong khi ướt $\varphi = 23^\circ 20'$, chỉ số SPT trung bình đạt 5 búa.
- Lớp 3 (Sét pha màu xám đen lẫn sạn sỏi): Trạng thái dẻo cứng, lực dính kết tiêu chuẩn $C = 0,24\text{ kG/cm}^2$, góc ma sát trong $\varphi = 12^\circ 23'$, độ ẩm tự nhiên $W = 24,6%$.
- Lớp 4 (Sét pha màu xám vàng lẫn sạn sỏi - Lớp đất nhạy cảm): Chiều dày biến thiên mạnh từ 2,1 m đến 18 m, phủ trực tiếp trên đá gốc carbonat. Lớp này có độ ẩm tự nhiên cao $W = 40,58%$, hệ số rỗng $e_0 = 1,2$, độ sệt $I_s = 1,7$ (trạng thái dẻo mềm đến dẻo chảy), lực dính kết rất thấp $C = 0,06\text{ kG/cm}^2$, góc ma sát trong $\varphi = 6^\circ$. Hệ số bất đồng nhất thành phần hạt dao động từ 30 đến 60, là tầng đất xung yếu nhất dễ bị xói ngầm.
- Lớp 5 (Đá vôi carbonat hệ tầng Bắc Sơn): Nứt nẻ mạnh, tồn tại các hang rỗng ngầm có chiều cao từ 2,5 m đến 8,0 m nằm ở độ sâu phổ biến từ -8 m đến -16 m.
Về quy luật xuất hiện, hơn 85% các vụ sụt đất xảy ra trong mùa mưa (từ tháng 5 đến tháng 9), với tháng 7 là đỉnh điểm (chiếm hơn 50% tổng số vụ). Phân tích chuỗi quan trắc cho thấy mực nước ngầm dao động theo chu kỳ 38-41 ngày. Tại thời điểm sụt đất Thủy Sơn ngày 29/7/2014, mực nước ngầm trước đó 1 ngày đã sụt giảm đột ngột từ 24 cm đến 31 cm, trùng khớp đợt mưa lớn đạt 77 mm/ngày và thời kỳ triều kiệt có biên độ triều thấp nhất trong tháng.
Thảo luận kết quả
Cơ chế gây sụt đất tại Cẩm Phả được xác định là quá trình biến dạng động lực phức hợp. Khi thủy triều rút kiệt kết hợp nước ngầm thoát nhanh ra biển sau mưa lớn, gradien thủy lực tăng đột ngột tạo dòng thấm xói ngầm cuốn trôi lớp sét pha xám vàng dẻo chảy (Lớp 4) vào các hang hốc đá vôi Lớp 5. Quá trình này hình thành các vòm rỗng ngầm phát triển dần lên phía trên. Dưới tác dụng của tải trọng tĩnh từ lớp đất đắp Lớp 1 dày tới 7,1 m và tải trọng động từ phương tiện giao thông trên Quốc lộ 18, vòm đất mất ổn định và sụp đổ tức thì dạng hình phễu hoặc hàm ếch.
Dữ liệu tính toán từ ma trận hệ số tương quan cặp đôi cho thấy 3 yếu tố đóng vai trò chi phối lớn nhất trong mô hình tích hợp gồm:
- Bề dày lớp đất nhạy cảm (tỷ trọng $g_i$ chiếm trên 28%).
- Độ cao bề mặt đá vôi carbonat (chiếm khoảng 22%).
- Khoảng cách đến đới đứt gãy kiến tạo hoạt động (chiếm khoảng 19%).
Bản đồ phân vùng nguy cơ sụt đất thể hiện rõ sự phân hóa không gian: vùng nguy cơ rất cao và cao phân bố tập trung thành dải kéo dài theo phương Tây Tây Bắc - Đông Đông Nam, chạy dọc chân đồi Núi Nhện qua các khu vực dân cư Thủy Sơn, Tây Nam Đá Mài và Cẩm Đông.
Đề xuất và khuyến nghị
Dựa trên kết quả phân vùng nguy cơ tai biến, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật và quản lý cụ thể:
- Thiết lập mạng lưới quan trắc cảnh báo sớm tự động: Ủy ban nhân dân thành phố Cẩm Phả phối hợp cùng Sở Tài nguyên và Môi trường tỉnh Quảng Ninh tiến hành lắp đặt bổ sung 8 trạm quan trắc tự động mực nước ngầm và áp lực nước lỗ rỗng tại các phường Cẩm Sơn, Cẩm Đông trước quý IV/2018. Hệ thống cần cài đặt ngưỡng cảnh báo tự động khi mực nước ngầm hạ thấp vượt quá 20 cm trong vòng 24 giờ vào mùa mưa bão.
- Kiểm soát quy hoạch tải trọng san lấp và cốt nền đô thị: Phòng Quản lý Đô thị thành phố Cẩm Phả cần giới hạn chiều dày lớp đất đắp mặt bằng không vượt quá 3,5 m tại các khu vực thuộc vùng nguy cơ cao. Đối với 100% công trình xây dựng mới trên nền đất có Lớp 4 dày trên 5 m, bắt buộc phải áp dụng giải pháp móng cọc xuyên qua hang karst ngầm tựa vào đá gốc ổn định hoặc phụt vữa gia cố nền đáy móng.
- Kiên cố hóa và chỉnh trị hệ thống tiêu thoát nước mặt: Ban Quản lý Dự án Đầu tư Xây dựng thành phố tiến hành bê tông hóa toàn bộ các tuyến mương dẫn dòng phương Bắc - Nam đổ ra biển, ngăn chặn hiện tượng rò rỉ nước mặt ngấm vào các đới nứt nẻ kiến tạo, hoàn thành trong giai đoạn 2018-2020 nhằm giảm 90% nguy cơ hình thành phễu xói ngầm nhân tạo.
- Ban hành quy chế quản lý xây dựng theo bản đồ phân vùng nguy cơ: Ủy ban nhân dân tỉnh Quảng Ninh tích hợp bản đồ phân vùng nguy cơ sụt đất tỷ lệ 1/25.000 vào quy hoạch chung xây dựng đô thị Cẩm Phả đến năm 2030, thiết lập hành lang an toàn cách mép các đứt gãy hoạt động tối thiểu 50 m.
Đối tượng nên tham khảo luận văn
Luận văn mang lại giá trị khoa học và ứng dụng thực tiễn cao cho nhiều nhóm đối tượng:
- Cơ quan quản lý nhà nước và quy hoạch đô thị: Sở Xây dựng Quảng Ninh, Ủy ban nhân dân thành phố Cẩm Phả sử dụng bản đồ phân vùng để phê duyệt các đồ án quy hoạch sử dụng đất, cấp phép xây dựng công trình ngầm và hạ tầng giao thông an toàn.
- Doanh nghiệp khai thác khoáng sản và chủ đầu tư bất động sản: Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV) và các chủ dự án khu đô thị ven biển tham khảo để kiểm soát việc đổ đất đá thải mỏ, hạn chế gia tải gây mất ổn định địa tầng.
- Kỹ sư địa kỹ thuật và đơn vị tư vấn thiết kế: Các kỹ sư công trình ứng dụng tập chỉ tiêu cơ lý 5 lớp đất (đặc biệt là chỉ số biến dạng $E_0=60\text{ kg/cm}^2$ và độ sệt $I_s=1,7$) để tính toán kết cấu móng thích ứng cho từng vị trí cụ thể.
- Học viên, nghiên cứu sinh và giảng viên chuyên ngành: Các nhà nghiên cứu thuộc lĩnh vực Địa chất công trình, Địa chất môi trường và Địa kỹ thuật sử dụng làm tài liệu tham khảo phương pháp luận kết hợp GIS với mô hình chỉ tiêu tích hợp trong đánh giá tai biến địa chất đô thị.
Câu hỏi thường gặp
1. Nguyên nhân cốt lõi gây sụt lún đất đột ngột tại thành phố Cẩm Phả là gì?
Sụt đất tại Cẩm Phả xảy ra do sự tương tác giữa cấu trúc hang ngầm trong đá vôi carbonat hệ tầng Bắc Sơn, sự tồn tại của lớp sét pha xám vàng dẻo chảy nhạy cảm (Lớp 4) và sự thay đổi đột ngột của mực nước dưới đất do triều kiệt kết hợp lượng mưa lớn.
2. Tại sao tai biến sụt đất thường tập trung chủ yếu vào tháng 7 hàng năm?
Tháng 7 là thời điểm mưa bão tập trung cao nhất trong năm (lượng mưa đo được đạt tới 77 mm/ngày). Khi mưa lớn xuất hiện trùng với thời kỳ biên độ thủy triều rút sâu nhất, gradien thủy lực ngầm tăng vọt làm gia tăng áp lực xói trôi đất phủ vào hang karst, kích hoạt sụt lún.
3. Lớp đất số 4 trong địa tầng khu vực có đặc điểm gì khiến nó trở nên nhạy cảm?
Lớp sét pha xám vàng có hệ số bất đồng nhất thành phần hạt cao từ 30 đến 60, độ ẩm tự nhiên $40,58%$, độ sệt $I_s = 1,7$ (trạng thái dẻo chảy) và lực dính kết rất yếu ($C = 0,06\text{ kG/cm}^2$), khiến đất nhanh chóng bị hóa lỏng và chảy trôi khi có dòng nước dịch chuyển.
4. Mô hình chỉ tiêu tích hợp có ưu thế vượt trội nào so với các phương pháp truyền thống?
Mô hình định lượng hóa chính xác tỷ trọng của từng yếu tố phát sinh tai biến thông qua ma trận tương quan nhiều chiều ($R > 0,7$), loại bỏ tính chủ quan của phương pháp chuyên gia truyền thống và cho phép phân vùng nguy cơ chi tiết theo từng ô lưới không gian GIS.
5. Người dân và chính quyền cần làm gì ngay khi phát hiện dấu hiệu nứt đất, nứt tường?
Cần khoanh vùng phong tỏa hiện trường trong bán kính tối thiểu 20 m, sơ tán người và tài sản khỏi vùng ảnh hưởng, thông báo ngay cho cơ quan chức năng để kiểm tra mực nước ngầm và áp dụng biện pháp phụt vữa xi măng bịt lỗ rỗng đáy móng, tuyệt đối không tự ý san lấp bằng đất bở rời.
Kết luận
- Luận văn đã làm sáng tỏ bản chất tai biến sụt đất ven biển Cẩm Phả là quá trình xói ngầm cơ học trong tầng sét pha dẻo chảy phủ trên nền đá vôi nứt nẻ chứa hang karst.
- Hệ thống hóa cấu trúc địa chất công trình 5 lớp, nhận diện chính xác "lớp đất nhạy cảm" (Lớp 4) có độ sệt $I_s = 1,7$ là tầng phát sinh tai biến trực tiếp.
- Xác lập mối liên hệ định lượng giữa dao động hạ mực nước ngầm từ 24 đến 31 cm với thời kỳ mưa lớn và triều kiệt kích hoạt sụt đất đột ngột.
- Ứng dụng thành công mô hình chỉ tiêu tích hợp trên nền tảng GIS để thành lập bản đồ phân vùng nguy cơ sụt đất tỷ lệ 1/25.000 với hệ số tương quan tin cậy cao.
- Đề xuất đồng bộ 4 nhóm giải pháp kỹ thuật công trình và quản lý quy hoạch đô thị giúp giảm thiểu thiệt hại và bảo đảm an toàn hạ tầng ven biển.
Để nâng cao hiệu quả phòng chống tai biến địa chất, các cơ quan ban ngành và đơn vị nghiên cứu cần tiếp tục mở rộng mạng lưới quan trắc tự động và ứng dụng công nghệ viễn thám radar giao thoa (InSAR) để theo dõi dịch chuyển bề mặt liên tục theo thời gian thực. Hãy tải toàn văn luận văn để tiếp cận toàn bộ cơ sở dữ liệu địa kỹ thuật chi tiết và bản đồ phân vùng nguy cơ sụt đất Cẩm Phả!