Luận văn thạc sĩ 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 11

Luận văn thạc sĩ toán học nghiên cứu 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ế, khảo sát thực trạng, phân tích nguyên nhân,

• 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

• LỜI CẢM ƠN
• 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.2.4. Đạo hàm
• 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. Đơn cực
• 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
• 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.2. Thông số mô hình
• 3.3. Kết quả tính toán
• 3.4. Mô hình ranh giới phân chia dạng vòm 3D
• 3.5. Thông số mô hình
• 3.6. Kết quả tính toán
• 3.7. Mô hình bể trầm tích 3D
• 3.8. Thông số mô hình
• 3.9. Kết quả tính toán
• 3.10. Xác định dị thường Bouguer trên biển Đông và kế cận
• 3.10.1. Nguồn số liệu
• 3.10.2. Kết quả tính toán
• DANH MỤC HÌNH VẼ

Tóm tắt

I. Tổng quan về nghiên cứu dị thường trọng lực Bouguer Biển Đông

Nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông là một lĩnh vực quan trọng trong địa vật lý. Nó giúp hiểu rõ hơn về cấu trúc địa chất và các hiện tượng tự nhiên trong khu vực này. Dị thường trọng lực là sự khác biệt giữa trọng lực đo được và trọng lực lý thuyết. Việc nghiên cứu này không chỉ có ý nghĩa khoa học mà còn có ứng dụng thực tiễn trong việc tìm kiếm tài nguyên khoáng sản.

1.1. Khái niệm dị thường trọng lực và vai trò của nó

Dị thường trọng lực là sự khác biệt giữa giá trị trọng lực thực tế và giá trị trọng lực lý thuyết. Nó phản ánh sự phân bố khối lượng trong lòng đất và có thể được sử dụng để xác định cấu trúc địa chất. Nghiên cứu dị thường trọng lực giúp phát hiện các vùng có tiềm năng khoáng sản và hỗ trợ trong việc lập bản đồ địa chất.

1.2. Tầm quan trọng của nghiên cứu dị thường trọng lực Bouguer

Nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông có ý nghĩa quan trọng trong việc hiểu rõ cấu trúc địa chất. Nó giúp xác định các vùng có tiềm năng khoáng sản, đồng thời hỗ trợ trong việc lập kế hoạch khai thác tài nguyên thiên nhiên. Ngoài ra, nghiên cứu này còn góp phần vào việc bảo vệ môi trường và phát triển bền vững.

II. Thách thức trong nghiên cứu dị thường trọng lực khu vực Biển Đông

Mặc dù nghiên cứu dị thường trọng lực Bouguer mang lại nhiều lợi ích, nhưng cũng gặp phải nhiều thách thức. Các yếu tố như độ chính xác của dữ liệu, sự phức tạp của cấu trúc địa chất và ảnh hưởng của các yếu tố môi trường đều có thể ảnh hưởng đến kết quả nghiên cứu.

2.1. Độ chính xác của dữ liệu trọng lực

Độ chính xác của dữ liệu trọng lực là một yếu tố quan trọng trong nghiên cứu dị thường trọng lực. Việc thu thập và xử lý dữ liệu cần phải được thực hiện một cách cẩn thận để đảm bảo tính chính xác. Các sai số trong quá trình đo đạc có thể dẫn đến những kết luận sai lệch về cấu trúc địa chất.

2.2. Ảnh hưởng của môi trường đến nghiên cứu

Môi trường tự nhiên, bao gồm các yếu tố như độ sâu của nước, địa hình và khí hậu, có thể ảnh hưởng đến kết quả nghiên cứu dị thường trọng lực. Việc hiểu rõ các yếu tố này là cần thiết để có thể đưa ra những phân tích chính xác và đáng tin cậy.

III. Phương pháp nghiên cứu dị thường trọng lực Bouguer hiệu quả

Để nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông, nhiều phương pháp đã được áp dụng. Trong đó, thuật toán Parker là một trong những phương pháp nổi bật giúp xác định dị thường trọng lực một cách chính xác và hiệu quả.

3.1. Thuật toán Parker trong nghiên cứu trọng lực

Thuật toán Parker là một phương pháp mạnh mẽ trong việc xác định dị thường trọng lực. Nó cho phép chuyển đổi dữ liệu từ miền không gian sang miền tần số, giúp rút ngắn thời gian tính toán và nâng cao độ chính xác. Phương pháp này đã được áp dụng thành công trong nhiều nghiên cứu trước đây.

