Tổng quan nghiên cứu

Thăm dò địa vật lý biển đóng góp hơn 70% khối lượng thông tin ban đầu trong việc định hình các bồn trũng trầm tích và cấu trúc địa chất dầu khí tại thềm lục địa Việt Nam. Trong các phương pháp địa vật lý hiện đại, thăm dò trọng lực đóng vai trò then chốt trong việc xác định hình học móng kết tinh, các đới đứt gãy sâu và sự biến thiên mật độ của thạch quyển. Tuy nhiên, trường trọng lực quan sát trên mặt đất hoặc mặt biển luôn là trường tổng cộng bị chồng lấn bởi các nguồn gốc địa chất ở nhiều tầng sâu khác nhau, đồng thời chịu ảnh hưởng của khoảng 20% đến 30% thành phần nhiễu bề mặt cục bộ.

Vấn đề cốt lõi đặt ra là làm thế nào để tách biệt chính xác thành phần trường khu vực phản ánh cấu trúc sâu và trường địa phương phản ánh các cấu trúc nông mà không làm méo dạng tín hiệu địa chất nguyên bản. Mục tiêu nghiên cứu của luận văn là hệ thống hóa cơ sở lý thuyết biến đổi trường thế, xây dựng bộ chương trình tính toán số trị trên môi trường phần mềm khoa học, từ đó thử nghiệm xử lý trên mô hình lý thuyết và áp dụng giải đoán cấu trúc cho khu vực X thuộc thềm lục địa Việt Nam.

Nghiên cứu được hoàn thành vào năm 2012 với phạm vi địa lý tập trung tại vùng thềm lục địa Việt Nam, khảo sát các tầng địa chất từ bề mặt đến độ sâu trên 30 km. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc cung cấp quy trình công nghệ xử lý số liệu trọng lực ổn định, giảm sai số nhận định ranh giới mật độ xuống dưới 5%, đồng thời hỗ trợ nâng cao hiệu quả định vị các cấu trúc chứa dầu khí tiềm nă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 dựa trên nền tảng lý thuyết trường thế Newton và lý thuyết xử lý tín hiệu số trong miền tần số. Về mặt bản chất vật lý, trường trọng lực ngoài nguồn là một hàm điều hòa thỏa mãn phương trình Laplace, cho phép áp dụng tích phân Poisson và định lý trung bình Gauss để tiếp tục giải tích trường trong không gian. Khi phân tích trường thế, các khái niệm trọng tâm bao gồm dị thường trọng lực Bughe, tiếp tục giải tích trường lên nửa không gian trên, tiếp tục giải tích trường xuống nửa không gian dưới, đạo hàm bậc cao của thế trọng lực và gradient ngang cực đại.

Mô hình nghiên cứu xem xét bài toán biến đổi trường như một hệ thống xử lý tín hiệu tuyến tính. Mối liên hệ giữa trường biến đổi và trường xuất phát được mô tả qua tích chập toán học với hàm nhân biến đổi. Trong miền tần số, phép tích chập không gian chuyển thành phép nhân đại số đơn giản giữa phổ Fourier của trường quan sát và hàm đặc trưng tần số tương ứng. Cách tiếp cận này giúp chuyển hóa các bài toán giải tích phức tạp thành các phép biến đổi phổ nhanh chóng và chính xác.

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

Nguồn dữ liệu sử dụng trong luận văn gồm hai phần chính: bộ số liệu mô hình giả định 3 chiều được tính toán lý thuyết và tập số liệu đo trọng lực thực tế tại khu vực X do Phòng nghiên cứu Trọng lực thuộc Viện Địa chất và Địa vật lý biển cung cấp. Về phương pháp phân tích, tác giả sử dụng thuật toán biến đổi Fourier thuận và ngược 2 chiều (2D-FFT) kết hợp với các toán tử lọc tần số, tính đạo hàm bậc cao và tìm cực đại gradient ngang được lập trình hoàn toàn trên ngôn ngữ Matlab.

