Tổng quan luận án

Đề tài luận án tiến sĩ: "Nghiên cứu ước lượng hiệu ứng nền đất khu vực nội thành thành phố Hà Nội và lân cận"
Tác giả: Nguyễn Tiến Hùng
Chuyên ngành: Vật lý địa cầu (Mã số: 9 44 01 11)
Người hướng dẫn khoa học: PGS. Nguyễn Hồng Phương và TS. Nguyễn Lê Minh
Cơ sở đào tạo: Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam (Hà Nội, 2023).

Tính cấp thiết và khoảng trống nghiên cứu

Thành phố Hà Nội là đô thị loại đặc biệt, trung tâm chính trị – hành chính quốc gia với quy mô dân số xấp xỉ 8,5 triệu người. Khu vực nội thành tập trung mật độ dân cư cao cùng nhiều công trình hạ tầng trọng điểm, nhà cao tầng và siêu cao tầng. Về bối cảnh địa chất và địa chấn, Hà Nội nằm trên đới đứt gãy Sông Hồng dài hơn 1000 km, gồm ba nhánh chính: đứt gãy Sông Hồng, đứt gãy Sông Chảy và đứt gãy Sông Lô. Trong đó, đứt gãy Sông Chảy cắt ngang trung tâm thành phố có tốc độ dịch chuyển ngang trung bình hiện nay đạt 2,7 ± 1,6 mm/năm, tiềm ẩn độ nguy hiểm địa chấn với cường độ chấn động cực đại có thể đạt cấp VIII theo thang MSK-64. Mặc dù lịch sử ghi nhận các trận động đất mạnh tại chỗ vào các năm 1277, 1278 và 1285 (cấp VII–VIII), trong giai đoạn hiện đại Hà Nội chủ yếu chịu tác động rung lắc từ các trận động đất xa. Cảm nhận rung chấn gia tăng mạnh tại các công trình cao tầng do cấu trúc nền đất mềm bồi tích Đệ tứ dày hàng chục đến hàng trăm mét và sự khuếch đại dao động trong môi trường các lớp trầm tích nằm ngang song song.

Khoảng trống nghiên cứu được xác định: Các nghiên cứu trước đây tại Việt Nam và Hà Nội chủ yếu tập trung vào việc xác định gia tốc đỉnh (PGA) trên nền đá cứng theo chu kỳ lặp lại, khảo sát tần số/chu kỳ dao động riêng lẻ hoặc mô phỏng phản ứng nền cho một số công trình đơn lẻ theo kịch bản cực đoan. Chưa có công trình nào đánh giá toàn diện, định lượng về mức độ khuếch đại dao động nền (KĐDĐN) do điều kiện nền đất địa phương gây ra, cũng như chưa xây dựng được sơ đồ phân loại nền đất chi tiết theo giá trị vận tốc sóng ngang trung bình trong 30 m trên cùng ($V_{S30}$) theo Tiêu chuẩn quốc gia TCVN 9386:2012 cho toàn bộ khu vực nội thành và các huyện lân cận.

Mục tiêu và nhiệm vụ nghiên cứu

  • Mục tiêu nghiên cứu: Làm rõ mức độ khuếch đại sóng động đất tại mặt đất (bề mặt lớp phủ) khu vực nội thành thành phố Hà Nội và lân cận, chủ yếu do điều kiện nền đất địa phương, nhằm phục vụ công tác lập quy hoạch và thiết kế công trình chịu động đất, góp phần giảm thiểu rủi ro cho cộng đồng đô thị.
  • Nhiệm vụ nghiên cứu (gồm 4 nhiệm vụ cụ thể):
    1. Thu thập, cập nhật và tổng hợp số liệu, tài liệu đo dao động vi địa chấn (DĐVĐC) và địa chất công trình (ĐCCT);
    2. Đo bổ sung số liệu DĐVĐC một trạm theo các tuyến mặt cắt ĐCCT, dải dọc sông Hồng – sông Đuống, các quận nội thành còn lại và vùng phía Nam;
    3. Nghiên cứu, hoàn thiện quy trình đánh giá hiệu ứng nền đất (HƯNĐ) bằng kỹ thuật đo DĐVĐC;
    4. Xác định hệ số khuếch đại dao động các loại nền khu vực nội thành thành phố Hà Nội và lân cận.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Nền đất địa phương khu vực nội thành thành phố Hà Nội và lân cận, cùng mối liên quan giữa cấu trúc nền đất với sự khuếch đại sóng động đất.
  • Phạm vi không gian: Tổng diện tích khảo sát khoảng 1000 $\text{km}^2$, bao gồm toàn bộ 12 quận nội thành (Ba Đình, Bắc Từ Liêm, Cầu Giấy, Đống Đa, Hai Bà Trưng, Hà Đông, Hoàn Kiếm, Hoàng Mai, Long Biên, Nam Từ Liêm, Tây Hồ, Thanh Xuân) và 4 huyện ngoại thành lân cận (Đông Anh, Gia Lâm, Sóc Sơn, Thanh Trì).
  • Phạm vi thời gian và nguồn dữ liệu: Kế thừa các đợt khảo sát từ năm 2009, 2012, 2013 và hoàn thiện tính toán, đo đạc bổ sung đến năm 2023.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Nghiên cứu về HƯNĐ bắt đầu từ quan sát của Milne (1898) về sự chênh lệch biên độ dao động giữa hai trạm đo đặt trên hai loại nền đất khác nhau ở cùng cự ly chấn tâm. Các công trình tiếp theo của Gutenberg (1957) và Seed (1968, 1969) chứng minh sóng địa chấn bị khuếch đại mạnh khi truyền qua các lớp phù sa bở rời. Hiện tượng này được ghi nhận rõ nét trong trận động đất Michoacan 1985 tại Mê-hi-cô (Anderson và nnk, 1986; Celebi và nnk, 1987), khi bồn trầm tích làm biên độ sóng tăng từ 8 đến 50 lần so với vùng đá cứng, và trận động đất Chi-Chi 1999 tại bồn trũng Đài Bắc (Wen và nnk, 1999).

