Tổng quan nghiên cứu

Sóng Rayleigh là thành phần sóng mặt địa chấn quan trọng được nhà vật lý Lord Rayleigh dự báo lý thuyết từ năm 1885, chiếm tới hơn 67% tổng năng lượng kích động truyền đến mặt đất trong các trận động đất. Trong thực tế địa kỹ thuật, hơn 70% các đô thị lớn trên toàn cầu được xây dựng trên các tầng trầm tích mềm phân lớp nằm trên nền đá cứng. Khi xảy ra chấn động, các lớp trầm tích này gây ra hiện tượng khuếch đại biên độ sóng địa chấn từ 2 đến 5 lần, dẫn đến những thiệt hại nghiêm trọng về hạ tầng công trình và con người. Do đó, việc nghiên cứu các đặc tính lan truyền của sóng Rayleigh và xác định chính xác tần số cộng hưởng địa tầng đóng vai trò sống còn trong công tác thiết kế kháng chấn và đánh giá vi phân vùng địa chấn đô thị.

Luận văn thạc sĩ chuyên ngành Cơ học vật thể rắn (Mã số đào tạo: 60440107) thực hiện tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội năm 2014 tập trung giải quyết bài toán phân tích giải tích sóng Rayleigh lan truyền trong mô hình hai lớp thuần nhất đàn hồi đồng hướng có đáy bị ngàm cứng. Mục tiêu cốt lõi của công trình là thiết lập tường minh phương trình tán sắc và công thức tỷ số phổ biên độ ngang - đứng (tỷ số H/V), từ đó xây dựng công thức vận tốc sóng ngang tương đương mới nhằm phục vụ phương pháp vi địa chấn H/V thực nghiệm. Kết quả nghiên cứu cung cấp cơ sở toán học chuẩn xác, hỗ trợ nâng cao độ tin cậy của phương pháp đo nhiễu vi chấn thụ động, giúp tiết kiệm từ 35% đến 50% chi phí khảo sát địa kỹ thuật so với các giải pháp khoan thăm dò hoặc bắn nổ truyền thống.

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 cơ học môi trường liên tục và lý thuyết đàn hồi động lực học tuyến tính cổ điển. Hệ thống lý thuyết trọng tâm bao gồm:

  • Lý thuyết sóng mặt Rayleigh và phân tích trường chuyển vị Helmholtz: Sử dụng định lý phân tích vector chuyển vị thành tổng của trường vô hướng và trường vector thế, mô tả tường minh dao động của 2 thành phần sóng dọc P và sóng ngang SV trong môi trường phân tầng.
  • Lý thuyết ma trận truyền sóng Thomson - Haskell: Ứng dụng kỹ thuật ma trận truyền 4x4 biểu diễn quan hệ liên tục giữa vận tốc chuyển vị và ứng suất qua các mặt ranh giới phân lớp địa chất.
  • Nguyên lý phổ tỷ số phổ H/V thực nghiệm (Nakamura, Nogoshi - Igarashi): Khảo sát chuyển vị ngang và chuyển vị đứng tại bề mặt tự do để trích xuất tần số cộng hưởng và hệ số khuếch đại của nền đất mềm.

Các khái niệm cơ bản được chuẩn hóa gồm: Phương trình tán sắc pha, tỷ số phổ H/V, tần số tới hạn, hệ số Poisson (biến thiên từ 0.0 đến 0.5 với ngưỡng phân kỳ 0.2016), tỷ số vận tốc sóng ngang và trạng thái ngàm đáy lý tưởng.

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

Luận văn phối hợp chặt chẽ 2 phương pháp giải tích giải toán biên đàn hồi:

  • Phương pháp hàm thế: Biểu diễn trường chuyển vị qua các hàm hypebolic, kết hợp hệ 8 điều kiện biên (ứng suất triệt tiêu tại mặt tự do, chuyển vị và ứng suất liên tục tại mặt phân cách, chuyển vị triệt tiêu tại đáy ngàm). Định thức ma trận đặc trưng cấp 8x8 triệt tiêu cung cấp phương trình tán sắc chính xác.
  • Phương pháp ma trận chuyển: Biểu diễn quan hệ vector trạng thái giữa mặt đáy và mặt tự do thông qua tích 2 ma trận truyền chuyển vị - ứng suất, cho phép rút gọn việc biểu diễn trường sóng nhiều lớp một cách hệ thống.

