Tổng quan nghiên cứu

Trong lĩnh vực địa vật lý thăm dò, phương pháp từ đóng vai trò thiết yếu khi chiếm hơn 70% khối lượng khảo sát ban đầu phục vụ tìm kiếm khoáng sản và nghiên cứu cấu trúc sâu của vỏ Trái Đất. Bài toán xác định các yếu tố thế nằm, ranh giới và hình dạng thực tế của các khối đá có từ tính mạnh luôn là thách thức cốt lõi. Trong tự nhiên, các thể địa chất thường có hình thái bất đối xứng và cấu trúc phức tạp, khiến các công thức giải tích cổ điển chỉ áp dụng cho dạng hình học đơn giản như hình cầu, hình trụ hay vỉa phẳng bộc lộ nhiều hạn chế về độ chuẩn xác.

Nghiên cứu này tập trung giải quyết bài toán ngược hai chiều xác định tọa độ các đỉnh của vật thể có tiết diện ngang dạng đa giác bất kỳ dựa trên tài liệu dị thường từ toàn phần. Mục tiêu trọng tâm là xây dựng quy trình thuật toán tự động hóa hoàn toàn trên môi trường tính toán số, giúp khôi phục chính xác ranh giới thực tế và các tham số phông khu vực. Phạm vi nghiên cứu được thực hiện với các mô hình hai chiều mô phỏng tuyến khảo sát gồm 32 điểm đo, độ sâu nguồn từ biến thiên trong khoảng từ 1.50 km đến 10.00 km, hoàn thành trong giai đoạn 2011 đến 2013 tại Khoa Vật lý, Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội. Ý nghĩa của công trình thể hiện rõ qua việc nâng cao độ chính xác định vị ranh giới cấu trúc ngầm với sai số bình phương trung bình cực tiểu đạt mức 0.0002, đồng thời cắt giảm hơn 45% thời gian xử lý thủ công so với các phương pháp thử dần truyền thố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 được xây dựng trên nền tảng vững chắc của lý thuyết trường thế và giải tích vector địa vật lý. Hai trụ cột lý thuyết chính gồm có:

Thứ nhất, định luật cảm ứng từ và phương trình Poisson liên hệ giữa thế từ và thế trọng trường. Thông qua việc biểu diễn thế từ dưới dạng tích phân thể tích của momen từ phân bố đồng nhất, phương trình Poisson cho phép chuyển đổi đạo hàm không gian của thế trọng trường sang các thành phần cường độ từ trường. Từ đó, mối liên hệ vi phân Laplace được áp dụng triệt để nhằm mô tả sự biến thiên của các thành phần trường thẳng đứng và nằm ngang.

Thứ hai, mô hình giải tích dị thường từ hai chiều cho vật thể có tiết diện ngang đa giác N cạnh bất kỳ. Mô hình tích hợp góc từ khuynh hiệu dụng và góc từ thiên của trường từ Trái Đất cũng như của véc-tơ từ hóa. Các khái niệm then chốt bao gồm: dị thường từ toàn phần delta T, véc-tơ từ hóa tổng hợp J, góc nghiêng từ hóa trong mặt phẳng thẳng đứng với hai trường hợp khảo sát đặc trưng là 45 độ và 90 độ, cùng các hệ số phông tuyến tính khu vực A và B đặc trưng cho trường khu vực biến đổi chậm.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng tổng hợp trên tuyến đo địa vật lý gồm 32 điểm quan sát phân bố đều đặn với bước điểm 0.50 km, tổng chiều dài tuyến đạt 15.00 km. Dữ liệu bao trùm toàn bộ dải biến thiên dị thường từ -1475.00 nT đến hơn 8119.00 nT. Cỡ mẫu 32 điểm đo này được lựa chọn dựa trên nguyên tắc lấy mẫu đối xứng qua tâm dị thường, đảm bảo phản ánh đầy đủ cả hai cánh âm dương và vùng cực trị của đường cong trường từ.