Về phân loại nền đất:

  • Sandra và Ralph (1989) nghiên cứu sự biến đổi biên độ gia tốc theo độ sâu, chỉ ra PGA tại mặt lớp phủ lớn gấp 5 lần trong đá cứng.
  • Anderson (1996) đề xuất dùng vận tốc sóng ngang ($V_S$) các lớp bề mặt để phân chia nền đất ($S_1, S_2, S_3, S_4$).
  • Dobry (2000) chuẩn hóa bảng phân loại nền đất dựa trên vận tốc sóng ngang trung bình 30 m trên cùng ($V_{S30}$) gồm 5 loại (A, B, C, D, E), trở thành nền tảng cho tiêu chuẩn NEHRP và Eurocode 8.

Về phương pháp xác định $V_S$:

  • Phương pháp trực tiếp trong lỗ khoan: Down-hole (Kanai, 1966; Woods, 1978), Up-hole (Gassmann, 1951; McDonal và nnk, 1958), Cross-hole (Stokoe và Woods, 1972; Hoar và Stokoe, 1978), carota PS-logging (Kitsunezaki, 1980; Nigbor và Imai, 1994).
  • Phương pháp tương quan thực nghiệm giữa $V_S$ và số nhát đập xuyên tiêu chuẩn $N_{SPT}$: Ohsaki và Iwasaki (1973), Imai (1977), Ohta và Goto (1978), Imai và Tonouchi (1982), Lee (1990), Yokota và nnk (1991), Athanasopoulos (1995), Jafari và nnk (2002), Hasancebi và Ulusay (2007), Hanumanthrao và Ramana (2008), Uma và nnk (2010), Tsiambaos và Sabatakakis (2011), Anbazhagan và nnk (2012), Fauzi và nnk (2014), Kirar và nnk (2016), Shukla và Solanki (2020).
  • Phương pháp sóng mặt và vi địa chấn: Đo mảng DĐVĐC bằng kỹ thuật F-K (Capon, 1969) và SPAC (Aki, 1957); phân tích tỷ số phổ H/V (HVSR) trên một trạm đo (Nakamura, 1989); đánh giá tương quan tần số trội $F_0$ và chiều dày lớp phủ $D$ (Seht và Wohlenberg, 1999); nghịch đảo cấu trúc vận tốc bằng thuật toán di truyền GA (Holland, 1975).

Các nghiên cứu tại Việt Nam và thành phố Hà Nội

  • Quy mô quốc gia: Nguyễn Đình Xuyên (1976, 2004) và Cao Đình Triều (2006) ghi nhận biến dạng đường đẳng chấn qua các trận động đất Điện Biên 1935, Vĩnh Phúc 1958, Bắc Giang 1961, Tuần Giáo 1983 và Điện Biên 2001 ($M_S = 5,3$). Các bản đồ gia tốc nền toàn quốc được xây dựng cho nền đá cứng chuẩn ($V_{S30} = 800\text{ m/s}$ hoặc $1100\text{ m/s}$) của Nguyễn Hồng Phương (2001, 2004, 2012).
  • Khu vực Hà Nội:
    • Giai đoạn 1964–1978: Phòng Vật lý địa cầu lập bản đồ phân vùng nhỏ động đất tỷ lệ 1:50.000 bằng phương pháp độ cứng địa chấn; hiệu chỉnh sau trận động đất Côn Minh (1970).
    • Năm 1990: Lê Tử Sơn và nnk đo 50 điểm DĐVĐC bằng máy $\Pi 06\text{-}12M$ và đầu đo CM-3, phân loại vùng đất Pleistocen ứng với cấp VII, Holocen ứng với cấp VIII (MSK-64).
    • Năm 1992–1994: Viện Kỹ thuật xây dựng Hà Nội và Viện Vật lý địa cầu sử dụng chương trình SHAKE tính toán thông số dao động cho 15 loại đất.
    • Năm 1997: Nguyễn Ngọc Thủy thành lập bản đồ phân vùng nhỏ động đất Hà Nội tỷ lệ 1:25.000 chia khu vực thành 3 vùng cấp chấn động VII, VIII và VIII–IX.
    • Năm 2003–2004: Nguyễn Huy Phương và nnk mở rộng khảo sát tại Sóc Sơn và bãi sông Hồng bằng máy ghi K2, đầu đo SS-1, xác định 58 giá trị chu kỳ trội $T_0$.
    • Năm 2004: Tuladhar và nnk đo 63 điểm vi địa chấn một trạm, xác định $T_0$ vùng đồng bằng từ 0,4 đến 1,2 s ($F_0 = 0,8\text{--}2,5\text{ Hz}$).
    • Năm 2012–2014: Hợp tác Nghị định thư Việt Nam – Đài Loan bổ sung 244 điểm đo máy Samtac, lập bản đồ phân vùng $V_{S30}$, $F_0$ và rủi ro động đất quận Đống Đa, Thanh Xuân.
    • Năm 2015: Sở Xây dựng Hà Nội và Viện Vật lý địa cầu khảo sát 162 điểm đo lặp tại quận Hà Đông và 6 tuyến đo khúc xạ sóng ngang bằng máy GEODE-24 kênh.
    • Năm 2017–2022: Bùi Thị Nhung (2017) đánh giá tiềm năng hóa lỏng nền (tiêu chuẩn NEHRP 1997); Phạm Thế Truyền (2020) đánh giá rủi ro động đất cho 5 quận trung tâm; Tran Ngoc Long và nnk (2021), Giang Kien Trung và nnk (2022) dùng phần mềm DEEPSOIL mô phỏng phổ phản ứng nền tại quảng trường Ba Đình và 18 lỗ khoan.

