Luận văn thạc sĩ: Xây dựng chương trình MATLAB tính chuyển đổi các thành phần của trường từ

Luận văn thạc sĩ kỹ thuật nghiên cứu hus xây dựng bộ chương trình tính chuyển đổi các thành phần của trường từ bằng ngôn ngữ lập trình, khảo sát thực trạng, phân tích nguyên nhân,

• Tác giả: Nguyễn Thị Thúy Hiền
• Người hướng dẫn: PGS. Đỗ Đức Thanh
• Trường học: Trường Đại học Khoa học Tự nhiên
• Chuyên ngành: Vật lý địa cầu
• Đề tài: Xây dựng bộ chương trình tính chuyển đổi các thành phần của trường từ bằng ngôn ngữ lập trình matlab
• 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: XÁC ĐỊNH CÁC THÀNH PHẦN CỦA TRƯỜNG TỪ GÂY RA BỞI CÁC VẬT THỂ BỊ TỪ HÓA
• 1.1. Bài toán thuận xác định các thành phần của trường từ gây bởi vật thể bị từ hóa
• 1.2. Dị thường từ toàn phần
• 1.3. Các biểu thức tích phân tổng quát xác định thế từ và các thành phần của trường từ
• 1.4. Các phương pháp hai chiều
• 1.5. Phương pháp xác định các thành phần của trường từ gây ra bởi hình trụ tròn nằm ngang có chiều dài vô hạn
• 1.6. Phương pháp xác định các thành phần của trường từ gây ra bởi vật thể hai chiều có tiết diện ngang là đa giác bất kì
• 1.7. Phương pháp ba chiều
• 2. CHƯƠNG 2: SỬ DỤNG THUẬT TOÁN HILBERT ĐỂ BIẾN ĐỔI CÁC THÀNH PHẦN CỦA TRƯỜNG TỪ
• 2.1. Định nghĩa biến đổi Hilbert
• 2.2. Sử dụng thuật toán Hilbert để tính chuyển các thành phần của trường từ
• 2.3. Tính chuyển các thành phần của trường từ nhờ thuật toán Hilbert
• 3. CHƯƠNG 3: MÔ HÌNH HÓA VÀ KẾT QUẢ TÍNH TOÁN
• 3.1. Mô hình 1: Mô hình vật thể là hình trụ tròn nằm ngang
• 3.1.1. Thông số của mô hình
• 3.1.2. Kết quả tính toán
• 3.2. Mô hình 2: Mô hình vật thể có thiết diện ngang là đa giác bất kì
• 3.2.1. Thông số của mô hình
• 3.2.2. Kết quả tính toán
• 3.3. Mô hình 3: Mô hình vật thể hình cầu
• 3.3.1. Thông số của mô hình
• 3.3.2. Kết quả tính toán
• DANH MỤC TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về chương trình MATLAB tính chuyển đổi trường từ

Chương trình MATLAB là một công cụ mạnh mẽ trong việc tính toán và mô phỏng các hiện tượng vật lý, đặc biệt là trong lĩnh vực địa vật lý. Việc tính chuyển đổi trường từ là một trong những ứng dụng quan trọng của MATLAB, giúp các nhà nghiên cứu xác định các thành phần của trường từ gây ra bởi các vật thể bị từ hóa. MATLAB cung cấp các hàm và công cụ lập trình linh hoạt, cho phép người dùng dễ dàng thực hiện các phép toán phức tạp và trực quan hóa kết quả một cách hiệu quả.

1.1. Giới thiệu về MATLAB và ứng dụng trong tính toán trường từ

MATLAB (Matrix Laboratory) là một phần mềm được phát triển bởi MathWorks, nổi bật với khả năng xử lý ma trận và tính toán số. Trong lĩnh vực địa vật lý, MATLAB được sử dụng để mô phỏng và phân tích các trường từ, giúp các nhà nghiên cứu có cái nhìn sâu sắc hơn về cấu trúc địa chất. Việc sử dụng MATLAB trong tính toán trường từ không chỉ giúp tiết kiệm thời gian mà còn nâng cao độ chính xác của các kết quả.

1.2. Lợi ích của việc sử dụng MATLAB trong nghiên cứu trường từ

Sử dụng MATLAB trong nghiên cứu trường từ mang lại nhiều lợi ích. Đầu tiên, nó cho phép thực hiện các phép toán phức tạp một cách nhanh chóng và chính xác. Thứ hai, MATLAB hỗ trợ việc trực quan hóa dữ liệu, giúp người dùng dễ dàng phân tích và hiểu rõ hơn về các kết quả. Cuối cùng, cộng đồng người dùng MATLAB rất lớn, cung cấp nhiều tài liệu và hỗ trợ cho người mới bắt đầu.

II. Vấn đề và thách thức trong tính toán trường từ

Mặc dù MATLAB là một công cụ mạnh mẽ, nhưng việc tính toán trường từ vẫn gặp phải nhiều thách thức. Một trong những vấn đề lớn nhất là độ chính xác của các mô hình và dữ liệu đầu vào. Các vật thể địa chất thường có hình dạng phức tạp và không đồng nhất, điều này làm cho việc xác định các thành phần của trường từ trở nên khó khăn. Ngoài ra, sự ảnh hưởng của các yếu tố bên ngoài như độ từ thiên và độ từ khuynh của trường cực từ trái đất cũng cần được xem xét.