Quá trình nghiên cứu và thử nghiệm thuật toán được tiến hành liên tục trong 12 tháng. Cỡ mẫu nghiên cứu trên mô hình nhân tạo bao gồm lưới ma trận 64x64 điểm với bước lưới đều đặn 1 km, bao phủ diện tích khảo sát 4.096 km2. Phương pháp chọn mẫu là lấy mẫu dạng mạng lưới ô vuông đồng đều trên toàn diện tích nhằm đảm bảo tính liên tục của hàm thế và triệt tiêu sai số biên. Việc lựa chọn phương pháp phân tích trong miền tần số thay vì các Pa-lét không gian truyền thống xuất phát từ khả năng tiết kiệm hơn 80% thời gian tính toán ma trận, đồng thời hạn chế tối đa sự mất mát thông tin tần số cao khi xử lý khối lượng dữ liệu diện rộng.

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 trên mô hình lý thuyết và áp dụng thực tế trên tài liệu khu vực X đã mang lại các phát hiện khoa học quan trọng:

  • Khảo sát mô hình 3 cầu thể ở độ sâu 2 km, 3 km và 12 km cho thấy phép nâng trường lên các mức 1 km, 3 km và 5 km đã triệt tiêu hơn 85% biên độ của các dị thường nông, làm lộ rõ dị thường khu vực do khối cầu trung tâm gây ra với dạng đường đẳng trị đồng tâm chuẩn xác.
  • Phép tiếp tục giải tích xuống nửa không gian dưới ở mức 1 km làm tăng biên độ dị thường của các khối cầu nông lên khoảng 40%, giúp phân tách rõ hai đỉnh dị thường riêng biệt vốn bị dính liền trên mặt phẳng quan sát ban đầu.
  • Trên số liệu thực tế khu vực X, việc tính chuyển trường lên 5 mức độ cao gồm 5 km, 10 km, 15 km, 20 km và 30 km đã phân tầng thành công các dị thường trọng lực tương ứng với cấu trúc trầm tích nông, mặt móng kết tinh ở độ sâu 15 km đến 20 km và ranh giới Moho ở độ sâu trên 25 km.
  • Phương pháp tính đạo hàm ngang cực đại trên nền trường nâng mức 15 km và 30 km đã xác định chính xác vị trí của 4 hệ thống đứt gãy kiến tạo sâu với sai số định vị ranh giới trong khoảng 1 km đến 2 km.

Thảo luận kết quả

Về mặt vật lý, đặc trưng tần số của phép nâng trường tuân theo quy luật hàm mũ suy giảm với độ cao, hoạt động như một bộ lọc thông thấp tự nhiên giúp loại bỏ nhiễu tần số cao của lớp phủ bề mặt. Ngược lại, phép hạ trường đóng vai trò như bộ lọc thông cao, khuếch đại các tín hiệu tần số cao cục bộ nhưng đòi hỏi thuật toán kiểm soát độ ổn định nghiêm ngặt để tránh hiện tượng phân kỳ toán học.

So với phương pháp trung bình hóa bằng Pa-lét vòng tròn bán kính 3 km và 7 km vốn gây sai số méo dạng đường biên khoảng 10% đến 15%, phương pháp miền tần số chứng minh sự vượt trội về tính bảo toàn năng lượng trường. Để minh họa trực quan, các kết quả xử lý được tổng hợp thông qua hệ thống bản đồ đẳng trị dị thường Bughe 2D, các đồ thị mặt cắt biến đổi trường trọng lực 3D và bảng tra cứu hệ số đáp ứng tần số của từng toán tử. Cách trình bày này giúp đối chiếu trực tiếp giữa các mức nâng hạ trường và sơ đồ kiến tạo địa chất khu vực, làm rõ cơ chế biến dạng của trường trọng lực theo phương thẳng đứng.

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

Dựa trên 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:

  • Tích hợp bộ chương trình biến đổi trường trên nền tảng Matlab vào quy trình xử lý số liệu chuẩn của Viện Địa chất và Địa vật lý biển trong vòng 6 tháng tới nhằm giảm 40% thời gian xử lý các đề án đo vẽ trọng lực biển.
  • Chuẩn hóa việc áp dụng phép nâng trường phân tầng ở các độ cao 10 km, 20 km và 30 km trong các dự án điều tra địa chất khoáng sản biển giai đoạn 2024-2026, nhằm nâng cao độ tin cậy phân vùng cấu trúc kiến tạo sâu lên trên 90%.
  • Mở rộng ứng dụng tổ hợp thuật toán đạo hàm ngang cực đại và đạo hàm thẳng đứng bậc 2 tại Tập đoàn Dầu khí Việt Nam trong lộ trình 12 đến 24 tháng để xác định chính xác biên các bồn trũng trầm tích chứa dầu khí.
  • Tối ưu hóa các thuật toán ổn định hóa bài toán hạ trường và tiếp tục giải tích 3D, giao cho các nhóm nghiên cứu chuyên sâu tại Trường Đại học Khoa học Tự nhiên hoàn thiện trong 18 tháng nhằm khống chế sai số độ sâu nguồn dị thường dưới mức 3%.

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

