Phân rã địa chấn khác gì so với lọc tần số địa chấn thông thường?

Lọc tần số địa chấn thông thường chỉ triệt tiêu hoặc giữ lại các dải tần số cố định trên toàn bộ đường ghi giả định tín hiệu là dừng. Ngược lại, phân rã địa chấn phân tách tín hiệu thành các thành phần thời gian tần số tức thời hoặc các nguyên tử sóng thích ứng, cho phép bảo toàn các đặc trưng phi dừng và định lượng các biến đổi địa tầng cục bộ.

Hiện tượng bóng tần số thấp trong phân rã địa chấn giúp ích gì cho tìm kiếm dầu khí?

Khi sóng địa chấn truyền qua tầng chứa khí tự nhiên, năng lượng tần số cao bị hấp thụ mạnh bởi tính nhớt đàn hồi của chất lưu, tạo ra vùng năng lượng dị thường tập trung ở dải tần số thấp từ 10 Hz đến 20 Hz ngay bên dưới vỉa chứa. Hiện tượng này đóng vai trò như một dấu hiệu chỉ thị hydrocarbon trực tiếp giúp các nhà địa chất khoanh định vỉa khí tiềm năng.

Giới hạn độ phân giải của phân rã địa chấn là bao nhiêu?

Độ phân giải của phân rã địa chấn bị chi phối bởi nguyên lý bất định thời gian tần số, trong đó không thể đồng thời đạt độ phân giải thời gian và độ phân giải tần số cực đại tùy ý. Mặc dù vậy, theo nghiên cứu của Hao và cộng sự (2023), các phương pháp phân rã thưa tiên tiến sử dụng từ điển sóng Ricker cho phép tách biệt các tập vỉa xen kẹp mỏng nằm sâu dưới giới hạn điều chỉnh cộng hưởng cổ điển.

Phân rã địa chấn (seismic decomposition) là gì?

Phân rã địa chấn là phương pháp toán học và xử lý tín hiệu địa vật lý nhằm phân tách một đường ghi địa chấn tổng hợp thành các thành phần sóng cơ bản, các nguyên tử thời gian tần số hoặc các dải tần số đơn lẻ. Phương pháp này cho phép các nhà khoa học Trái Đất và kỹ sư địa vật lý vượt qua giới hạn phân giải giao thoa truyền thống để định lượng các dị thường địa tầng mỏng, phân tích cấu trúc vỉa chứa và đánh giá hệ số hấp thụ năng lượng sóng đàn hồi dưới sâu. Bài viết trình bày cơ sở lý thuyết, các thuật toán phân rã chủ đạo, ứng dụng thực tiễn trong thăm dò dầu khí và địa chất công trình, cùng những giới hạn kỹ thuật cần lưu ý.

Bản chất khoa học và nguyên lý phân rã địa chấn

Trong địa chấn học thăm dò, tín hiệu phản xạ địa chấn ghi nhận trên mặt đất là sự tích chập giữa xung sóng nguồn phát và chuỗi hệ số phản xạ tại các ranh giới thạch học của các tầng đất đá trong lòng đất. Theo chuyên khảo kinh điển của Yilmaz (2001), xử lý số liệu địa chấn bao gồm các công đoạn lọc ngược, cộng sóng và dịch chuyển nhằm chuyển đổi chuỗi các đường ghi địa chấn thô thành mặt cắt cấu trúc địa chất tin cậy. Khi sóng đàn hồi truyền qua các cấu trúc địa chất phân lớp phức tạp, hiện tượng giao thoa sóng giữa nóc và đáy của các lớp đất đá có bề dày mỏng gây ra sự méo mó biên độ và biến dạng pha nghiêm trọng. Tình trạng này khiến việc nhận diện trực quan các ranh giới địa chất trên các mặt cắt thời gian thông thường gặp nhiều trở ngại.

Phân rã địa chấn giải quyết thách thức trên bằng cách biểu diễn đường ghi địa chấn dưới dạng tổ hợp tuyến tính của các dạng sóng thành phần có đặc tính thời gian và tần số xác định. Biểu thức tổng quát của phép phân rã một đường ghi địa chấn $s(t)$ theo thời gian $t$ được thiết lập như sau:

s(t) = \sum_{k} a_k \psi_k(t) + r(t)

