Luận văn: Áp dụng thuật toán Parker xác định dị thường trọng lực Bouguer Biển Đông

Luận văn nghiên cứu ứng dụng phương pháp CALUX trong phân tích, đánh giá mức độ ô nhiễm dioxin tại một số khu vực ô nhiễm nặng ở Việt Nam.

• Tác giả: Phạm Thành Luân
• Người hướng dẫn: PGS. Đỗ Đức Thanh
• Trường học: Đại học Quốc gia Hà Nội, Trường Đại học Khoa học Tự nhiên
• Chuyên ngành: Vật lý địa cầu
• Đề tài: Nghiên Cứu Áp Dụng Thuật Toán Parker Xác Định Dị Thường Trọng Lực Bouguer Khu Vực Biển Đông Và Kế Cận
• Loại tài liệu: Luận văn Thạc sĩ Khoa học
• Năm xuất bản: 2015
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỞ ĐẦU
• 1. CHƯƠNG 1: PHÉP BIẾN ĐỔI FOURIER
• 1.1. Khái niệm biến đổi Fourier
• 1.2. Một số tính chất của biến đổi Fourier
• 1.2.1. Tính chất đối xứng
• 1.2.2. Tính chất tuyến tính
• 1.2.3. Sự chuyển dịch
• 1.3. Biến đổi Fourier của một số dị thường đơn giản
• 1.3.1. Các nguồn ba chiều
• 1.3.2. Các nguồn hai chiều
• 2. CHƯƠNG 2: THUẬT TOÁN PARKER VÀ MỘT SỐ PHƯƠNG PHÁP KHÁC XÁC ĐỊNH DỊ THƯỜNG TRỌNG LỰC
• 2.1. Thuật toán Parker
• 2.2. Một số phương pháp khác xác định dị thường trọng lực
• 3. CHƯƠNG 3: MÔ HÌNH HÓA VÀ KẾT QUẢ ÁP DỤNG THUẬT TOÁN PARKER XÁC ĐỊNH DỊ THƯỜNG BOUGUER KHU VỰC BIỂN ĐÔNG VÀ KẾ CẬN
• 3.1. Mô hình ranh giới phân chia dạng vòm 2D
• 3.1.1. Thông số mô hình
• 3.1.2. Kết quả tính toán
• 3.2. Mô hình ranh giới phân chia dạng vòm 3D
• 3.2.1. Thông số mô hình
• 3.2.2. Kết quả tính toán
• 3.3. Mô hình bể trầm tích 3D
• 3.3.1. Thông số mô hình
• 3.3.2. Kết quả tính toán
• 3.4. Xác định dị thường Bouguer trên biển Đông và kế cận
• 3.4.1. Nguồn số liệu
• 3.4.2. Kết quả tính toán

Tóm tắt

I. Giới thiệu về Thuật toán Parker và Ứng dụng trong Địa Vật Lý

Thuật toán Parker là một phương pháp tiên tiến trong xử lý dữ liệu dị thường trọng lực Bouguer, đặc biệt quan trọng trong nghiên cứu cấu trúc địa chất khu vực Biển Đông. Phương pháp này được phát triển dựa trên biến đổi Fourier và cho phép xác định chính xác các dị thường trọng lực gây ra bởi sự thay đổi mật độ các lớp địa chất. Trong bối cảnh khảo sát địa vật lý biển Đông, thuật toán Parker đóng vai trò quan trọng trong việc phân tích cấu trúc vỏ Trái Đất, xác định các ranh giới địa chất, và tìm kiếm tài nguyên dầu khí. Các ứng dụng của phương pháp này không chỉ giới hạn ở việc xác định các dị thường mà còn hỗ trợ trong mô hình hóa ba chiều các cấu trúc địa chất phức tạp.

1.1. Khái niệm Cơ Bản về Biến Đổi Fourier

Biến đổi Fourier là nền tảng toán học cho thuật toán Parker, cho phép chuyển đổi dữ liệu từ miền không gian sang miền tần số. Phương pháp này giúp phân tách các thành phần dị thường khác nhau dựa trên tần số của chúng. Tính chất đối xứng và tuyến tính của biến đổi Fourier cho phép xử lý các dữ liệu trọng lực phức tạp, từ đó xác định mối quan hệ giữa địa hình và dị thường trọng lực.

