Tổng quan nghiên cứu

Thành phố Hà Nội nằm trên hệ thống đới đứt gãy kiến tạo Sông Hồng, Sông Chảy và Sông Lô đang hoạt động, tiềm ẩn nguy cơ động đất cực đại có thể đạt độ lớn Ms = 6,1 cùng cường độ chấn động cấp VII đến VIII theo thang MSK-64. Theo các số liệu quan trắc địa chấn lịch sử tại Việt Nam tính đến năm 2006, khu vực lãnh thổ và lân cận đã ghi nhận 146 trận động đất có độ lớn Ms từ 4,5 trở lên, trong đó khu vực Hà Nội từng chịu tác động của các trận động đất lịch sử năm 1278 và 1285 với độ lớn Ms = 5,1 gây chấn động cấp VII-VIII. Với tốc độ đô thị hóa nhanh chóng và mật độ dân cư cao, sự hiện diện của tầng trầm tích Đệ tứ bở rời có bề dày biến đổi phức tạp từ 0 m đến trên 90 m tạo nên hiệu ứng khuếch đại sóng địa chấn cục bộ rất nguy hiểm.

Nghiên cứu được triển khai trong phạm vi tọa độ từ 20º52’N đến 21º14’N và từ 105º42’E đến 106º02’E nhằm giải quyết bài toán đánh giá đặc tính dao động vi địa chấn, kiểm chứng bản chất vật lý của tỷ số phổ H/V, và thiết lập bản đồ vi phân vùng động đất phục vụ quy hoạch xây dựng thủ đô. Mục tiêu cụ thể là xác định chu kỳ trội của nền đất, phân loại nền đất theo tiêu chuẩn xây dựng, đồng thời xây dựng mô hình thực nghiệm định lượng mối quan hệ giữa chu kỳ trội và chiều dày lớp phủ trầm tích. Kết quả nghiên cứu mang lại giá trị thực tiễn to lớn khi giúp giảm thiểu khoảng 60% chi phí khảo sát so với các phương pháp khoan thăm dò truyền thống, cung cấp cơ sở dữ liệu nền móng chuẩn xác cho việc thiết kế kháng chấn theo tiêu chuẩn quốc gia, góp phần bảo vệ an toàn tính mạng và công trình hạ tầng đô thị.

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 ứng dụng lý thuyết phân tích tỷ số phổ biên độ dao động ngang trên dao động đứng (Horizontal-to-Vertical Spectral Ratio - HVSR) do Nakamura đề xuất năm 1989 và hoàn thiện trong giai đoạn 2000-2007. Khung lý thuyết này khẳng định tỷ số phổ H/V của dao động vi địa chấn tại một trạm quan trắc phản ánh trực tiếp hàm truyền sóng ngang qua tầng phủ bở rời và xác định chính xác tần số dao động riêng của nền đất.

Song song với đó, mô hình truyền sóng 2 lớp đơn giản và lý thuyết cộng hưởng một phần tư bước sóng được ứng dụng để mô hình hóa mối quan hệ phi tuyến giữa chiều dày lớp phủ và tần số dao động trội. Bốn khái niệm học thuật cốt lõi được sử dụng xuyên suốt gồm:

  • Dao động vi địa chấn: Những dao động liên tục có biên độ dịch chuyển cực nhỏ từ 0,1 đến 1 µm và biên độ vận tốc từ 0,001 đến 0,01 cm/s do nguồn gốc tự nhiên và nhân tạo sinh ra.
  • Chu kỳ trội (T0): Chu kỳ tương ứng với đỉnh khuếch đại phổ lớn nhất trên đường cong HVSR.
  • Tính elip của sóng Rayleigh: Đặc trưng phân cực của sóng mặt chuyển động theo quỹ đạo elip ngược chiều kim đồng hồ tại bề mặt.
  • Hiệu ứng khuếch đại nền địa phương: Hiện tượng sóng địa chấn bị phản xạ nhiều lần và tích lũy năng lượng khi truyền từ đá gốc cứng qua tầng trầm tích mềm hơn.

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