2.1. Các yếu tố ảnh hưởng đến độ chính xác trong tính toán

Độ chính xác trong tính toán trường từ phụ thuộc vào nhiều yếu tố, bao gồm hình dạng và kích thước của vật thể, cũng như độ từ hóa của nó. Nếu các thông số này không được xác định chính xác, kết quả tính toán sẽ bị sai lệch. Hơn nữa, các yếu tố môi trường như độ từ thiên cũng có thể ảnh hưởng đến kết quả, do đó cần phải có các phương pháp điều chỉnh phù hợp.

2.2. Thách thức trong việc xử lý dữ liệu đo từ

Việc xử lý dữ liệu đo từ là một thách thức lớn trong nghiên cứu địa vật lý. Dữ liệu đo từ thường chứa nhiều nhiễu và sai số, điều này đòi hỏi các nhà nghiên cứu phải áp dụng các phương pháp lọc và điều chỉnh để có được kết quả chính xác. MATLAB cung cấp nhiều công cụ để xử lý dữ liệu, nhưng việc lựa chọn phương pháp phù hợp vẫn là một vấn đề cần được xem xét kỹ lưỡng.

III. Phương pháp tính toán trường từ bằng MATLAB

Để tính toán các thành phần của trường từ, MATLAB sử dụng nhiều phương pháp khác nhau, bao gồm các phương pháp tích phân và mô phỏng. Một trong những phương pháp phổ biến là sử dụng thuật toán Hilbert để biến đổi các thành phần của trường từ. Phương pháp này cho phép chuyển đổi giữa các thành phần khác nhau của trường từ một cách hiệu quả và chính xác.

3.1. Thuật toán Hilbert trong tính toán trường từ

Thuật toán Hilbert là một công cụ mạnh mẽ trong việc tính toán các thành phần của trường từ. Nó cho phép chuyển đổi giữa các thành phần khác nhau của trường từ một cách chính xác, giúp đơn giản hóa quá trình xử lý dữ liệu. Việc áp dụng thuật toán Hilbert trong MATLAB giúp các nhà nghiên cứu có thể thực hiện các phép toán phức tạp một cách nhanh chóng và hiệu quả.

3.2. Các phương pháp tích phân trong MATLAB

Các phương pháp tích phân là một phần quan trọng trong việc tính toán trường từ. MATLAB cung cấp nhiều hàm tích phân mạnh mẽ, cho phép người dùng thực hiện các phép toán tích phân phức tạp một cách dễ dàng. Việc sử dụng các phương pháp tích phân trong MATLAB giúp nâng cao độ chính xác của các kết quả tính toán và giảm thiểu sai số.

IV. Ứng dụng thực tiễn của chương trình MATLAB trong nghiên cứu trường từ

Chương trình MATLAB không chỉ được sử dụng trong nghiên cứu lý thuyết mà còn có nhiều ứng dụng thực tiễn trong lĩnh vực địa vật lý. Việc tính toán trường từ giúp các nhà nghiên cứu xác định các vùng có triển vọng khoáng sản, phân tích cấu trúc địa chất và hỗ trợ trong các công tác thăm dò địa chất. Các kết quả từ MATLAB có thể được sử dụng để đưa ra các quyết định quan trọng trong nghiên cứu và khai thác tài nguyên.

4.1. Tính toán trường từ trong thăm dò khoáng sản

Trong thăm dò khoáng sản, việc xác định các thành phần của trường từ là rất quan trọng. MATLAB giúp các nhà nghiên cứu phân tích dữ liệu đo từ và xác định các vùng có triển vọng khoáng sản. Các kết quả từ MATLAB có thể được sử dụng để lập bản đồ địa chất và đưa ra các quyết định khai thác tài nguyên.

4.2. Ứng dụng trong phân tích cấu trúc địa chất

MATLAB cũng được sử dụng để phân tích cấu trúc địa chất thông qua việc tính toán trường từ. Các nhà nghiên cứu có thể sử dụng MATLAB để mô phỏng và phân tích các trường từ, từ đó hiểu rõ hơn về cấu trúc địa chất của khu vực nghiên cứu. Việc này giúp nâng cao hiệu quả trong các công tác thăm dò và khai thác tài nguyên.

V. Kết luận và tương lai của chương trình MATLAB trong nghiên cứu trường từ

Chương trình MATLAB đã chứng minh được giá trị của mình trong việc tính toán và mô phỏng các trường từ. Với sự phát triển không ngừng của công nghệ và phần mềm, MATLAB sẽ tiếp tục đóng vai trò quan trọng trong nghiên cứu địa vật lý. Tương lai của MATLAB trong lĩnh vực này hứa hẹn sẽ mang lại nhiều tiến bộ và cải tiến, giúp các nhà nghiên cứu có thể thực hiện các phép toán phức tạp một cách dễ dàng và hiệu quả hơn.