1.2. Tầm Quan Trọng trong Nghiên Cứu Biển Đông

Biển Đông là khu vực có cấu trúc địa chất phức tạp với nhiều ranh giới phân chia dạng vòm. Ứng dụng thuật toán Parker giúp nhà khoa học xác định độ sâu các ranh giới, phân tích cấu trúc bể trầm tích và đánh giá tiềm năng khám phá tài nguyên. Việc áp dụng phương pháp này trên dữ liệu biển Đông cung cấp thông tin quan trọng về hình học địa chất ba chiều.

II. Nguyên Lý và Phương Pháp Toán Học của Thuật Toán Parker

Thuật toán Parker hoạt động dựa trên mối liên hệ giữa dị thường trọng lực và mô hình mật độ của các lớp địa chất. Phương pháp này sử dụng biến đổi Fourier hai chiều để chuyển đổi dữ liệu từ miền không gian sang miền tần số, từ đó tính toán dị thường trọng lực với độ chính xác cao. So với các phương pháp truyền thống, thuật toán Parker cho phép xử lý dữ liệu lớn một cách hiệu quả hơn. Phương pháp này đặc biệt hữu ích khi làm việc với ranh giới phân chia phức tạp và các mô hình ba chiều của bể trầm tích. Độ chính xác của kết quả phụ thuộc vào chất lượng dữ liệu đầu vào và các tham số mô hình được sử dụng.

2.1. Quy Trình Tính Toán Chi Tiết

Quy trình áp dụng thuật toán Parker bao gồm các bước: (1) thu thập dữ liệu dị thường Free-air và thông tin địa hình; (2) chuẩn hóa dữ liệu và xác định tham số mô hình; (3) thực hiện biến đổi Fourier hai chiều; (4) tính toán dị thường Bouguer trong miền tần số; (5) biến đổi ngược sang miền không gian. Mỗi bước đòi hỏi sự kiểm soát chất lượng chặt chẽ để đảm bảo độ tin cậy của kết quả.

2.2. So Sánh với Các Phương Pháp Khác

Ngoài thuật toán Parker, các phương pháp như Bhaskara Rao trong miền không gian và Chai-Hinze trong miền tần số cũng được sử dụng. Thuật toán Parker nổi bật với tốc độ xử lý nhanh và độ chính xác cao, đặc biệt khi xử lý dữ liệu không đều hoặc lớn. Tuy nhiên, mỗi phương pháp có những ưu điểm riêng tùy theo bối cảnh ứng dụng cụ thể.

III. Ứng Dụng Mô Hình Hóa trong Khu Vực Biển Đông

Mô hình hóa dị thường Bouguer ở Biển Đông sử dụng thuật toán Parker giúp xác định chính xác các cấu trúc địa chất quan trọng. Các mô hình nghiên cứu bao gồm: ranh giới phân chia dạng vòm hai chiều, ranh giới ba chiều, và bể trầm tích ba chiều. Thông qua mô phỏng số, các nhà khoa học có thể so sánh kết quả tính toán với dữ liệu thực tế, từ đó hiệu chỉnh mô hình và nâng cao độ chính xác. Dữ liệu được sử dụng bao gồm bản đồ địa hình đáy biển, dị thường Free-air, và các hiệu chỉnh trọng lực. Kết quả mô hình cung cấp thông tin quý báu về cấu trúc địa chất sâu và tiềm năng khám phá tài nguyên.

3.1. Mô Hình Ranh Giới Phân Chia Dạng Vòm

Ranh giới phân chia dạng vòm là một đặc trưng phổ biến ở Biển Đông, được tạo thành bởi sự chuyển động của các địa khối. Mô hình hai chiều giúp xác định độ sâu và hình dạng ranh giới, trong khi mô hình ba chiều cung cấp hình ảnh chi tiết hơn về cấu trúc không gian. Thuật toán Parker cho phép tính toán dị thường trọng lực từ các mô hình này và so sánh với dữ liệu đo đạc thực tế.

3.2. Phân Tích Bể Trầm Tích Ba Chiều