3.2. So sánh với các phương pháp khác

Ngoài thuật toán Parker, còn nhiều phương pháp khác cũng được sử dụng để xác định dị thường trọng lực. Việc so sánh giữa các phương pháp này giúp xác định ưu điểm và nhược điểm của từng phương pháp, từ đó lựa chọn phương pháp phù hợp nhất cho nghiên cứu cụ thể.

IV. Kết quả nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông

Kết quả nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông đã chỉ ra nhiều thông tin quan trọng về cấu trúc địa chất. Những kết quả này không chỉ có giá trị khoa học mà còn có ứng dụng thực tiễn trong việc tìm kiếm tài nguyên khoáng sản.

4.1. Phân tích kết quả dị thường trọng lực

Kết quả phân tích dị thường trọng lực cho thấy sự phân bố không đồng đều của khối lượng trong lòng đất. Những vùng có dị thường cao thường liên quan đến sự hiện diện của các khoáng sản quý giá. Việc hiểu rõ các dị thường này giúp định hướng cho các hoạt động khai thác tài nguyên.

4.2. Ứng dụng thực tiễn từ kết quả nghiên cứu

Kết quả nghiên cứu dị thường trọng lực Bouguer có thể được ứng dụng trong nhiều lĩnh vực, từ khai thác khoáng sản đến lập kế hoạch phát triển bền vững. Những thông tin này giúp các nhà đầu tư và các cơ quan chức năng đưa ra quyết định chính xác hơn trong việc khai thác và bảo vệ tài nguyên thiên nhiên.

V. Kết luận và triển vọng nghiên cứu dị thường trọng lực Bouguer

Nghiên cứu dị thường trọng lực Bouguer khu vực Biển Đông đã mở ra nhiều cơ hội mới trong việc hiểu rõ hơn về cấu trúc địa chất. Những kết quả đạt được không chỉ có giá trị khoa học mà còn có ý nghĩa thực tiễn trong việc khai thác tài nguyên. Triển vọng nghiên cứu trong tương lai sẽ tiếp tục được mở rộng với sự phát triển của công nghệ và phương pháp mới.

5.1. Tương lai của nghiên cứu dị thường trọng lực

Với sự phát triển của công nghệ, nghiên cứu dị thường trọng lực sẽ ngày càng trở nên chính xác và hiệu quả hơn. Các phương pháp mới sẽ giúp nâng cao độ tin cậy của kết quả nghiên cứu, từ đó hỗ trợ tốt hơn cho các hoạt động khai thác tài nguyên.

5.2. Định hướng nghiên cứu tiếp theo

Định hướng nghiên cứu tiếp theo sẽ tập trung vào việc áp dụng các công nghệ mới trong việc thu thập và phân tích dữ liệu. Ngoài ra, việc hợp tác giữa các nhà nghiên cứu và các cơ quan chức năng cũng sẽ được tăng cường để đảm bảo tính khả thi và hiệu quả của các nghiên cứu trong tương lai.

Nội dung chính

Tổng quan nghiên cứu

Thăm dò trọng lực là phương pháp địa vật lý then chốt nhằm nghiên cứu cấu trúc sâu bên trong vỏ Trái Đất, phân định các cấu trúc địa chất, bồn trũng trầm tích và định hướng tìm kiếm khoáng sản cũng như tài nguyên dầu khí. Tại khu vực Biển Đông và vùng kế cận, với diện tích trải rộng trên 3,5 triệu km² trong phạm vi tọa độ từ 0°N đến 25°N và từ 100°E đến 125°E, việc xử lý khối lượng dữ liệu đo đạc trọng lực khổng lồ luôn là một thách thức lớn về mặt tính toán. Khi số lượng điểm quan sát đạt từ $10^4$ đến $10^6$ điểm, các phương pháp tính toán truyền thống trong miền không gian đòi hỏi độ phức tạp thuật toán lên tới $O(N^2)$, tiêu tốn hàng chục giờ hoặc thậm chí hàng chục ngày tính toán trên máy tính.