Nguồn dữ liệu thực nghiệm bao gồm 93 điểm đo dao động vi địa chấn được đo đạc bằng hệ thống máy ghi phân giải cao 24-bit SAMTAC-801H kết hợp đầu thu địa chấn 6 thành phần VSE315D và VSE355EV của hãng Tokyo Sokushin (dải tần 0,1 Hz đến 50 Hz). Mẫu nghiên cứu kết hợp đối sánh trực tiếp với 42 lỗ khoan địa chất công trình, trong đó có 33 lỗ khoan khoan sâu tới tầng đá gốc Đệ tam và 24 điểm đo vi địa chấn được bố trí trùng khớp tại vị trí lỗ khoan.

Phương pháp chọn mẫu dựa trên mạng lưới không gian phủ khắp các vùng địa hình bậc thang của Hà Nội, kết hợp lấy mẫu có chủ đích dọc theo 2 tuyến mặt cắt địa chất công trình tiêu biểu từ Bắc xuống Nam và từ Tây sang Đông. Lý do lựa chọn phương pháp phân tích tỷ số phổ H/V một trạm của Nakamura là tính ưu việt vượt trội trong môi trường đô thị: phương pháp không phá hủy, không cần sử dụng nguồn nổ nhân tạo, không cần trạm đá gốc tham chiếu cố định, và triệt tiêu hiệu quả ảnh hưởng của sóng mặt nhân tạo. Quy trình phân tích tín hiệu được thực hiện qua việc cắt băng ghi 18 phút thành các cửa sổ 20,48 giây (tương đương 4096 điểm số liệu với tần số lấy mẫu 200 Hz), áp dụng biến đổi Fourier nhanh (FFT), làm mịn phổ và tính tỷ số phổ trung bình để xác định đỉnh chu kỳ trội T0.

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

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

Thứ nhất, nghiên cứu đã kiểm chứng và làm sáng tỏ bản chất vật lý của dao động vi địa chấn thông qua phân tích quỹ đạo dao động hạt. Kết quả tại điểm đo H02 chứng minh rằng xung quanh dải tần số đỉnh trội, các hạt đất dao động chủ yếu theo phương ngang tương tự sóng cắt S. Ngược lại, tại dải tần số vùng lõm của đường cong HVSR, dao động hạt chuyển dịch chủ yếu theo phương thẳng đứng đặc trưng của sóng Rayleigh.

Thứ hai, chu kỳ trội T0 của nền đất Hà Nội biến thiên trong dải từ 0,51 giây đến 1,14 giây và có xu hướng tăng dần rõ rệt từ Bắc xuống Nam và từ Tây sang Đông. Dựa trên số liệu đo đạc, lãnh thổ thành phố được phân chia thành 4 vùng chu kỳ rõ rệt: Vùng I có T0 nhỏ hơn hoặc bằng 0,6 giây; Vùng II có T0 từ trên 0,6 đến 0,8 giây; Vùng III có T0 từ trên 0,8 đến 1,2 giây; và Vùng IV có T0 lớn hơn 1,2 giây.

Thứ ba, theo Tiêu chuẩn Xây dựng Việt Nam TCXDVN 375:2006 kết hợp tiêu chuẩn phân loại nền đất của Đài Loan, khu vực nghiên cứu được xác lập thành 2 loại nền chính:

  • Nền loại C (chu kỳ trội T0 nhỏ hơn hoặc bằng 0,6 giây): Chiếm diện tích khu vực phía Bắc và Tây Bắc thành phố (thuộc Sóc Sơn, Đông Anh).
  • Nền loại D (chu kỳ trội T0 lớn hơn 0,6 giây): Bao phủ toàn bộ khu vực nội thành trung tâm, các dải đất dọc sông Hồng, khu vực phía Đông và phía Nam.

Thứ tư, nghiên cứu đã thiết lập thành công phương trình tương quan thực nghiệm phi tuyến giữa chiều dày lớp phủ Đệ tứ (D, tính bằng mét) và chu kỳ trội vi địa chấn (T0, tính bằng giây) cho khu vực Hà Nội:

D = 82 · (T0)^1,12

Độ sâu tầng phủ tính toán dao động từ dưới 40 m ở phía Bắc đến trên 90 m tại vùng lõi trung tâm đô thị, cá biệt tại khu vực Long Biên chiều dày lớp phủ đạt từ 100 m đến 145 m.

Thảo luận kết quả

Nguyên nhân căn bản của sự phân hóa chu kỳ trội T0 bắt nguồn từ đặc điểm địa chất công trình và địa mạo khu vực. Phía Bắc thành phố nằm gần khối nâng Tam Đảo, tầng trầm tích mỏng và thành phần chủ yếu là sét cứng thuộc hệ tầng Vĩnh Phúc, tạo nên chu kỳ dao động riêng ngắn (T0 dưới 0,6 giây). Ngược lại, khu vực trung tâm và phía Nam được cấu thành bởi các lớp trầm tích sông, hồ, đầm lầy thuộc hệ tầng Hải Hưng và Thái Bình có bề dày lớn (trên 70-90 m) gồm bùn, than bùn, sét bột và cát mịn bở rời, khiến sóng địa chấn bị khuếch đại mạnh và kéo dài chu kỳ dao động lên mức 0,8 đến 1,2 giây.

So sánh với các nghiên cứu quốc tế, phương trình thực nghiệm tại Hà Nội có sự tương đồng cao về quy luật hàm mũ với kết quả tại bồn trũng Lower Rhine Embayment của Đức (với hệ số a = 96 và b = -1,388) và các kết quả nghiên cứu tại bồn trũng Đài Bắc của Đài Loan. So với các công bố trong nước của Nguyễn Ngọc Thủy (chu kỳ từ 0,6 đến 2,5 giây) và Tuladhar (chu kỳ từ 0,4 đến 1,2 giây), kết quả 93 điểm đo của luận văn có độ tập trung và độ tin cậy vượt trội nhờ số lượng mẫu lớn và đối chiếu chi tiết với lỗ khoan sâu.

Dữ liệu địa chấn và địa chất có thể được trực quan hóa thông qua biểu đồ so sánh chiều dày lớp phủ tính toán DT và chiều dày thực tế từ lỗ khoan DK dọc theo 2 tuyến mặt cắt địa chất công trình. Cụ thể, trên mặt cắt 1 (từ điểm H01 đến H59), chiều dày lớp phủ tăng từ 39,5 m lên 87 m tương ứng với T0 tăng từ 0,51 giây lên 0,82 giây. Trên mặt cắt 2 (từ điểm H60 đến H67), chiều dày biến thiên từ 35,1 m đến 94 m tương ứng với T0 từ 0,51 giây đến 1,14 giây. Sự trùng khớp giữa số liệu dự báo từ vi địa chấn và số liệu khoan đạt độ chính xác trên 85%, khẳng định phương pháp vi địa chấn hoàn toàn có thể thay thế hiệu quả mạng lưới khoan dày đặc.

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

Thứ nhất, chuẩn hóa và số hóa bản đồ vi phân vùng động đất: Sở Xây dựng Hà Nội phối hợp cùng Viện Vật lý Địa cầu cần tiến hành mở rộng mạng lưới quan trắc lên 300 đến 500 điểm đo vi địa chấn tại 100% các quận mới và khu đô thị vệ tinh trong giai đoạn 2024-2026, nhằm nâng cao độ phân giải của bản đồ phân vùng rủi ro địa chấn đô thị.

Thứ hai, tích hợp bắt buộc chỉ số chu kỳ trội T0 vào quy chuẩn cấp phép xây dựng: Các cơ quan quản lý quy hoạch cần yêu cầu các đơn vị tư vấn thiết kế kết cấu áp dụng nghiêm ngặt Tiêu chuẩn TCXDVN 375:2006. Đặc biệt, đối với các công trình có chiều cao từ 8 đến 25 tầng tại khu vực nền đất loại D (nội thành và dọc sông Hồng), bắt buộc phải tính toán tránh dải chu kỳ cộng hưởng từ 0,8 đến 1,2 giây, hoàn thành áp dụng ngay trong năm 2025.

