Tổng quan nghiên cứu

Đập vật liệu địa phương hiện chiếm hơn 75% tổng số công trình đập dâng giữ nước tại Việt Nam và giữ vai trò cốt lõi trong hệ thống hạ tầng thủy lợi, thủy điện quốc gia. Tuy nhiên, với vị trí địa lý nằm trong khu vực có hoạt động kiến tạo phức tạp cùng các hệ thống đứt gãy hoạt động mạnh như đứt gãy Sông Mã, đứt gãy Sơn La và đứt gãy Lai Châu - Điện Biên, nguy cơ xảy ra động đất là rất hiện hữu. Trong lịch sử, Việt Nam từng ghi nhận các trận động đất mạnh có độ lớn M = 6,75 tại Điện Biên năm 1935 và M = 6,8 tại Tuần Giáo năm 1983, gây sụt lở đất trên diện rộng và đe dọa nghiêm trọng đến an toàn công trình.

Trước đây, hầu hết các công trình đập đất đá tại Việt Nam đều được tính toán kháng chấn theo phương pháp giả tĩnh truyền thống. Phương pháp này chỉ đơn giản quy đổi tác động động đất thành một lực quán tính không đổi, bỏ qua tính chất dao động theo thời gian, hiện tượng tích tụ áp lực nước lỗ rỗng và sự suy giảm độ bền của đất. Bài học thực tế từ sự cố sạt trượt nghiêm trọng của mái thượng lưu đập Lower San Fernando cao hơn 40 m vào năm 1971 tại Mỹ đã chứng minh rằng những đập được thiết kế đạt hệ số an toàn tĩnh trên 1,3 vẫn có thể bị phá hủy hoàn toàn khi chịu tải trọng động đất thực tế.

Nghiên cứu tập trung giải quyết bài toán phân tích phản ứng động lực học của đập vật liệu địa phương dưới tác động của các giản đồ gia tốc động đất thực tế. Mục tiêu cụ thể là xác định chính xác trường chuyển vị, phân bố vận tốc, gia tốc khuếch đại, sự biến đổi của ứng suất hữu hiệu và nguy cơ hóa lỏng trong thân đập. Đề tài áp dụng tính toán thực chứng trên công trình đập chính Hàm Thuận có chiều cao gần 100 m, cung cấp các chỉ số an toàn định lượng với mức gia tốc đỉnh nền dao động từ 0,08g đến 0,20g, tạo cơ sở khoa học tin cậy cho công tác thiết kế và quản lý rủi ro đập dâng nước.

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 đồng thời ba khung lý thuyết cơ bản trong địa kỹ thuật động lực để mô phỏng ứng xử phi tuyến phức tạp của vật liệu đất đá:

  • Mô hình đàn hồi tuyến tính: Ứng dụng định luật Hooke để mô tả trạng thái ban đầu của các khối kết cấu cứng trong thân đập, thông qua mô đun đàn hồi E và hệ số Poisson.
  • Mô hình tuyến tính tương đương: Dựa trên các nghiên cứu kinh điển của Seed và Idriss năm 1970, mô hình sử dụng quan hệ trễ kiểu Masing để mô tả quá trình tiêu tán năng lượng khi đất chịu tải trọng chu kỳ. Trọng tâm của mô hình là xác định đường cong suy giảm mô đun cắt theo biến dạng cắt và đường cong biến thiên của tỷ số giảm chấn.
  • Mô hình phi tuyến Hyperbolic: Phát triển từ nghiên cứu của Hardin và Drnevich năm 1972 và cải tiến bởi Darendeli, mô tả chính xác lộ trình ứng suất - biến dạng phi tuyến thực tế, xét đến ảnh hưởng của áp lực nhốt hữu hiệu, chỉ số dẻo và tỷ số quá cố kết.

Bên cạnh đó, các khái niệm then chốt như Tỷ số ứng suất tuần hoàn (CSR), Tỷ số kháng cắt tuần hoàn (CRR), tỷ số áp lực nước lỗ rỗng dư và cơ chế hóa lỏng đất cát bão hòa được tích hợp đầy đủ để đánh giá khả năng phá hoại kết cấu.

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