Bể trầm tích Biển Đông chứa các lớp trầm tích có mật độ khác nhau. Mô hình ba chiều giúp xác định độ sâu đáy bể, cấu trúc nội bộ và phân bố mật độ. Kết quả phân tích cho thấy chênh lệch giữa dị thường tính toán và dữ liệu quan sát, từ đó điều chỉnh mô hình và cải thiện độ chính xác xác định các đặc trưng địa chất quan trọng.

IV. Kết Quả và Ý Nghĩa Khoa Học của Nghiên Cứu

Áp dụng thuật toán Parker xác định dị thường Bouguer khu vực Biển Đông và kế cận đã mang lại những kết quả quan trọng. Bản đồ dị thường Bouguer được tạo ra thể hiện rõ các cấu trúc địa chất lớn, từ các ranh giới phân chia tới các bể trầm tích sâu. Việc hiệu chỉnh trọng lực gây ra bởi địa hình và độ sâu đáy biển cho phép loại bỏ các tác động địa hình, từ đó làm nổi bật các dị thường gây bởi thay đổi mật độ sâu. Kết quả nghiên cứu cung cấp cơ sở dữ liệu vật lý quan trọng cho việc xây dựng mô hình cấu trúc địa chất, tìm kiếm tài nguyên dầu khí, và hiểu rõ hơn về lịch sử địa chất khu vực Biển Đông.

4.1. Các Phát Hiện Chính và Độ Chính Xác

Kết quả mô hình hóa cho thấy độ chính xác cao của thuật toán Parker trong xác định dị thường trọng lực. Chênh lệch giữa kết quả tính toán và dữ liệu thực tế được kiểm soát ở mức chấp nhận được. Bản đồ đồng mức dị thường Bouguer thể hiện các pattern rõ ràng liên quan đến cấu trúc địa chất, cho phép xác định vị trí và độ sâu các ranh giới quan trọng với độ tin cậy cao.

4.2. Ứng Dụng Thực Tiễn và Triển Vọng

Nghiên cứu này cung cấp công cụ hiệu quả cho các nhà khoa học trong việc phân tích dữ liệu địa vật lý biển. Kết quả có thể ứng dụng trong khám phá tài nguyên, xây dựng mô hình địa chất, và nghiên cứu cấu trúc vỏ Trái Đất. Thuật toán Parker có tiềm năng ứng dụng rộng rãi trong các khu vực biển khác, hỗ trợ phát triển bền vững và quản lý tài nguyên biển hiệu quả.

Nội dung chính

Tổng quan nghiên cứu

Thăm dò trọng lực đóng vai trò nền tảng trong việc làm sáng tỏ cấu trúc sâu vỏ Trái Đất, phân vùng kiến tạo và xác định các bồn trũng trầm tích có tiềm năng dầu khí. Trong các khảo sát địa vật lý biển quy mô lớn, việc xử lý khối lượng dữ liệu khổng lồ với hàng triệu điểm quan sát theo phương pháp tích phân trong miền không gian thường đòi hỏi chi phí tính toán rất lớn. Cụ thể, khi xử lý một mạng lưới chứa 1.000.000 điểm đo, thuật toán trong miền không gian tiêu tốn thời gian gấp khoảng 72.000 lần so với miền tần số, kéo dài thời gian xử lý thực tế từ 1 phút lên tới hơn 50 ngày làm việc liên tục.

Vấn đề cốt lõi mà nghiên cứu giải quyết là tối ưu hóa tốc độ và độ chính xác khi xác định dị thường trọng lực và giải bài toán ngược xác định ranh giới mật độ. Mục tiêu cụ thể của luận văn là nghiên cứu, làm chủ và ứng dụng thuật toán Parker trên nền tảng biến đổi Fourier nhanh (FFT) bằng ngôn ngữ lập trình chuyên dụng. Nghiên cứu tập trung đánh giá hiệu năng giải thuật trên các mô hình 2D và 3D lý thuyết, sau đó áp dụng tính toán dị thường trọng lực Bouguer hoàn chỉnh cho toàn bộ khu vực Biển Đông và kế cận.