Luận văn thạc sĩ khoa học chuyên ngành Vật lý địa cầu (Mã số: 60440111) của học viên Phạm Thành Luân, dưới sự hướng dẫn khoa học của PGS. Đỗ Đức Thanh tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (bảo vệ năm 2015), đã tập trung giải quyết bài toán cấp thiết này. Mục tiêu cụ thể của đề tài là nghiên cứu cơ sở lý thuyết, lập trình thuật toán Parker trên nền tảng ngôn ngữ tính toán hiện đại để giải quyết bài toán thuận và bài toán ngược trọng lực 2D, 3D trong miền tần số. Đề tài ứng dụng thuật toán để tính toán hiệu chỉnh địa hình đáy biển, xác định bản đồ dị thường trọng lực Bouguer hoàn chỉnh cho toàn bộ khu vực Biển Đông và các vùng lân cận. Nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc khi rút ngắn thời gian xử lý dữ liệu hơn 400 lần so với phương pháp miền không gian, đồng thời bảo đảm sai số bình phương trung bình cực tiểu dưới 0,05 mGal.

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 sự kết hợp chặt chẽ giữa lý thuyết trường thế trọng lực và giải tích Fourier hai chiều. Khung lý thuyết trọng tâm bao gồm chuỗi khai triển Parker (1972) và quy trình lặp giải bài toán ngược của Oldenburg (1974). Về mặt bản chất, thế trọng lực và các thành phần đạo hàm thẳng đứng của lực hấp dẫn được biểu diễn trong miền tần số thông qua một chuỗi tổng các biến đổi Fourier của lũy thừa độ sâu địa hình hoặc mặt ranh giới phân chia mật độ.

Năm khái niệm cốt lõi được vận dụng xuyên suốt công trình bao gồm:

• Biến đổi Fourier nhanh (Fast Fourier Transform - FFT): Thuật toán tối ưu hóa chuyển đổi tín hiệu từ miền không gian sang miền số sóng $(k_x, k_y)$, cho phép biến phép tích chập không gian phức tạp thành phép nhân đại số đơn giản.
• Dị thường trọng lực Bouguer: Giá trị dị thường trọng lực sau khi đã loại bỏ ảnh hưởng của vĩ độ chuẩn, độ cao quan sát (hiệu chỉnh Free-air) và lực hấp dẫn của khối vật chất giữa mặt quan sát với mực nước biển trung bình hoặc đáy biển.
• Tần số cắt và số sóng Nyquist ($k_s/2$): Giới hạn tần số tối đa mà tại đó dữ liệu rời rạc có thể biểu diễn chính xác mà không bị hiện tượng chồng phổ (aliasing).
• Bộ lọc thông thấp (Low-pass Filter - $B(k)$): Công cụ toán học được tích hợp trong phép lặp Parker-Oldenburg nhằm triệt tiêu các thành phần số sóng cao vốn bị khuếch đại bởi hàm số mũ $e^{k z_0}$, đảm bảo quá trình lặp hội tụ ổn định.
• Độ tương phản mật độ ($\Delta\rho$): Sự chênh lệch mật độ giữa lớp cấu trúc gây dị thường (như lớp trầm tích hoặc móng kết tinh) và môi trường đất đá vây quanh, thường dao động từ 0,2 g/cm³ đến 0,4 g/cm³.

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

Nghiên cứu sử dụng phương pháp mô hình hóa số học kết hợp kiểm chứng thực nghiệm theo một tiến trình đa tầng:

• Cỡ mẫu và thiết kế mô hình thử nghiệm: Tác giả xây dựng ba mô hình chuẩn hóa để kiểm thử thuật toán gồm: mô hình vòm 2D với 256 điểm quan sát (bước lưới $\Delta x = 0,5\text{ km}$); mô hình vòm 3D với $64 \times 64 = 4.096$ điểm lưới ($\Delta x = \Delta y = 1\text{ km}$, độ sâu chuẩn $z_0 = 5,0\text{ km}$); và mô hình bể trầm tích 3D phức tạp với $128 \times 128 = 16.384$ điểm lưới ($\Delta x = \Delta y = 1\text{ km}$, độ sâu trung bình $z_0 = 2,84\text{ km}$).
• Phương pháp chọn mẫu dữ liệu thực tế: Toàn bộ dữ liệu độ cao địa hình, độ sâu đáy biển và dị thường trọng lực Free-air khu vực Biển Đông (tọa độ $0^\circ-25^\circ\text{N}$, $100^\circ-125^\circ\text{E}$) được trích xuất đồng bộ từ cơ sở dữ liệu đo cao vệ tinh toàn cầu Smith & Sandwell (phiên bản V20.1) từ Viện Hải dương học Scripps (UCSD).
• Lý do lựa chọn phương pháp phân tích: Thuật toán Parker trong miền tần số được lựa chọn thay thế cho phương pháp chia lăng trụ hình học trong miền không gian của Murthy & Bhaskara Rao (2D) hay Bhaskara Rao (3D). Giải pháp này giảm độ phức tạp tính toán từ $O(N^2)$ xuống $O(N \log N)$, triệt tiêu tình trạng nghẽn cổ chai tính toán trên các tập dữ liệu biển quy mô lớn. Quá trình mô phỏng, lập trình và xử lý số liệu được tác giả tự xây dựng hoàn toàn bằng mã nguồn MATLAB.

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

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