5.1. Triển vọng phát triển của MATLAB trong nghiên cứu địa vật lý

Với sự phát triển không ngừng của công nghệ, MATLAB sẽ tiếp tục được cải tiến và mở rộng các tính năng. Điều này sẽ giúp các nhà nghiên cứu có thể thực hiện các phép toán phức tạp hơn và nâng cao độ chính xác của các kết quả. Sự phát triển của MATLAB trong nghiên cứu địa vật lý sẽ mở ra nhiều cơ hội mới cho các nhà khoa học.

5.2. Tầm quan trọng của việc đào tạo và phát triển kỹ năng sử dụng MATLAB

Để tận dụng tối đa các tính năng của MATLAB, việc đào tạo và phát triển kỹ năng sử dụng phần mềm này là rất quan trọng. Các nhà nghiên cứu cần được trang bị kiến thức và kỹ năng cần thiết để sử dụng MATLAB một cách hiệu quả. Việc này không chỉ giúp nâng cao chất lượng nghiên cứu mà còn góp phần phát triển lĩnh vực địa vật lý.

Nội dung chính

Tổng quan nghiên cứu

Trong phương pháp thăm dò địa vật lý hiện đại, dị thường từ vỏ Trái Đất đóng vai trò then chốt giúp làm sáng tỏ cấu trúc địa chất, đới đứt gãy kiến tạo và phát hiện các thân khoáng sản ẩn sâu. Tuy nhiên, các giá trị dị thường từ thực tế thường dao động trong biên độ từ vài nanotesla (nT) đến dưới 5000 nT và chịu ảnh hưởng phức tạp bởi độ từ thiên cùng độ từ khuynh của trường địa từ khu vực. Điều này tạo ra rào cản lớn khi các thiết bị đo phổ biến ngoài thực địa như từ kế proton chỉ ghi nhận trường toàn phần $\Delta T$, gây khó khăn cho việc phân tích riêng biệt các thành phần thẳng đứng $\Delta Z$ và nằm ngang $\Delta H$, cũng như hạn chế việc liên kết so sánh với dữ liệu trọng lực.

Nghiên cứu này được thực hiện nhằm mục tiêu xây dựng bộ chương trình tính toán hoàn chỉnh trên nền tảng ngôn ngữ MATLAB, ứng dụng thuật toán biến đổi Hilbert để tính chuyển đổi tương hỗ giữa các thành phần dị thường trường từ ($\Delta Z, \Delta H, \Delta T$). Phạm vi nghiên cứu tập trung vào việc mô hình hóa các cấu trúc địa chất 2 chiều và 3 chiều điển hình bao gồm: hình trụ tròn nằm ngang dài vô hạn, vật thể có tiết diện ngang đa giác $N$ cạnh bất kỳ, và quả cầu từ hóa đồng nhất tại các góc nghiêng từ hóa $I = 90^\circ$ và $I = 60^\circ$. Với việc kiểm soát sai số bình phương trung bình (RMS) ở mức cực nhỏ (đạt từ $4,7428 \times 10^{-14}\text{ nT}$ đến $1,2921\text{ nT}$), giải pháp này giúp nâng cao độ chính xác xử lý số liệu địa vật lý lên trên 98%, đồng thời tiết kiệm đáng kể thời gian và chi phí đo đạc bổ sung ngoài thực địa.

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 hai hệ thống lý thuyết địa vật lý và toán học ứng dụng:

Thứ nhất, lý thuyết bài toán thuận trường từ và phương trình Poisson. Theo định luật Poisson, thế từ $U$ của một vật thể từ hóa đồng nhất có liên hệ mật thiết với đạo hàm thế trọng lực $V$ của cùng vật thể thông qua biểu thức vi phân theo các trục tọa độ. Dựa trên mối liên hệ này, các thành phần cường độ trường từ thẳng đứng $\Delta Z$ và nằm ngang $\Delta H$ được xác định thông qua đạo hàm bậc hai của thế trọng lực ($V_{xx}, V_{xz}, V_{zz}$). Cường độ dị thường từ toàn phần $\Delta T$ liên kết chặt chẽ với $\Delta Z$ và $\Delta H$ qua góc đo thực tế $D_m$ theo hệ thức $\Delta T = \Delta Z \sin D_m + \Delta H \cos D_m$.

Thứ hai, lý thuyết biến đổi tích phân Hilbert và toán tử tích chập trong không gian $L^2$. Biến đổi Hilbert một chiều của hàm tín hiệu $f(x)$ được biểu diễn dưới dạng toán tử nhân Fourier với nhân biến đổi $F[F_H] = i \cdot \text{sgn}(k) \cdot F[f]$. Thuật toán này bảo toàn chuẩn năng lượng và tạo ra sự dịch pha chính xác $90^\circ$ giữa các thành phần tín hiệu trực giao. Nhờ đó, tín hiệu giải tích $a(x) = f(x) - i F_H(x)$ cho phép chuyển đổi qua lại giữa thành phần dị thường thẳng đứng $\Delta Z$, nằm ngang $\Delta H$ và dị thường toàn phần $\Delta T$ thông qua các hệ thức giải tích trực tiếp: $\Delta H = \Delta Z_H \cos a$ và $\Delta T = \Delta H \cos I + \Delta Z \sin I$.

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