Phạm vi không gian nghiên cứu bao trùm toàn bộ bồn trũng Biển Đông cùng các cấu trúc chuyển tiếp lục địa - đại dương lân cận, hoàn thành tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội vào năm 2015. Kết quả đạt được ghi nhận tốc độ xử lý nhanh hơn từ 17 lần đối với mô hình 2 chiều đến hơn 375 lần đối với mô hình 3 chiều phức tạp, trong khi sai số bình phương trung bình luôn duy trì ở mức dưới 0,01 mGal, mang lại giải pháp công nghệ hiệu quả cho công tác minh giải địa vật lý biển Việt Nam.

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 xây dựng trên nền tảng lý thuyết thế trường trọng lực và giải tích Fourier nâng cao. Trọng tâm của khung lý thuyết bao gồm:

• Thuật toán Parker (1972): Khai triển thế trọng lực của lớp vật liệu có ranh giới bất kỳ và mật độ biến thiên thành chuỗi biến đổi Fourier vô hạn. Bằng cách chuyển đổi phép tính tích phân phức tạp sang chuỗi lũy thừa của số sóng trong miền tần số, dị thường trọng lực thành phần thẳng đứng được xác định thông qua phép biến đổi Fourier ngược của tổng chuỗi Taylor.
• Phương pháp nghịch đảo Parker - Oldenburg (1974): Thuật toán giải bài toán ngược xác định hình thái ranh giới phân chia mật độ dựa trên quá trình lặp trong miền tần số từ trường dị thường quan sát.
• Lý thuyết lọc số và tần số cắt Nyquist: Ứng dụng nguyên lý triệt tiêu hiện tượng chồng phổ (aliasing) và thiết lập bộ lọc thông thấp (Low-Pass Filter) dạng cosin để ổn định quá trình lặp, ngăn chặn sự phân kỳ do số hạng khuếch đại tần số cao gây ra.

Các khái niệm then chốt trong nghiên cứu gồm: dị thường Free-air, hiệu chỉnh địa hình đáy biển 3D, dị thường Bouguer hoàn chỉnh, mật độ dư tương phản và mức độ sâu tham chiếu.

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

Phương pháp nghiên cứu kết hợp chặt chẽ giữa mô hình hóa toán học, thử nghiệm mô phỏng số và xử lý số liệu thực tế:

• Nguồn dữ liệu: Nghiên cứu sử dụng 3 bộ dữ liệu mô phỏng chuẩn hóa (mô hình vòm 2D với 256 điểm, mô hình vòm 3D với 4.096 điểm, mô hình bể trầm tích 3D với 16.384 điểm) và bộ dữ liệu độ cao địa hình, độ sâu đáy biển cùng dị thường Free-air thực tế khu vực Biển Đông từ các nguồn đo đạc vệ tinh và hải dương học quốc tế.
• Cỡ mẫu và phương pháp chọn mẫu: Lưới điểm tính toán được thiết lập đồng nhất theo quy chuẩn không gian. Cỡ mẫu mô hình 2D gồm 256 điểm với khoảng cách trạm 0,5 km; mô hình 3D dạng vòm gồm mạng lưới $64 \times 64$ điểm đo (khoảng cách 1 km); mô hình bể trầm tích gồm mạng lưới $128 \times 128$ điểm đo (khoảng cách 1 km). Phương pháp lấy mẫu cách đều đảm bảo đáp ứng đầy đủ điều kiện lấy mẫu Nyquist-Shannon.
• Lý do lựa chọn phương pháp: Thuật toán Parker chuyển toàn bộ phép tính tích phân chập phức tạp $O(N^2)$ trong miền không gian sang phép nhân ma trận $O(N \log N)$ thông qua thuật toán FFT, giúp rút ngắn thời gian xử lý xuống hàng trăm lần mà vẫn bảo toàn độ chính xác của các đường đồng mức dị thường. Toàn bộ quy trình được lập trình và thực thi trên môi trường phần mềm tính toán khoa học trong thời gian 12 thá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 mô hình và áp dụng số liệu thực tế đã mang lại các kết quả định lượng cụ thể:

• Mô hình ranh giới vòm 2D ($N = 256$ điểm): Thuật toán Parker hoàn thành tính toán bài toán thuận chỉ trong 0,0138 giây, nhanh hơn 17,3 lần so với phương pháp tích phân đa giác của Murthy và Bhaskara Rao trong miền không gian (0,2399 giây). Sai số bình phương trung bình (RMS) giữa hai phương pháp đạt mức cực nhỏ, chỉ xấp xỉ 0,0001 mGal. Bài toán ngược hội tụ hoàn toàn sau 4 vòng lặp với sai số độ sâu RMS đạt 0,002 km.
• Mô hình ranh giới vòm 3D ($64 \times 64 = 4.096$ điểm): Thuật toán Parker xác định dị thường trọng lực trong 0,2724 giây, vượt trội hơn 83,8 lần so với phương pháp lăng trụ chữ nhật của Bhaskara Rao trong miền không gian (22,8400 giây). Độ lệch dị thường giữa hai phương pháp có RMS chỉ 0,0100 mGal. Bài toán ngược kết hợp bộ lọc thông thấp với thông số $WH = 0,3$ hội tụ sau 5 vòng lặp, đạt sai số độ sâu RMS là 0,0100 km so với mô hình gốc.
• Mô hình bể trầm tích 3D ($128 \times 128 = 16.384$ điểm): Khi số điểm đo tăng lên, thuật toán Parker thể hiện ưu thế vượt trội khi chỉ mất 0,9800 giây để hoàn thành, trong khi phương pháp trong miền không gian tiêu tốn 367,5159 giây (tốc độ xử lý tăng hơn 375 lần, giảm 99,7% thời gian chạy máy). So sánh với phương pháp miền tần số của Chai và Hinze (có bước trượt mẫu 0,26), thuật toán Parker vẫn cho thời gian tính toán nhanh hơn đáng kể với sai số RMS tương đồng dưới 0,0200 mGal. Bài toán ngược đạt độ hội tụ chính xác sau 7 vòng lặp tại độ sâu tham chiếu $z_0 = 2,84\text{ km}$.
• Xác định dị thường Bouguer Biển Đông: Xây dựng thành công sơ đồ hiệu chỉnh địa hình đáy biển 3D và thành lập bản đồ đồng mức dị thường trọng lực Bouguer dạng 2D và 3D cho toàn bộ Biển Đông, phản ánh rõ nét ranh giới móng kết tinh, các trục đứt gãy sâu và sự biến thiên chiều dày vỏ Trái Đất.

Thảo luận kết quả

Kết quả tính toán dị thường trọng lực và độ sâu ranh giới có thể được biểu diễn trực quan thông qua bảng so sánh hiệu năng và hệ thống đồ thị đối sánh:

Mô hình thử nghiệm	Số điểm tính toán	Thời gian Parker (giây)	Thời gian miền không gian (giây)	Tốc độ cải thiện (lần)	Sai số RMS dị thường (mGal)	Số vòng lặp giải ngược
Vòm 2D	256	0,0138	0,2399	17,3	< 0,0001	4
Vòm 3D	4.096 ($64 \times 64$)	0,2724	22,8400	83,8	0,0100	5
Bể trầm tích 3D	16.384 ($128 \times 128$)	0,9800	367,5159	375,0	0,0200	7

Trên đồ thị lát cắt 2D và bản đồ đồng mức 3D, đường cong dị thường tính theo thuật toán Parker hoàn toàn trùng khít với đường cong tính theo phương pháp giải tích không gian truyền thống. Sự chênh lệch biên độ dị thường tập trung chủ yếu ở các rìa biên mô hình do ảnh hưởng biên của biến đổi Fourier, tuy nhiên biên độ sai lệch cực đại không vượt quá 0,05 mGal.

Nguyên nhân cốt lõi giúp thuật toán Parker đạt hiệu suất vượt bậc là việc khai thác tối đa cấu trúc nhân chập trong miền tần số, cho phép tính toán đồng thời tác động của toàn bộ ma trận số liệu thay vì phải cộng dồn từng cặp điểm quan sát - yếu tố thể tích. Việc lựa chọn mức tham chiếu $z_0$ phù hợp kết hợp với bộ lọc thông thấp (với tần số cắt thích hợp từ 0,04 đến 0,3) đóng vai trò then chốt trong việc kiểm soát các sóng ngắn nhiễu bậc cao, đảm bảo quá trình giải bài toán ngược ổn định và không bị phân kỳ.

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