Khoảng trống luận án giải quyết

Luận án giải quyết hạn chế của việc phân loại nền đất dựa trên định tính thạch học hoặc số liệu thưa thớt, bằng cách kết hợp đo đạc DĐVĐC một trạm mật độ dày (834 điểm), kỹ thuật tỷ số phổ H/V, nghịch đảo sóng cắt SH bằng thuật toán di truyền (GA) để xác định đồng bộ 816 cấu trúc $V_{S1D}$, tính toán $V_{S30}$ và định lượng hệ số khuếch đại dao động nền ($K$) theo các kịch bản động đất cho toàn bộ khu vực nội thành Hà Nội và 4 huyện lân cận.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khái niệm và phương pháp đo dao động vi địa chấn

DĐVĐC là trường rung động tự nhiên liên tục của mặt đất có biên độ dao động từ $10^{-7}\text{ m}$ đến $10^{-5}\text{ m}$, sinh ra từ các nguồn tự nhiên (gió, sóng biển, dòng chảy) và nhân tạo (giao thông, máy móc công nghiệp).

  • Kỹ thuật đo DĐVĐC điểm tựa: Sử dụng 2 máy đo đồng thời (1 máy cố định trên đá gốc, 1 máy di động tại điểm khảo sát), xác định chu kỳ trội $T_0$ qua tỷ số biên độ dao động. Kỹ thuật này chịu ảnh hưởng lớn bởi nhiễu địa phương.
  • Kỹ thuật đo DĐVĐC một trạm (Nakamura, 1989): Đặt 1 máy đo 3 thành phần (NS, EW, Z) tại điểm khảo sát. Tỷ số phổ H/V ($HVSR$) được tính bằng tỷ số giữa phổ Fourier tổng hợp hai thành phần ngang và phổ Fourier thành phần đứng:

$$HVSR = \frac{\sqrt{F_{NS}^2 + F_{EW}^2}}{\sqrt{2} \cdot F_Z}$$

Trong đó $F_{NS}, F_{EW}, F_Z$ lần lượt là phổ Fourier của các thành phần dao động Bắc–Nam, Đông–Tây và thẳng đứng.

  • Kỹ thuật đo mảng DĐVĐC: Bố trí các máy đo thành hình tam giác đều nội tiếp các vòng tròn đồng tâm bán kính $R_1, R_2$ ($R_2 = 2R_1$). Sử dụng phương pháp phổ mật độ năng lượng tần số – số sóng F-K (Capon, 1969) để xác định vận tốc pha của sóng mặt:

$$V(f) = \frac{2\pi f}{\sqrt{k_{x0}^2 + k_{y0}^2}}$$

Mô hình tần số trội trong môi trường phân lớp

Dựa trên hiện tượng sóng dừng tạo bởi giao thoa giữa sóng tới và sóng phản xạ từ mặt ranh giới phân lớp, tần số trội của sóng dao động trong môi trường lớp phủ được xác định theo công thức:

$$f = \frac{(2k + 1)V}{4h} \quad (k = 0, 1, 2, \dots)$$

Khi $k = 0$, tần số cơ bản là:

$$F_0 = \frac{V_S}{4D}$$

Mối quan hệ thực nghiệm giữa tần số trội $F_0$ và chiều dày lớp phủ nông $D$ (Seht và Wohlenberg, 1999) có dạng:

$$D = a \cdot F_0^b$$

Trong đó các hệ số $a$ và $b$ phụ thuộc vào vận tốc sóng ngang tại mặt đất ($V_0$) và hệ số biến thiên vận tốc theo độ sâu ($x$).