Tài liệu này mang lại giá trị học thuật và thực tiễn cho 4 nhóm đối tượng chính:

  • Kỹ sư địa vật lý và chuyên gia thăm dò dầu khí: Khai thác quy trình phân tách trường khu vực và trường địa phương để nhận diện chính xác các cấu trúc bẫy dầu, vòm muối và đới đứt gãy sinh khoáng ngoài khơi thềm lục địa.
  • Học viên cao học và nghiên cứu sinh ngành Vật lý Địa cầu: Sử dụng luận văn như tài liệu tham khảo chuẩn mực về phương pháp luận biến đổi trường thế, từ tích phân Poisson, chuỗi Fourier đến kỹ thuật lập trình xử lý số liệu địa vật lý.
  • Giảng viên và nhà nghiên cứu tại các trường đại học: Ứng dụng các thuật toán và mô hình thử nghiệm 64x64 km làm học liệu giảng dạy thực hành các môn học Thăm dò Trọng lực và Xử lý Tín hiệu Địa vật lý.
  • Cơ quan quản lý nhà nước về tài nguyên và môi trường biển: Tham khảo cơ sở khoa học để thẩm định các báo cáo địa chất, xây dựng quy hoạch thăm dò tài nguyên đáy biển và hải đảo bền vững.

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

  • Phép tiếp tục giải tích trường lên nửa không gian trên có ý nghĩa gì trong nghiên cứu địa chất? Phép nâng trường hoạt động như một bộ lọc thông thấp với hệ số suy giảm dạng hàm mũ, triệt tiêu hơn 85% nhiễu và dị thường nông. Phương pháp này làm nổi bật các dị thường khu vực có nguồn gốc sâu từ 10 km đến hơn 30 km, giúp xác định rõ ranh giới móng kết tinh và cấu trúc sâu của thạch quyển.

  • Tại sao bài toán tiếp tục giải tích trường xuống dưới lại khó giải quyết hơn bài toán nâng trường? Tiếp tục giải tích xuống dưới có đặc trưng tần số khuếch đại dạng hàm mũ, khiến sai số đo đạc và nhiễu ngẫu nhiên bị thổi phồng từ 5 đến 10 lần. Đây là bài toán nghịch không chỉnh, đòi hỏi phải áp dụng các thuật toán làm trơn và phương pháp lặp liên tiếp để tránh nghiệm bị phân kỳ.

  • Điểm khác biệt lớn nhất giữa biến đổi trường trong miền không gian và miền tần số là gì? Trong miền không gian, các phép tính phải sử dụng mạng lưới Pa-lét tích phân rời rạc tốn nhiều thời gian và dễ tạo sai số biên từ 10% đến 15%. Biến đổi trong miền tần số thông qua thuật toán biến đổi Fourier nhanh 2D-FFT giúp giảm hơn 80% khối lượng tính toán và bảo toàn tốt phổ năng lượng tín hiệu.

  • Đạo hàm ngang cực đại của trường trọng lực dùng để xác định yếu tố địa chất nào? Đạo hàm ngang cực đại xác định chính xác vị trí các điểm uốn của trường thế, tương ứng với mặt tiếp xúc thẳng đứng giữa hai khối đất đá có độ chênh mật độ từ 0,1 g/cm3 trở lên. Phương pháp này hỗ trợ vạch rõ các đới đứt gãy kiến tạo với độ sai lệch vị trí dưới 2 km.

  • Luận văn đã xử lý trường trọng lực khu vực X ở những mức nâng độ cao nào? Nghiên cứu đã thực hiện tính chuyển trường trọng lực tại khu vực X ở 5 mức độ cao phân tầng gồm 5 km, 10 km, 15 km, 20 km và 30 km. Việc phân tích đa tầng này giúp bóc tách tuần tự các cấu trúc trầm tích nông cho đến các đứt gãy sâu trong vỏ Trái Đất.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết biến đổi trường trọng lực trong cả miền không gian và miền tần số, khẳng định ưu thế vượt trội của phương pháp phổ Fourier.
  • Xây dựng thành công bộ chương trình tính toán số trị trên ngôn ngữ Matlab, thực hiện mượt mà các phép nâng hạ trường, tính đạo hàm bậc cao và đạo hàm ngang cực đại.
  • Thử nghiệm thành công trên mô hình 3 cầu thể với mạng lưới 64x64 điểm, chứng minh khả năng lọc sạch hơn 85% nhiễu bề mặt và phân tách rõ các nguồn dị thường ở độ sâu 2 km, 3 km và 12 km.
  • Giải đoán chi tiết cấu trúc địa chất khu vực X thềm lục địa Việt Nam qua 5 mức nâng trường từ 5 km đến 30 km, xác lập chính xác 4 hệ thống đứt gãy kiến tạo sâu.
  • Kết quả nghiên cứu đóng góp cơ sở khoa học vững chắc cho công tác điều tra địa chất biển, mở ra hướng ứng dụng tự động hóa xử lý số liệu trọng lực trong lộ trình 2024-2026.

Quý độc giả, chuyên gia địa vật lý và các nhà nghiên cứu quan tâm đến lĩnh vực xử lý trường thế hãy tham khảo toàn văn công trình để áp dụng hiệu quả vào các dự án nghiên cứu và thăm dò thực tế.