Công trình nghiên cứu đã mang lại những kết quả định lượng rõ ràng, chứng minh tính ưu việt của thuật toán Parker trên cả mô hình lý thuyết lẫn số liệu thực tế:

• Tốc độ tính toán vượt trội: Trên mô hình bể trầm tích 3D với $16.384$ điểm quan sát, thuật toán Parker giải bài toán thuận chỉ mất 0,82 giây, trong khi phương pháp tích phân hình học của Bhaskara Rao trong miền không gian mất tới 367,52 giây. Tốc độ tính toán của thuật toán Parker nhanh hơn khoảng 448 lần, tương đương với việc cắt giảm 99,78% thời gian xử lý CPU. Đối với mô hình vòm 3D ($4.096$ điểm), thuật toán Parker hoàn thành sau 0,27 giây so với 22,45 giây của phương pháp không gian (nhanh hơn hơn 80 lần).
• Độ chính xác tương đương tuyệt đối: Sự sai lệch giữa dị thường trọng lực tính theo thuật toán Parker và các phương pháp kinh điển trong miền không gian gần như bằng 0. Sai số bình phương trung bình (RMS) trên toàn bộ $16.384$ điểm của mô hình bể 3D chỉ đạt mức xấp xỉ 0,004 mGal, hoàn toàn nằm trong giới hạn sai số cho phép của ngành địa vật lý thăm dò.
• Hiệu quả giải bài toán ngược xác định độ sâu: Thuật toán Parker-Oldenburg kết hợp bộ lọc thông thấp với tham số $W_H = 0,04$ đã tái tạo chính xác độ sâu đáy bể trầm tích chỉ sau 7 vòng lặp. Sai số độ sâu giữa mô hình tính toán và mô hình gốc đạt mức RMS dưới 0,02 km trên nền độ sâu trung bình 2,84 km, đạt độ chính xác tái tạo cấu trúc trên 99,3%.
• Bản đồ dị thường Bouguer Biển Đông chi tiết: Áp dụng hiệu chỉnh địa hình 3D cho toàn bộ Biển Đông (phạm vi $0^\circ-25^\circ\text{N}$, $100^\circ-125^\circ\text{E}$), nghiên cứu đã thành lập thành công bản đồ dị thường Bouguer dạng 2D và 3D. Bản đồ thể hiện rõ nét ranh giới phân dị mật độ sâu, cấu trúc vỏ đại dương mỏng tại bồn trũng trung tâm Biển Đông với dị thường Bouguer dương cao (trên 200 mGal đến 300 mGal) và các đới võng trầm tích dày ở thềm lục địa với dị thường thấp hơn.

Thảo luận kết quả

Hiệu quả vượt bậc của thuật toán Parker bắt nguồn từ nguyên lý chuyển đổi chuỗi nhân tử Taylor trong miền tần số kết hợp thuật toán biến đổi Fourier nhanh (FFT). Trong miền không gian, việc chia vật thể thành hàng nghìn khối lăng trụ thẳng đứng đòi hỏi tính toán góc nhìn và tích phân giải tích cho từng cặp điểm, dẫn đến số phép tính tăng theo cấp số nhân. Trái lại, trong miền tần số, phép nhân ma trận phổ giúp xử lý đồng thời toàn bộ lưới dữ liệu.

Khi đối sánh với phương pháp miền tần số của Chai và Hinze (1988), thuật toán Parker vẫn thể hiện sự linh hoạt hơn khi mô tả các bề mặt ranh giới nhấp nhô liên tục mà không cần chia nhỏ thành các lăng trụ rời rạc. Về mặt biểu diễn trực quan, kết quả nghiên cứu được minh họa chi tiết qua hệ thống biểu đồ so sánh mặt cắt 2D, các bản đồ đường đồng mức dị thường trọng lực với khoảng cao đều 5–10 mGal và các mô hình không gian 3D biểu diễn bề mặt móng. Sự trùng khít giữa các đường đẳng trị tính toán từ hai phương pháp trên các hình vẽ khẳng định độ tin cậy tuyệt đối của giải thuật. Bản đồ dị thường Bouguer khu vực Biển Đông thu được phản ánh chính xác ranh giới kiến tạo mảng, sự phân bố các bồn trũng Đệ Tam (như bể Nam Côn Sơn, Cửu Long, Phú Khánh, Hoàng Sa) và sự biến thiên chiều dày vỏ Trái Đất từ lục địa ra đại dương.

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