Kỹ thuật chuyển đổi sóng SH và thuật toán di truyền (GA)

  • Mô hình chuyển đổi sóng SH (Haskell, 1953, 1960): Mô hình hóa môi trường gồm $n$ lớp song song, đẳng hướng và thuần nhất. Hệ phương trình lan truyền sóng cắt ngang SH qua các lớp được thiết lập dưới dạng ma trận chuyển tiếp từ mặt đá cứng lên bề mặt tự do.
  • Thuật toán di truyền (GA): Dùng để giải bài toán nghịch đảo tìm mô hình $V_{S1D}$ tối ưu từ đường cong HVSR thực nghiệm. Hàm thích nghi (fitness function) kết hợp 80% mức độ tương quan dạng phổ và 20% mức độ tương quan tần số đỉnh trội:

$$Fitness = 0{,}8 \cdot r + 0{,}2 \cdot \left(1 - \frac{|F_{SH} - F_{HV}|}{F_{HV}}\right)$$

Trong đó $r$ là hệ số tương quan tuyến tính giữa biểu đồ phổ sóng SH lý thuyết và biểu đồ HVSR thực nghiệm; $F_{SH}$ và $F_{HV}$ là các tần số trội tương ứng.

Hàm suy giảm chấn động và hệ số khuếch đại dao động nền

Sử dụng mô hình suy giảm chấn động của Campbell và Bozorgnia (2008 – CB08) để tính toán dao động nền nằm ngang ($GM$) theo độ lớn động đất ($M_W$), khoảng cách đến nguồn phát sinh và tham số nền đất:

$$K = \frac{GMS}{GMR}$$

Trong đó:

  • $GMR$: Dao động nền ngang tại mặt đá cứng chuẩn ($V_{S30} = 800\text{ m/s}$).
  • $GMS$: Dao động nền ngang tại bề mặt lớp phủ ứng với giá trị $V_{S30}$ thực tế.
  • $K$: Hệ số khuếch đại dao động nền.

Bảng tương quan thực nghiệm $V_S$ và chỉ số $N_{SPT}$

Luận án tổng hợp các tương quan thực nghiệm đã công bố để phục vụ xây dựng mô hình khởi tạo $V_S$ từ số liệu lỗ khoan ĐCCT:

Tác giả công bố Hàm tương quan thực nghiệm ($V_S$ tính bằng m/s) Điều kiện nền đất áp dụng
Ohsaki và Iwasaki (1973) $V_S = 82 \cdot N^{0{,}39}$ Tất cả các loại đất
Imai (1977) $V_S = 128 \cdot N^{0{,}337}$
$V_S = 80{,}6 \cdot N^{0{,}331}$
$V_S = 102 \cdot N^{0{,}292}$		 Tất cả các loại đất
Đất cát
Đất sét
Ohta và Goto (1978) $V_S = 85{,}35 \cdot N^{0{,}348}$ Tất cả các loại đất
Imai và Tonouchi (1982) $V_S = 97 \cdot N^{0{,}314}$ Tất cả các loại đất
Lee (1990) $V_S = 57{,}4 \cdot N^{0{,}49}$
$V_S = 114{,}43 \cdot N^{0{,}31}$		 Tất cả các loại đất
Đất sét
Yokota và nnk (1991) $V_S = 121 \cdot N^{0{,}27}$ Tất cả các loại đất
Athanasopoulos (1995) $V_S = 107{,}6 \cdot N^{0{,}36}$ Tất cả các loại đất
Jafari và nnk (2002) $V_S = 121 \cdot N^{0{,}27}$
$V_S = 80{,}0 \cdot N^{0{,}33}$
$V_S = 100 \cdot N^{0{,}33}$		 Tất cả các loại đất
Đất cát
Đất sét
Hasancebi và Ulusay (2007) $V_S = 90 \cdot N^{0{,}309}$
$V_S = 90{,}8 \cdot N^{0{,}319}$
$V_S = 97{,}9 \cdot N^{0{,}269}$		 Tất cả các loại đất
Đất cát
Đất sét
Hanumanthrao và Ramana (2008) $V_S = 82{,}6 \cdot N^{0{,}43}$
$V_S = 79 \cdot N^{0{,}434}$
$V_S = 86 \cdot N^{0{,}42}$		 Tất cả các loại đất
Đất cát
Đất sét
Uma và nnk (2010) $V_S = 95{,}64 \cdot N^{0{,}301}$
$V_S = 100{,}53 \cdot N^{0{,}265}$
$V_S = 89{,}31 \cdot N^{0{,}358}$		 Tất cả các loại đất
Đất cát
Đất sét
Tsiambaos và Sabatakakis (2011) $V_S = 105{,}7 \cdot N^{0{,}301}$ Tất cả các loại đất
Anbazhagan và nnk (2012) $V_S = 68{,}96 \cdot N^{0{,}51}$ Tất cả các loại đất
Fauzi và nnk (2014) $V_S = 105{,}03 \cdot N^{0{,}268}$ Tất cả các loại đất
Kirar và nnk (2016) $V_S = 99{,}5 \cdot N^{0{,}345}$ Tất cả các loại đất
Shukla và Solanki (2020) $V_S = 77{,}448 \cdot N^{0{,}41}$ Tất cả các loại đất