Trong công thức trên, đại lượng $s(t)$ biểu thị tín hiệu địa chấn ghi nhận tại một trạm đo, đại lượng $a_k$ là trọng số biên độ tương ứng với thành phần thứ $k$ , đại lượng $\psi_k(t)$ là hàm sóng cơ bản hoặc nguyên tử thời gian tần số được trích xuất từ một từ điển toán học, và đại lượng $r(t)$ đại diện cho tín hiệu phần dư hoặc nhiễu không tương quan sau quá trình phân rã.

Khác biệt căn bản giữa phân rã địa chấn và các phép lọc tần số tuyến tính cổ điển nằm ở khả năng bảo toàn thông tin phi dừng. Các đặc trưng địa chấn của tầng đất đá biến đổi liên tục theo cả chiều sâu chôn vùi và khoảng cách đo đạc. Phân rã địa chấn không giả định tín hiệu là dừng trên toàn bộ chiều dài đường ghi, mà tập trung định vị các biến thiên năng lượng tức thời tại từng vị trí phản xạ cục bộ.

Các phương pháp phân rã địa chấn chủ đạo

Sự phát triển của lý thuyết xử lý tín hiệu số đã hình thành ba nhóm kỹ thuật phân rã địa chấn chủ yếu, bao gồm phân rã phổ thời gian tần số, phân rã theo từ điển sóng nguyên tử và phân rã thích ứng dựa trên chế độ dữ liệu.

1. Phân rã phổ dựa trên các biến đổi tích phân

Năm 1999, Partyka và các cộng sự đã thiết lập một bước tiến nền tảng khi công bố phương pháp phân rã phổ địa chấn ứng dụng trong minh giải vỉa chứa. Bằng cách sử dụng biến đổi Fourier cửa sổ ngắn hoặc phép biến đổi rời rạc, tín hiệu địa chấn được chuyển đổi sang miền tần số để sinh ra các lát cắt biên độ và pha tại các tần số đơn lẻ. Phương pháp này cho phép nhận diện các kênh cát cổ bị vùi lấp và các bẫy địa tầng có bề dày nằm dưới giới hạn phân giải Rayleigh.

Năm 1996, Stockwell và các cộng sự đã xây dựng biến đổi S nhằm khắc phục nhược điểm cửa sổ cố định của biến đổi Fourier cửa sổ ngắn. Biến đổi S của một tín hiệu địa chấn $s(t)$ được định nghĩa thông qua tích phân sau:

S(\tau, f) = \int_{-\infty}^{\infty} s(t) \frac{|f|}{\sqrt{2\pi}} \exp\left(-\frac{(t-\tau)^2 f^2}{2}\right) \exp(-i 2\pi f t) dt

Trong đó, biến số $\tau$ biểu diễn vị trí thời gian của tâm cửa sổ Gauss, biến số $f$ là tần số phân tích, số vô tỷ $\pi$ là hằng số Archimedes, số $e$ là cơ số logarit tự nhiên, và ký hiệu $i$ là đơn vị ảo thỏa mãn bình phương bằng trừ một. Hàm cửa sổ chuẩn Gauss trong biến đổi S có độ rộng tỷ lệ nghịch với tần số, tạo ra độ phân giải thời gian tối ưu ở vùng tần số cao và độ phân giải tần số vượt trội ở vùng tần số thấp, đồng thời giữ nguyên liên hệ pha trực tiếp với phổ Fourier cổ điển.

2. Phân rã từ điển sóng bằng giải thuật truy đuổi khớp

Năm 1993, Mallat và Zhang đã đề xuất thuật toán truy đuổi khớp phục vụ phân rã tín hiệu trong các từ điển thời gian tần số dư thừa. Trong xử lý dữ liệu địa chấn, từ điển sóng thường được xây dựng từ tập hợp các xung sóng Ricker hoặc sóng Morlet với các tham số biến thiên về tần số trung tâm, vị trí thời gian và pha ban đầu.

Thuật toán hoạt động theo nguyên tắc xấp xỉ tham lam lặp từng bước. Tại mỗi chu trình tính toán, thuật toán quét qua toàn bộ từ điển để chọn ra nguyên tử sóng có tích vô hướng cực đại với tín hiệu hiện tại, trừ thành phần đó ra khỏi dữ liệu để thu được phần dư mới, và lặp lại quá trình cho đến khi năng lượng phần dư giảm xuống dưới ngưỡng quy định. Phương pháp này triệt tiêu hoàn toàn các thành phần giao thoa giả mạo vốn thường xuất hiện trong các biến đổi song tuyến tính như phân bố Wigner-Ville.

3. Phân rã chế độ thực nghiệm thích ứng

