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

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ả.

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. Năng lượng tổng cộng của f(x) là: +∞ 𝐸= 𝐹 𝑘 2𝑑𝑥 −∞ và 𝐹 𝑘 2đượ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) 2 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 Đ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 +∞ 1 2𝜋𝑗 𝐹𝐷 𝑘 = 𝐹 𝑘– (1.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 𝑘𝑠 ∆𝑥 � 6 ∆𝑥 đượ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. 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ì 8 𝑑𝑛 𝑓 𝑥 ↔ 𝑖𝑘 𝑛 F 𝑘 (1.9) 𝑑𝑥 𝑛 Nếu hàm phụ thuộc vào hai biến và nếu f(x,y)↔F(𝑘𝑥 , 𝑘𝑦 ), thì: 𝜕𝑛 𝜕 𝑓 𝑥, 𝑦 = 𝑖𝑘𝑥 𝑛 𝑖𝑘𝑦 𝐹 𝑘𝑥 , 𝑘𝑦 � (1.