Tổng quan nghiên cứu

Khu vực tỉnh Cao Bằng sở hữu vị trí địa chính trị và kinh tế đặc biệt quan trọng tại vùng Đông Bắc Việt Nam với tổng diện tích tự nhiên đạt 6.724,6 km² cùng đường biên giới dài 314 km giáp ranh với Trung Quốc. Đây là địa bàn giàu tiềm năng tài nguyên khi tập trung 142 mỏ và điểm quặng với 22 loại khoáng sản khác nhau, nổi bật là trữ lượng quặng sắt đạt khoảng 50-60 triệu tấn và quặng bauxit ước tính 200 triệu tấn. Về mặt địa chất động lực, Cao Bằng nằm trọn trong khối kiến tạo Cao Bằng – Móng Cái, chịu tác động trực tiếp từ trạng thái ứng suất siết ép mạnh mẽ từ khối lục địa phía Bắc, dẫn đến sự hình thành nhiều hệ thống đứt gãy sâu đang hoạt động.

Vấn đề nghiên cứu trọng tâm xuất phát từ thực trạng các tai biến địa chất như nứt đất, sạt lở và đặc biệt là động đất thường xuyên xảy ra với mức độ phức tạp ngày càng tăng. Điển hình là trận động đất có độ lớn 5,4 độ Richter xảy ra vào ngày 25/11/2019 tại khu vực Trùng Khánh đã gây nứt sụt đất nghiêm trọng, đe dọa an toàn dân sinh của 28 dân tộc trong tỉnh. Mục tiêu cụ thể của công trình là xác định chi tiết các đặc điểm cấu trúc sâu của vỏ Trái đất, bao gồm độ sâu bề mặt móng kết tinh, bề mặt Moho và bản chất hình học của hệ thống đứt gãy thông qua việc xử lý định lượng số liệu dị thường trọng lực.

Phạm vi không gian nghiên cứu được giới hạn trong khung tọa độ địa lý từ 105°00’ đến 107°00’ kinh độ Đông và từ 22°00’ đến 23°30’ vĩ độ Bắc (tương ứng tọa độ UTM X: 500000 – 704229, Y: 2432827 – 2600293). Nghiên cứu kế thừa dữ liệu đo đạc trọng lực tỷ lệ 1:500.000 kết hợp tỷ lệ 1:200.000 hoàn thiện vào năm 2022. Ý nghĩa thực tiễn của công trình thể hiện qua việc nâng cao độ chính xác giải đoán cấu trúc địa chất, kiểm soát sai số mô hình hóa dưới 5%, cung cấp luận cứ khoa học tin cậy phục vụ công tác dự báo vùng nguồn động đất cực đại và quy hoạch khai thác khoáng sản 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

Nghiên cứu ứng dụng tổ hợp lý thuyết trường thế địa vật lý hiện đại nhằm phân tích định lượng trường trọng lực. Lý thuyết giải chập Euler (Euler Deconvolution) dựa trên phương trình thuần nhất Euler được vận dụng để xác định vị trí biên và tọa độ nguồn gây dị thường (x₀, y₀, z₀) với chỉ số cấu trúc N biến thiên từ 0 đến 2 cho trường trọng lực. Phương pháp này cho phép nhận dạng nhanh các ranh giới mật độ mà không đòi hỏi giả thiết khắt khe về hình học nguồn.

Lý thuyết nghịch đảo mật độ có giới hạn theo phương pháp Guillen và Menichetti (1984) được triển khai nhằm chia cắt môi trường hai chiều thành hệ thống lăng trụ chữ nhật, giải bài toán tối ưu hóa ma trận nhân với hệ số giảm trừ để xác định mật độ dư của các khối đất đá. Bên cạnh đó, lý thuyết nghịch đảo ba chiều miền tần số của Oldenburg (1974) trên nền tảng thuật toán Parker thiết lập mối quan hệ phi tuyến giữa biến đổi Fourier của dị thường trọng lực với chuỗi lũy thừa độ sâu ranh giới mật độ. Các khái niệm then chốt được chuẩn hóa bao gồm: dị thường trọng lực Bouguer, gradient ngang cực đại, góc nghiêng trường thế (Tilt Angle), mặt móng kết tinh (ranh giới giữa tầng trầm tích và móng đá kết tinh) và mặt ranh giới Moho (ranh giới phân chia vỏ Trái đất và Manti trên).

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