Thứ ba, thiết lập mạng lưới quan trắc dao động mạnh thời gian thực: Ủy ban Nhân dân Thành phố Hà Nội cần đầu tư lắp đặt tối thiểu 15 đến 20 trạm đo gia tốc nền liên tục tại các phân vùng nền C và D trọng điểm trong giai đoạn 2025-2027, nhằm thu thập dữ liệu gia tốc đỉnh thực tế phục vụ cảnh báo sớm và hiệu chỉnh quy chuẩn xây dựng.

Thứ tư, xây dựng cơ sở dữ liệu số địa kỹ thuật và địa vật lý tích hợp: Sở Quy hoạch - Kiến trúc Hà Nội chủ trì số hóa toàn diện dữ liệu từ hơn 1000 hố khoan địa chất hiện có, ứng dụng công thức thực nghiệm D = 82 · (T0)^1,12 để mô hình hóa không gian 3D tầng đá gốc, phục vụ trực tiếp cho công tác quy hoạch không gian ngầm và metro trong giai đoạn 2024-2025.

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

Thứ nhất, các kỹ sư kết cấu và đơn vị tư vấn thiết kế xây dựng: Luận văn cung cấp bảng thông số chi tiết về chu kỳ dao động riêng và phân loại nền đất theo từng khu vực quận huyện, giúp kỹ sư dễ dàng tra cứu để lựa chọn giải pháp móng và kết cấu chịu lực tối ưu, tránh hiện tượng cộng hưởng phá hủy công trình.

Thứ hai, các nhà quản lý đô thị và chuyên gia quy hoạch không gian: Cán bộ thuộc Sở Xây dựng, Sở Quy hoạch - Kiến trúc có thể sử dụng bản đồ vi phân vùng động đất làm căn cứ khoa học để phân bổ mật độ xây dựng, định hướng phát triển các khu đô thị mới và quy hoạch mạng lưới công trình ngầm an toàn.

Thứ ba, các nhà khoa học, nghiên cứu sinh và giảng viên chuyên ngành Địa vật lý, Địa chất công trình: Tài liệu đóng vai trò là một tài liệu tham khảo chuẩn mực về quy trình xử lý tín hiệu tỷ số phổ H/V, phương pháp phân tích dao động hạt và mô hình hóa nghịch đảo trong địa chấn công trình.

Thứ tư, các ban quản lý dự án hạ tầng giao thông và công trình trọng điểm: Cung cấp giải pháp kỹ thuật khảo sát địa tầng nhanh chóng, tiết kiệm kinh phí để đánh giá độ sâu tầng đá gốc và rủi ro địa chấn tại các tuyến đường sắt đô thị, cầu vượt sông và các tổ hợp cao ốc.

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

Phương pháp tỷ số phổ HVSR một trạm của Nakamura có ưu điểm gì vượt trội so với phương pháp khoan thăm dò truyền thống? Phương pháp HVSR một trạm sử dụng nguồn vi địa chấn tự nhiên có sẵn, hoàn toàn không phá hủy môi trường, không cần sử dụng nguồn nổ nhân tạo và không đòi hỏi trạm đá gốc tham chiếu. Kỹ thuật này giúp rút ngắn trên 70% thời gian khảo sát thực địa, tiết kiệm kinh phí đáng kể và triển khai dễ dàng tại các khu vực đô thị có mật độ dân cư và giao thông đông đúc.

Tại sao khu vực nội thành Hà Nội lại có nguy cơ chịu rủi ro chấn động địa chấn cao hơn khu vực phía Bắc? Khu vực nội thành Hà Nội sở hữu tầng trầm tích Đệ tứ bở rời rất dày (từ 70 m đến trên 90 m) gồm bùn, than bùn và sét mềm, thuộc loại nền đất D với chu kỳ trội T0 kéo dài từ 0,8 đến 1,2 giây. Lớp phủ dày này làm tăng hệ số khuếch đại biên độ sóng địa chấn khi xảy ra động đất, kết hợp với mật độ xây dựng cao tạo nên nguy cơ tổn thất lớn.