Phương pháp phân tích trung tâm là thuật toán tối ưu hóa phi tuyến Marquardt. Thuật toán này kết hợp ưu điểm của phương pháp Gradient Descent ở các bước lặp đầu khi điểm khởi tạo còn xa nghiệm và phương pháp Gauss-Newton ở các bước lặp sau nhằm đạt tốc độ hội tụ bậc hai. Việc lựa chọn Marquardt giúp giải quyết triệt để sự bất ổn định của ma trận Hessian, tránh hiện tượng nghiệm bị phân kỳ hoặc rơi vào cực tiểu địa phương. Quá trình tính toán vi phân từng phần của hàm mục tiêu theo tọa độ các đỉnh và hệ số phông được tự động hóa bằng ngôn ngữ Matlab. Toàn bộ quy trình nghiên cứu, xây dựng thuật toán và thử nghiệm mô hình được thực hiện liên tục trong khung thời gian từ tháng 09/2011 đến tháng 12/2013.

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

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

Quá trình thử nghiệm tính toán trên các mô hình địa chất nhân tạo đã mang lại 4 phát hiện quan trọng:

Thứ nhất, đối với mô hình vật thể có tiết diện ngang dạng đẳng thước có góc nghiêng từ hóa 90 độ, thuật toán Marquardt đạt trạng thái hội tụ tuyệt đối sau 20 bước lặp. Độ lệch bình phương trung bình cuối cùng giảm mạnh từ hàng trăm đơn vị xuống còn 0.0002, tọa độ các đỉnh của đa giác tính toán trùng khớp hoàn toàn với mô hình thực tế với sai số vị trí dưới 0.0001 km.

Thứ hai, khi đưa thêm thành phần phông tuyến tính khu vực với các hệ số A = 200 và B = 1000 nT, chương trình vẫn tự động tách biệt và xác định chính xác 100% giá trị của A và B sau 22 bước lặp. Độ lệch bình phương trung bình ở lần lặp cuối đạt mức 0.0002, minh chứng cho khả năng tự khử phông mà không cần bước tiền xử lý lọc trường phức tạp.

Thứ ba, đối với mô hình vật thể có tiết diện ngang dạng kéo dài tương đương dạng tường hoặc vỉa quặng cắm nghiêng, tính phi tuyến của bài toán tăng lên đáng kể. Thuật toán cần 51 bước lặp để hội tụ với độ lệch bình phương trung bình 0.0007, thể hiện mức độ phức tạp hình học làm tăng số chu kỳ lặp thêm hơn 120% so với mô hình đẳng thước.

Thứ tư, khi góc nghiêng từ hóa thay đổi từ 90 độ sang 45 độ, hình dáng đường cong dị thường chuyển từ đối xứng sang dạng bất đối xứng mạnh với sự dịch chuyển rõ rệt của các cực đại và cực tiểu. Tuy nhiên, thuật toán vẫn khôi phục chính xác các đỉnh nằm sâu từ 1.50 km đến 8.00 km với độ lệch tọa độ cuối cùng không vượt quá 0.001 km.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp thuật toán đạt độ chính xác cao là nhờ cơ chế điều chỉnh động của tham số suy giảm trong phương pháp Marquardt. Khi sai số còn lớn, thuật toán tăng cường tính ổn định của ma trận nghịch đảo, sau đó chuyển nhanh sang hướng giải tích để tinh chỉnh vị trí các đỉnh đa giác.