Mô hình nghiên cứu sử dụng phương pháp phần tử hữu hạn (FEM) hai chiều thông qua bộ công cụ Geostudio 2007, kết hợp chặt chẽ giữa hai mô đun SEEP/W và QUAKE/W:

  • Dữ liệu và lưới phần tử: Mặt cắt đập Hàm Thuận được chia lưới với quy mô 1.293 nút và 1.240 phần tử phẳng tứ giác và tam giác bậc hai. Dữ liệu địa chất đầu vào bao gồm chỉ số thí nghiệm xuyên tiêu chuẩn SPT dao động từ 15 đến 45 búa, hệ số rỗng từ 0,45 đến 0,75 và dung trọng khô từ 1,65 đến 2,15 tấn trên mét khối.
  • Dữ liệu động đất: Lựa chọn giản đồ gia tốc thực đo của trận động đất Irpinia năm 1980 với phổ phản ứng tương thích với phổ thiết kế của nền loại B theo tiêu chuẩn TCVN 375-2006. Tải trọng động được mô phỏng liên tục trong 36 giây với bước tích phân thời gian là 0,02 giây, tương đương 1.800 bước tính toán chi tiết.
  • Lý do lựa chọn phương pháp: Quy trình liên hợp cho phép SEEP/W xác định chính xác trạng thái ứng suất hữu hiệu và mạng thấm tĩnh ban đầu trước khi QUAKE/W tiếp nhận dữ liệu để mô phỏng chính xác sự gia tăng áp lực nước lỗ rỗng và biến dạng động theo từng phần trăm giây.

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

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

Phân tích số liệu từ 1.800 bước thời gian cho thấy các phản ứng cơ học nổi bật của đập vật liệu địa phương:

  • Hiện tượng khuếch đại dao động tại đỉnh đập: Gia tốc ngang cực đại tại đỉnh đập ghi nhận mức khuếch đại gấp 1,8 đến 2,4 lần so với gia tốc đỉnh của nền đất. Chuyển vị ngang tích lũy cực đại của đỉnh đập đạt 0,26 m và tập trung chủ yếu vào khoảng thời gian dao động mạnh nhất từ giây thứ 12 đến giây thứ 16 của trận động đất.
  • Sự phân bố lại và tập trung ứng suất cắt: Ứng suất cắt cực đại tại vùng tiếp giáp giữa lõi sét chống thấm và khối đá gia tải hai phía tăng mạnh từ mức tĩnh 145 kPa lên giá trị đỉnh 210 kPa, tương ứng với mức tăng 44,8%. Hiện tượng này đi kèm với sự sụt giảm từ 30% đến 55% ứng suất hữu hiệu theo phương đứng tại các vị trí xung yếu.
  • Gia tăng áp lực nước lỗ rỗng và vùng nguy cơ hóa lỏng: Tại vùng chân đập và phần dưới của mái thượng lưu bão hòa nước, áp lực nước lỗ rỗng tăng vọt từ 65 kPa lên vượt ngưỡng 180 kPa. Tỷ số áp lực nước lỗ rỗng dư tại vùng này đạt từ 0,65 đến 0,78 ở độ sâu từ 10 m đến 25 m, làm xuất hiện các vùng có nguy cơ hóa lỏng cục bộ nếu gia tốc nền vượt quá 0,15g.

Thảo luận kết quả

Khi chịu tác động của tải trọng động chu kỳ, các hạt đất trong thân đập bị sắp xếp lại làm gia tăng áp lực nước kẽ rỗng tức thời. Do tính chất thoát nước chậm của đất đắp, ứng suất hữu hiệu bị suy giảm nghiêm trọng, kéo theo sự giảm sút của mô đun cắt tiếp tuyến xuống chỉ còn khoảng 15% đến 30% so với giá trị mô đun cắt cực đại ban đầu.

So sánh giữa hai cách tiếp cận, mô hình phi tuyến Hyperbolic cho kết quả phản ánh chân thực hơn về biến dạng dẻo vĩnh cửu và sự tập trung trượt cục bộ so với mô hình tuyến tính tương đương. Mô hình tuyến tính tương đương có xu hướng đánh giá thấp độ lún đỉnh đập khoảng 20% đến 35% do không xét đến quá trình tích lũy biến dạng không phục hồi sau từng chu kỳ tải.

Trong thực tế tính toán và báo cáo kỹ thuật, toàn bộ diễn biến động lực học này được trực quan hóa rất hiệu quả thông qua các đồ thị lịch sử thời gian tại các nút quan trắc đại diện như nút 4 tại đỉnh đập, nút 45 và nút 51 tại sườn đập, kết hợp với các biểu đồ đường đẳng mức thể hiện trường ứng suất và vùng hóa lỏng 2D trên toàn bộ mặt cắt công trình.

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