Hiện tượng cộng hưởng địa chấn giữa nền đất và công trình xây dựng tại Hà Nội dễ xảy ra với loại công trình nào? Hiện tượng cộng hưởng địa chấn nguy hiểm nhất xuất hiện khi chu kỳ dao động riêng của công trình trùng với chu kỳ trội của nền đất. Tại khu vực nội thành có T0 từ 0,8 đến 1,2 giây, các tòa nhà có chiều cao phổ biến từ 8 đến 25 tầng (mỗi tầng tương đương khoảng 0,1 giây chu kỳ dao động riêng) có nguy cơ xảy ra hiện tượng cộng hưởng rất cao, cần được tính toán kháng chấn kỹ lưỡng.

Công thức thực nghiệm D = 82 · (T0)^1,12 có thể áp dụng trực tiếp cho các địa phương khác ngoài Hà Nội không? Công thức thực nghiệm này được xây dựng dựa trên đặc tính vận tốc sóng ngang và điều kiện địa tầng trầm tích Đệ tứ đặc thù của thành phố Hà Nội. Khi áp dụng cho các bồn trũng hoặc đô thị khác như Thành phố Hồ Chí Minh hay Điện Biên, các nhà nghiên cứu cần thu thập số liệu lỗ khoan cục bộ để hiệu chỉnh lại các hệ số tương quan a và b cho phù hợp.

Dao động vi địa chấn tại dải chu kỳ trội T0 có bản chất vật lý là sóng mặt hay sóng khối? Kết quả phân tích dao động hạt và mô phỏng hàm truyền sóng trong luận văn khẳng định: Xung quanh dải tần số đỉnh trội T0, năng lượng dao động theo phương ngang chiếm ưu thế tuyệt đối và có bản chất là sóng cắt phân cực ngang SH (sóng khối). Sóng mặt Rayleigh chỉ chiếm ưu thế ở dải tần số cao hơn và tại vùng lõm của đường cong phổ.

Kết luận

  • Luận văn đã phân tích toàn diện 93 điểm đo dao động vi địa chấn và 42 lỗ khoan địa chất, kiểm chứng thành công bản chất sóng cắt SH chiếm ưu thế tại đỉnh chu kỳ trội T0 và sóng Rayleigh tại vùng lõm phổ H/V.
  • Thành lập bản đồ phân vùng chu kỳ trội nền đất Hà Nội với 4 vùng đặc trưng, chỉ ra sự biến thiên tăng dần của T0 từ Bắc xuống Nam (từ dưới 0,6 giây đến trên 1,2 giây).
  • Xác lập bản đồ phân loại nền đất đô thị gồm 2 vùng chính: Vùng nền C ở phía Bắc/Tây Bắc và Vùng nền D có độ rủi ro khuếch đại địa chấn cao tại khu vực nội thành trung tâm, phía Đông và Nam.
  • Thiết lập thành công công thức thực nghiệm đầu tiên cho Hà Nội liên hệ giữa chiều dày trầm tích và chu kỳ trội: D = 82 · (T0)^1,12 với độ chính xác cao khi kiểm chứng trên 2 tuyến mặt cắt địa chất.
  • Đóng góp giải pháp khoa học mang tính đột phá, hỗ trợ đắc lực cho công tác thiết kế kháng chấn theo tiêu chuẩn TCXDVN 375:2006 và định hướng mở rộng mô hình hóa 3D vận tốc sóng ngang trong giai đoạn 2026-2028.

Các kỹ sư xây dựng, nhà quy hoạch đô thị và cơ quan quản lý hãy tích hợp ngay dữ liệu vi phân vùng địa chấn vào thực tiễn thiết kế và cấp phép công trình để kiến tạo một thủ đô Hà Nội an toàn và phát triển bền vững trước hiểm họa thiên tai.