Nội dung chính theo từng chương

CHƯƠNG 1. TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU VÀ ĐIỀU KIỆN NỀN ĐẤT KHU VỰC NGHIÊN CỨU

Chương 1 hệ thống hóa các kết quả nghiên cứu hiệu ứng nền đất trên phạm vi quốc tế, quy mô quốc gia và khu vực Hà Nội; đồng thời mô tả các đặc trưng địa hình, địa tầng và địa chất công trình vùng nghiên cứu:

  • Địa hình: Phân chia thành 3 khu vực rõ rệt:
    1. Vùng núi thấp: Phía Bắc huyện Sóc Sơn (độ cao từ 270 đến 374 m, đỉnh An Lom 462 m), cấu tạo bởi đá lục nguyên phun trào phong hóa mạnh.
    2. Vùng gò đồi chuyển tiếp: Dải hẹp chân núi Sóc Sơn.
    3. Vùng đồng bằng: Chiếm 90% diện tích, nghiêng về phía Đông Nam, gồm đồng bằng cao (độ cao 6–15 m, tập trung ở Đông Anh, một phần Sóc Sơn và Cổ Nhuế – Xuân Đỉnh) và đồng bằng thấp (độ cao 2–6 m, thuộc các quận trung tâm và phía Nam).
  • Địa tầng Đệ tứ và địa chất công trình: Gồm 5 hệ tầng trầm tích phủ trên nền đá gốc (đá gốc lộ ở phía Bắc Sóc Sơn và chìm sâu trên 100 m tại Hoàng Mai, Long Biên):
    1. Hệ tầng Thái Bình ($Q_2^{tb}$): Trầm tích trẻ nhất gồm sét bột, cát lẫn mùn hữu cơ, cát hạt nhỏ.
    2. Hệ tầng Hải Hưng ($Q_1^{3hh}$ / $Q_2^{1-2hh}$): Bùn, than bùn, sét bột xám đen chứa hữu cơ – là lớp đất yếu nhất.
    3. Hệ tầng Vĩnh Phúc ($Q_1^{3vp}$): Sét bột, cát lẫn sét, cát vàng xây dựng.
    4. Hệ tầng Hà Nội ($Q_1^{2-3hn}$): Cát thô, sỏi sạn, cuội tảng.
    5. Hệ tầng Lệ Chi ($Q_1^{1lc}$): Nằm ở độ sâu 45 đến 69,5 m, thành phần gồm cuội, sỏi, cát, ít sét.
  • Trong giới hạn 60 m trên cùng, các loại đất được phân chia chi tiết thành 26 loại đất địa chất công trình dựa trên đặc điểm thạch học và trạng thái vật lý.

CHƯƠNG 2. PHƯƠNG PHÁP NGHIÊN CỨU HIỆU ỨNG NỀN ĐẤT

Chương 2 trình bày toàn bộ cơ sở toán – lý và phương pháp luận áp dụng trong luận án:

  • Nguyên lý kỹ thuật đo DĐVĐC một trạm và mảng trạm F-K/SPAC.
  • Lý thuyết lan truyền sóng dừng trong môi trường đàn hồi phân lớp xác định tần số trội cơ bản $F_0 = V_S / (4D)$ và các họa ba bậc lẻ ($3F_0, 5F_0$).
  • Cơ chế phân tích tỷ số phổ Fourier $H/V$ theo Nakamura, giả định sóng ngang bị khuếch đại bởi sóng S phản xạ trong khi thành phần đứng đại diện cho sóng P và sóng Rayleigh trong môi trường đá rắn.
  • Công thức giải tích Haskell cho chùm sóng cắt SH truyền trong môi trường phân lớp song song đẳng hướng.
  • Thuật toán di truyền (GA) với toán tử lai ghép, đột biến để giải bài toán tối ưu hàm mục tiêu tương quan phổ $HVSR$.
  • Quy chuẩn phân loại nền đất dựa trên tiêu chuẩn TCVN 9386:2012 (Eurocode 8) chia thành các loại nền A, B, C, D, E theo giá trị $V_{S30}$:
Loại nền đất (TCVN 9386:2012) Mô tả địa tầng / Trạng thái cơ lý Vận tốc sóng ngang trung bình $V_{S30}$ (m/s)
A Đá hoặc địa tầng địa chất tương tự đá $> 800$
B Cát, cuội sỏi rất chặt hoặc sét rất cứng $360 - 800$
C Cát, sỏi chặt/chặt vừa hoặc sét cứng $180 - 360$
D Đất dính từ mềm đến cứng vừa hoặc đất rời rời $< 180$
E Lớp phù sa bồi tích bề mặt có chiều dày $\ge 5\text{ m}$ thuộc loại D Biến thiên

CHƯƠNG 3. NGUỒN SỐ LIỆU VÀ CÁC QUY TRÌNH MINH GIẢI