Dựa trên các kết quả mô phỏng động lực học, nghiên cứu đưa ra 4 nhóm giải pháp cụ thể nhằm đảm bảo an toàn tuyệt đối cho các đập vật liệu địa phương:

  • Lắp đặt hệ thống quan trắc gia tốc và áp lực nước lỗ rỗng tự động: Trang bị mạng lưới cảm biến đo gia tốc đa chiều và đầu đo áp lực nước lỗ rỗng điện tử tại đỉnh đập, thân đập và chân đập với tần số lấy mẫu tối thiểu 100 Hz. Hệ thống phải đảm bảo truyền dữ liệu thời gian thực và kích hoạt cảnh báo sự cố trong vòng dưới 3 giây khi rung chấn vượt ngưỡng 0,05g. Thời gian triển khai hoàn thành trong 6 đến 12 tháng do Ban quản lý công trình và các đơn vị vận hành hồ chứa chủ trì.
  • Gia cố kết cấu thoát nước và nâng cao độ ổn định mái dốc: Thiết kế bổ sung tầng đệm lọc ngược bằng vật liệu đá dăm chọn lọc có hệ số thấm lớn hơn 0,01 cm/s và góc ma sát trong đạt từ 38 độ trở lên tại vùng chân đập thượng và hạ lưu. Giải pháp này giúp triệt tiêu nhanh áp lực nước lỗ rỗng dư, nâng hệ số an toàn trượt động từ mức 1,05 lên trên 1,30. Thời gian thực hiện dự kiến 18 tháng do đơn vị tư vấn và nhà thầu xây lắp đảm nhiệm.
  • Cập nhật quy chuẩn thiết kế kháng chấn cho đập dâng nước: Thay thế dần phương pháp giả tĩnh bằng phương pháp phân tích lịch sử thời gian sử dụng phần tử hữu hạn phi tuyến đối với tất cả các đập cấp I và cấp đặc biệt có chiều cao trên 30 m hoặc dung tích trên 50 triệu mét khối. Lộ trình thực hiện trong giai đoạn từ năm 2024 đến năm 2026 dưới sự chỉ đạo của Bộ Xây dựng và Bộ Nông nghiệp và Phát triển Nông thôn.
  • Định kỳ đánh giá lại an toàn địa chấn công trình: Tổ chức kiểm toán định kỳ 5 năm một lần đối với hiện trạng hồ đập, ứng dụng mô hình số cập nhật theo dữ liệu đo đạc thực tế để kiểm soát độ lún dư đỉnh đập luôn dưới ngưỡng 0,5% đến 1,0% tổng chiều cao đập. Nhiệm vụ do các cơ quan kiểm định an toàn đập độc lập và các viện nghiên cứu chuyên ngành thực hiện.

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

Nội dung và phương pháp nghiên cứu trong luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư tư vấn và thiết kế công trình thủy lợi - thủy điện: Tiếp cận quy trình mô hình hóa liên hợp giữa SEEP/W và QUAKE/W để tối ưu hóa hình học mặt cắt đập, lựa chọn thông số vật liệu và bố trí tầng lọc kháng chấn cho các dự án xây mới chịu động đất từ cấp 7 đến cấp 9.
  • Đơn vị quản lý và vận hành hồ chứa thủy điện: Nắm rõ các quy luật tập trung ứng suất và các vị trí xung yếu ở cao trình từ 580 m đến 620 m để xây dựng quy trình quan trắc rung chấn, kiểm tra thân đập và lập phương án ứng phó khẩn cấp khi có động đất xảy ra.
  • Giảng viên, nghiên cứu sinh và học viên cao học: Sử dụng tài liệu làm tài liệu tham khảo chuyên khảo về cơ học đất động lực, mô hình trễ Masing, cơ chế hóa lỏng đất và phương pháp số giải bài toán tương tác đất - công trình dưới tải trọng chu kỳ.
  • Cơ quan quản lý nhà nước về an toàn đập: Có thêm cơ sở khoa học vững chắc để xây dựng các văn bản hướng dẫn kỹ thuật, tiêu chí kiểm toán an toàn và chính sách quản lý rủi ro thiên tai cho hơn 7.000 hồ chứa trên toàn quốc.

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