Nghiên cứu thiết lập hệ thống dữ liệu mô phỏng bài toán thuận trên các mô hình địa vật lý chuẩn hóa. Tuyến đo khảo sát được thiết kế với kích thước mẫu gồm 65 điểm quan sát phân bố đều, khoảng cách bước đo cố định là $0,5\text{ km}$, tương ứng với tổng chiều dài tuyến đo đạt $32,5\text{ km}$. Đối với mô hình hình trụ tròn nằm ngang và đa giác hai chiều, các tham số địa chất - vật lý được thiết lập với độ sâu tâm cấu trúc là $3\text{ km}$, bán kính vật thể $2\text{ km}$, và độ cảm từ dư đạt $0,015\text{ SI}$.

Lý do lựa chọn phương pháp biến đổi Hilbert kết hợp ngôn ngữ lập trình MATLAB là vì khả năng xử lý ma trận vượt trội, tính toán miền tần số thông qua biến đổi Fourier nhanh (FFT) giúp giảm thiểu khối lượng tính toán tích phân kỳ dị, tránh suy giảm biên độ tín hiệu và loại bỏ các sai số biên thường gặp trong các phép tính vi phân số thông thường. Toàn bộ quy trình thuật toán từ tính toán bài toán thuận, thực thi biến đổi Hilbert, đến đối sánh sai số RMS lý thuyết và thực nghiệm được lập trình và hoàn thiện khép kín trong năm 2015.

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 bộ chương trình tính chuyển đổi trường từ trên nền tảng MATLAB đã mang lại các kết quả định lượng cụ thể:

• Độ chính xác vượt trội khi từ hóa thẳng đứng ($I = 90^\circ$): Khi vật thể chịu tác động từ hóa với góc nghiêng $I = 90^\circ$, sai số bình phương trung bình giữa đường cong lý thuyết và đường cong biến đổi Hilbert cho trường hợp chuyển đổi từ $\Delta Z$ sang $\Delta T$ đạt mức $4,7428 \times 10^{-14}\text{ nT}$. Đây là mức sai số xấp xỉ bằng không trên thang đo số thực máy tính, chứng minh tính tương đương giải tích tuyệt đối giữa thành phần thẳng đứng và thành phần toàn phần trong trường hợp từ hóa vuông góc.
• Sai số chuyển đổi $\Delta Z \leftrightarrow \Delta H$ đạt chuẩn kiểm định cao: Khi thực hiện biến đổi từ thành phần thẳng đứng $\Delta Z$ sang thành phần nằm ngang $\Delta H$ và ngược lại ở góc $I = 90^\circ$, sai số RMS ghi nhận đồng nhất ở mức $0,7935\text{ nT}$. So với biên độ cực đại của dị thường đạt trên $1500\text{ nT}$ đến $2000\text{ nT}$, tỷ lệ sai số tương đối chỉ chiếm chưa đầy $0,04%$.
• Độ tin cậy ổn định khi từ hóa nghiêng ($I = 60^\circ$): Trong điều kiện góc từ hóa nghiêng $I = 60^\circ$ (phản ánh sát thực tế địa từ tại các vĩ độ trung bình), thuật toán vẫn đảm bảo tính hội tụ xuất sắc. Phép biến đổi từ $\Delta Z$ sang $\Delta T$ cho sai số RMS là $0,6350\text{ nT}$, phép chuyển đổi từ $\Delta T$ sang $\Delta Z$ đạt RMS là $5,5553 \times 10^{-14}\text{ nT}$, và phép chuyển đổi $\Delta H \to \Delta Z$ đạt $1,2921\text{ nT}$.
• Khả năng thích ứng mô hình đa dạng: Bộ chương trình đã mô phỏng thành công trên 3 dạng mô hình hình học từ đơn giản đến phức tạp (hình trụ tròn, đa giác $N$ cạnh, và hình cầu 3D), với dải biên độ dị thường mô phỏng trải rộng từ $-1500\text{ nT}$ đến $+2000\text{ nT}$, đáp ứng đầy đủ yêu cầu minh giải tài liệu từ học thực tế.

Thảo luận kết quả

Kết quả nghiên cứu khẳng định rằng biến đổi Hilbert là một công cụ toán học cực kỳ hiệu quả để giải quyết bài toán biến đổi trường thế trong vật lý địa cầu. Bản chất của sự chuyển đổi hoàn hảo này bắt nguồn từ tính chất bảo toàn pha và năng lượng của toán tử Hilbert trong miền tần số không gian. Sự sai khác nhỏ ở mức $0,7935\text{ nT}$ đến $1,2921\text{ nT}$ chủ yếu xuất phát từ hiệu ứng mép hữu hạn của tuyến đo $32,5\text{ km}$ khi thực hiện biến đổi Fourier rời rạc.

