Tổng quan nghiên cứu

Theo thống kê của Hội Đập lớn Thế giới (ICOLD), toàn cầu hiện có hơn 45.000 đập lớn, trong đó Việt Nam sở hữu khoảng 460 công trình, đứng thứ 16 trên thế giới. Cùng với sự phát triển của công nghệ bê tông đầm lăn (RCC), hàng loạt công trình thủy điện quy mô lớn đã được xây dựng như Thủy điện Sơn La (cao 138,1 m), Thủy điện Lai Châu (cao 137 m) và Thủy điện Bản Vẽ (cao 137 m). Tại các đập bê tông trọng lực này, hệ thống tràn xả lũ đóng vai trò quyết định đến an toàn công trình đầu mối, đòi hỏi các trụ pin phân khoang có chiều cao vượt bậc từ 15 m đến 30 m để bố trí dầm cầu giao thông, cầu công tác và hệ thống máy đóng mở cửa van cung hoặc van phẳng.

Tuy nhiên, đặc điểm kết cấu của các trụ đập tràn cao thường có bề rộng mỏng, dẫn tới tỷ lệ độ mảnh hình học rất lớn. Khi chịu tác động của tải trọng động đất, hiện tượng rung lắc theo phương ngang dọc trục đập sẽ tạo ra lực quán tính xô ngang rất mạnh. Dao động này gây ra hai nguy cơ phá hoại nghiêm trọng: tạo chuyển vị ngang lớn làm biến dạng khe van gây kẹt cứng cửa van xả lũ, và phát sinh ứng suất kéo uốn cục bộ cực đại tại chân trụ vượt quá khả năng chịu lực của bê tông, dẫn đến nứt gãy sụp đổ công trình. Thực tế đã ghi nhận các sự cố nguy hiểm như tại Thủy điện Hố Hô năm 2010 khi cửa van kẹt khiến nước tràn qua đỉnh đập 1 m làm phá hỏng nhà máy, hay chuỗi động đất kích thích tại Thủy điện Sông Tranh 2 năm 2012 làm mở rộng 10 khe nhiệt từ 2 mm đến 6 mm với lưu lượng thấm lên tới 75 lít/giây.

Trước thực trạng các quy trình thiết kế truyền thống chủ yếu tập trung vào ổn định trượt lật của thân đập mà chưa đánh giá đầy đủ dao động ngang của trụ pin, luận văn tập trung nghiên cứu định lượng ảnh hưởng của động đất (cấp 7, 8 và 9) đến ứng suất và biến dạng ngang của trụ pin đập tràn cao. Nghiên cứu thực hiện phân tích trên 27 kịch bản không gian tổng quát và áp dụng tính toán thực chứng cho công trình Thủy điện Lai Châu. Kết quả nghiên cứu xác lập quan hệ độ mảnh giới hạn và định lượng hiệu quả của giải pháp dầm giằng ngang, giúp kiểm soát chuyển vị đỉnh trụ dưới 5,5 cm và giảm thiểu ứng suất kéo nguy hiểm tại chân trụ tới 67,0%, đảm bảo an toàn vận hành tuyệt đối cho hệ thống xả lũ.

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 dựa trên ba trụ cột lý thuyết cơ học và tiêu chuẩn kháng chấn công trình thủy:

Thứ nhất là lý thuyết cơ học đàn hồi và độ bền kết cấu bê tông khối lớn. Khối trụ pin đập tràn được mô hình hóa như một kết cấu dạng console thẳng đứng ngàm vào thân đập tràn, làm việc dưới trạng thái nén uốn lệch tâm không gian phức tạp. Khi có gia tốc địa chấn theo phương trục đập, lực xô ngang tác động làm phát sinh ứng suất kéo chính lớn nhất tại mép chân trụ.

Thứ hai là lý thuyết trạng thái giới hạn chống nứt bê tông cốt thép. Điều kiện ổn định chống nứt của trụ pin chịu uốn ngang được xác định thông qua hệ số an toàn Kn = sigma_1n / S1 >= Kcp. Trong đó, S1 là ứng suất kéo lớn nhất sinh ra tại mép chân trụ khi có động đất; sigma_1n là ứng suất kéo giới hạn gây nứt của bê tông, được xác định qua công thức sigma_1n = mh * eta * Rkc. Hệ số biến dạng dẻo của bê tông miền kéo lấy eta = 1,75; hệ số xét đến chiều dày trụ mh = 0,9 + 0,1D (với D là chiều dày trụ tính bằng mét); Rkc là cường độ chịu kéo tiêu chuẩn của các cấp mác bê tông M300, M350 và M400; Kcp là hệ số an toàn cho phép tương ứng với tổ hợp tải trọng đặc biệt.