Phương pháp phân tích lịch sử thời gian vượt trội hơn phương pháp giả tĩnh như thế nào?
Phương pháp giả tĩnh chỉ áp dụng một lực quán tính không đổi từ 0,05g đến 0,15g, không phản ánh được tính chất dao động và sự tích tụ biến dạng dẻo. Ngược lại, phân tích lịch sử thời gian mô phỏng chi tiết 1.800 bước tải, xác định chính xác thời điểm đỉnh đập đạt chuyển vị cực đại 0,26 m và mô tả sự suy giảm sức kháng cắt của đất.

Tại sao cần kết hợp đồng thời hai mô đun SEEP/W và QUAKE/W trong tính toán?
Sự kết hợp này giải quyết trọn vẹn bài toán liên hợp thủy lực - cơ học. SEEP/W thiết lập mạng thấm ổn định và trường áp lực nước lỗ rỗng ban đầu trên 1.200 phần tử. Dữ liệu này là điều kiện biên bắt buộc để QUAKE/W tính toán chính xác biến dạng động và dự báo vùng hóa lỏng trong suốt 36 giây động đất.

Vị trí nào trong đập vật liệu địa phương dễ bị tổn thương nhất khi có động đất?
Đỉnh đập và vùng chuyển tiếp mái thượng lưu là hai khu vực chịu ảnh hưởng nặng nề nhất. Đỉnh đập chịu hiện tượng khuếch đại gia tốc nền lên gấp 1,8 đến 2,4 lần gây nứt nẻ và lún sụt, trong khi mái thượng lưu bão hòa nước có áp lực nước lỗ rỗng tăng từ 40% đến 55%, dễ dẫn đến trượt mái dốc.

Hiện tượng hóa lỏng trong thân đập xảy ra dưới những điều kiện nào?
Hóa lỏng xuất hiện khi đất cát bão hòa hoặc vật liệu hạt mịn chịu tải trọng rung động lặp lại từ 10 đến 30 chu kỳ. Áp lực nước kẽ rỗng gia tăng làm tỷ số áp lực dư đạt mức từ 0,7 đến 1,0, triệt tiêu hoàn toàn ứng suất hữu hiệu và khiến đất mất đi khả năng chịu lực cắt.

Cách thức lựa chọn giản đồ gia tốc phù hợp cho công trình đập tại Việt Nam?
Kỹ sư căn cứ vào bản đồ phân vùng gia tốc nền theo TCVN 375-2006 với giá trị từ 0,08g đến 0,16g và phân loại nền đất từ A đến E. Sau đó, sử dụng các phần mềm chuyên dụng như PGMD để tìm kiếm các sóng động đất thực đo có phổ phản ứng gia tốc tương thích trong dải chu kỳ từ 0,1 giây đến 2,0 giây.

Kết luận

  • Làm sáng tỏ bản chất ứng xử động lực học phi tuyến và phản ứng biến dạng của đập vật liệu địa phương dưới tác động của chấn động địa chấn thực tế.
  • Xác định rõ quy luật khuếch đại gia tốc đỉnh tại đỉnh đập đạt từ 1,8 đến 2,4 lần và chuyển vị ngang cực đại đạt mức 0,26 m.
  • Định lượng chính xác sự gia tăng áp lực nước lỗ rỗng và khoanh vùng nguy cơ hóa lỏng cục bộ tại mái thượng lưu với tỷ số áp lực dư từ 0,65 đến 0,78.
  • Xây dựng hoàn chỉnh quy trình tính toán số liên hợp giữa SEEP/W và QUAKE/W phục vụ kiểm toán an toàn đập dâng quy mô lớn.
  • Đề xuất 4 nhóm giải pháp công trình và quan trắc giúp nâng cao hệ số an toàn kháng chấn cho công trình thêm từ 20% đến 30%.

Đóng góp lớn nhất của luận văn là đã cung cấp phương pháp luận tính toán động lực học thực chứng, tạo bước chuyển đổi cần thiết từ phương pháp giả tĩnh truyền thống sang phân tích lịch sử thời gian hiện đại cho ngành xây dựng công trình thủy tại Việt Nam.

Trong giai đoạn từ năm 2024 đến năm 2026, nghiên cứu định hướng mở rộng mô hình phân tích không gian ba chiều và xét đến tương tác phức tạp giữa đập, hồ chứa và nền đá sâu. Các đơn vị quản lý, tư vấn thiết kế và nghiên cứu chuyên ngành được khuyến khích áp dụng ngay quy trình mô phỏng số này vào công tác thẩm định và kiểm định an toàn đập để chủ động bảo vệ nguồn tài nguyên và an toàn hạ du.