Chương 3 chi tiết hóa việc thu thập dữ liệu, kiểm tra độ tin cậy và xây dựng các quy trình xử lý chuẩn hóa:

  • Dữ liệu sử dụng: 834 điểm đo DĐVĐC một trạm (máy ghi Samtac VSE315D, VSE355EV, K2/SS-1), các điểm đo mảng DĐVĐC, cùng cơ sở dữ liệu hàng trăm lỗ khoan ĐCCT sâu đến 60 m.
  • Kiểm nghiệm độ tin cậy của số liệu đo:
    • Đo lặp lại theo thời gian: Kiểm tra độ ổn định của giá trị $F_0$ và biên độ đỉnh $A_0$ tại cùng vị trí ở các thời điểm khác nhau.
    • Đo đồng thời bằng nhiều thiết bị: Thực hiện 84 phép đo đồng thời bằng 7 bộ máy Samtac tại 12 điểm đo trong các năm 2009, 2012, 2013; kết quả cho thấy sai số tương đối giữa các thiết bị nằm trong giới hạn cho phép.
  • Các quy trình minh giải:
    1. Quy trình xây dựng biểu đồ tỷ số phổ H/V: Chọn lọc đoạn tín hiệu ổn định, khử xu hướng, lọc dải tần, làm trơn phổ (Konno-Ohmachi), tính tỷ số trung bình hình học $H/V$.
    2. Quy trình xác lập tương quan $F_0 - D$: Kết hợp tần số $F_0$ với chiều dày lớp phủ $D$ tại các vị trí lỗ khoan ĐCCT.
    3. Quy trình mô phỏng $V_{S1D}$ tối ưu bằng GA: Xây dựng mô hình khởi tạo các lớp nông (từ 26 loại đất ĐCCT) và các lớp đá sâu (từ đo mảng DĐVĐC); chạy thuật toán GA để tìm phân bố vận tốc sóng ngang $V_S$ theo độ sâu.
    4. Quy trình thành lập bản đồ phân loại nền và tính hệ số $K$: Tính toán $V_{S30}$, gán mã loại nền theo TCVN 9386:2012 và xác định $K$ qua hàm suy giảm CB08.

CHƯƠNG 4. ĐẶC ĐIỂM DAO ĐỘNG VI ĐỊA CHẤN KHU VỰC NỘI THÀNH THÀNH PHỐ HÀ NỘI VÀ LÂN CẬN

Chương 4 công bố kết quả phân tích trường dao động vi địa chấn từ 834 điểm đo:

  • Hình dạng biểu đồ tỷ số phổ H/V: Phân thành các dạng đặc trưng gồm dạng 1 đỉnh rõ nét, dạng 2 đỉnh (thể hiện 2 mặt ranh giới tương phản âm học), dạng đa đỉnh và dạng đường nằm ngang bằng phẳng (khu vực núi đá lộ hoặc lớp phủ rất mỏng).
  • Phân bố tần số trội $F_0$:
    • Khu vực núi phía Bắc Sóc Sơn không xuất hiện đỉnh trội rõ rệt ($F_0 > 10\text{ Hz}$).
    • Vùng đồng bằng cao (Đông Anh, Bắc Sóc Sơn): $F_0$ dao động từ $1{,}5\text{ Hz}$ đến $5{,}0\text{ Hz}$.
    • Vùng đồng bằng thấp (12 quận nội thành, Gia Lâm, Thanh Trì): $F_0$ giảm mạnh, phổ biến từ $0{,}4\text{ Hz}$ đến $1{,}5\text{ Hz}$.
  • Hàm tương quan thực nghiệm $F_0 - D$: Dựa trên số liệu khớp giữa đo đạc vi địa chấn và lỗ khoan địa chất công trình tại khu vực Hà Nội, thiết lập phương trình tương quan thực nghiệm:

$$D = 84{,}735 \cdot F_0^{-1{,}239}$$

  • Sơ đồ phân vùng chiều dày lớp phủ nông ($D$): Chiều dày $D$ tăng dần từ Bắc xuống Nam và từ Tây sang Đông; đạt giá trị nhỏ ($< 10\text{ m}$) ở Sóc Sơn, trung bình ($20\text{--}40\text{ m}$) ở Đông Anh, Bắc Từ Liêm và lớn nhất ($> 70\text{--}100\text{ m}$) tại Hoàng Mai, Thanh Trì, Gia Lâm.

CHƯƠNG 5. HIỆU ỨNG NỀN ĐẤT KHU VỰC NỘI THÀNH THÀNH PHỐ HÀ NỘI VÀ LÂN CẬN