Thứ ba là lý thuyết động lực học công trình và phương pháp phổ phản ứng đàn hồi. Phổ phản ứng gia tốc thiết kế được thiết lập với hệ số động lực cực đại alpha_max = 2,0, chu kỳ dao động góc Tg = 0,2 giây và gia tốc nền tương ứng với các cấp động đất 7, 8, 9 theo phân vùng địa chấn khu vực miền Bắc Việt Nam (hệ số địa chấn k dao động từ 0,1 đến 0,4).

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

Nghiên cứu ứng dụng Phương pháp phần tử hữu hạn (FEM) không gian 3 chiều trên nền tảng phần mềm ANSYS 16.0 kết hợp xây dựng mô hình hình học chi tiết từ AutoCAD. Mô hình bài toán bao gồm hệ tương tác liên hợp "Khoang mẫu đập tràn – Trụ pin – Nền đá cứng". Vùng nền đá được mô phỏng mở rộng về phía thượng lưu và hạ lưu một khoảng bằng bề rộng đáy đập B, chiều sâu xuống đáy nền bằng chiều cao đập H nhằm mô tả chính xác sự phân bố biến dạng đàn hồi tại vùng tiếp giáp.

Cỡ mẫu nghiên cứu tổng quát gồm 27 kịch bản tính toán độc lập, được lựa chọn theo phương pháp phân tầng tổ hợp đầy đủ giữa các thông số biến thiên:

  • Chiều cao trụ pin trên đỉnh tràn: H1 = 20 m, 25 m, 30 m (chiều cao tường phân dòng H2 = 0,4 * H1).
  • Chiều dày trụ D được xác định theo tỷ lệ độ mảnh lambda = H/D = 5, 6, 7 (tương ứng chiều dày từ 4 m đến 6 m).
  • Cấp chấn động địa chấn: cấp 7, cấp 8 và cấp 9.
  • Cầu công tác được mô hình liên kết ngàm liên tục và cầu giao thông liên kết khớp gối tự do lên đỉnh trụ.

Dữ liệu kiểm chứng thực tế được thu thập từ hồ sơ thiết kế và tài liệu quan trắc địa kỹ thuật của dự án Thủy điện Lai Châu (công trình cấp đặc biệt trên đới đứt gãy sông Đà có nguy cơ chấn động lên đến 6,8 độ Richter). Việc lựa chọn phương pháp phân tích phổ và kiểm tra bằng phương pháp lịch sử thời gian cho phép xác định chính xác vị trí tập trung ứng suất cục bộ và giá trị chuyển vị tức thời lớn nhất, khắc phục hoàn toàn nhược điểm xem dao động mọi điểm là đồng nhất của phương pháp tĩnh lực học cổ điển.

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

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

Quá trình mô phỏng số trên 27 kịch bản tổng quát và bài toán ứng dụng Thủy điện Lai Châu đã mang lại những phát hiện quan trọng:

Thứ nhất, mức độ gia tăng chuyển vị ngang tỷ lệ thuận phi tuyến với độ mảnh và cấp động đất. Khi chiều cao trụ pin tăng từ 20 m lên 30 m với độ mảnh lambda = 7, chuyển vị ngang đỉnh trụ Uz dưới tác động của động đất cấp 9 tăng tới 165% so với trường hợp động đất cấp 7. Biên độ dao động ngang lớn ở đỉnh trụ là nguyên nhân trực tiếp đe dọa làm biến dạng các dầm đỡ cửa van và gây kẹt khe van có kích thước tiêu chuẩn 0,5 m x 0,5 m.