Cơ sở dữ liệu đầu vào của luận văn được tổng hợp từ nguồn bản đồ dị thường trọng lực Bouguer tỷ lệ 1:500.000 của Tổng cục Địa chất và Khoáng sản Việt Nam (2011), bổ sung tài liệu đo trọng lực chi tiết tỷ lệ 1:200.000 từ đề tài nghiên cứu địa chấn mã số VAST05.06/21-22. Mô hình độ cao số độ phân giải cao ASTER-GDEM 30 m được tích hợp để hiệu chỉnh ảnh hưởng địa hình. Dữ liệu địa chấn gồm vị trí và độ lớn của các trận chấn động lịch sử do Viện Vật lý Địa cầu quan trắc đến năm 2020 được dùng làm thông tin tiên nghiệm.

Về phương pháp chọn mẫu và quy mô phân tích, nghiên cứu thiết lập 7 tuyến khảo sát 2D cắt vuông góc với phương cấu trúc địa chất chính Tây Bắc – Đông Nam, sử dụng bước lấy mẫu đồng đều 2.000 m trên mỗi tuyến, đảm bảo bao quát toàn diện các đới cấu trúc dị thường trên diện tích hơn 6.700 km².

Quy trình phân tích được thực hiện theo chuỗi công nghệ khép kín:

  • Tính chuyển nâng trường trọng lực lên các độ cao 10 km và 20 km nhằm triệt tiêu trường cục bộ nông, làm nổi bật trường khu vực từ các cấu trúc sâu.
  • Kết hợp phương pháp gradient ngang cực đại và phương pháp góc nghiêng để truy vết chính xác các đới đứt gãy kiến tạo.
  • Xây dựng mô hình 2D ban đầu từ kết quả giải chập Euler và nghịch đảo mật độ, sau đó tiến hành giải ngược 2D trên mô-đun GM-SYS thuộc phần mềm chuyên dụng Oasis Montaj (Geosoft).
  • Chuẩn hóa kết quả 2D để xây dựng lưới không gian 3D, tiếp tục giải bài toán ngược 3D bằng thuật toán miền tần số để xác định chính xác bình đồ độ sâu mặt móng kết tinh và mặt Moho với thời gian tính toán tối ưu và triệt tiêu tính đa trị.

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

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

Thứ nhất, nghiên cứu đã khoanh định và chuẩn hóa hệ thống đứt gãy sâu khu vực tỉnh Cao Bằng và lân cận theo hai hệ thống phương chính: phương Tây Bắc – Đông Nam và phương cánh cung Đông Bắc. Đới đứt gãy Cao Bằng – Tiên Yên thể hiện quy mô lớn nhất với chiều dài vượt 250 km trên lãnh thổ Việt Nam, bề rộng đới phá hủy kiến tạo dao động từ 3-5 km và mở rộng lên đến 15 km ở đoạn giữa Cao Bằng và Thất Khê. Mặt trượt đứt gãy đoạn Hà Quảng – Lạng Sơn cắm về phía Tây Nam với góc dốc lớn từ 60° đến 70°, trong khi đoạn Lạng Sơn – Tiên Yên cắm về phía Đông Bắc với góc dốc 80° đến 85°. Các đới đứt gãy khác như Yên Minh – Phú Lương (dài trên 400 km) và Văn Sơn – Hà Giang (dài 250 km) cũng được làm sáng tỏ về mặt ranh giới phân bố không gian.

Thứ hai, sơ đồ độ sâu bề mặt móng kết tinh xác định được sự phân dị mạnh mẽ của địa hình móng cổ. Độ sâu mặt móng biến thiên trong khoảng từ 0,5 km tại các khối nâng uốn nếp đến trên 4,5 km tại các bồn trũng sụt lún Kainozoi. Bồn trũng Cao Bằng được xác định là cấu trúc dạng kéo toạc (pull-apart) kéo dài gần 20 km theo phương Tây Bắc – Đông Nam với bề dày tầng trầm tích đạt từ 900 m đến 1.000 m.