Chương 5 là kết quả tính toán định lượng hiệu ứng nền đất:

  • Biểu đồ cấu trúc $V_{S1D}$ và giá trị $V_{S30}$: Mô phỏng thành công 816 cấu trúc $V_{S1D}$ tối ưu bằng thuật toán GA với hệ số tương quan giữa phổ SH lý thuyết và phổ HVSR thực nghiệm đạt trên 80%.
  • Sơ đồ phân loại nền đất theo TCVN 9386:2012:
    • Nền loại B ($V_{S30} = 360\text{--}800\text{ m/s}$): Phân bố ở khu vực đồi núi và chân núi phía Bắc huyện Sóc Sơn.
    • Nền loại C ($V_{S30} = 180\text{--}360\text{ m/s}$): Phân bố tại khu vực đồng bằng cao gồm phần lớn huyện Đông Anh, Sóc Sơn, Bắc Từ Liêm, Cầu Giấy, Tây Hồ.
    • Nền loại D ($V_{S30} < 180\text{ m/s}$): Chiếm diện tích lớn ở vùng trũng ngập, ven sông Hồng, sông Đuống và các quận trung tâm phía Nam như Đống Đa, Thanh Xuân, Hai Bà Trưng, Hoàng Mai, Thanh Trì, Gia Lâm.
  • Mô phỏng dao động nền theo 3 kịch bản động đất: Tính toán chấn động trên mặt đá cứng ($GMR$) và mặt lớp phủ ($GMS$) qua 3 kịch bản động đất giả định với các độ sâu chấn tiêu khác nhau trên hệ đứt gãy Sông Hồng/Sông Chảy.
  • Hệ số khuếch đại dao động nền ($K$): Xác lập giá trị hệ số khuếch đại trung bình theo từng loại nền đất khu vực nghiên cứu:
    • Nền loại B: Hệ số $K$ dao động từ $0{,}9$ đến $1{,}2$ (ít khuếch đại).
    • Nền loại C: Hệ số $K$ dao động từ $1{,}3$ đến $1{,}8$.
    • Nền loại D: Hệ số $K$ dao động từ $1{,}8$ đến $2{,}5$ (khuếch đại mạnh nhất tại các khu vực bùn sét dày).

Kết quả và những đóng góp mới

Những đóng góp mới về mặt khoa học và phương pháp

  1. Làm rõ mối quan hệ định lượng giữa đặc điểm trường dao động vi địa chấn (tần số trội $F_0$, tỷ số phổ $HVSR$) với cấu trúc địa chất công trình và điều kiện nền đất địa phương tại khu vực nội thành Hà Nội và lân cận.
  2. Xây dựng hoàn chỉnh quy trình tích hợp mô phỏng chuyển đổi sóng cắt ngang SH và thuật toán di truyền (GA) để nghịch đảo cấu trúc vận tốc sóng ngang $V_{S1D}$ từ dữ liệu đo DĐVĐC một trạm.
  3. Thiết lập hàm tương quan thực nghiệm địa phương giữa tần số trội vi địa chấn $F_0$ và chiều dày lớp phủ trầm tích $D$ dành riêng cho cấu trúc địa tầng khu vực Hà Nội: $D = 84{,}735 \cdot F_0^{-1{,}239}$.

Những đóng góp mới về mặt thực tiễn

  1. Thành lập bộ bản đồ chuyên đề tỷ lệ lớn cho khu vực nội thành Hà Nội và 4 huyện phụ cận ($1000\text{ km}^2$), bao gồm:
    • Bản đồ phân vùng tần số trội $F_0$ và chiều dày lớp phủ nông $D$.
    • Bản đồ phân vùng giá trị vận tốc sóng ngang $V_{S30}$.
    • Sơ đồ phân loại nền đất theo Tiêu chuẩn thiết kế công trình chịu động đất TCVN 9386:2012 trên cơ sở 834 điểm đo thực tế.
  2. Xác định định lượng hệ số khuếch đại dao động nền ($K$) cho các loại nền B, C, D trên cơ sở 816 giá trị $V_{S30}$ tính toán bằng thuật toán GA.
  3. Cung cấp cơ sở dữ liệu số về đặc trưng động lực học nền đất phục vụ trực tiếp công tác quy hoạch không gian đô thị, quản lý rủi ro thiên tai địa chấn và tính toán kháng chấn cho các công trình xây dựng tại thủ đô Hà Nội.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế:
    • Nghiên cứu tập trung vào mô hình cấu trúc 1 chiều (1D) theo phương thẳng đứng, chưa xét đến hiệu ứng bồn trũng 2D/3D phức tạp tại các vùng ranh giới uốn lượn ven sông.
    • Dải tần số khảo sát phụ thuộc vào giới hạn kỹ thuật của đầu đo vi địa chấn (VSE315D, VSE355EV, SS-1); việc đánh giá lớp đá gốc sâu trên 100 m chủ yếu ngoại suy qua đo mảng và số liệu địa chất sâu.
    • Các kịch bản tính toán gia tốc dao động nền sử dụng hàm suy giảm CB08 toàn cầu, cần thêm số liệu ghi nhận động đất mạnh tại chỗ để hiệu chỉnh.
  • Hướng nghiên cứu tiếp:
    • Phát triển mô hình lan truyền sóng 2D và 3D nhằm đánh giá chi tiết hiệu ứng tập trung năng lượng sóng tại các thung lũng, lòng sông cổ và đới đứt gãy.
    • Mở rộng mạng lưới trạm quan trắc gia tốc kế chuyển động mạnh (Strong Motion) trong khu vực đô thị để ghi nhận trực tiếp các trận động đất thực tế, từ đó kiểm chứng và hoàn thiện hệ số khuếch đại nền địa phương.