Thứ hai, hiện tượng tập trung ứng suất kéo nguy hiểm tại vị trí ngàm chân trụ. Khi xét trụ có chiều cao H = 30 m, bề dày D = 4 m (độ mảnh lambda = 7,5), dưới tác động của động đất cấp 8, ứng suất kéo chính S1 tại mép chân trụ vượt mức 500 T/m2 (tương đương 5,0 MPa). Giá trị này vượt xa cường độ chịu kéo tiêu chuẩn của bê tông mác M300 (khoảng 140 T/m2 hay 1,4 MPa), dẫn tới việc hình thành các vết nứt uốn xuyên suốt chân trụ và đẩy kết cấu vào trạng thái phá hủy giòn.

Thứ ba, hiệu quả vượt trội của giải pháp dầm chống ngang tại Thủy điện Lai Châu. Khi chưa xảy ra động đất (tổ hợp tải trọng cơ bản), chuyển vị ngang đỉnh trụ là 4,9 cm. Khi xuất hiện động đất cấp 8 trên mô hình không có giằng chống, chuyển vị Uz tăng lên 6,6 cm (tăng 34,7%) và ứng suất kéo tại chân trụ đạt tới 540 T/m2, không đảm bảo điều kiện an toàn chống nứt. Tuy nhiên, khi bổ sung dầm chống ngang BTCT liên kết giữa các trụ pin, chuyển vị ngang giảm xuống còn 5,3 cm (giảm 19,7%) và ứng suất kéo cực đại S1 giảm mạnh từ 540 T/m2 xuống 178 T/m2 (giảm tới 67,0%), đưa hệ số an toàn chống nứt Kn đạt mức 1,25, thỏa mãn tiêu chuẩn an toàn cho phép Kcp.

Thảo luận kết quả

Cơ chế cơ học giải thích cho sự suy giảm ổn định của trụ pin đập tràn cao nằm ở sự cộng hưởng dao động. Trụ pin thanh mảnh có chu kỳ dao động riêng nằm trong dải tần số đặc trưng của sóng địa chấn kiến tạo vùng Tây Bắc (chu kỳ 0,1 đến 0,3 giây), làm khuếch đại lực quán tính rung lắc ngang. Khi các trụ rung lắc độc lập, tải trọng động đất tác động đồng thời với áp lực thủy tĩnh truyền từ cửa van tạo ra mô men uốn tập trung cực lớn tại tiết diện chân trụ tiếp giáp đỉnh đập tràn.

Khi so sánh với các báo cáo khảo sát sự cố rò rỉ nước tại Thủy điện Sông Tranh 2 và kẹt cửa van tại Thủy điện Hố Hô, kết quả tính toán của luận văn khẳng định tính đúng đắn về sự tồn tại của biến dạng lệch ngang gây hở khe lún và kẹt cơ khí. Khi bổ sung dầm giằng ngang, hệ thống các trụ pin độc lập chuyển hóa thành hệ khung không gian làm việc đồng thời, giúp phân phối lại nội lực và tăng độ cứng chống uốn tổng thể theo phương ngang trục đập.

Các kết quả nghiên cứu đã được trực quan hóa thông qua hệ thống biểu đồ quan hệ giữa chuyển vị ngang Uz và độ mảnh lambda, biểu đồ quan hệ giữa ứng suất kéo S1 và chiều dày D theo từng cấp mác bê tông M300, M350, M400. Đặc biệt, luận văn đã xây dựng thành công ma trận bảng tra độ mảnh cho phép lambda_cp: đối với động đất cấp 7, độ mảnh cho phép của bê tông M300 là 6,2; nhưng khi động đất nâng lên cấp 9, công trình bắt buộc phải sử dụng bê tông M400 và khống chế độ mảnh lambda dưới 4,8 để ngăn chặn hoàn toàn nguy cơ nứt gãy.

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