Thứ ba, bình đồ độ sâu mặt ranh giới Moho khu vực nghiên cứu được xác lập nằm trong khoảng từ 28,0 km đến 34,5 km, đạt giá trị trung bình khoảng 31,5 km. Chiều dày vỏ Trái đất có xu hướng mỏng dần tại các đới đứt gãy phân đới sâu và dày lên về phía các khối uốn nếp nâng cao ở khu vực Tây Bắc.

Thứ tư, thông số mật độ vật chất của các lớp vỏ được chuẩn hóa thông qua giải bài toán ngược: lớp vỏ trầm tích có mật độ trung bình đạt 2,62 g/cm³ (biến thiên 2,54 – 2,65 g/cm³), lớp vỏ kết tinh đạt 2,83 g/cm³ (biến thiên 2,80 – 2,85 g/cm³) và tầng Manti trên đạt 3,15 g/cm³ (biến thiên 3,10 – 3,23 g/cm³).

Thảo luận kết quả

Sự tương quan giữa các phát hiện địa vật lý và đặc điểm kiến tạo cho thấy tính phân dị sâu sắc của vỏ Trái đất Cao Bằng hình thành do lịch sử tiến hóa địa chất đa kỳ, từ pha tạo núi Indosini đến các pha biến dạng kéo toạc Kainozoi. Dọc đới đứt gãy Cao Bằng – Tiên Yên, tốc độ dịch chuyển ngang thời kỳ Tân kiến tạo ghi nhận từ 0,1-0,7 mm/năm và dịch chuyển đứng đạt 0,1-0,5 mm/năm.

Trong phân tích dữ liệu, các kết quả mô hình hóa 2D trên 7 mặt cắt tuyến được thể hiện trực quan qua các biểu đồ so sánh đường cong dị thường trọng lực thực tế và đường cong tính toán. Mức sai số bình phương trung bình trên toàn bộ 7 tuyến đều duy trì ở mức nhỏ hơn hoặc bằng 5% so với biên độ dị thường quan sát, đáp ứng tiêu chuẩn độ tin cậy theo các quy chuẩn địa vật lý quốc tế. Kết quả giải ngược 3D dạng bản đồ đẳng sâu mặt móng kết tinh và mặt Moho phản ánh tính liên tục không gian vượt trội so với các lát cắt 2D đơn lẻ.

