Lý thuyết cơ bản và ứng dụng của chụp ảnh địa chấn
Chụp ảnh địa chấn là một phương pháp địa vật lý được sử dụng để "chụp ảnh" bên trong Trái đất bằng sóng địa chấn. Tương tự như chụp CT trong y học, tạo ra hình ảnh ba chiều của cơ thể từ tia X, chụp ảnh địa chấn mô hình hóa các cấu trúc dưới lòng đất bằng cách sử dụng dữ liệu về thời gian truyền, biên độ hoặc dạng sóng được ghi lại bởi mạng lưới các máy đo địa chấn. Phương pháp này rất quan trọng vì bên trong Trái đất không thể quan sát trực tiếp, trong khi các động lực như kiến tạo mảng, hoạt động núi lửa và động đất lại chịu ảnh hưởng lớn bởi sự biến đổi các tính chất vật lý của đá ở độ sâu.
1. Khái niệm cơ bản về sóng địa chấn
Sóng địa chấn là sóng đàn hồi lan truyền qua đá. Sóng địa chấn thường được chia thành:
1. Sóng cơ thể
– Sóng P (sóng sơ cấp/sóng nén): lan truyền nhanh nhất, có thể xuyên qua chất rắn và chất lỏng, nhạy cảm với sự thay đổi về tốc độ nén và mật độ.
– Sóng S (sóng thứ cấp/sóng cắt): chậm hơn sóng P, không thể lan truyền trong chất lỏng, nhạy cảm với độ cứng (mô đun cắt) của đá.
2. Sóng bề mặt
Ví dụ, sóng Rayleigh và sóng Love, thường chiếm ưu thế trong các bản ghi động đất ở xa, có sự phân tán và cung cấp thông tin rất hữu ích về cấu trúc của thạch quyển xuống đến tầng mềm nông.
Tốc độ sóng địa chấn phụ thuộc vào các thông số đàn hồi (modulus) và mật độ. Sự thay đổi về nhiệt độ, thành phần khoáng chất, áp suất, độ rỗng, các vết nứt và sự hiện diện của chất lỏng hoặc các chất nóng chảy một phần sẽ ảnh hưởng đến tốc độ. Đây là cơ sở vật lý để lập bản đồ sự biến đổi tốc độ nhằm cung cấp manh mối về điều kiện địa chất dưới lòng đất.
2. Nguyên lý chụp cắt lớp: Vấn đề phía trước và vấn đề phía sau
Chụp ảnh địa chấn dựa trên hai khái niệm tính toán:
a) Bài toán thuận (forward problem)
Với một mô hình Trái Đất (ví dụ: phân bố vận tốc sóng), ta có thể tính toán dự đoán dữ liệu địa chấn: thời gian truyền sóng từ nguồn (động đất hoặc nguồn nhân tạo) đến trạm đo. Trong phương pháp lý thuyết tia đơn giản, đường đi của sóng được giả định là một tia tuân theo nguyên lý Fermat: sóng chọn con đường có thời gian truyền ngắn nhất.
Về mặt toán học, thời gian di chuyển \(T\) có thể được viết như sau:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
trong đó \(v(\mathbf{x})\) là vận tốc sóng tại vị trí \(\mathbf{x}\), và \(ds\) là phần tử chiều dài đường đi.
b) Bài toán nghịch đảo
Điều ngược lại cũng đúng: chúng ta có dữ liệu quan sát (thời gian truyền sóng, sai số thời gian truyền sóng, độ phân tán sóng bề mặt hoặc dạng sóng) và muốn ước tính mô hình vận tốc phù hợp nhất. Bài toán nghịch đảo thường không duy nhất và không ổn định: nhiều mô hình có thể phù hợp như nhau, và dữ liệu chứa đựng sự không chắc chắn/nhiễu. Do đó, cần phải có sự điều chỉnh, chẳng hạn như làm mịn, làm giảm hoặc các ràng buộc địa chất.
Trên thực tế, phép nghịch đảo thường được thực hiện theo phương pháp tuyến tính: mô hình ban đầu được làm nhiễu nhẹ, sau đó sự khác biệt về thời gian truyền được tính như sau:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
trong đó \(\delta s = \delta(1/v)\) là nhiễu loạn độ chậm. Các phương trình này được sắp xếp thành một hệ phương trình tuyến tính lớn \( \mathbf{d} = \mathbf{Gm} \), sau đó được giải bằng phương pháp bình phương tối thiểu có giảm chấn hoặc các biến thể của nó.
3. Các loại chụp cắt lớp địa chấn
1) Chụp cắt lớp theo thời gian di chuyển
Phương pháp phổ biến nhất sử dụng tín hiệu sóng P và sóng S từ nhiều trận động đất. Dữ liệu là phần dư thời gian truyền sóng so với mô hình tham chiếu (ví dụ: IASP91 hoặc ak135 trên quy mô toàn cầu). Thích hợp cho mô hình vận tốc 3D trong lớp vỏ và lớp phủ trên, đặc biệt nếu mạng lưới trạm đo dày đặc.
2) Chụp cắt lớp sóng bề mặt
Sử dụng sự phân tán (thay đổi pha/vận tốc nhóm theo chu kỳ) của sóng Rayleigh/Love. Nhạy cảm với các cấu trúc từ nông đến trung bình, phương pháp chụp cắt lớp này rất hữu ích để lập bản đồ độ dày thạch quyển, các vùng vận tốc thấp (LVZ) và sự biến đổi nhiệt độ.
3) Chụp ảnh địa chấn từ xa (teleseismic tomography)
Sử dụng động đất từ xa (địa chấn viễn thám) có sóng truyền qua lớp phủ Trái đất và được ghi lại bởi mạng lưới địa phương. Ưu điểm: chúng có nhiều nguồn và đến từ nhiều hướng khác nhau, do đó giúp "làm sáng tỏ" thể tích dưới lòng đất của khu vực nghiên cứu, ví dụ như bên dưới núi lửa hoặc các vùng hút chìm.
4) Chụp ảnh cắt lớp dạng sóng / đảo ngược toàn bộ dạng sóng (FWI)
Về mặt lý thuyết, việc sử dụng toàn bộ dạng sóng, chứ không chỉ thời gian đến, sẽ mang lại độ phân giải cao vì nó tận dụng tối đa thông tin về biên độ và pha, nhưng đòi hỏi tài nguyên tính toán lớn và một mô hình ban đầu tốt để tránh bẫy cực tiểu cục bộ.
4. Các giai đoạn chung của quá trình xử lý và đảo ngược
1. Thu thập dữ liệu
Tập hợp các bản ghi từ các máy đo địa chấn cố định hoặc tạm thời. Ở quy mô thăm dò, nguồn có thể là một vụ nổ hoặc rung chấn; ở quy mô khu vực-toàn cầu, nguồn chính là một trận động đất.
2. Nhận dạng và lựa chọn pha
Xác định thời gian đến của sóng P, sóng S hoặc sóng bề mặt. Chất lượng của việc chọn điểm quyết định chất lượng của mô hình.
3. Hiệu chỉnh và mô hình hóa ban đầu
Hiệu chỉnh thời gian (sai lệch đồng hồ), hiệu chỉnh độ cao và lựa chọn mô hình tham chiếu. Mô hình ban đầu có thể là 1D hoặc 3D đơn giản.
4. Mô phỏng dò tia / mô phỏng sóng
Tính toán đường đi của các tia hoặc dạng sóng tổng hợp để xây dựng ma trận độ nhạy.
5. Phép đảo ngược và phép chuẩn hóa
Giải hệ phương trình để thu được mô hình vận tốc. Phương pháp điều chỉnh được lựa chọn dựa trên mục tiêu: liệu có nên làm nổi bật các điểm bất thường sắc nét hay làm mượt các xu hướng.
6. Đánh giá độ phân giải và kiểm tra độ tin cậy
Ví dụ, có thể sử dụng phép thử bàn cờ, phép thử đỉnh nhọn, hoặc phân tích hàm hiệp phương sai/phân tán điểm để xem phần nào thực sự được dữ liệu phân giải rõ ràng.
5. Giải thích mô hình vận tốc
Kết quả của phép chụp cắt lớp thường là bản đồ về sự bất thường của vận tốc tương đối: các vùng có vận tốc cao thường được hiểu là đá lạnh hơn, đặc hơn hoặc cứng hơn (ví dụ: các mảng kiến tạo đang hút chìm). Các vùng có vận tốc thấp có thể cho thấy nhiệt độ cao, đá bị biến đổi, các vùng đứt gãy chứa đầy chất lỏng hoặc sự nóng chảy một phần - thường liên quan đến núi lửa đang hoạt động hoặc tầng quyển mềm nóng.
Tuy nhiên, việc giải thích cần phải thận trọng vì vận tốc bị ảnh hưởng bởi nhiều yếu tố. Lý tưởng nhất là nên kết hợp chụp ảnh địa chấn với các dữ liệu khác như trọng lực, điện từ trường (MT), trắc địa, địa chất bề mặt và thạch học.
6. Ứng dụng chụp cắt lớp địa chấn
a) Nghiên cứu về các vùng hút chìm và động lực học mảng kiến tạo
Kỹ thuật chụp cắt lớp toàn cầu và khu vực có thể lập bản đồ các mảng kiến tạo đang hút chìm trong lớp phủ, bao gồm hình dạng, độ sâu và sự phân đoạn của chúng. Thông tin này rất quan trọng để hiểu được nguồn gốc của các trận động đất lớn, cơ chế liên kết mảng kiến tạo và sự tiến hóa kiến tạo của một khu vực.
b) Hệ thống núi lửa và giảm thiểu thiên tai
Tại các khu vực núi lửa, chụp cắt lớp có thể xác định các vùng có vận tốc thấp liên quan đến buồng magma, đường dẫn chất lỏng dâng lên hoặc đá bị biến đổi. Với việc giám sát thường xuyên, sự thay đổi vận tốc (chụp cắt lớp theo thời gian) có thể cung cấp các chỉ số về chuyển động của chất lỏng/magma, rất quan trọng cho việc cảnh báo sớm phun trào.
c) Thăm dò năng lượng và tài nguyên
Ở quy mô vỏ Trái Đất nông, phương pháp chụp cắt lớp được sử dụng trong thăm dò hydrocarbon, địa nhiệt và khai thác mỏ. Sự biến đổi vận tốc giúp lập bản đồ địa tầng, cấu trúc đứt gãy, vùng nứt gãy và các tầng chứa. Trong địa nhiệt, sự kết hợp giữa chụp cắt lớp và chụp ảnh vi mô thường rất hiệu quả: chụp cắt lớp cung cấp thông tin về tính đàn hồi, trong khi chụp ảnh vi mô cung cấp thông tin về độ dẫn của chất lỏng.
d) Đặc điểm của các đứt gãy hoạt động và nguy cơ động đất
Chụp cắt lớp cục bộ có thể phát hiện các vùng yếu, vùng hư hỏng và sự không đồng nhất xung quanh các đứt gãy. Điều này giúp hiểu rõ hơn về sự phân đoạn đứt gãy, khả năng khóa chặt và sự biến đổi vận tốc ảnh hưởng đến sự khuếch đại sóng (hiệu ứng vị trí).
e) Cấu trúc vỏ Trái Đất và thạch quyển khu vực
Bằng phương pháp chụp ảnh cắt lớp sóng bề mặt, các nhà nghiên cứu có thể lập bản đồ độ dày thạch quyển, ranh giới Moho và các dị thường tầng mềm. Kết quả này tạo cơ sở cho các mô hình địa động lực học, bao gồm sự hình thành bồn địa, quá trình tạo núi và sự tiến hóa lục địa.
7. Những hạn chế và thách thức
Phương pháp chụp ảnh địa chấn phụ thuộc rất nhiều vào sự phân bố của các nguồn và trạm đo. Các khu vực có ít động đất hoặc trạm đo thưa thớt sẽ có độ chiếu sáng kém, dẫn đến độ phân giải thấp. Hơn nữa, các giả định của lý thuyết tia có thể kém chính xác hơn đối với sự không đồng nhất mạnh hoặc một số tần số nhất định. Nhiễu, lỗi chọn điểm và sự không chắc chắn về vị trí tâm chấn cũng có thể ảnh hưởng đến kết quả đảo ngược. Do đó, việc đánh giá độ phân giải và tích hợp đa phương pháp là rất quan trọng để tránh diễn giải quá mức.
Đóng cửa
Chụp ảnh địa chấn là một công cụ quan trọng để nghiên cứu cấu trúc bên trong Trái đất từ quy mô địa phương đến toàn cầu. Cơ sở lý thuyết của nó dựa trên mối quan hệ giữa các đặc tính đàn hồi của đá và sự lan truyền sóng, cũng như việc giải quyết các bài toán nghịch đảo đòi hỏi sự chuẩn hóa và kiểm tra độ phân giải. Ứng dụng của nó rất rộng: từ việc lập bản đồ các mảng kiến tạo đang hút chìm và hệ thống magma núi lửa đến thăm dò năng lượng và đánh giá rủi ro động đất. Với những tiến bộ trong mạng lưới địa chấn, điện toán và các phương pháp nghịch đảo (bao gồm cả nghịch đảo toàn sóng), chụp ảnh địa chấn ngày càng có khả năng tạo ra hình ảnh sắc nét và nhiều thông tin hơn cho cả mục đích khoa học và giảm thiểu thiên tai.
Nếu muốn, tôi có thể chỉnh sửa bài viết này thành một phiên bản mang tính học thuật hơn (có trích dẫn và thư mục tham khảo), hoặc tập trung vào một trong các ứng dụng (ví dụ: chụp cắt lớp cho núi lửa hoặc năng lượng địa nhiệt).