Dựa trên các phân tích định lượng, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể nhằm nâng cao độ ổn định hướng ngang của trụ pin đập tràn cao:

  1. Chuẩn hóa giới hạn độ mảnh hình học trong thiết kế sơ bộ: Các đơn vị tư vấn thiết kế công trình thủy cần khống chế nghiêm ngặt tỷ lệ độ mảnh lambda = H/D không vượt quá 5,0 đối với các công trình xây dựng tại vùng động đất cấp 8 và không quá 4,5 tại vùng động đất cấp 9. Quy định này cần được đưa vào chỉ dẫn thiết kế cơ sở trong giai đoạn lập báo cáo khả thi để loại bỏ nguy cơ quá tải ứng suất uốn ngay từ bước tạo hình kết cấu.

  2. Bố trí kết cấu dầm chống ngang liên kết giữa các trụ pin: Ban quản lý dự án và các kỹ sư thiết kế cần tích hợp hệ dầm giằng bê tông cốt thép tại cao trình từ 1/2 đến 2/3 chiều cao trụ pin. Giải pháp này giúp triệt tiêu tới 67,0% ứng suất kéo tại chân trụ và giảm chuyển vị ngang xuống dưới mức 5,5 cm. Thời gian triển khai thiết kế và thi công bổ sung dầm chống nên hoàn thành trong vòng 6 đến 12 tháng đối với các công trình đang xây dựng hoặc đang vận hành có nguy cơ chấn động địa chấn cao.

  3. Nâng cao phẩm cấp vật liệu bê tông và mật độ cốt thép chân trụ: Các nhà thầu xây lắp cần sử dụng bê tông cường độ cao từ mác M350 đến M400 cho toàn bộ khối trụ pin đập tràn cao trên 15 m, đồng thời tăng cường mật độ cốt thép chịu uốn cục bộ với hàm lượng tối thiểu mu_min đạt từ 0,5% đến 0,8% tại đoạn 3 m tính từ chân ngàm. Giải pháp này giúp nâng cao cường độ chịu kéo giới hạn của bê tông lên trên 1,5 MPa, đáp ứng điều kiện chống nứt bền vững dưới tổ hợp tải trọng đặc biệt.

  4. Lắp đặt hệ thống quan trắc chấn động và chuyển vị tự động: Đơn vị quản lý vận hành các nhà máy thủy điện lớn cần trang bị hệ thống cảm biến đo gia tốc nền 3 chiều và máy đo chuyển vị laser tự động tại đỉnh trụ pin với tần số lấy mẫu 100 Hz. Hệ thống phải được hoàn thiện lắp đặt và hiệu chuẩn trong vòng 3 tháng trước thời điểm tích nước hồ chứa, nhằm cung cấp dữ liệu cảnh báo sớm nguy cơ kẹt cửa van khi xuất hiện động đất tự nhiên hoặc động đất kích thích.

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

Luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

  1. Kỹ sư kết cấu và chuyên gia thiết kế công trình thủy điện, thủy lợi: Luận văn cung cấp quy trình mô phỏng phần tử hữu hạn 3D trên ANSYS và bảng tra cứu nhanh độ mảnh cho phép lambda_cp ứng với từng cấp mác bê tông, giúp rút ngắn thời gian tính toán và tối ưu hóa kích thước hình học thân trụ pin đập tràn cao.

  2. Đơn vị tư vấn thẩm tra và cơ quan quản lý an toàn hồ đập: Cung cấp hệ thống tiêu chí định lượng về ứng suất kéo giới hạn và chuyển vị cho phép để đánh giá tính chuẩn xác của các hồ sơ thiết kế kỹ thuật, phục vụ công tác nghiệm thu công trình xây dựng tại các khu vực nhạy cảm địa chấn.

  3. Ban quản lý dự án và các tập đoàn năng lượng: Nắm vững các cơ sở khoa học của giải pháp bố trí dầm chống ngang để chủ động lập phương án gia cường, cải tạo cho các công trình đập tràn hiện hữu đang gặp sự cố nứt chân trụ hoặc rung lắc kẹt cửa van xả lũ.

  4. Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Xây dựng Công trình thủy: Sử dụng tài liệu như một cẩm nang nghiên cứu chuyên sâu về phương pháp phân tích động lực học kết cấu bê tông khối lớn tương tác với nền đá và tải trọng địa chấn.

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

  1. Tại sao độ mảnh của trụ pin đập tràn lại quyết định trực tiếp đến nguy cơ kẹt cửa van khi xảy ra động đất? Khi tỷ lệ độ mảnh lambda = H/D vượt quá 6,0, độ cứng chống uốn ngang của trụ pin giảm mạnh. Dưới tác động của lực quán tính động đất cấp 8 hoặc 9, chuyển vị ngang đỉnh trụ có thể vượt quá 6,6 cm. Dao động lệch phương này làm biến dạng hình học của khe van và gối đỡ cơ khí, gây kẹt cứng cửa van cung hoặc van phẳng và làm tê liệt khả năng điều tiết xả lũ của hồ chứa.

  2. Việc bổ sung dầm chống ngang giữa các trụ pin mang lại hiệu quả cơ học cụ thể như thế nào? Dầm chống ngang biến các trụ pin đơn lẻ thành hệ khung không gian liên hợp có độ cứng uốn ngang tăng vọt. Kết quả mô phỏng tại Thủy điện Lai Châu chứng minh dầm chống ngang giúp giảm 19,7% chuyển vị ngang đỉnh trụ (từ 6,6 cm xuống 5,3 cm) và giảm tới 67,0% ứng suất kéo uốn tại chân trụ (từ 540 T/m2 xuống 178 T/m2), đưa hệ số an toàn chống nứt đạt yêu cầu tiêu chuẩn.

  3. Động đất kích thích tại các hồ chứa thủy điện lớn như Sông Tranh 2 ảnh hưởng ra sao đến trụ pin? Động đất kích thích phát sinh do áp lực của khối nước hồ chứa tích sâu hàng trăm mét tác động lên các đứt gãy địa chất. Dù cường độ thường ở mức 4,5 đến 5,5 độ Richter, nhưng tần suất rung chấn liên tục trong 2 đến 5 năm đầu tích nước sẽ làm gia tăng hiện tượng mỏi kết cấu, mở rộng các khe biến dạng từ 2 mm đến 6 mm và gây nứt cục bộ chân trụ pin nếu không được giằng chống ngang.

  4. Phương pháp phổ phản ứng có ưu điểm gì vượt trội so với phương pháp tĩnh lực học truyền thống? Phương pháp tĩnh lực học xem gia tốc rung động tại mọi điểm trên công trình là như nhau và bỏ qua tính chất đàn hồi. Ngược lại, phương pháp phổ phản ứng đưa vào hệ số động lực alpha_max = 2,0 và chu kỳ dao động riêng của từng mode, phản ánh chính xác hiện tượng khuếch đại dao động uốn ở đỉnh trụ và sự tập trung ứng suất kéo ở chân ngàm.

  5. Khi thiết kế trụ pin đập tràn cao 25 m tại vùng địa chấn cấp 8, nên lựa chọn mác bê tông và bề dày trụ bao nhiêu? Theo bảng tra cứu độ mảnh của luận văn, tại vùng động đất cấp 8, công trình nên sử dụng bê tông mác tối thiểu M350 hoặc M400. Với chiều cao H = 25 m, độ mảnh cho phép khống chế ở mức lambda <= 5,0, tương ứng bề dày trụ pin D tối thiểu phải đạt từ 5,0 m trở lên kết hợp bố trí dầm giằng ngang để đảm bảo ứng suất kéo không vượt quá cường độ chịu kéo tiêu chuẩn 1,45 MPa.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán cấp thiết về đánh giá ổn định hướng ngang của trụ pin đập tràn cao từ 15 m đến 30 m chịu tác động của tải trọng động đất.
  • Xây dựng thành công quy trình mô phỏng số phần tử hữu hạn 3D trên ANSYS 16.0 với 27 kịch bản không gian, làm sáng tỏ quy luật phân bố chuyển vị ngang và ứng suất kéo uốn tập trung tại chân trụ.
  • Chứng minh bằng số liệu thực chứng tại Thủy điện Lai Châu rằng giải pháp dầm chống ngang giúp giảm 19,7% chuyển vị đỉnh và triệt tiêu 67,0% ứng suất kéo chân trụ, bảo đảm an toàn vận hành cửa van xả lũ.
  • Thiết lập hệ thống bảng tra cứu độ mảnh cho phép lambda_cp tương ứng với các cấp mác bê tông M300, M350, M400 và các cấp động đất 7, 8, 9 phục vụ trực tiếp công tác thiết kế thực tế.
  • Khuyến nghị các đơn vị tư vấn và quản lý hồ đập nhanh chóng chuẩn hóa quy chuẩn kháng chấn cho trụ pin, tích hợp hệ giằng ngang và hệ thống quan trắc tự động trong giai đoạn tới. Hãy liên hệ các chuyên gia công trình thủy để nhận bộ file mô hình tính toán chi tiết và áp dụng giải pháp bảo vệ an toàn đập tràn ngay hôm nay.