So với các nghiên cứu trước đây chỉ sử dụng phương pháp thống kê khu vực hoặc phân tích hàm nhận địa chấn tại một số trạm thưa thớt (như mạng lưới 23 điểm đo trên toàn miền Bắc), công trình này đã khắc phục hoàn toàn hiện tượng phản ánh giá trị trung bình hóa. Việc kết hợp giải chập Euler và giải ngược mật độ để tạo mô hình xuất phát cho bài toán 3D đã triệt tiêu hiệu quả tính đa nghiệm, làm sáng tỏ mối liên hệ mật thiết giữa các chấn tiêu động đất có độ lớn M = 5,0 – 5,4 với các nút giao cắt của đứt gãy bậc cao.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp và khuyến nghị ứng dụng thực tiễn:

  1. Ứng dụng sơ đồ đứt gãy sâu vào phân vùng vi địa chấn và quy hoạch hạ tầng giao thông: Sở Xây dựng và Sở Giao thông Vận tải tỉnh Cao Bằng cần đưa bản đồ hệ thống đứt gãy tỷ lệ 1:200.000 vào quy chuẩn thẩm định các công trình công cộng cấp I và cấp II. Mục tiêu đặt ra là kiểm soát 100% các dự án hạ tầng lớn, tuyến đường huyết mạch (như Quốc lộ 4 và các tuyến cao tốc kết nối cửa khẩu) tránh xây dựng trực tiếp trên đới phá hủy động lực rộng 3-15 km của đứt gãy Cao Bằng – Tiên Yên, hoàn thành cập nhật quy hoạch trước năm 2025.

  2. Nâng cấp mạng lưới trạm quan trắc địa chấn địa phương: Viện Vật lý Địa cầu phối hợp với Ủy ban Nhân dân tỉnh Cao Bằng cần xúc tiến lắp đặt thêm 5 đến 8 trạm địa chấn kỹ thuật số băng rộng tại các khu vực ghi nhận hoạt động kiến tạo mạnh như huyện Trùng Khánh, Hà Quảng và Nguyên Bình. Mục tiêu nâng cao khả năng phát hiện các vi chấn có độ lớn M ≥ 1,5, hoàn thành triển khai trong giai đoạn 2024-2026 nhằm giám sát liên tục tích tụ ứng suất kiến tạo.

  3. Định hướng công tác thăm dò khoáng sản ẩn sâu: Tổng cục Địa chất và Khoáng sản cùng các doanh nghiệp khai khoáng cần áp dụng mô hình ranh giới móng kết tinh và các thể xâm nhập magma (như khối granitoit Pia Oắc diện tích 20 km², khối Phia Bioc) để tái định vị mục tiêu khoan sâu. Tập trung thăm dò các thân quặng thiếc – wolfram, mangan và bauxit ở độ sâu từ 500 m đến 2.000 m, phấn đấu gia tăng 20-30% trữ lượng khoáng sản xác minh trong kế hoạch giai đoạn 2024-2030.

  4. Xây dựng hệ thống cảnh báo sớm trượt lở và nứt đất theo thời gian thực: Ban Chỉ huy Phòng chống Thiên tai và Tìm kiếm Cứu nạn tỉnh Cao Bằng cần tích hợp dữ liệu độ dốc ASTER-GDEM 30 m và mật độ đứt gãy nông để xây dựng bản đồ số cảnh báo sạt lở đất. Đảm bảo cung cấp cảnh báo trước 24 đến 48 giờ trong mùa mưa bão (từ tháng 5 đến tháng 8 hàng năm) cho các cộng đồng dân cư sinh sống tại các thung lũng đá vôi và sườn núi dốc.

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

  1. Các nhà nghiên cứu và chuyên gia Địa vật lý, Địa chất cấu trúc: Cung cấp quy trình tích hợp phương pháp giải chập Euler, nghịch đảo mật độ 2D và mô hình hóa 3D trên phần mềm Geosoft Oasis Montaj, ứng dụng trực tiếp cho các đề tài nghiên cứu cấu trúc sâu thạch quyển tại các khu vực đồi núi phức tạp.

  2. Cơ quan quản lý nhà nước và nhà hoạch định chính sách đô thị: Cung cấp cơ sở dữ liệu địa chất tin cậy cho các Sở Tài nguyên và Môi trường, Sở Xây dựng khu vực Đông Bắc Bộ nhằm đánh giá rủi ro tai biến tự nhiên, cấp phép xây dựng các công trình dân sinh và hồ đập thủy điện an toàn.

  3. Kỹ sư thăm dò địa chất và doanh nghiệp khai thác mỏ: Hỗ trợ các kỹ sư khai thác khoáng sản nắm bắt quy luật phân bố của các đới đứt gãy dẫn magma, khoanh vùng chính xác các tích tụ quặng kim loại quý hiếm, tối ưu hóa chi phí và vị trí các mũi khoan thăm dò.

  4. Giảng viên, học viên cao học và sinh viên ngành Kỹ thuật Địa vật lý: Nguồn tài liệu học thuật giá trị phục vụ công tác giảng dạy và học tập các học phần: Phương pháp Trọng lực – Từ, Minh giải tài liệu địa vật lý, Xử lý số liệu trường thế và Mô hình hóa cấu trúc vỏ Trái đất.

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

Nghiên cứu cấu trúc vỏ Trái đất khu vực tỉnh Cao Bằng có tính cấp thiết như thế nào? Cao Bằng nằm tại vị trí biên giới trọng yếu dài 314 km, tập trung 142 mỏ quặng và chịu ứng suất siết ép kiến tạo mạnh. Trận động đất 5,4 độ Richter năm 2019 tại Trùng Khánh gây sụt nứt đất nghiêm trọng đã chỉ rõ tính cấp bách của việc xác định các cấu trúc đứt gãy sâu và tầng sinh chấn để bảo vệ tính mạng dân cư và công trình hạ tầng.