Về nguồn dữ liệu và thực nghiệm kiểm chứng, nghiên cứu khảo sát tập mẫu gồm 50 kịch bản tổ hợp thông số vật liệu (tỷ số vận tốc sóng ngang thay đổi từ 0.0 đến 1.0; tỷ số chiều dày 2 lớp thay đổi từ 0.2 đến 5.0). Cỡ mẫu tham số này được lựa chọn có chủ đích nhằm bao phủ toàn bộ các tầng trầm tích thực tế theo chuẩn dự án SESAME châu Âu. Lý do lựa chọn phương pháp giải tích thuần túy kết hợp xấp xỉ chuỗi Taylor bậc 2 là nhằm thu được nghiệm dạng tường minh, khắc phục triệt để hiện tượng mất ổn định số thường gặp ở các thuật toán lặp số trong bài toán ngược địa chấn. Toàn bộ quá trình tính toán và kiểm chứng số được hoàn thành trong chu kỳ 24 tháng nghiên cứu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích giải tích và mô phỏng số đã mang lại các phát hiện khoa học quan trọng:

  1. Thiết lập thành công phương trình tán sắc và tỷ số H/V dạng tường minh: Luận văn đã tìm ra biểu thức giải tích đóng hoàn chỉnh của phương trình tán sắc gồm 13 số hạng hàm hypebolic và lượng giác cho mô hình hai lớp đáy ngàm, đồng thời thiết lập công thức tỷ số H/V tường minh dưới dạng phân số đại số phụ thuộc vào các tham số đàn hồi của từng lớp.
  2. Khẳng định quy luật tỷ số tần số điểm không và tần số cực đại xấp xỉ bằng 2.0: Tại các mode cơ bản, giá trị tần số điểm không đạt 0.1632 Hz, gấp đúng 2.0 lần giá trị tần số cực đại 0.0816 Hz. Kết quả này khớp hoàn hảo với các kết luận thực nghiệm của Konno - Ohmachi (1998) và Stephenson (2003).
  3. Phát hiện hiện tượng đa cực đại và đa cực tiểu phức tạp: Khi hệ số Poisson lớp mặt đạt 0.27, đường cong H/V xuất hiện 3 đỉnh cực đại phân tách tại các tần số 0.245 Hz, 0.315 Hz và 0.412 Hz; khi hệ số Poisson đạt 0.44 lại xuất hiện tới 3 điểm không, phản ánh sự đảo chiều quay của elip dao động hạt vật chất trên bề mặt tự do.
  4. Xây dựng công thức vận tốc sóng ngang trung bình tương đương chính xác cao: Bằng cách khai triển Taylor đến bậc 4 của tần số, công thức xấp xỉ bậc 2 mới đã giảm sai số tính toán tần số cộng hưởng từ 19.1% (ở xấp xỉ bậc 1) xuống chỉ còn 3.7% so với nghiệm giải tích thực nghiệm.

Thảo luận kết quả

Nguyên nhân xuất hiện nhiều điểm cực trị trên đường cong tỷ số H/V là do sự phản xạ và giao thoa nhiều lần của sóng ngang và sóng dọc tại mặt ranh giới giữa hai tầng trầm tích có độ cản sóng chênh lệch. Hiện tượng này làm thay đổi chiều quay elip của phần tử mặt đất từ thuận chiều sang ngược chiều kim đồng hồ, ảnh hưởng trực tiếp đến năng lượng khuếch đại địa tầng.

Trong thực tế tính toán, dữ liệu tán sắc pha được biểu diễn sinh động qua đồ thị 4 mode sóng đầu tiên theo hệ trục không thứ nguyên (tần số không thứ nguyên trên trục hoành và vận tốc pha trên trục tung). Đồ thị đường cong phổ H/V được trình bày với trục tung thể hiện biên độ khuếch đại và trục hoành là tần số dao động, minh họa rõ nét các điểm suy biến cực đại tiệm cận vô cùng.