Năm 2013, Han và van der Baan đã nghiên cứu ứng dụng phương pháp phân rã chế độ thực nghiệm trong phân tích thời gian tần số dữ liệu địa chấn. Điểm đặc thù của phương pháp này là không yêu cầu lựa chọn bất kỳ hàm cơ sở toán học định sẵn nào trước khi tính toán. Thay vào đó, thuật toán trực tiếp trích xuất các hàm chế độ nội tại từ chính dữ liệu ghi nhận thông qua quá trình đãi lọc sóng dựa trên các điểm cực trị cục bộ.

Mỗi hàm chế độ nội tại đại diện cho một thành phần dao động đơn giản với tần số tức thời có ý nghĩa vật lý. Sau khi tách các hàm chế độ nội tại, phép biến đổi Hilbert được áp dụng để thiết lập phổ Hilbert-Huang, cung cấp độ tập trung năng lượng cao trên mặt phẳng thời gian và tần số, hỗ trợ đắc lực cho việc phát hiện các biến đổi thạch học phi tuyến.

Bảng so sánh các phương pháp phân rã địa chấn

Mỗi thuật toán phân rã sở hữu những ưu thế toán học và hạn chế kỹ thuật riêng biệt khi áp dụng cho dữ liệu thực tế:

Phương pháp phân rã	Cơ chế toán học chủ đạo	Ưu điểm nổi bật	Hạn chế chính
Biến đổi Fourier cửa sổ ngắn	Tích chập với cửa sổ thời gian có độ dài cố định	Tốc độ tính toán nhanh, dễ chuẩn hoá quy trình công nghiệp	Độ phân giải thời gian và tần số bị ràng buộc cố định trên toàn dải phổ
Biến đổi S	Cửa sổ Gauss biến thiên bề rộng nghịch đảo theo tần số	Bảo toàn liên hệ pha tuyệt đối, độ phân giải thích ứng theo tần số	Chi phí tính toán cao hơn Fourier, nhạy cảm với biên độ nhiễu ngẫu nhiên
Truy đuổi khớp	Phân rã tham lam lặp từng bước theo từ điển sóng dư thừa	Độ phân giải cực cao, không tạo thành phần giao thoa giả mạo	Khối lượng tính toán rất lớn, kết quả phụ thuộc vào thiết kế từ điển sóng
Phân rã chế độ thực nghiệm	Đãi lọc thích ứng trực tiếp từ các cực trị của đường ghi	Hoàn toàn phi tham số, thích ứng hoàn hảo với tín hiệu phi dừng	Dễ gặp hiện tượng trộn chế độ khi tín hiệu chứa các xung gián đoạn

Ứng dụng thực tiễn trong khoa học Trái Đất và dầu khí

Phân rã địa chấn đóng vai trò trung tâm trong nhiều quy trình minh giải địa chấn hiện đại, mang lại giá trị gia tăng rõ rệt so với các thuộc tính biên độ truyền thống.

1. Nhận diện các thể cát vỉa chứa mỏng và bẫy địa tầng

Theo tiêu chí phân giải Rayleigh, hai mặt ranh giới phản xạ chỉ có thể phân biệt độc lập khi khoảng cách giữa chúng không nhỏ hơn một phần tư bước sóng địa chấn chiếm ưu thế. Khi chiều dày vỉa cát mỏng hơn ngưỡng này, các xung sóng phản xạ giao thoa và tạo ra hiệu ứng điều chỉnh cộng hưởng. Partyka và các cộng sự (1999) đã chứng minh rằng tần số cộng hưởng điều chỉnh tỷ lệ nghịch với chiều dày vỉa đất đá.

Bằng cách phân rã địa chấn thành các lát cắt tần số đơn lẻ, các nhà địa chất có thể lập bản đồ phân bố các thân cát hình dẻ quạt dưới đáy biển hoặc các kênh lòng sông cổ bị vùi lấp với độ chính xác cao. Theo nghiên cứu của Hao và cộng sự (2023), việc tích hợp thuật toán phân rã địa chấn thưa sử dụng từ điển sóng Ricker cho phép bóc tách hiệu quả các tập trầm tích mỏng nằm sâu dưới giới hạn điều chỉnh cộng hưởng truyền thống.

2. Phát hiện trực tiếp bẫy khí bằng hiện tượng bóng tần số thấp