• Tích hợp thuật toán vào hệ thống xử lý dữ liệu địa vật lý quốc gia: Bộ Tài nguyên và Môi trường cùng các cơ quan điều tra địa chất biển cần đưa mã nguồn thuật toán Parker vào phần mềm xử lý chuyên dụng, hướng tới mục tiêu rút ngắn 95% thời gian xử lý các bộ dữ liệu trọng lực biển quy mô trên 1.000.000 trạm đo trong giai đoạn 2026-2028.
• Xây dựng quy chuẩn xác lập thông số lọc thông thấp tự động: Các nhóm nghiên cứu chuyên ngành tại các trường đại học cần nghiên cứu thuật toán tự động tối ưu hóa hai tham số $WH$ và $SH$ của bộ lọc cho từng bồn trũng cấu trúc cụ thể, nhằm duy trì sai số bài toán ngược dưới mức 1,5% độ sâu thực tế, hoàn thành thử nghiệm trong vòng 18 tháng.
• Ứng dụng minh giải chiều sâu móng các bể trầm tích dầu khí: Tập đoàn Dầu khí Việt Nam nên áp dụng mô hình giải ngược Parker - Oldenburg 3D để tái tạo cấu trúc móng và bề mặt Moho tại các khu vực nước sâu thuộc bồn trũng Nam Côn Sơn, Cửu Long và Tư Chính - Vũng Mây, hoàn thành báo cáo đánh giá cấu trúc trước quý IV/2027.
• Mở rộng giải thuật cho dị thường từ và mô hình mật độ phi tuyến: Các viện nghiên cứu địa vật lý cần tiếp tục phát triển mở rộng thuật toán Parker cho tài liệu trường từ hàng không - biển và các mô hình mật độ thay đổi phi tuyến theo chiều sâu, triển khai thử nghiệm trên diện tích tối thiểu 50.000 $\text{km}^2$ khảo sát thực địa trong 2 năm tới.

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

• Học viên cao học và nghiên cứu sinh ngành Vật lý địa cầu, Địa chất: Tài liệu cung cấp cơ sở toán học chặt chẽ về lý thuyết thế trường, kỹ thuật biến đổi Fourier 2D/3D và thuật toán giải bài toán ngược với các ví dụ code tính toán cụ thể.
• Kỹ sư xử lý dữ liệu thăm dò tại các doanh nghiệp dầu khí và khoáng sản biển: Hỗ trợ quy trình tự động hóa phân tích dị thường trọng lực trên các mảng dữ liệu lớn với hơn 10.000 điểm đo, nâng cao năng suất xử lý bản đồ cấu trúc đáy biển.
• Chuyên viên địa chính và trắc địa biển tại các cơ quan quản lý nhà nước: Ứng dụng phương pháp hiệu chỉnh địa hình 3D và tính toán dị thường Bouguer phục vụ xây dựng mô hình Geoid cục bộ với độ phân giải cao tại thềm lục địa Việt Nam.
• Giảng viên và nhà nghiên cứu thuật toán xử lý tín hiệu số trong khoa học Trái Đất: Tham khảo cấu trúc thuật toán FFT, kỹ thuật lọc chống nhiễu số sóng cao và mô hình lặp phi tuyến để phát triển các bài giảng và công trình nghiên cứu chuyên sâu.

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

Thuật toán Parker vượt trội hơn các phương pháp trong miền không gian ở điểm nào?
Thuật toán Parker chuyển đổi phép tích phân thể tích phức tạp thành chuỗi biến đổi Fourier nhanh (FFT), giúp giảm độ phức tạp tính toán từ $O(N^2)$ xuống $O(N \log N)$. Khi xử lý 16.384 điểm đo, thuật toán chỉ mất 0,98 giây thay vì 367,52 giây, tăng tốc độ xử lý hơn 375 lần mà vẫn đảm bảo sai số RMS dưới 0,02 mGal.

Tại sao cần sử dụng bộ lọc thông thấp trong bài toán ngược Parker - Oldenburg?
Trong công thức giải ngược, thành phần lũy thừa $e^{|k|z_0}$ sẽ khuếch đại mạnh mẽ các tần số cao (bước sóng ngắn). Nếu không sử dụng bộ lọc thông thấp, nhiễu bề mặt và sai số tính toán biên sẽ làm quá trình lặp bị phân kỳ. Bộ lọc dạng cosin giúp triệt tiêu tần số lớn hơn ngưỡng $SH$, đảm bảo quá trình lặp hội tụ ổn định sau 4-7 chu kỳ.