Giá trị tham khảo

  • Đối với các nhà nghiên cứu, nghiên cứu sinh và học viên cao học: Luận án là tài liệu tham khảo trong các chuyên ngành Địa vật lý, Địa chất công trình, Địa kỹ thuật xây dựng và Địa chấn công trình. Phần phương pháp luận về nghịch đảo phổ H/V bằng thuật toán di truyền GA (Chương 2 và Chương 3) là tài liệu mẫu mực về ứng dụng tin học tính toán trong khoa học Trái đất.
  • Đối với các cơ quan quản lý và quy hoạch đô thị: Cung cấp cơ sở dữ liệu khoa học để Viện Quy hoạch Xây dựng Hà Nội, Sở Xây dựng Hà Nội tham khảo khi cấp phép xây dựng, quy hoạch mạng lưới hạ tầng kỹ thuật và bố trí các công trình cao tầng, công trình ngầm tránh vùng nền đất yếu có nguy cơ cộng hưởng dao động cao.
  • Đối với các kỹ sư tư vấn thiết kế công trình: Bản đồ phân loại nền đất theo TCVN 9386:2012 và bảng hệ số khuếch đại $K$ theo từng loại nền đất (Chương 5) cung cấp thông số đầu vào trực tiếp để xác định phổ phản ứng thiết kế và hệ số nền ($S$) khi tính toán tải trọng động đất cho công trình.

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

1. Luận án sử dụng phương pháp nào để xác định cấu trúc vận tốc sóng ngang $V_S$ tại các điểm đo?

Luận án sử dụng phương pháp đo dao động vi địa chấn một trạm để thu thập dữ liệu dao động 3 thành phần, tính toán đường cong tỷ số phổ $HVSR$. Sau đó, áp dụng thuật toán di truyền (GA) để nghịch đảo mô hình truyền sóng cắt SH (theo lý thuyết ma trận Haskell), tìm ra biểu đồ cấu trúc vận tốc sóng ngang 1 chiều ($V_{S1D}$) tối ưu khớp với phổ thực nghiệm.

2. Số lượng điểm đo và khối lượng dữ liệu thực tế phục vụ cho kết quả luận án là bao nhiêu?

Luận án đã thu thập và xử lý tổng cộng 834 điểm đo dao động vi địa chấn một trạm (kế thừa và đo bổ sung phủ kín 12 quận và 4 huyện ngoại thành), thực hiện 84 phép đo kiểm tra đồng thời bằng 7 bộ máy Samtac, phân tích các điểm đo mảng F-K, đồng thời tính toán thành công 816 mô hình $V_{S1D}$ và giá trị $V_{S30}$ chi tiết.

3. Phương trình tương quan thực nghiệm giữa chiều dày lớp phủ $D$ và tần số trội $F_0$ được xác lập như thế nào?

Dựa trên việc đối sánh tần số trội $F_0$ từ phổ vi địa chấn và chiều dày $D$ xác định từ các lỗ khoan địa chất công trình, luận án xác lập hàm tương quan thực nghiệm: $D = 84{,}735 \cdot F_0^{-1{,}239}$. Phương trình này cho phép ước tính nhanh chiều dày lớp đất mềm bồi tích từ số liệu đo vi địa chấn trên mặt đất tại khu vực Hà Nội.

4. Kết quả phân loại nền đất khu vực nội thành Hà Nội theo TCVN 9386:2012 phân bố ra sao?

Khu vực nghiên cứu phân bố chủ yếu 3 loại nền đất:

  • Nền loại B: Phân bố hẹp tại vùng núi đá và đồi thấp huyện Sóc Sơn ($V_{S30} = 360\text{--}800\text{ m/s}$).
  • Nền loại C: Chiếm phần lớn diện tích đồng bằng cao như Đông Anh, Bắc Từ Liêm, Tây Hồ ($V_{S30} = 180\text{--}360\text{ m/s}$).
  • Nền loại D: Phân bố tại các vùng đồng bằng thấp, trũng ngập, dải ven sông Hồng, Hoàng Mai, Thanh Trì, Gia Lâm, Đống Đa, Thanh Xuân ($V_{S30} < 180\text{ m/s}$), đây là khu vực có nền đất yếu nhất và có hệ số khuếch đại dao động lớn nhất ($K$ đạt từ $1{,}8$ đến $2{,}5$).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Tiến Hùng đã giải quyết bài toán ước lượng hiệu ứng nền đất tại khu vực nội thành thành phố Hà Nội và phụ cận bằng phương pháp tiếp cận định lượng hiện đại, kết hợp chặt chẽ giữa đo đạc dao động vi địa chấn mật độ cao với số liệu địa chất công trình và kỹ thuật nghịch đảo thuật toán di truyền. Công trình đã xác lập thành công bản đồ phân loại nền đất theo tiêu chuẩn TCVN 9386:2012 và định lượng hệ số khuếch đại dao động nền cho các loại nền B, C, D trên diện tích $1000\text{ km}^2$. Kết quả nghiên cứu là tài liệu khoa học có độ tin cậy cao, đóng góp cơ sở dữ liệu quan trọng cho công tác quy hoạch đô thị và thiết kế kháng chấn nhằm giảm thiểu rủi ro động đất tại thủ đô Hà Nội.