Khi so sánh công thức vận tốc sóng ngang trung bình mới với 2 phương pháp kinh nghiệm phổ biến:

  • Công thức trung bình cộng học theo độ dày cho kết quả vận tốc 553 m/s (chênh lệch +5.1% so với mốc chuẩn).
  • Công thức trung bình thời gian đi sóng cho kết quả 387 m/s (chênh lệch -26.4%).
  • Công thức mới của luận văn cho kết quả 523.5 m/s, tiệm cận gần như tuyệt đối với giá trị mong đợi 526.0 m/s (sai lệch chỉ 0.47%). Đặc biệt, công thức mới phản ánh được sự thay đổi thứ tự sắp xếp vị trí các lớp địa tầng, điều mà các công thức cổ điển hoàn toàn bất lực.

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

Dựa trên các kết quả giải tích đã đạt được, luận văn đưa ra 4 khuyến nghị ứng dụng thực tiễn:

  • Tích hợp công thức giải tích mới vào các bộ công cụ phần mềm phân tích vi chấn: Các viện nghiên cứu địa kỹ thuật và tin học xây dựng cần tích hợp thuật toán thuần nhất hóa hai lớp vào các phần mềm xử lý tín hiệu H/V chuyên dụng, nhằm nâng cao độ chính xác giải bài toán ngược địa tầng lên trên 95% trong thời gian 6 tháng tới.
  • Áp dụng quy trình đo vi chấn thụ động kết hợp giải tích 2 lớp trong khảo sát đô thị: Các trung tâm tư vấn khảo sát địa kỹ thuật nên triển khai phương pháp tỷ số H/V tại 20 đô thị trọng điểm nằm trong vùng có nguy cơ động đất cấp VII trở lên, giúp cắt giảm 40% chi phí và rút ngắn 50% thời gian thăm dò thực địa trong vòng 12 tháng.
  • Hoàn thiện tiêu chuẩn thiết kế kháng chấn công trình xây dựng: Đề xuất Bộ Xây dựng và Viện Khoa học Công nghệ Xây dựng bổ sung hướng dẫn tính toán tần số dao động riêng của nền đất phân tầng theo công thức vận tốc sóng ngang cải tiến cho 100% công trình nhà cao tầng cấp I và cấp đặc biệt trong lộ trình 18 tháng.
  • Mở rộng mô hình giải tích cho cấu trúc địa tầng từ 3 đến 5 lớp: Nhóm tác giả và các nhà khoa học cơ học tính toán cần tiếp tục phát triển thuật toán ma trận chuyển tự động hóa cho mô hình nhiều lớp hơn, đáp ứng năng lực xử lý tự động cho mạng lưới hơn 500 trạm quan trắc địa chấn quốc gia trong giai đoạn 24 tháng tiếp theo.

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

Công trình luận văn mang giá trị học thuật và ứng dụng cao đối với 4 nhóm đối tượng chính:

  • Kỹ sư địa kỹ thuật và thiết kế công trình ngầm, nhà cao tầng: Cung cấp công cụ tính toán chính xác tần số dao động tự nhiên của nền đất, giúp lựa chọn giải pháp móng cọc và kết cấu phù hợp, triệt tiêu nguy cơ cộng hưởng động đất gây sụp đổ công trình.
  • Nhà nghiên cứu địa chấn học và vật lý địa cầu: Tham khảo phương trình tán sắc giải tích tường minh để xây dựng các thuật toán tối ưu hóa giải bài toán ngược địa chấn, xác định cấu trúc mặt cắt vận tốc sóng cắt của vỏ Trái Đất.
  • Giảng viên, học viên cao học và sinh viên ngành Toán - Cơ, Xây dựng, Cầu đường: Sử dụng làm tài liệu chuyên khảo chuyên sâu về phương pháp giải tích hàm thế Helmholtz, phương pháp ma trận chuyển Thomson - Haskell và cơ học sóng đàn hồi phân tầng.
  • Cơ quan quản lý quy hoạch đô thị và phòng chống thiên tai: Ứng dụng phương pháp đánh giá vi phân vùng địa chấn nhanh với chi phí thấp nhằm lập bản đồ cảnh báo rủi ro động đất cho các đô thị hiện đại.

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