Ưu điểm chính của việc kết hợp giải bài toán ngược 2D và 3D tài liệu trọng lực là gì? Phương pháp trọng lực vốn có tính đa trị cao. Việc giải ngược 2D trên 7 tuyến khảo sát chuẩn giúp xác lập mô hình ban đầu chính xác, sau đó áp dụng giải ngược 3D miền tần số giúp chuẩn hóa độ sâu móng kết tinh và Moho trên toàn diện tích, khống chế sai số dưới 5% so với tài liệu thực tế.

Bề mặt móng kết tinh và mặt Moho tại Cao Bằng phân bố ở độ sâu bao nhiêu? Theo kết quả mô hình hóa 3D, bề mặt móng kết tinh biến thiên từ 0,5 km tại các khối nâng uốn nếp đến trên 4,5 km tại các bồn trũng Kainozoi (như trũng Cao Bằng dày 900-1000 m). Bề mặt Moho dao động trong khoảng từ 28,0 km đến 34,5 km, phản ánh đặc trưng vỏ lục địa uốn nếp cổ.

Làm thế nào để xác định chính xác các đới đứt gãy sâu ẩn giấu dưới lớp phủ trầm tích? Nghiên cứu phối hợp phép nâng trường trọng lực lên độ cao 10 km và 20 km để lọc nhiễu nông, kết hợp thuật toán gradient ngang cực đại ($G_{max}$) và phương pháp góc nghiêng (Tilt Angle), sau đó đối chiếu trực tiếp với tọa độ các chấn tiêu động đất thực tế để xác định phương và góc dốc của đứt gãy.

Kết quả của luận văn đóng góp gì cho công tác thăm dò khoáng sản địa phương? Mô hình mật độ đã bóc tách rõ ranh giới các khối magma xâm nhập như khối granitoit Pia Oắc (diện tích 20 km²) và Phia Bioc. Điều này giúp các nhà địa chất xác định chính xác các đới tiếp xúc khoáng hóa thiếc, wolfram, vàng và bauxit, nâng cao hiệu quả kinh tế cho 142 điểm quặng trên địa bàn tỉnh.

Kết luận

  • Luận văn đã xây dựng thành công quy trình xử lý định lượng kết hợp giữa giải chập Euler, giải ngược mật độ 2D và giải ngược 3D trên tài liệu trọng lực Bouguer với sai số mô hình hóa được kiểm soát tối ưu dưới mức 5%.
  • Xác định và làm rõ đặc điểm hình học của các đới đứt gãy sâu chủ đạo, nổi bật là đới đứt gãy Cao Bằng – Tiên Yên (dài trên 250 km, rộng 3-15 km, mặt trượt dốc 60°-70°) và đứt gãy Yên Minh – Phú Lương.
  • Xây dựng hoàn chỉnh sơ đồ độ sâu bề mặt móng kết tinh (dao động từ 0,5 km đến trên 4,5 km) và bề mặt Moho (dao động từ 28,0 km đến 34,5 km) cho toàn bộ khu vực tỉnh Cao Bằng và lân cận.
  • Chuẩn hóa hệ thống thông số mật độ vật chất cho 3 tầng thạch quyển chính: vỏ trầm tích (2,62 g/cm³), vỏ kết tinh (2,83 g/cm³) và Manti trên (3,15 g/cm³), tạo cơ sở dữ liệu đồng bộ cho các nghiên cứu địa động lực tiếp theo.
  • Cung cấp luận cứ khoa học tin cậy hỗ trợ phân vùng hiểm họa động đất, giảm thiểu tai biến sụt nứt đất và định hướng thăm dò chiến lược cho 22 loại tài nguyên khoáng sản trong vùng.

Kế hoạch tiếp theo cần tập trung mở rộng mạng lưới đo đạc địa vật lý từ Tellua và địa chấn sâu trong giai đoạn 2024-2026 để kiểm chứng chi tiết các dị thường mật độ cục bộ. Các cơ quan quản lý và đơn vị chuyên môn hãy tích cực tham khảo và áp dụng các kết quả nghiên cứu này vào công tác quy hoạch phát triển kinh tế – xã hội và phòng chống thiên tai ngay hôm nay.