I. Giới thiệu về phương pháp ma trận chuyển
Phương pháp ma trận chuyển là một kỹ thuật toán học quan trọng được phát triển bởi Thomson (1950) và Haskell (1953) để nghiên cứu sóng Rayleigh trong môi trường phân lớp. Phương pháp này biểu diễn mối liên hệ giữa chuyển vị và ứng suất tại hai mặt của mỗi lớp thông qua ma trận chuyển của lớp. Việc ứng dụng phương pháp này đã được cải tiến bởi các nhà khoa học như Knopoff (1964), Dunkin (1965), Kennett (1983) và Chen (1993). Phương pháp ma trận chuyển vẫn được sử dụng rộng rãi trong địa chất học và kỹ thuật địa chất hiện nay để phân tích sóng lan truyền trong các cấu trúc phân lớp phức tạp.
1.1. Lịch sử phát triển phương pháp ma trận chuyển
Năm 1885, Rayleigh phát hiện ra sóng mặt Rayleigh lan truyền tương tự sóng mặt nước. Từ đó đến nay, nhiều nghiên cứu về sóng mặt Rayleigh trong bán không gian phân lớp nhận được sự quan tâm lớn trong lĩnh vực địa chất học. Phương pháp ma trận chuyển được Thomson đề xuất năm 1950, sau đó Haskell cải tiến năm 1953, tạo nền tảng cho các nghiên cứu hiện đại.
1.2. Nguyên lý cơ bản của phương pháp
Ma trận chuyển liên kết chuyển vị và ứng suất giữa hai mặt của lớp. Tích các ma trận này xác định mối quan hệ giữa chuyển vị ở mặt tự do lớp trên cùng và mặt trên của bán không gian. Các điều kiện biên được áp dụng để tìm phương trình tán sắc cho sóng mặt, giúp xác định tần số lan truyền của sóng.
II. Ứng dụng trong môi trường phân lớp trực hướng
Khi vật liệu cấu thành các lớp không còn là đang hướng mà là vật liệu trực hướng, bài toán trở nên phức tạp hơn. Số sóng theo phương thẳng đứng của sóng phẳng lan truyền trong môi trường phân lớp trực hướng trở thành giá trị phức. Điều này dẫn đến ma trận chuyển có thể nhận giá trị phức, gây khó khăn trong tính toán số. Dạng hiển của ma trận chuyển không còn đơn giản như trong trường hợp đang hướng. Phương pháp Thomson-Haskell cổ điển cần được cải tiến để xử lý các điều kiện biên phức tạp của sóng mặt trong bán không gian có tính chất đặc biệt.
2.1. Khó khăn khi tính toán với vật liệu trực hướng
Ma trận chuyển của lớp trực hướng có các phần tử nhận giá trị phức, làm phức tạp quá trình tính toán số. Dạng hiển của ma trận không còn rõ ràng như trong vật liệu đang hướng. Các công trình của Crampin (1970) và Crampin-Taylor (1971) đã chỉ ra những thách thức này trong phân tích sóng P-SV.
2.2. Biểu diễn ma trận matrizant
Ma trận matrizant là ma trận nghiệm cơ bản của hệ phương trình vi phân chuyển động của sóng P-SV. Cách tiếp cận mới sử dụng biểu diễn hiển của ma trận matrizant dưới dạng cụ thể, được Solyanik (1977) và Rokhlin-Wang (1992) đề xuất, giúp tính toán số thuận tiện hơn.
III. Phương trình tán sắc và tỷ số H V
Phương trình tán sắc được thiết lập bằng cách áp dụng điều kiện biên tại mặt tự do của lớp trên cùng và điều kiện liên tục tại các mặt tiếp xúc giữa các lớp. Tỷ số H/V (Horizontal to Vertical) là tỷ số giữa thành phần ngang và thành phần thẳng đứng của sóng, được sử dụng rộng rãi trong địa chất kỹ thuật để xác định tính chất của các lớp địa chất. Ảnh hưởng của tính chất không đẳng hướng của vật liệu lên tỷ số H/V là một yếu tố quan trọng cần được nghiên cứu. Phương pháp ma trận chuyển cải tiến cho phép tính toán chính xác tỷ số H/V trong các môi trường phân lớp phức tạp với vật liệu trực hướng.
3.1. Xây dựng phương trình tán sắc
Phương trình tán sắc được xác định từ tích các ma trận chuyển của từng lớp, áp dụng điều kiện ứng suất tiếp bằng không tại mặt tự do. Determinant của ma trận tổng hợp bằng không cho ra các giá trị tần số tương ứng với các chế độ sóng mặt Rayleigh.
3.2. Ảnh hưởng của tính trực hướng lên tỷ số H V
Vật liệu trực hướng làm thay đổi mối quan hệ giữa các thành phần chuyển vị trong sóng mặt, ảnh hưởng trực tiếp đến tỷ số H/V. Tính chất không đẳng hướng tạo ra các chế độ sóng phức tạp, cần được tính toán số chính xác bằng phương pháp ma trận cải tiến.
IV. Kết quả và ứng dụng thực tiễn
Phương pháp ma trận chuyển được cải tiến cho môi trường phân lớp trực hướng đã cho phép giải quyết nhiều bài toán phức tạp trong địa chất học ứng dụng. Các kết quả minh họa số cho thấy tác động rõ rệt của tính chất không đẳng hướng lên phương trình tán sắc và tỷ số H/V. Phương pháp này có ứng dụng quan trọng trong đánh giá nguy hiểm động đất, phân loại vật liệu địa chất và dự báo phản ứng của đất trong công trình xây dựng. Việc sử dụng ma trận matrizant giúp cải thiện độ ổn định số và chính xác của tính toán so với phương pháp cổ điển, đáp ứng nhu cầu phân tích các mô hình địa chất thực tế phức tạp.
4.1. Kết quả minh họa số
Các mô phỏng số cho thấy ảnh hưởng của tính trực hướng lên phương trình tán sắc là đáng kể, đặc biệt ở các tần số cao. Tỷ số H/V biến đổi theo mức độ trực hướng của vật liệu, cung cấp thông tin quý giá cho việc xác định cấu trúc địa chất.
4.2. Ứng dụng trong địa chất kỹ thuật
Phương pháp này được ứng dụng trong phân tích sóng địa chấn, đánh giá nguy hiểm động đất, và thiết kế công trình xây dựng chống chấn động. Ma trận chuyển cải tiến cho phép dự báo chính xác phản ứng của từng lớp địa chất trước các tải động.