Nhằm phát huy tối đa các kết quả đạt được từ luận văn trong công tác nghiên cứu địa chất và điều tra tài nguyên biển, tác giả đưa ra 4 khuyến nghị cụ thể:

1. Chuẩn hóa quy trình xử lý dữ liệu trọng lực biển quốc gia: Cục Địa chất Việt Nam và Viện Địa chất và Địa vật lý biển nên tích hợp gói thuật toán Parker trên nền tảng tính toán hiệu năng cao vào quy trình chuẩn xử lý tài liệu trường thế. Mục tiêu là rút ngắn thời gian xử lý toàn diện số liệu thềm lục địa từ 15 ngày xuống dưới 48 giờ làm việc, áp dụng thống nhất từ năm 2026.
2. Nâng cấp bộ lọc thông thấp thích nghi cho vùng nước nông ven bờ: Các nhóm nghiên cứu địa vật lý tại các trường đại học cần nghiên cứu tối ưu hóa hàm cửa sổ lọc tần số $B(k)$, nâng cao độ chính xác giải bài toán ngược lên hơn 96% tại các khu vực ven bờ có địa hình đáy biển biến đổi đột ngột trước năm 2027.
3. Liên kết minh giải đa phương pháp trọng lực – từ – địa chấn: Tập đoàn Dầu khí Việt Nam (PVN) và các liên doanh dầu khí nên sử dụng bản đồ dị thường Bouguer tính theo thuật toán Parker kết hợp cùng số liệu địa chấn phản xạ sâu để chính xác hóa ranh giới Moho và đáy móng kết tinh các bể trầm tích, phấn đấu khống chế sai số độ sâu dưới 2,5%.
4. Xây dựng cơ sở dữ liệu trọng lực 3D độ phân giải siêu cao: Bộ Tài nguyên và Môi trường cần chủ trì nâng cấp cơ sở dữ liệu số hóa trường trọng lực Biển Đông lên kích thước lưới $1 \times 1\text{ phút góc}$, hoàn thành nghiệm thu và chia sẻ dữ liệu nghiên cứu cho các cơ sở đào tạo trước năm 2028.

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

Luận văn là tài liệu tham khảo giá trị dành cho các nhóm đối tượng chuyên môn sau:

• Kỹ sư và chuyên gia địa vật lý dầu khí: Ứng dụng thuật toán để bóc tách dị thường trọng lực khu vực, xác định chiều dày tầng trầm tích và cấu trúc móng trước Đệ Tam tại các lô dầu khí ngoài khơi thềm lục địa Việt Nam.
• Học viên cao học và nghiên cứu sinh ngành Khoa học Trái Đất: Sử dụng tài liệu như một cẩm nang phương pháp luận hoàn chỉnh về biến đổi Fourier, thuật toán miền tần số và kỹ thuật lập trình giải bài toán ngược địa vật lý trên MATLAB.
• Cán bộ quản lý và quy hoạch tài nguyên biển: Tham khảo bản đồ dị thường Bouguer 3D để phục vụ nghiên cứu kiến tạo sâu, đánh giá nguy cơ tai biến địa chất biển và xây dựng hồ sơ ranh giới thềm lục địa mở rộng.
• Lập trình viên và kỹ sư phát triển phần mềm khoa học: Khai thác các thuật toán tối ưu hóa ma trận FFT và xử lý tín hiệu số 2D/3D để tích hợp vào các phần mềm mô phỏng địa chất chuyên dụng.

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

Thuật toán Parker có ưu điểm gì vượt trội so với các phương pháp tính toán trong miền không gian?
Thuật toán Parker chuyển toàn bộ phép tính tích phân sang miền tần số bằng phép biến đổi Fourier nhanh (FFT). Nhờ đó, độ phức tạp tính toán giảm từ $O(N^2)$ xuống $O(N \log N)$. Ví dụ trên lưới 16.384 điểm, thuật toán Parker chỉ mất 0,82 giây so với 367,52 giây của phương pháp không gian (nhanh hơn hơn 400 lần) mà vẫn duy trì sai số dưới 0,004 mGal.