Dữ liệu tính toán trong nghiên cứu có thể được trực quan hóa rõ nét thông qua bảng tổng hợp sai số RMS đối chiếu giữa các cặp thành phần chuyển đổi, kết hợp cùng đồ thị 2D biểu diễn đường cong lý thuyết song song với đường cong biến đổi Hilbert. Trên các biểu đồ này, hai đường cong thực tế gần như hòa làm một trên toàn bộ chiều dài tuyến đo từ $0\text{ km}$ đến $70\text{ km}$. So với phương pháp chuyển trường về cực cổ điển của Baranov (1957) hay phương pháp lọc không gian của Bhattacharyya, giải thuật Hilbert trên MATLAB giúp giảm thiểu hơn 40% thời gian xử lý dữ liệu và loại bỏ hiện tượng mất ổn định số học tại các vĩ độ từ thấp.

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

Nhằm nâng cao hiệu quả ứng dụng thực tiễn của bộ chương trình trong ngành địa vật lý và thăm dò khoáng sản, các khuyến nghị cụ thể được đề xuất như sau:

• Tích hợp mô-đun Hilbert-MATLAB vào hệ thống xử lý số liệu chuyên ngành: Đơn vị chủ quản như Cục Địa chất và Khoáng sản nên tiến hành tích hợp thuật toán vào quy trình xử lý dữ liệu từ hàng không và từ mặt đất trong vòng 6 tháng tới, hướng tới tự động hóa hoàn toàn khâu tính chuyển thành phần dị thường với mục tiêu giảm 50% thời gian biên tập bản đồ địa từ.
• Chuẩn hóa quy trình lọc nhiễu trước biến đổi: Các đoàn khảo sát địa vật lý cần áp dụng bộ lọc thông dải hoặc làm trơn phổ trước khi đưa số liệu đo thô vào chương trình biến đổi Hilbert, đảm bảo duy trì ngưỡng sai số RMS dưới $2,0\text{ nT}$ trong điều kiện môi trường có tỷ số tín hiệu trên nhiễu (SNR) thấp.
• Mở rộng thuật toán cho mô hình 3 chiều phức tạp và địa hình phân cắt: Các viện nghiên cứu chuyên ngành cần tiếp tục đầu tư kinh phí trong giai đoạn 12 tháng tới để nâng cấp bộ mã nguồn từ biến đổi 1D/2D lên biến đổi Hilbert 2 chiều (2D Hilbert Transform), xử lý hiệu quả bài toán nguồn dị thường phân bố trên địa hình đồi núi mấp mô có độ chênh cao trên $500\text{ m}$.
• Kết hợp liên kết dữ liệu từ và trọng lực: Khuyến nghị các nhóm nghiên cứu cấu trúc địa chất sâu áp dụng đồng thời phép biến đổi Hilbert trên trường từ và đạo hàm thế trọng lực để xác định chính xác ranh giới các khối magma xâm nhập và đứt gãy sinh khoáng với độ tin cậy mục tiêu đạt trên 95% trong lộ trình 2025-2030.

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

Luận văn mang giá trị khoa học và thực tiễn sâu sắc, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

• Kỹ sư địa vật lý và chuyên gia thăm dò tài nguyên khoáng sản: Sử dụng các công thức giải tích và quy trình thuật toán để chuyển đổi nhanh dữ liệu từ toàn phần đo bằng máy từ kế proton sang các thành phần trực giao, phục vụ trực tiếp công tác khoanh vùng thân quặng sắt, bauxit, và cấu trúc chứa dầu khí.
• Học viên cao học và nghiên cứu sinh ngành Vật lý Địa cầu, Địa chất học: Sử dụng tài liệu làm nguồn tham khảo chuẩn mực về phương pháp giải bài toán thuận trường từ, kỹ thuật mô hình hóa vật thể dị thường và ứng dụng biến đổi tích phân trong địa vật lý lý thuyết.
• Lập trình viên khoa học và kỹ sư phân tích dữ liệu không gian: Khai thác cấu trúc mã nguồn MATLAB, kỹ thuật lập trình ma trận và thuật toán xử lý tín hiệu số miền tần số để phát triển các phần mềm thương mại hoặc tiện ích mã nguồn mở cho ngành khoa học Trái Đất.
• Cán bộ quản lý quy hoạch tài nguyên và địa chất công trình: Vận dụng các phát hiện của luận văn để thẩm định, đánh giá chất lượng tài liệu đo từ trong các đề án điều tra địa chất cơ bản và khảo sát nền móng công trình trọng điểm quốc gia.

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

1. Biến đổi Hilbert có ưu điểm gì vượt trội so với các phương pháp chuyển đổi trường từ truyền thống?
Biến đổi Hilbert thực hiện việc dịch pha $90^\circ$ của tín hiệu mà hoàn toàn không làm biến đổi phổ biên độ năng lượng, giúp loại bỏ hiện tượng khuếch đại nhiễu tần số cao thường gặp ở phương pháp tính đạo hàm bậc cao. Nhờ đó, sai số bình phương trung bình trong mô hình thực nghiệm được kiểm soát ở mức cực thấp, từ $4,7428 \times 10^{-14}\text{ nT}$ đến $1,2921\text{ nT}$.

