Nghiên cứu cấu trúc nền địa chất khu vực xảy ra sụt lún đất tại xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa bằng phương pháp đo ảnh điện đa cực

Bài viết Nghiên cứu cấu trúc nền địa chất khu vực xảy ra sụt lún đất tại xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa bằng phương pháp đo ảnh điện đa cực trình bày việc sử dụng phương pháp thăm đò địa vật lý điện để xác định vị trí, kích thước các hang ngầm tại khu vực này, tạo cơ sở khoa học giúp các nhà quản lý hoạch định chính sách quản lý, quy hoạch hợp lý, nhằm giảm thiểu thiệt hại về người và của cho người dân trong tương lai.

NGHIÊN CỨU CẤU TRÚC NỀN ĐỊA CHẤT KHU VỰC XẢY RA SỤT LÚN ĐẤT TẠI XÃ VĨNH THÀNH, HUYỆN VĨNH LỘC, TỈNH THANH HÓA BẰNG PHƯƠNG PHÁP ĐO ẢNH ĐIỆN ĐA CỰC Đỗ Trọng Quốc1, Phạm Nguyễn Hà Vũ Hồ Tiến Chung2, Trịnh Thị Thúy2, Nguyễn Văn Đông Nguyễn Bá Duẩn3 TÓM TẮT Tai biến sụt lún đất xảy nước ta ngày tăng số lượng hậu tác động Kết nghiên cứu sụt lún đất khu vực xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa phương pháp địa vật lý đo ảnh điện đa cực xác định: 1) Khu vực nghiên cứu nằm địa hình karst với hệ thống hang ngầm lớn liên thông 2) Có thể phân chia địa chất khu vực nghiên cứu thành lớp cấu trúc có dị thường điện trở suất khác nhau, gồm lớp trầm tích bở rời (ρ = 13 ÷ 30 Ωm); lớp phong hóa mạnh (ρ = 30 ÷ 60 Ωm); lớp phong hóa ( = 60 ÷ 150 Ωm) lớp đá tươi (ρ> 150 Ohm.m.) 3) Có tầng hang karst khu vực nghiên cứu gồm: Tầng hang karst mở dạng hàm ếch phổ biến độ sâu 20 ÷ 30 m bị lấp nhét vật liệu trầm tích bở rời; Tầng hang karst độ sâu từ 40 ÷ 60 m, gồm hang kín chứa nước bị lấp nhét vật liệu ngậm nước; Tầng hang kín độ sâu lớn 90 m có kích thước lên đến hàng chục mét Từ khóa: Sụt lún, phương pháp đo ảnh điện đa cực, karst Nhận bài: 26/6/2022; Sửa chữa: 28/6/2022; Duyệt đăng: 30/6/2022 Đặt vấn đề Biến đổi khí hậu gia tăng dân số kéo theo nhu cầu sử dụng nước gia tăng Mực nước ngầm bị hạ thấp khai thác mức dẫn đến nhiều hậu quả, có tai biến, sụt lún đất Tại Việt Nam, năm gần đây, tai biến, sụt lún đất gia tăng số lượng hậu tác động Sụt lún xảy vùng có đa chất yếu, khai thác nước ngầm mức đồng sông Cửu Long [9], đồng sông Hồng tỉnh miền núi, nơi nằm đá vôi Thanh Hóa, Hịa Bình, Sơn La, Hà Giang, Cao Bằng, Lạng Sơn, Nghệ An… Chu trình nước ngầm đá vơi hịa tan CaCO3, tạo nên hang ngầm gọi trình karst [2] Nước ngầm giống chống trần hang, mực nước cao, mực nước ngầm suy giảm tượng sụt lún vị trí hang gần bề mặt có trần yếu xảy Việc xác định vùng có hang ngầm, kích thước hang phân bố chúng vùng núi đá vôi quan trọng công tác dự báo tai biến, giúp người dân giảm thiểu thiệt hại người của, đồng thời, giúp nhà quản lý đề xuất phương án quy hoạch hợp lý Khu vực huyện Vĩnh Lộc Trường Đại học Khoa học tự nhiên - ĐHQGHN Viện Khoa học địa chất khoáng sản Viện Vật lý địa cầu 66 Chuyên đề II, tháng năm 2022 tỉnh Thanh Hóa nằm chủ yếu địa hình karst, từ năm 1983 đến nay, tài trợ Chương trình nước UNICEF xây dựng khoảng 600 giếng khoan công suất lớn khai thác nước ngầm [10] gia tăng đáng kể giếng khoan hộ dân dẫn đến gia tăng tai biến, sụt lún đất thời gian gần đây, đặc biệt khu vực xã Vĩnh Thành Trong nghiên cứu này, chúng tơi sử dụng phương pháp thăm đị địa vật lý điện để xác định vị trí, kích thước hang ngầm khu vực này, tạo sở khoa học giúp nhà quản lý hoạch định sách quản lý, quy hoạch hợp lý, nhằm giảm thiểu thiệt hại người cho người dân tương lai Đối tượng phương pháp nghiên cứu 2.1 Phạm vi đối tượng nghiên cứu Khu vực lựa chọn nghiên cứu thuộc xã Vĩnh Thành, huyện Vĩnh lộc, tỉnh Thanh Hóa, nơi ghi nhận xuất nhiều hố sụt đất liên tiếp xảy khu dân cư cánh đồng [3,5] Khu vực phẳng phong hóa đá vơi nằm kẹp hai sơng lớn sơng Mã phía KẾT QUẢ NGHIÊN CỨU KHOA HỌC VÀ ỨNG DỤNG CÔNG NGHỆ Tây Bắc sơng Bưởi phía Đơng (Hình 1) Theo tài liệu địa chất, khu vực nghiên cứu nằm hệ tầng Đồng Giao (T2a đg) hệ tầng Cò Nòi (T1 cn) đá vôi xám, xám sẫm phân lớp bị dập vỡ mạnh [7] Hệ thống khe nứt, đứt gãy có xu hướng bị căng kéo hướng Tây Bắc Hệ thống hang động thường phát triển dọc theo đới dập vỡ, vậy, tuyến địa vật lý thiết kế cho bắt gặp nhiều đới cấu trúc, nhiều tầng hang tài liệu địa chất, kiến tạo sở cho việc thiết kế tuyến, hay nói cách khác tuyến địa vật lý thiết kế chủ yếu vng góc với hệ thống khe nứt, đứt gãy Theo đó, có tổng cộng 10 tuyến đo địa vật lý điện thiết kế, bao gồm tuyến làng tuyến đồng xung quanh hố sụt cũ Trong đó, tuyến số 1, 1b, 2, 4, 4b, 5a, 5b: Tiến hành đo song song theo hướng Đông Bắc - Tây Nam, chiều dài điện cực a = m a = m; tuyến số 03a, 03b BS: Tiến hành đo theo hướng cắt gần vng góc với tuyến chạy theo hướng Tây Bắc - Đông Nam, chiều dài điện cực a = m (Hình 1) Nguyên lý phương pháp nghiên cứu thay đổi tính chất dẫn điện loại đất đá theo phương thẳng đứng, độ sâu khác theo chiều từ xuống Mỗi loại đất đá có tính chất dẫn điện khác nhau, tính chất xác định qua giá trị điện trở suất (ký hiệu ρ có thứ ngun Ωm) Trong mơi trường địa chất tự nhiên, hang karst tồn hai trạng thái: Hang khơ (hang khơng khí) hang bị lấp nhét (bởi vật liệu bùn sét nước) Trường hợp hang trạng thái khô rỗng, điện trở suất đo cao Ngược lại, hang trạng thái bị lấp nhét, điện trở suất đo nhỏ nhiều so với môi trường xung quanh Chính khác biệt giá trị điện trở suất làm sở cho phương pháp thăm dò địa vật lý điện phát dị thường cục liên quan đến đới nứt nẻ, vị trí hổng hốc hang karst có mặt môi trường địa chất Hình là sơ đồ minh họa các điểm số liệu đo được bằng hệ cực Wenner với hệ thống 28 điện cực Kết quả đo đạc được biểu diễn dưới dạng mặt cắt điện trở suất biểu kiến đo đạc: Trục nằm ngang là khoảng cách các điện cực tuyến; trục thẳng đứng là độ sâu khảo sát tương ứng với cự ly điện cực, độ dài tuyến và loại hệ cực được sử dụng Số liệu sau đã thu thập (mặt cắt điện trở suất biểu kiến đo đạc) cùng thông tin về bề mặt địa hình sẽ được đưa vào phần mềm xử lý 2D chuyên dụng để tính toán sự phân bố điện trở suất thực của môi trường cần nghiên cứu (mặt cắt địa điện) cuối mặt cắt minh giải (mặt cắt địa chất - địa vật lý) ▲Hình Vị trí nghiên cứu sơ đồ tuyến địa vật lý 2.2 Phương pháp nghiên cứu Phương pháp đo ảnh điện đa cực Cho đến nay, phương pháp đo sâu (nghiên cứu sự thay đổi điện trở suất theo phương thẳng đứng) mặt cắt điện (nghiên cứu sự thay đổi điện trở suất theo phương nằm ngang) với cực thường sử dụng rộng rãi để giải toán địa chất chiều (1D) Tuy nhiên, thực tế cần giải toán địa chất hai chiều (2D) phức tạp, nghiên cứu thay đổi điện trở suất môi trường theo hai hướng thẳng đứng nằm ngang, phương pháp ảnh điện đời kết hợp phương pháp đo sâu, phương pháp mặt cắt điện sử dụng cho khảo sát cần độ phân giải cao đặc biệt nghiên cứu tìm kiếm hang động ngầm [4,8] Trong nghiên cứu này, dùng phương pháp sử dụng hệ thống thiết bị SUPERSTING R8/IP + 56, hãng AGI Mỹ trường Đại học Khoa học Tự nhiên ▲Hình Sơ đồ bố trí các điện cực phương pháp ảnh điện và vị trí các điểm số liệu thu mặt cắt điện trở suất biểu kiến đo đạc Phương pháp thu thập tài liệu thực địa Trong khu vực khảo sát bố trí tuyến đo địa vật lý theo mặt cắt dọc ngang Trên tuyến đo bố trí điểm đo sâu điện Tại điểm đo, tiến hành đo máy giá trị cường độ dòng điện (I) qua lưỡng cực phát AB hiệu điện (∆U) lưỡng cực thu MN, sau dịch chuyển tăng kích thước hệ cực với khoảng cách AB = MN = a Các hang độ sâu không lớn khơng chứa nước có khả gây sụt lún Chuyên đề II, tháng năm 2022 67 lại hang độ sâu lớn thường mức độ gây sụt không cao Do vậy, độ sâu nghiên cứu thiết kế để tìm kiếm hang độ sâu khoảng 100 m đổ lại, tùy vị trí, độ mở cực lựa chọn tương ứng a = m a = m Các tuyến khảo sát lấy chi tiết tọa độ cực cắm sở để luận giải, xác định xác vị trí dị thường điện trở suất (hang ngầm) Số liệu đo lưu trữ tự động nhớ máy trút máy tính thực địa để tiến hành bước xử lý kết ▲Hình Sơ đồ đo “cuốn chiếu” thiết bị SuperSting R8/ IP+56 thực địa Phương pháp xử lý số liệu Tài liệu thu tuyến đo mảng số liệu, bao gồm tập hợp giá trị điện trở suất biểu kiến điểm đo, độ sâu khác tuyến qua q trình xử lý, phân tích mà thơng số như: Điện trở suất thực, bề dày lớp đất đá xác định Quá trình xử lý thực sở xây dựng lát cắt đẳng trị điện trở suất Trên tuyến đo trục hoành vị trí điểm đo, trục tung chiều sâu khảo sát Giá trị điện trở suất ghi lát cắt cộng tương hỗ để khử dị thường ảo đưa vị trí thẳng đứng tương ứng với vị trí thực đối tượng theo cơng thức: U Z (Zn ) >U ki = n  U ki  n = n @ Ở đây: i - vị trí điểm đo sâu tuyến đo n - thứ tự dịch chuyển cự ly đo sâu điểm đo ( n1) a Zn = chiều sâu khảo sát Quá trình xử lý phân tích, định lượng theo mơ hình hai chiều (2D) [1,6] Mục đích xác định tham số phần tử toàn mặt cắt cho trường quan sát mặt cắt trùng với trường tính tốn cách tốt nhất, thực chất cực tiểu hóa phiếm hàm độ lệch bình phương trung bình giá trị điện trở suất đo toàn tuyến đo với giá trị tính mơ hình lý thuyết Tồn q trình thực phần mềm xử lý chuyên dụng Res2Dinvi cho hình ảnh lát cắt điện trở suất (ảnh điện) [4] Trên thể độ sâu điện trở suất đối tượng nằm bên dưới, từ cho phép thành lập mặt cắt địa vật lý luận giải, nhận dạng đối tượng 68 Chuyên đề II, tháng năm 2022 Chất lượng số liệu đo đạc đánh giá qua sai số phép đo Sai số chênh lệch lần đo phép đo Từ số liệu đo đạc (mặt cắt điện trở suất biểu kiến đo đạc), dựa vào loại hệ cực đo hệ số thiết bị, chương trình tính điện trở suất tính tốn (mặt cắt điện trở suất biểu kiến tính tốn) Từ số liệu tính tốn, số liệu địa hình độ sâu khảo sát loại hệ cực đo, chương trình tạo mơ hình ban đầu cách phân chia mơi trường thành khối chữ nhật mà vị trí, kích thước chúng dựa vị trí điểm số liệu đo đạc cự ly chúng Sau số vòng lặp so sánh, điều chỉnh (nghịch đảo) lại mơ hình ban đầu, chương trình chọn mơ hình mà giả mặt cắt điện trở suất chúng tạo (tính tốn) trùng với số liệu đo đạc thực tế (mặt cắt nghịch đảo) Từ mặt cắt địa điện, với kiến thức địa chất, kinh nghiệm, thu mặt cắt địa chất tuyến khảo sát Chất lượng trình xử lý đánh giá qua hai loại sai số: Sai số RMS chuẩn L2 (xem mặt cắt) Hai loại sai số 10% hoàn toàn đáng tin cậy để minh giải 10 mặt cắt thu thỏa mãn yêu cầu Để tiện cho việc đối sánh, phân tích tuyến, số liệu xử lý theo lưu đồ, với tham số hiển thị với thang màu điện trở suất Kết nghiên cứu 3.1 Phân chia cấu trúc Dựa vào giá trị thu mặt cắt nghịch đảo, đưa thang điện trở suất (ρ) nằm khoảng từ 13 ÷ 300 Ωm Giá trị phơng chung đại diện cho điện trở suất mặt cắt khu vực nghiên cứu sở để nhận diện đối tượng có giá trị điện trở suất Từ kết thu 10 mặt cắt phân chia địa chất khu vực nghiên cứu thành lớp đất đá sau: (I) Lớp lớp vật liệu bở rời điện trở suất thấp, ρ nằm khoảng 13 ÷ 30 Ohm.m Lớp phủ liên quan đến trầm tích Đệ tứ, có bề dày thay đổi theo tuyến trung bình dày 10 - 12 m (line 1b, 2, 3a, 4a, 5a, 5b) Ngoài ra, số tuyến vật liệu trầm tích cịn lấp nhét vào hang hốc karst hang karst mở nên bề dày lớp lớn, lên đến 20 m line 1a, 3b, BS (Hình8,9,10,11,12) (II) Lớp đá phong hóa mạnh Có giá trị điện trở suất từ 30 ÷ 60 Ohm.m (III) Lớp đá phong hóa vừa, phong hóa Có giá trị điện trở suất từ 60 ÷ 150 Ohm.m (IV) Lớp đá tươi Có giá trị điện trở suất lớn 150 Ohm.m Để kiểm chứng kết luận giải, nhóm tác giả khảo sát hộ dân khoan giếng quanh khu vực KẾT QUẢ NGHIÊN CỨU KHOA HỌC VÀ ỨNG DỤNG CÔNG NGHỆ khảo sát cho thấy, lớp (I) lớp đất bở rời; lớp (II) lớp đá sột, đá non có màu xám đến nâu đỏ; lớp (III) lớp đá cứng có kẽ, người dân thường khoan giếng đến tầng gặp nước; lớp (IV) lớp đá xanh (đá vôi) cứng nguyên khối, khoan nên lấy thỏi đá Các kết khoan trùng với luận giải địa vật lý 3.2 Kết nhận dạng dị thường hang Từ mặt cắt địa điện 10 tuyến đo (Hình1), nhóm tác giả lựa chọn tuyến có phát dị thường điện trở suất (phát hang) bao gồm tuyến 5a, 5b, 4, 3a Cụ thể: Tuyến 3a Tuyến chạy theo hướng Tây Bắc Đơng Nam nằm phía cánh đồng (Hình1) Tuyến có chiều dài 555 m với khoảng cách điện cực m đạt độ sâu tối đa 62 m Trên mặt cắt địa điện (Hình4) nhận dạng đới dị thường giá trị điện trở suất thấp Các dị thường nằm lớp đá vơi phong hóa nứt nẻ mạnh (lớp II III) Cụ thể: + (1) Dị thường thứ xuất từ mét thứ 70 ÷ 120 Vị trí cấu trúc liên quan đến cấu trúc hang độ sâu từ 25 - 33 m rộng khoảng 15 m; + (2) Dị thường thứ từ mét thứ 90 105 Vị trí đới phong hóa mạnh nằm lớp đá phong hóa nên đới dập vỡ chứa nước phát triển trở thành hang; + (3) Dị thường thứ từ mét 385 - 520 mặt cắt địa điện điển hình cho đới bị phong hóa dập vỡ mạnh hình thành nên cấu trúc hai hang sâu từ 22 - 30 m rộng khoảng 30 m ▲Hình Mặt cắt địa điện tuyến line 3A Tuyến line Tuyến line dài 888 m, chạy theo hướng Đông Bắc - Tây Nam với 112 điện cực khoảng cách điện cực m, độ sâu tối đa 99 m Số liệu đo, sau xử lý thu mặt cắt địa vật lý tuyến Line - Hình số Trên mặt cắt nghịch đảo thấy 02 đới dị thường gồm: + (1) Dị thường thứ xuất từ mét thứ 56 224 Trên mặt cắt cho thấy hang mở với cấu trúc dạng hàm ếch đạt độ sâu đến 55 m, gây nên dị thường điện trở suất thấp; đới cấu trúc có khả gây sụt lún cao + (2) Dị thường mét thứ 353 - 367 sâu tới khoảng 25 xuống tới 75 m, đới dập vỡ bị phong hóa mạnh, nằm lớp đá gốc rắn với kích thước 14 x 50 m Tại vị trí hình thành hang tương lai Tuyến 5a Tuyến chạy dọc phía sau khu dân cư ven QL45 theo hướng Đông Bắc Tây Nam, bắt đầu phía trung tâm thị trấn kết thúc cổng làng Hà Lương Tuyến có chiều dài 664 m với khoảng cách điện cực m, đạt độ sâu tối đa 99 m Trên mặt cắt địa điện (Hình6) xác định đới dị thường, cụ thể: ▲Hình Mặt cắt địa điện tuyến line + Dị thường thứ nằm lớp đá vôi phong hóa nứt nẻ mạnh (lớp II III) Dị thường xuất từ mét thứ 115 - 256 Đây đới bị sụt lún, bị lấp nhét vật liệu phía trên, độ sâu đạt đến 60 m; + (2) Dị thường nằm mét thứ 512 - 624 m mặt cắt Vị trí liên quan đến dị thường hang độ sâu từ 57 - 99 m rộng 20 m Phần trần hang đới phong hóa mạnh nên có nguy xảy sụt lún cao tương lai ▲Hình Mặt cắt địa điện tuyến 5a Tuyến 5b Tuyến chạy dọc phía sau khu dân cư ven QL45, phía bên sơng Bưởi theo hướng Đơng Bắc Tây Nam, bắt đầu phía gần ngã tư đường trục làng với QL45, cuối tuyến gần làng Giáng Tuyến có chiều dài 664 m, khoảng cách điện cực m độ sâu 99 m Trên mặt cắt địa điện tuyến Line 5b - Hình số cho thấy có đới dị thường (hang) Các dị thường nằm lớp đá vơi phong hóa nứt nẻ mạnh (lớp II III) Cụ thể: + (1) Dị thường thứ xuất từ mét thứ 223 250 Trên mặt cắt địa điện cho thấy dị thường liên quan đến cấu trúc hang độ sâu tới 90 m; + (2) Dị thường nằm mét thứ 296 - 416 m mặt cắt Vị trí liên quan đến dị thường hang độ sâu từ 52 - 80 m Hang nằm lớp đá tươi có phần trần hang phần bị phong hóa mạnh nên có nguy sụt lún tương lai ▲Hình Mặt cắt địa điện tuyến line 5b Chuyên đề II, tháng năm 2022 69 Một số mặt cắt khác: ▲Hình Mặt cắt địa điện tuyến line ▲Hình Mặt cắt địa điện tuyến line 1b ▲Hình 10 Mặt cắt địa điện tuyến line ▲Hình 11 Mặt cắt địa điện tuyến line 3b ▲Hình 12 Mặt cắt địa điện tuyến line BS 3.3 Kết thảo luận Như vậy, khu vực nghiên cứu phát nhiều hang karst với kích thước lớn bị sập trần, tạo địa hình dạng hàm ếch chưa bị sập (đây đối tượng dễ sinh sụt lún tương lai) Các dị thường (hang Karst) xuất chủ yếu tuyến, gồm line 3a, 4, 5a 5b Dựa vào vị trí xuất hang kích thước chúng mối quan hệ với đới xung quanh, phân thành tầng hang khác khu vực nghiên cứu sau: + Tầng 1: Các hang karst mở dạng hàm ếch phổ biến độ sâu 20 ÷ 30 m (Line 1, 3a, 4), bị lấp nhét vật liệu bên trên; 70 Chuyên đề II, tháng năm 2022 + Tầng 2: Các hang karst độ sâu từ 40 ÷ 60 m (Line 1, 1b, 4, 5a, 5b), bao gồm hang lớn, hang kín hang hở; + Tầng 3: Các hang lớn độ sâu 90 m (line 5a, 5b) bắt gặp tuyến khảo sát Chiều rộng hang karst dao động từ vài mét đến 10 m (Line 4) Về vật chất lấp nhét hang dự kiến lấp nhét phần bùn, sét, hữu trạng thái ngậm nước Căn theo vị trí xuất dị thường (biểu hang karst) tuyến đo, theo quy luật phát triển karst thơng thường hang phát khu vực nghiên cứu có phát triển liên thông với Tuy nhiên, mật độ điểm đo, khoảng cách tuyến đo lớn, nên phát triển liên thông với (ranh giới tầng hang này) cịn mang tính chất dự kiến, ngoại suy Trong thực tế có sai khác (phụ thuộc vào khoảng cách tuyến điểm đo) Đá gốc vùng đá vôi lẫn sét vơi, thành phần có chứa vật chất khác sulphur, thạch anh, than Than sulphur trạng thái ngậm nước vật chất dẫn điện tốt (điện trở suất thấp) Sự có mặt khống vật đối tượng tạo dị thường giả, phần ảnh hưởng đến kết đo địa vật lý, gây nên sai số định Ngồi ra, việc liên thơng tầng hang thực tế cịn có biểu sau: + Tụt nước ngầm giếng khoan hộ dân; + Đục nước giếng xung quanh khoan giếng gần đó; Các tượng vận động thủy văn trên, dịng chảy liên thơng hang karst (nếu có) làm trôi rửa phần vật chất lấp nhét (sang lưu vực khác) hang, gây lượng lún sụt mặt đất khu vực xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa Kết luận - Về phân chia địa tầng: Có thể phân chia địa tầng khu vực nghiên cứu làm lớp tương ứng với giá trị điện trở suất khác nhau, gồm: Lớp đất phủ (lớp I) với điện trở suất thấp (ρ = 13 ÷ 30 Ωm) Đá vơi phong hóa nứt nẻ mạnh (lớp II) vừa (lớp III): Các lớp nằm lớp phủ, có mặt toàn tuyến đo với độ dày mỏng, m Giá trị điện trở suất hai lớp thấp, có giá trị ρ = 30 60 Ωm ρ = 60 ÷ 150 Ωm Đới đá tươi (lớp IV): Có giá trị điện trở suất lớn 150 Ohm.m, nằm sâu phía dưới, nhiên đơi lộ bề mặt tuyến line KẾT QUẢ NGHIÊN CỨU KHOA HỌC VÀ ỨNG DỤNG CÔNG NGHỆ - Về phát triển hang karst khu vực: 1) Khu vực nghiên cứu nằm địa hình karst có hệ thống hang ngầm lớn liên thơng với 2) Có tầng hang karst khu vực nghiên cứu gồm: Tầng 1: Các hang karst mở dạng hàm ếch phổ biến độ sâu 20 ÷ 30 m, bị lấp nhét vật liệu trầm tích bở rời; Tầng 2: Các hang karst độ sâu từ 40 ÷ 60 m, gồm hang kín chứa nước bị lấp nhét vật liệu ngậm nước; Tầng 3: Các hang kín độ sâu lớn 90 m có kích thước lên đến hàng chục mét Kết phát hang karst có mặt vị trí mà tuyến đo cắt qua Mặc dù khoảng cách tuyến điểm đo phù hợp với tiêu chuẩn hành khối lượng cơng việc cịn khiêm tốn so với mức độ phức tạp phát triển karst mặt đất Do đó, thực tế có thêm hang khác nằm ngồi (hoặc nằm giữa) vị trí tuyến đo điểm đo Vì vậy, cần bổ sung thêm khối lượng địa vật lý (tuyến đo điện) cho khu vực nghiên cứu Khu vực xã Vĩnh Thành nói riêng xã lân cận thuộc huyện Vĩnh Lộc nói chung có mức độ phát triển phức tạp hang karst nằm đá vôi hệ tầng Đồng Giao (T2a đg) hệ tầng Cò Nòi (T1 cn) nên tượng sụt lún nước không tránh khỏi Vì vậy, cần thiết phải có nghiên cứu chi tiết rộng (ngoài xã Vĩnh Thành) địa chất - địa vật lý để có sở khoa học phục vụ cảnh báo c tai biến liên quan hoạch định sách quy hoạch gắn với phát triển bền vững địa phương■ TÀI LIỆU THAM KHẢO Advanced Geosciences, Inc., 2009 Intruction Manual of AGI EarthImager 2D, p139 "What is Karst?"  (PDF).  Environmental Science Institute The University of Texas at Austin May 16, 2006 Retrieved  25 December 2020 Đoàn Mỏ - Địa chất Thanh Hóa, 2020 Nhiệm vụ: Điều tra đánh giá tổng thể giải pháp xử lý khu vực có nguy sụt, lún đất xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa Geotomo Software (2007) RES2DINV Rapid 2D resistivity and IP inversion using the least-squares methods User’s manual, p138 UBND huyện Vĩnh Lộc, 2020 Quy hoạch sử dụng đất đến năm 2030 kế hoạch sử dụng đất năm 2021 Jackson, Julia A., ed (1997) "Karst".  Glossary of geology (Fourth ed.) Alexandria, Viriginia: American Geological Institute. ISBN 0922152349 Lê Duy Bách, Đặng Trần Quân, Địa chất khoáng sản tỷ lệ 1: 50.000 nhóm tờ Thanh Hóa ( E-48-IV) Cục Địa chất Khoáng sản Việt Nam , 2001 Loke M.H., 2004 Electrical imaging surveys for environmental and engineering studies: a practical guide to 2-D and 3-D surveys, p61 https://sites ualbe rta.ca/~unswo rth/UA-class es/223/lokecours e-notes pdf Trần Văn Tỷ, Huỳnh Văn Hiệp, 2017 Hiện trạng khai thác nước đất mối tương quan hạ thấp cao độ mực nước sụt lún đất STUDYING THE GEOLOGICAL STRUCTURE IN VINH THANH COMMUNE, VINH LOC DISTRICT, THANH HOA PROVINCE WHERE LAND SUBSIDENCE OCCURRED USING MULTI-ELECTRODE RESISTIVITY IMAGING METHOD Do Trong Quoc, Pham Nguyen Ha Vu VNU University of Science Ho Tien Chung, Trinh Thi Thuy, Nguyen Van Dong Vietnam Institute of Geosciences and Mineral Resources Nguyen Ba Duan Institute of Geophysics ABSTRACT Land subsidence hazards occurring in our country are increasing both in number and impact consequences The results of land subsidence study in Vinh Thanh commune, Vinh Loc district, Thanh Hoa province by multi-electrode resistivity imaging method have determined: 1) The study area is located on karst terrain with large underground caves connected into a system 2) The geology of the study area can be divided into layers with different resistivity anomalies, including unconsolidated sedimentary layer (ρ = 13 ÷ 30 Ωm), strong belt of weathering (ρ = 30 ÷ 60 Ωm), less weathered layer (ρ = 60 ÷ 150 Ωm), and bed rock layer (ρ >150 Ωm); 3) There are floors of karst caves in the study area including: 1st floor: karst caves with collapsed ceilings are common at depths from 20 ÷ 30 m and are filled with unconsolidated sedimentary materials; 2nd floor: Karst caves at a depth from 40 60 m including closed caves containing water or filled with hydrated materials; 3rd floor: including caves at great depths over 90m and up to tens of meters in size Key words: Land subsidence, multi-electrode resistivity imaging method, karst Chuyên đề II, tháng năm 2022 71 ... vật chất lấp nhét (sang lưu vực khác) hang, gây lượng lún sụt mặt đất khu vực xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa Kết luận - Về phân chia địa tầng: Có thể phân chia địa tầng khu vực nghiên. .. 25 December 2020 Đo? ?n Mỏ - Địa chất Thanh Hóa, 2020 Nhiệm vụ: Điều tra đánh giá tổng thể giải pháp xử lý khu vực có nguy sụt, lún đất xã Vĩnh Thành, huyện Vĩnh Lộc, tỉnh Thanh Hóa Geotomo Software... cần nghiên cứu (mặt cắt địa điện) cuối mặt cắt minh giải (mặt cắt địa chất - địa vật lý) ▲Hình Vị trí nghiên cứu sơ đồ tuyến địa vật lý 2.2 Phương pháp nghiên cứu Phương pháp đo ảnh điện đa cực