So với các phương pháp nghịch đảo xấp xỉ tuyến tính cổ điển vốn thường mắc sai số tọa độ từ 5% đến 10% ở các góc nhọn của thân quặng, phương pháp này duy trì sai số tương đối dưới 0.05% trên mọi góc đỉnh. Kết quả nghiên cứu có thể được trình bày trực quan thông qua bảng đối sánh tọa độ đỉnh đa giác và biểu đồ đường cong dị thường. Trên đồ thị, đường cong dị thường tính toán lặp cuối bám khít hoàn toàn vào đường cong quan sát thực tế, trong khi mặt cắt hai chiều thể hiện sự trùng khít giữa đa giác tính toán và đa giác thực thể. Đồ thị bán logarit của độ hội tụ biểu diễn đường dốc suy giảm sai số thẳng đứng sau 5 bước lặp đầu tiên, phản ánh hiệu năng vượt trội của thuật toán.

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

Dựa trên kết quả nghiên cứu, 4 đề xuất cụ thể được đưa ra nhằm hoàn thiện và ứng dụng thuật toán trong thực tiễn:

Thứ nhất, tích hợp mô-đun lọc nhiễu ngẫu nhiên tiền xử lý. Các đơn vị địa vật lý cần ứng dụng bộ lọc biến đổi Wavelet kết hợp biến đổi Fourier để triệt tiêu từ 15% đến 20% biên độ nhiễu trắng trên tài liệu đo thực địa trước khi đưa vào giải bài toán ngược. Thời gian triển khai dự kiến trong 6 tháng tới.

Thứ hai, mở rộng thuật toán từ không gian hai chiều sang không gian hai chiều rưỡi và ba chiều. Các viện nghiên cứu chuyên ngành cần chuyển đổi mã nguồn Matlab sang ngôn ngữ C++ hoặc Python kết hợp tính toán song song trên GPU, nhằm mục tiêu giảm hơn 60% thời gian xử lý các mạng lưới đo dị thường quy mô lớn, hoàn thành trong giai đoạn 2024 đến 2025.

Thứ ba, ban hành quy chuẩn xác lập mô hình tiên nghiệm. Tổng cục Địa chất và Khoáng sản nên xây dựng cẩm nang kỹ thuật kết hợp dữ liệu địa chấn nông và trọng lực nhằm khoanh vùng tọa độ ban đầu, giúp giảm số bước lặp từ 50 xuống dưới 15 bước lặp cho mỗi cấu trúc.

Thứ tư, áp dụng thử nghiệm quy trình nghịch đảo vào thăm dò khoáng sản từ tính cao. Các liên đoàn địa chất nên triển khai thuật toán trên các thân quặng sắt và đa kim ở độ sâu từ 500 m đến 2000 m tại các mỏ khoáng sản trọng điểm, hướng tới mục tiêu nâng tỷ lệ chính xác của các lỗ khoan thăm dò lên trên 95% trước năm 2026.

Đố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 mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

Thứ nhất, các kỹ sư địa vật lý thăm dò. Nghiên cứu cung cấp công cụ tính toán chính xác để minh giải các tuyến đo từ kế mặt đất và hàng không, giúp xác định chuẩn xác bề dày, góc cắm và độ sâu đáy của các vỉa quặng ngầm.

Thứ hai, học viên cao học và nghiên cứu sinh chuyên ngành Vật lý Địa cầu. Tài liệu là nguồn tham khảo mẫu mực về quy trình toán học giải bài toán ngược trường thế, kỹ thuật lấy đạo hàm giải tích theo tọa độ đỉnh đa giác và phương pháp lập trình tối ưu trên Matlab.

Thứ ba, các nhà quản lý và chuyên gia quy hoạch mỏ địa chất. Việc ứng dụng thuật toán giúp xây dựng mô hình ranh giới quặng tin cậy ở độ sâu từ 1.00 km đến 3.00 km, tối ưu hóa mạng lưới lỗ khoan đắt đỏ và tiết kiệm đáng kể ngân sách thăm dò.

Thứ tư, các kỹ sư phát triển phần mềm địa tin học. Hệ thống công thức vi phân và thuật toán Marquardt trong luận văn là cơ sở thuật toán vững chắc để tích hợp vào các mô-đun phần mềm địa vật lý thương mại nội địa.

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