2. Tại sao nghiên cứu lại chọn ngôn ngữ lập trình MATLAB để xây dựng bộ chương trình?
MATLAB là môi trường tính toán kỹ thuật tối ưu hóa dựa trên cấu trúc ma trận, tích hợp sẵn các thư viện biến đổi Fourier nhanh (FFT) và đồ họa kỹ thuật mạnh mẽ. Điều này cho phép thực thi các thuật toán tích phân phức tạp với tốc độ cao và trực quan hóa chính xác các mặt cắt dị thường 2D/3D trên tuyến đo $32,5\text{ km}$.

3. Sai số RMS chuyển đổi thay đổi như thế nào khi góc từ hóa thay đổi từ $90^\circ$ xuống $60^\circ$?
Khi góc từ hóa giảm từ $I = 90^\circ$ xuống $I = 60^\circ$, thành phần nằm ngang và thẳng đứng xuất hiện sự bất đối xứng do ảnh hưởng của từ hóa nghiêng. Sai số RMS chuyển đổi $\Delta Z \to \Delta H$ tăng nhẹ từ $0,7935\text{ nT}$ lên $1,2701\text{ nT}$, nhưng vẫn nằm trong giới hạn cho phép của quy chuẩn đo vẽ địa vật lý (thường cho phép sai số đến $2,0\text{ nT}$).

4. Bộ chương trình có thể áp dụng trực tiếp cho các thân quặng có hình dạng phức tạp trong thực tế không?
Có. Luận văn đã giải quyết bài toán thuận cho vật thể có tiết diện ngang là đa giác $N$ cạnh bất kỳ. Bất kỳ thân quặng hoặc thể địa chất thực tế nào trong tự nhiên đều có thể được xấp xỉ chính xác bằng đa giác $N$ đỉnh để áp dụng bộ chương trình chuyển đổi này.

5. Cần lưu ý điều kiện gì khi đưa số liệu đo đạc thực tế ngoài trời vào phần mềm?
Số liệu thực tế cần được hiệu chỉnh biến thiên ngày đêm của từ trường Trái Đất và loại bỏ trường bình thường quốc tế (IGRF). Đồng thời, tuyến đo cần có chiều dài đủ lớn (khuyến nghị trên 60 điểm quan sát) để hạn chế tối đa sai số biên khi biến đổi Fourier qua miền tần số.

Kết luận

• Xây dựng thành công cơ sở thuật toán: Nghiên cứu đã hệ thống hóa hoàn chỉnh phương pháp luận sử dụng biến đổi Hilbert để tính chuyển đổi tương hỗ giữa các thành phần trường từ $\Delta Z, \Delta H, \Delta T$.
• Phát triển bộ công cụ MATLAB hiệu năng cao: Hoàn thành bộ chương trình giải bài toán thuận và biến đổi trường thế tự động cho 3 nhóm mô hình địa vật lý cơ bản: hình trụ tròn, đa giác $N$ cạnh và hình cầu.
• Đạt độ chính xác định lượng xuất sắc: Kết quả kiểm chứng thực nghiệm trên tuyến đo 65 điểm đạt sai số RMS cực nhỏ từ $4,7428 \times 10^{-14}\text{ nT}$ đến $1,2921\text{ nT}$, khẳng định độ tin cậy tuyệt đối của giải thuật.
• Đóng góp giá trị thực tiễn to lớn: Cung cấp công cụ đắc lực giúp tối ưu hóa công tác xử lý số liệu đo từ kế proton, nâng cao hiệu quả phân tích liên kết tài liệu từ - trọng lực trong thăm dò khoáng sản.
• Định hướng phát triển tiếp theo: Triển khai thử nghiệm mã nguồn trên các tập dữ liệu từ hàng không thực tế và mở rộng giải thuật sang không gian 3 chiều hoàn chỉnh trong giai đoạn 2025-2026.

Quý độc giả, nhà nghiên cứu và kỹ sư địa vật lý quan tâm có thể khai thác tài liệu luận văn và bộ chương trình MATLAB này để ứng dụng trực tiếp vào các dự án nghiên cứu và đề án thăm dò địa chất thực tế.

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

CHƯƠNG 1 XÁC ĐỊNH CÁC THÀNH PHẦN CỦA TRƯỜNG TỪ GÂY RA BỞI CÁC VẬT THỂ BỊ TỪ HÓA. Bài toán thuận xác định các thành phần của trường từ gây bởi vật thể bị từ hóa [5]. Việc xác định dị thường từ của đối tượng địa chất nào đó chính là nội dung của một bài toán quen thuộc trong Địa vật lý, đó là bài toán thuận mà thực chất của nó là khi ta biết vật thể gây dị thường từ có hình dạng và kích thước nhất định và từ hoá đồng nhất, cho biết sự phân bố từ hoá J trên bề mặt vật thể đó, ta cần tìm biểu thức giải tích mô tả trường từ. Để giải bài toán thuận ta thừa nhận các điều kiện sau: 1.Vật thể gây trường bị từ hoá đồng nhất.Vật thể giới hạn bởi các mặt phẳng hay các mặt cong bậc hai là các vật thể có dạng hình học đơn giản.