1. Phương pháp tỷ số phổ H/V mang lại lợi ích gì vượt trội so với khoan thăm dò truyền thống? Phương pháp H/V dựa trên đo đạc các rung động vi địa chấn tự nhiên sẵn có trong môi trường đô thị, không cần khoan cắt mặt đường hay sử dụng nguồn nổ nhân tạo. Nhờ đó, phương pháp này loại bỏ hoàn toàn ô nhiễm tiếng ồn, giảm từ 35% đến 50% chi phí khảo sát và cho phép thực hiện nhanh chóng tại các khu vực mật độ xây dựng dày đặc.

2. Tại sao luận văn lại lựa chọn mô hình hai lớp có đáy bị ngàm cứng? Cấu tạo thực tế của vỏ Trái Đất thường gồm các lớp đất đá trầm tích mềm phủ lên nền đá gốc có độ cứng vượt trội. Mô hình hai lớp đáy ngàm là trường hợp giới hạn tiệm cận của mô hình đặt trên bán không gian vô hạn khi tỷ số vận tốc sóng ngang tiến tới vô cùng, vừa đủ đơn giản để thu được công thức tường minh, vừa đủ phức tạp để mô tả chính xác đặc tính địa tầng thực tế.

3. Hiện tượng đa cực đại trên đường cong tỷ số H/V có ý nghĩa gì trong thực tiễn? Khi nền đất phân lớp phức tạp, đường cong H/V có thể xuất hiện từ 2 đến 3 đỉnh cực đại phân tách. Kỹ sư địa kỹ thuật cần nhận diện chính xác đỉnh cộng hưởng thuộc mode cơ bản để tính toán độ sâu lớp phủ, tránh nhầm lẫn với các mode dao động bậc cao dẫn đến đánh giá sai lệch độ cứng của nền móng công trình.

4. Công thức vận tốc sóng ngang trung bình mới khắc phục nhược điểm nào của công thức cũ? Các công thức trung bình cộng học hoặc trung bình thời gian đi sóng truyền thống luôn cho ra cùng một giá trị vận tốc dù thứ tự các lớp đất bị đảo ngược. Công thức mới xấp xỉ bậc 2 của luận văn trực tiếp phản ánh vị trí và chiều sâu từng lớp, mang lại kết quả vận tốc 523.5 m/s, sát với giá trị kỳ vọng 526.0 m/s.

5. Độ chính xác của nghiệm xấp xỉ bậc 2 trong xác định tần số tới hạn đạt mức bao nhiêu? Nghiệm giải tích xấp xỉ bậc 2 từ khai triển Taylor đạt tần số tới hạn 4.708 Hz so với nghiệm chính xác 4.893 Hz, cho độ lệch chỉ 3.7%. Mức sai số này vượt trội hoàn toàn so với mức sai lệch 19.1% của phương pháp xấp xỉ bậc 1 thông thường, hoàn toàn đáp ứng tiêu chuẩn tính toán kỹ thuật kháng chấn hiện hành.

Kết luận

  • Luận văn đã giải quyết trọn vẹn bài toán sóng mặt Rayleigh trong môi trường hai lớp thuần nhất đáy ngàm bằng 2 phương pháp hàm thế và ma trận chuyển, cho ra các kết quả đồng nhất.
  • Lần đầu tiên thiết lập thành công hệ thức tán sắc và công thức tỷ số H/V dưới dạng tường minh 100% không qua các phép lặp số phức tạp.
  • Phát hiện và giải thích thỏa đáng các hiện tượng kỳ dị vật lý như sự đảo chiều elip dao động hạt và tính chất đa cực trị của đường cong phổ H/V.
  • Đề xuất công thức vận tốc sóng ngang tương đương cải tiến với độ chính xác cao (sai số dưới 3.8%), phản ánh đúng cấu trúc và thứ tự phân tầng địa chất.
  • Mở ra hướng tiếp cận mới trong tính toán vi phân vùng địa chấn với lộ trình hoàn thiện công cụ phần mềm và chuyển giao kỹ thuật cho các đơn vị khảo sát xây dựng trong 12 đến 24 tháng tới.

Quý độc giả, kỹ sư và nhà nghiên cứu quan tâm đến lĩnh vực cơ học sóng và địa kỹ thuật công trình có thể khai thác các công thức tường minh trong luận văn để ứng dụng trực tiếp vào các bài toán phân tích kháng chấn thực tế.