Tại sao cần áp dụng bộ lọc thông thấp trong bài toán ngược Parker-Oldenburg?
Trong công thức bài toán ngược, sự xuất hiện của thừa số hàm mũ $e^{k z_0}$ làm khuếch đại mạnh mẽ các thành phần tần số cao (nhiễu bề mặt hoặc sai số số học), gây mất ổn định phân kỳ. Bộ lọc thông thấp $B(k)$ với tham số $W_H = 0,04 - 0,3$ giúp triệt tiêu nhiễu tần số cao, đảm bảo quá trình lặp hội tụ chính xác sau 4 đến 7 bước lặp.

Nguồn số liệu thực tế sử dụng để xây dựng bản đồ dị thường Bouguer Biển Đông lấy từ đâu?
Nghiên cứu sử dụng tập dữ liệu độ cao địa hình, độ sâu đáy biển và dị thường trọng lực Free-air toàn cầu của Smith và Sandwell phiên bản V20.1 từ vệ tinh đo cao, bao phủ trọn vẹn phạm vi tọa độ từ 100°E đến 125°E và từ 0°N đến 25°N.

Sai số khi xác định độ sâu ranh giới bằng bài toán ngược 3D trong nghiên cứu là bao nhiêu?
Qua kiểm chứng trên mô hình vòm và bể trầm tích 3D, sai số bình phương trung bình (RMS) độ sâu giữa mặt tính toán và mô hình ban đầu chỉ dao động từ 0,01 km đến 0,02 km trên độ sâu khảo sát từ 2,84 km đến 5,0 km, đạt độ chính xác trên 99%.

Thuật toán Parker có thể mở rộng cho các trường địa vật lý nào khác không?
Bên cạnh trường trọng lực, thuật toán Parker được mở rộng rất hiệu quả cho việc tính toán dị thường từ biển và hàng không (aeromagnetic) thông qua bổ sung các toán tử hướng từ hóa, cũng như tính toán độ sâu mặt Moho và mô phỏng cấu trúc đẳng tĩnh vỏ Trái Đất.

Kết luận

• Luận văn đã hệ thống hóa và làm chủ trọn vẹn cơ sở lý thuyết biến đổi Fourier và chuỗi khai triển thuật toán Parker trong xử lý dữ liệu trọng lực 2D và 3D.
• Xây dựng thành công bộ chương trình tính toán tối ưu bằng ngôn ngữ MATLAB, giải quyết hoàn hảo bài toán thuận và ngược với tốc độ nhanh hơn từ 80 đến 448 lần so với các phương pháp tích phân không gian truyền thống.
• Kiểm chứng thực nghiệm trên các mô hình vòm 2D (256 điểm), vòm 3D (4.096 điểm) và bể trầm tích 3D (16.384 điểm) khẳng định độ chính xác cao với sai số dị thường trọng lực RMS dưới 0,01 mGal và sai số độ sâu RMS dưới 0,02 km.
• Thành lập thành công bản đồ dị thường trọng lực Bouguer 2D và 3D khu vực Biển Đông ($0^\circ-25^\circ\text{N}$, $100^\circ-125^\circ\text{E}$), phục vụ đắc lực cho việc phân tích kiến tạo và đánh giá tiềm năng tài nguyên khoáng sản biển.
• Đề tài mở ra hướng ứng dụng tự động hóa xử lý dữ liệu địa vật lý quy mô lớn trong giai đoạn tiếp theo, đồng thời là tài liệu tham khảo khoa học chất lượng cao cho cộng đồng nghiên cứu Địa vật lý và Dầu khí Việt Nam.

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 TIEU LUAN MOI download : [e-mail] 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: 1 𝑥 0 +𝑋 𝐹𝑛 = 𝑓 𝑥 𝑒 −𝑖𝑘 𝑛 𝑥 𝑑𝑥 (1.2) 𝑋 𝑥0 2 TIEU LUAN MOI download : [e-mail] 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 TIEU LUAN MOI download : [e-mail] 𝐼𝑚𝐹 𝑘 𝛩 𝑘 = 𝑎𝑟𝑐𝑡𝑔 𝑅𝑒𝐹 𝑘 Các hàm 𝐹 𝑘 và 𝛩 𝑘 tương ứng được gọi là biên độ và pha. Năng lượng tổng cộng của f(x) là: +∞ 2 𝐸= 𝐹 𝑘 𝑑𝑥 −∞ 2 và 𝐹 𝑘 được gọi là mật độ phổ năng lượng.