Độ sâu tham chiếu $z_0$ có vai trò gì trong quá trình tính toán?
Độ sâu tham chiếu $z_0$ đóng vai trò là mặt phẳng chuẩn của lớp Bouguer tương đương, giúp cải thiện sự hội tụ của chuỗi Taylor trong thuật toán Parker. Việc lựa chọn $z_0$ bằng độ sâu trung bình của cấu trúc (ví dụ $z_0 = 5\text{ km}$ cho ranh giới vòm hoặc $z_0 = 2,84\text{ km}$ cho bể trầm tích) giúp triệt tiêu hiện tượng chồng phổ hiệu quả.

Thuật toán Parker có thể áp dụng cho tài liệu dị thường từ không?
Thuật toán Parker hoàn toàn có thể mở rộng cho trường từ thông qua việc nhân thêm các toán tử định hướng vectơ từ hóa và vectơ trường khu vực. Nguyên lý khai triển chuỗi Fourier vẫn được giữ nguyên để tính toán dị thường từ gây bởi các ranh giới từ tính phức tạp trong vỏ Trái Đất.

Ý nghĩa của việc thành lập bản đồ dị thường Bouguer khu vực Biển Đông là gì?
Bản đồ dị thường Bouguer loại bỏ hoàn toàn ảnh hưởng trọng trường của cột nước biển và địa hình đáy biển dạng 3D, phản ánh trung thực sự bất đồng nhất mật độ sâu trong vỏ và manti. Đây là tài liệu địa vật lý nền tảng để nghiên cứu ranh giới manti - vỏ (Moho), cơ chế tách giãn đáy biển và đánh giá tiềm năng bồn trũng dầu khí.

Kết luận

• Thuật toán Parker chứng minh tính ưu việt vượt trội về thời gian tính toán so với các phương pháp trong miền không gian, ghi nhận tốc độ xử lý nhanh hơn từ 17,3 lần (mô hình 2D) đến 375 lần (mô hình bể trầm tích 3D 16.384 điểm).
• Độ chính xác của thuật toán được kiểm chứng tin cậy với sai số bình phương trung bình cực nhỏ (RMS $\le 0,02\text{ mGal}$ cho bài toán thuận và RMS $\le 0,01\text{ km}$ cho bài toán ngược).
• Ứng dụng thành công quy trình tính toán trên môi trường số để xây dựng bản đồ dị thường Bouguer 2D và 3D khu vực Biển Đông và kế cận, phản ánh chính xác cấu trúc sâu vỏ Trái Đất.
• Việc kết hợp bộ lọc thông thấp dạng cosin và lựa chọn mặt phẳng tham chiếu $z_0$ hợp lý đã giải quyết triệt để bài toán phân kỳ trong nghịch đảo địa vật lý.
• Đề xuất kế hoạch hoàn thiện bộ phần mềm xử lý dữ liệu trọng lực biển tự động hóa cao trong giai đoạn 2026-2028 nhằm phục vụ công tác điều tra tài nguyên và nghiên cứu hải dương học.

Quý độc giả, nhà nghiên cứu và kỹ sư địa vật lý quan tâm có thể khai thác trực tiếp phương pháp luận và các kết quả mô hình số trong toàn văn luận văn để áp dụng vào các dự án nghiên cứu cấu trúc địa chất và thăm dò khoáng sản biển.

Trích đoạn nội dung tài liệu

Chương 1: Phép biến đổi Fourier - Chương 2: Thuật toán Parker và một số phương pháp khác xác định dị thường trọng lực. - Chương 3: Mô hình hóa và kết quả áp dụng thuật toán Parker xác định dị thường Bouguer khu vực biển Đông và kế cận. 1 CHƢƠNG 1: PHÉP BIẾN ĐỔI FOURIER Trong chương này, chúng tôi trình bày cơ sở lý thuyết của phép biến đổi Fourier. Blakely [17] đã khái quát những nội dung quan trọng của phép biến đổi Fourier và áp dụng chúng để phân tích, xử lý tài liệu từ và trọng lực.