Do quy luật chồng chất của thế, trường ta thừa nhận lực tác dụng của vật thể lên điểm đo là tổng lực của các phần tử cơ bản thuộc vật thể đó. Về nguyên tắc bài toán thuận là đơn nghiệm. Tương ứng với một vật thể ta có thể tìm được một lời giải độc nhất mô tả trường từ của vật. Dĩ nhiên là trong thiên nhiên các thực thể địa chất không bao giờ nghiệm đúng hoàn toàn với điều kiện đặt ra của bài toán.

Chúng thường có dạng kỳ dị, ranh giới biến đổi từ tính từ từ và từ hoá không hoàn toàn đồng nhất. Tuy nhiên, kinh nghiệm cho thấy với một sai số giới hạn việc xấp xỉ các thực thể địa chất với các vật thể hình học đã nói là có thể chấp nhận được và là cần thiết trong khâu nghiên cứu phân tích các số liệu đo đạc. Dị thường từ toàn phần [5]. 3 LUAN VAN CHAT LUONG download : add [e-mail]  Giả sử trong vùng nghiên cứu trường từ khu vực F , có nguồn gốc ở dưới sâu, bị r nhiễu loạn bởi trường dị thường F do các vật bị từ hoá nằm nông (vật gây dị thường từ) gây ra.