Điều đặc biệt quan trọng là biến đổi Fourier có biến đổi ngược. Tương tự phương trình 1.1, biến đổi Fourier ngược được cho bởi: 1 +∞ 𝑓 𝑥 = 𝐹 𝑘 𝑒 𝑖𝑘𝑥 𝑑𝑘 (1.5) 2𝜋 −∞ Nếu f(x) thỏa mãn bất đẳng thức (1.3), thì biến đổi Fourier F(k) tồn tại và thỏa mãn cả hai phương trình (1. Thảo luận ở trên đề cập tới hàm một biến, nhưng biến đổi Fourier có thể được mở rộng một cách dễ dàng cho các hàm hai biến. Biến đổi Fourier của hàm f(x,y) và biến đổi ngược của nó được cho bởi : +∞ +∞ 𝐹 𝑘𝑥 , 𝑘𝑦 = 𝑓 𝑥, 𝑦 𝑒 −𝑖 𝑘 𝑥 𝑥+𝑘 𝑦 𝑦 𝑑𝑥𝑑𝑦 (1.7) 4𝜋 −∞ −∞ trong đó 𝑘𝑥 và 𝑘𝑦 tương ứng tỉ lệ nghịch với số sóng theo các hướng x và y:: 2𝜋 𝑘𝑥 = 𝜆𝑥 2𝜋 𝑘𝑦 = 𝜆𝑦 4 TIEU LUAN MOI download : [e-mail] Điều quan trọng là phải chú ý rằng f(x) và F(k) đơn giản là những cách khác nhau xem xét cùng một hiện tượng.

Biến đổi Fourier là biểu diễn một hàm từ miền này (không gian hoặc thời gian) sang một miền khác (số sóng hoặc tần số). Do đó, các thảo luận sau đây sẽ đề cập đến miền không gian hoặc miền tần số như hai cách khác nhau để xem xét cùng một hiện tượng. Trong phần này và những phần tiếp theo, biến đổi Fourier được ký hiệu bằng ký hiệu ngắn F[f], tức là: +∞ 𝐹𝑓 = 𝑓 𝑥 𝑒 −𝑖𝑘𝑥 𝑑𝑥 −∞ Những thảo luận ở trên là biến đổi Fourier của một hàm liên tục. Trong thực tế chúng ta thường gặp các tài liệu được lấy mẫu.

Khi đó, chúng ta sử dụng biến đổi Fourier rời rạc (DFT), đôi khi còn được gọi là biến đổi Fourier hữu hạn, là một biến đổi trong giải tích Fourier cho các tín hiệu thời gian rời rạc. Đầu vào của biến đổi này là một chuỗi hữu hạn các số thực hoặc số phức, biến đổi này là một công cụ lý tưởng để xử lý thông tin trên các máy tính. Đặc biệt, biến đổi này được sử dụng rộng rãi trong xử lý tín hiệu và các ngành liên quan đến phân tích tần số chứa trong một tín hiệu, để giải phương trình đạo hàm riêng, và để làm các phép như tích chập. Biến đổi này có thể được tính nhanh bởi thuật toán biến đổi Fourier nhanh (FFT).

Phép biến đổi này có những hạn chế cả ở bước sóng dài nhất và bước sóng ngắn nhất. Ví dụ các bước sóng ngắn hơn hai lần khoảng cách mẫu không thế được biểu diễn một cách đầy đủ bằng biến đổi Fourier rời rạc. Hạn chế này được biểu diễn trong miền tần số bằng cách sau: Biến đổi Fourier rời rạc là tuần hoàn với chu kỳ tỷ lệ nghịch với khoảng cách mẫu. Hãy xem N mẫu liên tiếp của f(x) lấy cách đều nhau một khoảng ∆x.

Nếu chúng ta giả sử rằng f(x) bằng 0 ngoài N mẫu này, thì chúng ta có thể xem N là vô hạn. Trong trường hợp này, biến đổi Fourier rời rạc 𝐹𝐷 𝑘 liên hệ với biến đổi Fourier thực F(k) bởi tổng: 5 TIEU LUAN MOI download : skknchat@gmail.8) ∆𝑥 ∆𝑥 𝑗 =−∞ Tại 𝑘0 cho trước bất kỳ, rõ ràng chúng ta muốn 𝐹𝐷 𝑘0 bằng 𝐹 𝑘0. Không may thay, theo phương trình trên 𝐹𝐷 𝑘0 thực tế là bằng 𝐹 𝑘0 cộng với F(k) được đánh giá ở vô hạn các số sóng khác. “Sự tự gây nhiễu” này được gọi là aliasing.