Về quá trình phát triển, có thể nói Tsuboi và Fuchida [20, 21] là những người đầu tiên áp dụng các phép biến đổi Fourier vào việc minh giải các dị thường trường thế. Họ sử dụng chuỗi Fourier để chỉ ra mối liên hệ giữa các dị thường trọng lực và các phân bố khối lượng, trong cả hai trường hợp hai và ba chiều, được giới hạn trong các mặt phẳng ngang. Vào những năm 60, nhiều tác giả đã sử dụng biến đổi Fourier trong việc minh giải các dị thường từ biển, tiêu biểu là Gudmundson [9], Heirtzler và Le Pichon [13]. Tiếp đó, Harrison [12] đã đưa ra cái nhìn khái quát về chủ đề này.

Cùng thời gian đó, Bhattacharryya [5] đã công bố một số bài báo quan trọng về biến đổi Fourier của các dị thường từ và trọng lực. Có lẽ đóng góp có ý nghĩa nhất của ông là ông đã nhận ra rằng nhiều phép biến đổi, chẳng hạn như nâng trường, hạ trường, chuyển trường về cực, … dễ dàng được thực hiện trong miền tần số. Sau đây, chúng tôi sẽ tóm lược lại phép biến đổi này. Khái niệm biến đổi Fourier Một hàm tuần hoàn có thể được tổng hợp bằng tổng vô hạn các hàm sin có trọng số, trong đó các trọng số của các hàm sin được xác định qua phân tích hàm tuần hoàn đó.

Nếu f(x) là một hàm tuần hoàn với chu kỳ X, nó có thể được biểu diễn bằng: ∞ (1. Các trọng số 𝐹 trong tổng này là các số phức và được xác 𝑛 𝑋 𝑛 định bằng tích phân: 𝑥 0+𝑋 𝐹𝑛 = 1 𝑓 𝑥 𝑒 −𝑖 𝑘 𝑛 𝑥 𝑑𝑥 (1.2) 𝑋 𝑥0 2 Bây giờ chúng ta giả sử rằng f(x) không tuần hoàn trên một khoảng hữu hạn của trục x. Thay vì thế, chúng ta đòi hỏi rằng f(x) có dáng vẻ hợp lý và có biến thiên giam hãm trong một khoảng hữu hạn của trục x. Nói một cách khác, chúng ta đòi hỏi rằng: ∞ 𝑓 𝑥 𝑑𝑥 < ∞ (1.3) −∞ Các dị thường từ và trọng lực thoả mãn đặc trưng này nếu phạm vi đo đạc lớn hơn nhiều kích thước ngang của tất cả các vật thể gây dị thường.

Việc cho X → ∞ trong phương trình (1.2) tạo ra biến đổi Fourier của một hàm không chu kỳ f(x).4) −∞ Biến k trong phương trình (1.4) được gọi là số sóng và có đơn vị là nghịch đảo khoảng cách; tương tự với tần số góc trong biến đổi Fourier miền thời gian, số sóng có đơn vị nghịch đảo thời gian. Số sóng tỷ lệ nghịch với bước sóng λ, tức là: 2𝜋 𝑘= 𝜆 Chú ý từ phương trình (1.4) rằng biến đổi Fourier của hàm f(x) được ước lượng ở k = 0 đơn giản là trung bình của f(x) trên toàn trục x, tức là: +∞ 𝐹 0 = 𝑓 𝑥 𝑑𝑥 −∞ Biến đổi Fourier F(k) nói chung là một hàm phức với các phần thực và phần ảo, F(k) = ReF(k) + iImF(k). Ta cũng có thể biểu diễn: 𝐹 𝑘 = 𝐹 𝑘 𝑒 𝑖𝛩 𝑘 trong đó: 2 2 𝐹 𝑘 = 𝑅𝑒𝐹 𝑘 + 𝐼𝑚𝐹 𝑘 , 3 𝐼𝑚𝐹 𝑘 𝛩 𝑘 = 𝑎𝑟𝑐𝑡𝑔 𝑅𝑒𝐹 𝑘 Các hàm 𝐹 𝑘 và 𝛩 𝑘 tương ứng được gọi là biên độ và pha.