Năm 2003, Castagna và các đồng nghiệp đã làm sáng tỏ hiện tượng bóng tần số thấp xuất hiện ngay phía dưới các vỉa chứa hydrocarbon. Khi sóng nén địa chấn lan truyền qua tầng đá rỗng bão hoà chất lưu dạng khí, tính nhớt đàn hồi của môi trường gây ra sự suy giảm năng lượng nghiêm trọng ở các thành phần tần số cao.

Kết quả phân tích phổ tức thời của Castagna và cộng sự (2003) ghi nhận sự suy giảm mạnh mẽ của các thành phần tần số cao, trong khi năng lượng phản xạ tại dải tần số thấp từ 10 Hz đến 20 Hz được duy trì hoặc gia tăng tương đối. Hiện tượng dị thường bóng tần số thấp trở thành một dấu hiệu chỉ thị hydrocarbon trực tiếp có giá trị cao trong đánh giá rủi ro địa chất. Bên cạnh đó, mô hình thực nghiệm của Hao và cộng sự (2023) với sóng Ricker có tần số trung tâm 40 Hz cũng khẳng định sự dịch chuyển tần số trọng tâm về phía dải tần thấp khi truyền qua tầng đất đá có độ hấp thụ năng lượng cao.

3. Ước lượng hệ số phẩm chất hấp thụ Q

Năm 2023, Hao và cộng sự đã công bố giải pháp phân rã địa chấn sử dụng từ điển sóng Ricker để ước lượng hệ số phẩm chất địa chất $Q$ . Hệ số $Q$ là đại lượng không thứ nguyên phản ánh mức độ suy hao năng lượng của môi trường thạch học khi sóng truyền qua. Sự dịch chuyển tần số trọng tâm của các sóng phản xạ thu được từ quá trình phân rã cho phép thiết lập lát cắt thuộc tính $Q$ chi tiết, hỗ trợ nhận diện các đới nứt nẻ và biến đổi tướng đá chứa chất lưu.

4. Thực tiễn ứng dụng tại các bồn trũng trầm tích Việt Nam

Tại thềm lục địa Việt Nam, kỹ thuật phân rã địa chấn được quan tâm và triển khai trong các đề án xử lý và minh giải số liệu địa chấn tại các bồn trũng trầm tích như bồn trũng Sông Hồng, Cửu Long và Nam Côn Sơn. Việc phân tách các trường sóng địa chấn thành các thành phần tần số đơn lẻ hỗ trợ các nhà địa chất trong nước nâng cao độ tin cậy khi nhận diện các bẫy địa tầng phi cấu tạo, làm rõ ranh giới các thể cát kết dạng phân nhánh và đánh giá tiềm năng vỉa chứa trong các điều kiện địa chất phức tạp.

Hạn chế kỹ thuật và những vấn đề khoa học còn mở

Mặc dù là công cụ phân tích địa chấn mạnh mẽ, phân rã địa chấn tồn tại một số rào cản vật lý và toán học cần được cân nhắc kỹ lưỡng:

Nguyên lý bất định thời gian tần số: Tương tự như hệ thức bất định Heisenberg trong vật lý lượng tử, tích số giữa độ phân giải thời gian và độ phân giải tần số của các phép biến đổi tuyến tính luôn bị chặn dưới. Tăng cường độ phân giải ở một miền đại lượng tất yếu dẫn đến sự suy giảm độ phân giải ở miền còn lại.
Tính không duy nhất của nghiệm phân rã: Việc xấp xỉ một đường ghi địa chấn bằng từ điển sóng dư thừa là bài toán tối ưu hoá phi lồi. Các thuật toán tham lam như truy đuổi khớp có thể rơi vào các cực tiểu cục bộ, dẫn đến việc cùng một đường ghi địa chấn có thể phân rã thành các tổ hợp sóng khác nhau tuỳ thuộc vào điều kiện dừng và thiết kế từ điển.
Nhiễu địa chấn và sóng bội: Tín hiệu nhiễu ngẫu nhiên và sóng phản xạ nhiều lần có thể bị thuật toán phân rã nhầm lẫn thành các đơn vị sóng phản xạ địa chất thực tế, tạo ra các dị thường phổ nhân tạo gây sai lệch kết quả minh giải.
Tác động chồng lấn giữa thạch học và chất lưu: Cả sự biến đổi chiều dày vỉa mỏng lẫn sự hiện diện của khí tự nhiên đều có thể làm thay đổi phổ tần số của tín hiệu. Do đó, việc kết luận sự tồn tại của bẫy chứa dầu khí luôn đòi hỏi sự kiểm chứng đồng thời từ số liệu đo carota giếng khoan và mô hình hoá địa chấn tổng hợp.