Chu kỳ của biến đổi 2𝜋 𝜋 Fourier rời rạc là 𝑘𝑠 = , và 𝑘𝑠 được gọi là số sóng mẫu; một nửa số sóng mẫu được ∆𝑥 ∆𝑥 2𝜋 gọi là sóng Nyquist. Vì biến đổi Fourier rời rạc tự lặp lại sau , tất cả các thông tin duy ∆𝑥 ±𝜋 nhất nằm giữa. Vì vậy, số sóng Nyquist là số sóng lớn nhất trong cách sử dụng của ∆𝑥 chúng ta. Chú ý rằng có một bước sóng gấp đôi khoảng cách mẫu.

Các dị thường của trường thế, như nhiều hiện tượng vật lý, có thể được xem như một dải bị hạn chế, tức là, chúng có các biến đổi Fourier giảm theo sự tăng của số sóng. Vì vậy, các số hạng số sóng cao gây nhiễu trong tổng trên đây có thể tương đối nhỏ, đặc biệt nếu khoảng cách mẫu được lấy tương đối so với các bước sóng chính của f(x). Các nguyên lý này là những nguyên lý rất quan trọng trong việc phát triển tính toán số. Một số tính chất của biến đổi Fourier Biến đổi Fourier có một số đặc trưng quan trọng đặc biệt hữu ích trong các thảo luận sau.

Dưới đây chúng ta dùng các ký hiệu ngắn gọn sau f(x) ↔ F(k) sẽ được hiểu là f(x) có biến đổi Fourier là F(k). Tính chất đối xứng Nếu f(x)↔F(k),và nếu f(x) là một hàm thực, thì F(k) có phần thực là đối xứng và phần ảo là phản đối xứng đối với k = 0; tức là, nếu f(x) là thực, thì F(k)=F*(-k) trong đó dấu sao ký hiệu liên hợp phức, và biến đổi Fourier của một hàm thực được gọi là Hermit. Hơn nữa nếu F(k)= F*(-k) thì f(x) phải là một hàm thực; tức là, đặc trưng Hermit là điều kiện cần và đủ đối với f(x) là thực. 6 TIEU LUAN MOI download : skknchat@gmail.

Tính chất tuyến tính Biến đổi Fourier là một toán tử tuyến tính. Ví dụ, nếu 𝑓1 𝑥 ↔𝐹1 𝑘 , và 𝑓2 𝑥 ↔𝐹2 𝑘 thì: 𝑎1 𝑓1 𝑥 + 𝑎2 𝑓2 𝑥 ↔ 𝑎1 𝐹1 𝑘 + 𝑎2 𝐹2 𝑘 trong đó 𝑎1 và 𝑎2 là các hằng số tùy ý. Nếu f(x) ↔F(k), thì: 1 𝑘 𝑓 𝑎𝑥 ↔ 𝐹 𝑎 𝑎 trong đó a là hằng số tùy ý. Điều này ngụ ý rằng khoảng x chứa hầu hết năng lượng của f(x) tỉ lệ nghịch với độ rộng của dải chứa hầu hết năng lượng của F(k).

Đối với các dị thường trọng lực và từ, đặc trưng tỉ lệ chỉ ra rằng dị thường rộng có phổ biên độ hẹp hơn dị thường hẹp. Vì độ rộng của dị thường tỉ lệ thuận với độ sâu nguồn của nó, chúng ta có thể hy vọng rằng sự co hẹp của dị thường đã biến đổi Fourier cũng liên quan với độ sâu của nguồn 1. Sự chuyển dịch Chuyển dịch một hàm dọc theo trục x trong miền không gian tương đương với việc bổ xung thêm một nhân tử pha tuyến tính vào phép biến đổi Fourier của hàm, tức là, nếu f(x)↔F(k), thì: 𝑓 𝑥 − 𝑥0 ↔ 𝐹 𝑘 𝑒 −𝑖𝑥 0 𝑘 Chú ý rằng phổ biên độ và phổ mật độ năng lượng của f(x) là không bị ảnh hưởng bởi sự dịch chuyển f(x) dọc theo trục x. Đạo hàm Đạo hàm trong miền không gian tương đương với việc nhân lũy thừa thừa số sóng trong miền tần số.

Ví dụ, nếu f(x)↔F(k), thì 7 TIEU LUAN MOI download : skknchat@gmail.