Bài toán ngược dị thường từ hai chiều có gặp phải hiện tượng đa nghiệm không? Có. Hiện tượng đa nghiệm là bản chất tự nhiên của bài toán ngược trường thế khi nhiều cấu trúc hình học khác nhau có thể tạo ra cùng một trường dị thường trên mặt đất. Nghiên cứu giải quyết vấn đề này bằng cách đưa vào thông tin tiên nghiệm từ tài liệu địa chất và ràng buộc chặt chẽ số đỉnh đa giác ban đầu.

Tại sao nghiên cứu lại ưu tiên thuật toán Marquardt thay vì phương pháp Gradient thông thường? Phương pháp Gradient truyền thống thường hội tụ rất chậm khi tiến gần đến đáy phễu sai số và dễ mắc kẹt tại các cực tiểu địa phương. Thuật toán Marquardt bổ sung tham số suy giảm linh hoạt, giúp duy trì tốc độ hội tụ nhanh và đạt sai số cực tiểu 0.0002 chỉ sau 20 đến 23 bước lặp.

Góc nghiêng từ hóa ảnh hưởng như thế nào đến kết quả định hình thân thể địa chất? Khi góc nghiêng từ hóa chuyển từ 90 độ sang 45 độ, trường dị thường mất tính đối xứng, biên độ cực trị âm dương dịch chuyển mạnh. Thuật toán đã giải quyết triệt để sự thay đổi này thông qua hệ công thức góc từ khuynh hiệu dụng, đảm bảo khôi phục đúng hình học vật thể.

Làm thế nào để thuật toán tự động tách bỏ phông khu vực tuyến tính? Nghiên cứu tích hợp trực tiếp hai hệ số phông tuyến tính A và B vào hệ véc-tơ tham số cần tối ưu hóa. Qua mỗi chu kỳ lặp, các tham số phông được cập nhật đồng thời với tọa độ đỉnh, giúp khôi phục chính xác 100% giá trị phông mà không làm méo mó hình dạng vật thể.

Quy trình tính toán trong luận văn có thể áp dụng trực tiếp cho số liệu thực địa không? Hoàn toàn có thể. Quy trình được thiết kế tổng quát cho tuyến đo 32 điểm hoặc nhiều hơn. Khi áp dụng vào thực địa, người dùng chỉ cần bổ sung bước lọc nhiễu tần số cao và nạp dữ liệu địa chất khu vực để thiết lập mô hình đa giác ban đầu.

Kết luận

Năm đóng góp cốt lõi của công trình nghiên cứu bao gồm:

  • Hoàn thiện hệ thống biểu thức giải tích tính dị thường từ toàn phần và vi phân từng phần cho vật thể hai chiều có tiết diện ngang đa giác N cạnh bất kỳ.
  • Ứng dụng thành công thuật toán tối ưu hóa phi tuyến Marquardt trên môi trường Matlab, đạt độ chính xác cao với sai số bình phương trung bình chỉ 0.0002.
  • Giải quyết triệt để bài toán đồng thời xác định ranh giới hình học vật thể và tự động khử phông tuyến tính khu vực.
  • Kiểm chứng thành công độ tin cậy và tốc độ hội tụ trên các mô hình địa chất phức tạp gồm vật thể đẳng thước, vật thể dạng tường kéo dài và mô hình móng từ.
  • Cung cấp giải pháp tính toán định lượng tin cậy phục vụ công tác giải đoán tài liệu thăm dò từ trường tại Việt Nam.

Trong kế hoạch 12 tháng tới, hướng phát triển trọng tâm là thử nghiệm thuật toán trên tập dữ liệu đo từ thực tế tại các bể trầm tích và mở rộng sang bài toán ba chiều. Các chuyên gia và nhà nghiên cứu địa vật lý được khuyến khích áp dụng khung giải thuật này để nâng cao hiệu quả thăm dò cấu trúc ngầm ngay hôm nay.