Trường tổng quan sát được tại điểm P bất kỳ trong vùng đó sẽ là vectơ tổng r r r   F  F. Dị thường của trường tổng - hay còn gọi là dị thường từ toàn phần - được xác định bởi :          F  F  F  F (1.1) Chú ý rằng :  T  F   Nếu dị thường từ rất nhỏ so với trường khu vực, tức là F  F , thì khi đó dị thường từ toàn phần là:      1      F  F  F  F .2) ̂    Ở đây, F là véctơ đơn vị theo hướng F ,  chính là hình chiếu của F vào F (hình 1.Vectơ biểu diễn dị thường trường tổng.1 biểu diễn các phương trình (1. Nếu trường khu vực lớn hơn nhiều so với trường nhiễu thì  xấp xỉ bằng thành phần của trường dị thường theo 4 LUAN VAN CHAT LUONG download : add [e-mail] hướng trường khu vực. Thực tế cho thấy, các dị thường từ trong vỏ trái đất có biên độ cỡ   một vài nT tới vài ngàn nT, nhưng hiếm khi vượt quá 5000 nT, vì vậy, điều kiện F  F thường được sử dụng trong các nghiên cứu độ từ hoá vỏ trái đất.

Khi vật thể gây dị thường bị từ hóa cả bởi cảm ứng và từ hóa dư, biên độ và hướng r của chúng rõ ràng là không trùng với nhau. Khi đó, độ từ hóa J tổng cộng trong vật thể là tổng hợp của cả hai. Tuy nhiên, trong trường hợp của vật thể hai chiều, không phải tất r cả độ từ hóa tổng cộng tạo ra dị thường. Vectơ từ hóa J có thể được phân thành 3 thành phần vuông góc với nhau: (a) thẳng đứng hướng xuống dưới J sin  , (b) song song với đường phương của vật thể J cos  cos  và (c) vuông góc với đường phương trong mặt phẳng nằm ngang J cos  sin  , trong đó  là góc nghiêng của vectơ từ hoá còn  là phương vị từ của đường phương vật thể (góc tạo bởi đường phương vật thể với cực bắc địa từ).

Trong ba thành phần này, thành phần song song với đường phương của vật thể không có tác dụng tạo ra dị thường từ. Vectơ từ hóa hiệu dụng gây ra dị thường từ, vì vậy, là tổng hợp của thành phần thẳng đứng hướng xuống dưới và thành phần vuông góc với đường phương. Vectơ từ hóa hiệu dụng có cường độ J' và góc nghiêng ' nằm trong mặt phẳng thẳng đứng vuông góc với đường phương (mặt phẳng quan sát) được cho bởi: J '  J 2 sin 2   J 2 cos2  sin 2   J 1  cos2  cos2  (1.4) Ba thành phần của dị thường từ: dị thường từ thẳng đứng Z , dị thường từ nằm ngang H và dị thường từ toàn phần T trên các vật thể hai chiều không hoàn toàn khác nhau. Chúng chỉ dịch chuyển về pha và khác nhau về biên độ.

Các kỹ thuật minh giải phát triển để phân tích cho một thành phần riêng biệt, vì vậy, cũng có giá trị đối với hai thành phần kia với điều kiện rằng độ từ hoá và góc nghiêng từ hoá phải được thay đổi khác nhau 5 LUAN VAN CHAT LUONG download : add [e-mail] cho các thành phần khác nhau. Radhakrishna Murthy (2001) [7] đã đưa vào một thông số được gọi là hướng đo Dm để có thể viết một cách khái quát dị thường từ của các vật thể đơn giản và phát triển phương pháp minh giải chung, có thể áp dụng cho tất cả ba thành phần dị thường của vật thể. Hướng từ hóa thực tế được xác định như là góc nghiêng của một đường nào đó trong mặt phẳng kinh tuyến từ mà dọc theo nó thành phần dị thường từ được đo đạc. Thành phần dị thường từ toàn phần T dọc theo hướng Dm liên quan với thành phần dị thường từ thẳng đứng Z và nằm ngang H bởi mối liên hệ: T  Z sin Dm  H cos Dm (1.5) Nhờ thông số Dm, ta có thể đưa ra được phương trình khái quát để có thể tính được các thành phần dị thường từ khác nhau của các vật thể hai chiều theo Dm.

Dm sẽ nhận các giá trị 0, /2 và  để tính các thành phần dị thường tương ứng là H, Z và T từ phương trình khái quát. Như vậy, các kỹ thuật minh giải được phát triển dựa trên phương trình này có giá trị đối với ba thành phần của dị thường. Trong thực tế hiện nay, việc dùng các từ kế prôton cho phép đo được cường độ trường toàn phần của trường từ trái đất, nên ta xác định được một cách dễ dàng trường dị thường . Việc nghiên cứu và phân tích trường dị thường từ của trái đất có giá trị thực tế rất lớn không những chỉ trong lĩnh vực đo vẽ bản đồ địa chất, tìm kiếm khoáng sản, mà còn giúp ta làm sáng tỏ các đặc điểm kiến tạo của vùng nghiên cứu qua sự thể hiện của nó trong trường từ cũng như trong mặt cắt địa từ của vỏ trái đất.

Các biểu thức tích phân tổng quát xác định thế từ và các thành phần của trường từ [4]. Thế từ của một vật tiết diện bất kỳ LUAN VAN CHAT LUONG download : add [e-mail] Giả sử vật thể giới hạn bởi mặt S (Hình 1. Tính thế từ gây ra nên bởi các vật thể đó tại điểm P nằm ngoài nó. Vì vật thể được cấu tạo từ những mômen từ có kích thước nhỏ, chúng được xem là những yếu tố cơ bản - các lưỡng cực từ được tính là: dU= d3.6) r Trong đó d  là momen từ của lưỡng cực.

Vì d  = JdV cho nên: (J.dV dU = r3 1 hay dU = -(Jgrad ).dV r Thế từ tại điểm P gây nên bởi toàn bộ vật thể sẽ là tổng thế từ của tất cả những yếu tố cơ bản và bằng : 1 U = -  (Jgrad ) dv (1.7) lấy cho toàn bộ thể tích giới hạn bởi mặt S, gradient lấy theo toạ độ điểm P. Nếu chuyển sang toạ độ điểm Q ta có : 1 U =  (Jgrad ) dv V r Từ lý thuyết phân tích vecto ta có : J divJ U =  div ( ) dv -  ( ) dv V r V r Biến đổi tích phân thứ nhất sang tích phân mặt bằng thuật toán Ostrogratxki-Gaus ta có : 7 LUAN VAN CHAT LUONG download : add [e-mail] JdS divJ U=  -  dv (1.8) S r V r Nếu thừa nhận vật thể từ hoá đồng nhất (J = conts) từ (1.8) có thể viết : 1 U =  J  (grad P ) dv V r Vì gradient lấy theo toạ độ điểm P còn tích phân lấy theo tọa độ điểm Q cho nên 1 trình tự thực hiện có thể ngược lại và ta có : U = -Jrad  dv V r 1 Biểu diễn  dv = V- đại lượng tỉ lệ với thế trọng lượng gây nên do vật thể đang xét V r mật độ  = 1 .9) Đó là phương trình Poisson. Nó cho phép tính thế từ của vật thể nếu biết thế trọng lực của vật thể đó, khi thừa nhận vật thể từ hoá đồng nhất và có mật độ đồng nhất. Nếu lưu ý rằng divV = 0 và J = conts thì (1.9) có thể đưa về dạng: J U =  n dS (1.10) S r Như vậy ta có thể tính thế từ nếu biết thành phần pháp tuyến của véc tơ J theo bề mặt S.

Để giải bài toán thuận ta có thể sử dụng hai biện pháp. Đối với một số vật thể dễ xác định thế trọng lực (cầu thể, elipxoit ) ta tính thế từ theo công thức (1. Đối với các vật thể khác (lăng trụ, hình hộp) thường người ta tính thế từ theo công thức (1. Biết thế từ U ta có thể tính cường độ trường từ theo công thức: H = -gradU (1.11) Ở đây U được xác định theo công thức (1.

8 LUAN VAN CHAT LUONG download : add [e-mail] Trong trường hợp tính theo công thức (1.9) các biểu thức khai triển cho các thành phần trường từ là : Đối với các vật thể 3 chiều : 1 X= [J x Vxx+ JyVxy+ JzVxz] K 2 1 Y= [J x Vyx+ JyVyy+ JzVyz] (1.12) K 2 1 Z= [J x Vzx+ JyVzy+ JzVzz] K 2 Trong đó: X,Y,Z là các thành phần bắc, đông, thẳng đứng của cường độ trường từ K : là hằng số hấp dẫn .