Tổng quan nghiên cứu

Trong hệ thống hạ tầng năng lượng và thủy lợi toàn cầu, đập vật liệu địa phương chiếm hơn 70% tổng số các công trình đập dâng nước đã được xây dựng nhờ khả năng tận dụng tối đa nguồn vật liệu tại chỗ và thích ứng tốt với nhiều điều kiện địa chất phức tạp. Trên thế giới, nhiều công trình đập đất đá đã đạt chiều cao kỷ lục như đập Rogun cao 335 mét, đập Nurek cao 315 mét, và tại Việt Nam có đập thủy điện Hòa Bình cao 128 mét. Tuy nhiên, việc đánh giá độ an toàn của các công trình đập dâng có chiều cao lớn đòi hỏi các phương pháp tính toán cơ học đất ngày càng chuẩn xác.

Vấn đề cốt lõi đặt ra trong thiết kế đập đất đá hiện nay là các phương pháp giải tích truyền thống hoặc mô hình đàn hồi tuyến tính một chiều thường đơn giản hóa quá mức các ứng xử cơ học của vật liệu. Đất đá trên thực tế thể hiện tính phi tuyến rất mạnh ngay từ khi chịu cấp tải trọng nhỏ. Việc áp dụng mô hình tuyến tính dễ dẫn đến đánh giá sai lệch trường ứng suất, bỏ qua hiện tượng tập trung ứng suất kéo hoặc biến dạng cục bộ, từ đó làm gia tăng nguy cơ nứt gãy thủy lực và lún không đều trong thân đập.

Mục tiêu trọng tâm của nghiên cứu là phân tích chuyên sâu trạng thái ứng suất - biến dạng của đập vật liệu địa phương bằng mô hình phi tuyến đàn dẻo tăng cứng, áp dụng trực tiếp cho dự án đập thủy điện Khao Mang Thượng tại tỉnh Yên Bái. Công trình có công suất lắp máy 24.5 MW, chiều cao đập lớn nhất đạt 68.0 mét trên nền đá cứng và dung tích hồ chứa tương ứng mực nước dâng bình thường 915.00 mét là 13.37 triệu mét khối.

Phạm vi nghiên cứu bao trùm toàn bộ các mặt cắt ngang xung yếu và mặt cắt dọc trục đập trong các giai đoạn thi công đắp lớp từng phần cũng như giai đoạn tích nước vận hành đầu tiên. Kết quả tính toán đem lại giá trị thực tiễn cao, giúp dự báo chính xác độ lún, xác định vị trí tiềm ẩn nguy cơ nứt nẻ và tối ưu hóa giải pháp khối gia tải hạ lưu, nâng cao hệ số an toàn công trình theo Quy chuẩn kỹ thuật Quốc gia QCVN 04-05:2011/BNNPTNT.

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 nền tảng cơ học đất phi tuyến hiện đại, kết hợp đồng thời lý thuyết ứng suất hiệu quả của Terzaghi và lý thuyết đàn dẻo tăng cứng của Schanz.

Theo nguyên lý Terzaghi, ứng suất tổng trong khối đất đá bão hòa được phân tách thành ứng suất hiệu quả chịu lực bởi khung hạt khoáng và áp lực nước lỗ rỗng sinh ra trong các khe rỗng. Để mô tả chính xác sự thay đổi mô đun biến dạng theo độ sâu và trạng thái tải trọng, mô hình đất tăng cứng Hardening Soil được sử dụng thay thế cho mô hình đàn dẻo lý tưởng Mohr-Coulomb hay mô hình Hyperbolic của Duncan-Chang.

Mô hình Hardening Soil tích hợp hai cơ chế biến dạng dẻo độc lập: tăng cứng cắt mô tả biến dạng trượt không thuận nghịch do độ lệch ứng suất gây ra, và tăng cứng nén mô tả độ nén chặt vĩnh cửu của cấu trúc hạt dưới tác dụng của ứng suất nén đẳng hướng hoặc nén cố kết một trục. Độ cứng của đất đá biến thiên phi tuyến theo áp lực buồng hiệu quả thông qua hàm số mũ với số mũ biến dạng dao động trong khoảng từ 0.5 đến 1.0.

Các khái niệm cơ bản được chuẩn hóa trong nghiên cứu bao gồm: mô đun cát tuyến tham chiếu ứng với 50% độ lệch ứng suất cực đại, mô đun tiếp tuyến cố kết một trục tham chiếu, mô đun dỡ tải và nạp tải lại tham chiếu (thường có giá trị xấp xỉ gấp 3 lần mô đun cát tuyến), cùng hệ số phá hoại tiếp tuyến luôn nhỏ hơn 1.0. Bên cạnh đó, hiện tượng nứt gãy thủy lực được đánh giá qua điều kiện so sánh giữa ứng suất nén hiệu quả chính nhỏ nhất với áp lực nước lỗ rỗng tại từng cao trình.

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

Nguồn dữ liệu đầu vào của luận văn được tổng hợp từ hồ sơ khảo sát địa chất công trình và kết quả thí nghiệm cơ lý đất đá phục vụ thiết kế kỹ thuật của dự án thủy điện Khao Mang Thượng.

Quy mô mẫu thí nghiệm bao gồm hệ thống thí nghiệm nén 3 trục cố kết thoát nước và không thoát nước trên các mẫu đất sét chống thấm đầm nén, kết hợp thí nghiệm nén cố kết một trục Oedometer để xác định chính xác các mô đun biến dạng tham chiếu và các chỉ tiêu cường độ chống cắt. Đối với vật liệu đá đắp thượng hạ lưu có đường kính hạt thực tế ngoài hiện trường lên đến 1200 mm, phương pháp đường cấp phối song song với cỡ hạt khống chế tối đa 50 mm trong phòng thí nghiệm đã được áp dụng nhằm suy giảm kích thước mẫu nhưng vẫn bảo toàn trọn vẹn đặc trưng ma sát góc trong và độ cứng khung hạt.

Phương pháp phân tích chủ đạo là mô hình hóa số bằng phương pháp phần tử hữu hạn hai chiều và ba chiều phi tuyến. Lý do lựa chọn phương pháp này xuất phát từ khả năng mô phỏng liên tục sự thay đổi hình học và trạng thái ứng suất theo từng giai đoạn đắp đập thực tế. Toàn bộ tiến trình xây dựng được chia thành hơn 8 cấp đắp tải trọng theo sơ đồ kích hoạt phần tử từng lớp, kết hợp quy trình Nobari-Duncan gồm 3 bước với hơn 4 cấp dâng nước nhằm mô phỏng hiện tượng mềm hóa vật liệu khi đất đá chuyển từ trạng thái khô sang bão hòa nước.

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ố trị theo mô hình Hardening Soil trên các mặt cắt đặc trưng của đập Khao Mang Thượng đã mang lại nhiều phát hiện quan trọng:

Thứ nhất, trường ứng suất hiệu quả theo phương đứng và phương ngang trong thân đập phân bố phi tuyến và thể hiện rõ hiệu ứng vòm chuyển tải. Do lõi sét có tính biến dạng cao hơn các khối đá đắp hai bên, ứng suất tiếp xúc tại mặt phân chia vật liệu đã truyền bớt một phần tải trọng từ lõi sang khối đá thượng và hạ lưu. Ứng suất nén hiệu quả lớn nhất trong lõi sét ở cao trình 904.00 mét duy trì ổn định, loại trừ hoàn toàn sự xuất hiện của vùng ứng suất kéo gây mất ổn định cục bộ.

Thứ hai, quy luật phân bố chuyển vị đứng trong giai đoạn thi công phản ánh đúng bản chất vật lý khi đắp đập nhiều lớp. Chuyển vị lún lớn nhất không xuất hiện tại đỉnh đập như các tính toán đàn hồi tuyến tính tức thời quy ước, mà tập trung ở khu vực giữa thân đập từ cao độ 880.00 mét đến 895.00 mét. Trên mặt cắt lòng suối 0+100, độ lún tích lũy lớn nhất đạt giá trị nằm trong ngưỡng an toàn cho phép của kết cấu đập đá đổ lõi giữa.

Thứ ba, nghiên cứu chứng minh vai trò then chốt của khối đá gia tải chân đập hạ lưu tại mặt cắt xung yếu 0+75. Khi không có khối gia tải, chuyển vị ngang lớn nhất ghi nhận lên tới 14.8 cm tại sườn dốc; khi bố trí khối gia tải hạ lưu hợp lý, chuyển vị ngang cực đại giảm mạnh xuống chỉ còn 5.0 cm, tương đương mức giảm chênh lệch biến dạng hơn 66.2%.

Thứ tư, khi hồ chứa tích nước đến mực nước dâng bình thường 915.00 mét, dưới tác dụng của áp lực thấm và hiện tượng bão hòa nước, thân đập có xu hướng chuyển vị dịch nhẹ về phía hạ lưu. Tuy nhiên, hệ số an toàn chống nứt gãy thủy lực theo mặt đứng và mặt ngang tại mọi vị trí trong lõi sét đều đạt giá trị lớn hơn 1.0, khẳng định tuyến đập không xảy ra rủi ro nứt nẻ xuyên tâm.

Thảo luận kết quả

Các kết quả phân tích số liệu có thể được hệ thống hóa rõ ràng thông qua bảng tổng hợp mô đun đàn hồi trung bình của từng khối đắp và biểu đồ đường cong đẳng trị ứng suất chính:

Trong giai đoạn đắp khô, mô đun biến dạng trung bình của khối đá đắp đạt giá trị cao giúp giữ vững khung chịu lực. Khi chuyển sang giai đoạn tích nước bình thường, do hiện tượng ngâm nước làm suy giảm ma sát liên kết giữa các hạt đá, mô đun đàn hồi của vùng bão hòa giảm khoảng 15% đến 25% so với trạng thái ban đầu. Điều này giải thích nguyên nhân làm gia tăng độ lún bổ sung sau khi tích nước.

So sánh với các nghiên cứu tính toán theo mô hình tuyến tính đơn giản trước đây, việc áp dụng mô hình tăng cứng phi tuyến giúp triệt tiêu sai số dự báo chuyển vị lún lên tới 20% đến 30%. Các biểu đồ quan hệ ứng suất - biến dạng trích xuất từ mô hình phần tử hữu hạn hoàn toàn trùng khớp với quy luật biến dạng cố kết và nén hông thực tế, chứng minh tính ưu việt của việc mô phỏng quá trình thi công đắp lớp từng bước so với phương pháp chất tải tức thời một lần.

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

Dựa trên các phân tích chuyên sâu về ứng suất và chuyển vị của đập thủy điện Khao Mang Thượng, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm đảm bảo an toàn tuyệt đối cho công trình:

Thứ nhất, kiểm soát chặt chẽ độ ẩm và công tác đầm nén vật liệu lõi sét tại hiện trường. Đơn vị thi công và Tư vấn giám sát cần duy trì độ ẩm thi công trong giới hạn sai số không vượt quá 2% so với độ ẩm đầm nén tối ưu, đồng thời đảm bảo độ chặt đầm nén đạt từ 95% đến 100% dung trọng khô lớn nhất theo tiêu chuẩn Proctor tiêu chuẩn trong suốt thời gian thi công đắp đập.

Thứ hai, thi công hoàn thiện khối gia tải hạ lưu tại các mặt cắt sườn dốc hẹp như mặt cắt 0+75 theo đúng thiết kế trước khi tích nước. Ban Quản lý dự án và Nhà thầu xây lắp phải ưu tiên bố trí khối đá gia tải nhằm khống chế chuyển vị ngang dưới mức 5.0 cm, ngăn ngừa triệt để nguy cơ trượt mái đập và nứt tách tại vị trí tiếp giáp vai đập có độ dốc thay đổi đột ngột.

Thứ ba, thiết lập mạng lưới thiết bị quan trắc địa kỹ thuật đa chỉ tiêu gồm tối thiểu 8 mặt cắt quan trắc lún, chuyển vị ngang và áp lực nước lỗ rỗng tự động trong thân đập và nền móng. Chủ đầu tư phối hợp với Đơn vị tư vấn địa kỹ thuật tiến hành phân tích ngược từ các chu kỳ đo đạc thực tế trong 2 năm đầu vận hành để hiệu chỉnh chính xác các tham số cơ học của mô hình với độ tin cậy trên 90%.

Thứ tư, ban hành quy trình tích nước hồ chứa ban đầu theo lộ trình tối thiểu 4 giai đoạn kiểm soát mực nước dâng. Hội đồng nghiệm thu và Đơn vị quản lý vận hành cần duy trì tốc độ dâng nước không vượt quá 0.5 đến 1.0 mét mỗi ngày nhằm cho phép áp lực nước lỗ rỗng trong lõi sét tiêu tán nhịp nhàng, triệt tiêu nguy cơ phát sinh hiện tượng nứt gãy thủy lực do sốc áp suất.

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

Nội dung và phương pháp luận của luận văn mang giá trị ứng dụng thực tiễn và học thuật phong phú cho nhiều nhóm đối tượng chuyên ngành:

Nhóm thứ nhất là các kỹ sư thiết kế công trình thủy lợi và thủy điện. Tài liệu cung cấp quy trình chi tiết về việc ứng dụng mô hình phi tuyến Hardening Soil, cách thức xác định các tham số tham chiếu từ thí nghiệm phòng và phương pháp thiết lập bài toán phần tử hữu hạn phân tích ứng suất - biến dạng cho các đập dâng cao trên 50 mét.

Nhóm thứ hai là các chuyên gia tư vấn giám sát và kỹ sư công nghệ thi công. Luận văn cung cấp cơ sở khoa học để đánh giá mức độ ảnh hưởng của trình tự đắp lớp, tầm quan trọng của hệ số đầm nén và giải pháp bố trí khối gia tải chân đập nhằm kiểm soát độ lún lệch trong quá trình xây lắp.

Nhóm thứ ba là các ban quản lý dự án, chủ đầu tư các nhà máy thủy điện vừa và nhỏ (công suất từ 10 MW đến 30 MW). Luận văn hỗ trợ xây dựng kế hoạch tích nước an toàn, tối ưu hóa khối lượng đắp đập và thiết kế mạng lưới quan trắc biến dạng công trình đạt chuẩn kỹ thuật quốc gia.

Nhóm thứ tư là các giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Xây dựng công trình thủy, Địa kỹ thuật công trình và Cơ học đất. Đây là tài liệu tham khảo chất lượng cao về việc kết hợp giữa lý thuyết đàn dẻo phi tuyến và bài toán ứng dụng thực tế trên một công trình thủy điện trọng điểm tại vùng núi Tây Bắc.

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

Mô hình Hardening Soil có ưu điểm vượt trội gì so với mô hình Mohr-Coulomb khi tính toán đập đất đá? Mô hình Hardening Soil phản ánh chính xác sự phụ thuộc phi tuyến của độ cứng vật liệu vào trạng thái ứng suất thông qua hàm số mũ, đồng thời phân tách rõ ràng giữa biến dạng dẻo tăng cứng cắt và tăng cứng nén. Trong khi đó, mô hình Mohr-Coulomb chỉ sử dụng một mô đun đàn hồi không đổi, dễ dẫn đến ước tính sai lệch độ lún tới hơn 25% trong các đập đá đổ có chiều cao lớn trên 50 mét.

Khối gia tải hạ lưu đóng vai trò quan trọng như thế nào tại mặt cắt 0+75 của đập Khao Mang Thượng? Khối gia tải hạ lưu có tác dụng giữ ổn định mái dốc và cân bằng trường ứng suất tại các vị trí địa hình sườn dốc không đối xứng. Tại mặt cắt 0+75, việc đắp thêm khối gia tải đã giúp giảm chuyển vị ngang từ 14.8 cm xuống còn 5.0 cm, triệt tiêu hơn 66% biên độ lún lệch và loại bỏ hoàn toàn nguy cơ nứt tách thân đập tiếp giáp nền.

Nguyên nhân cốt lõi gây ra hiện tượng nứt gãy thủy lực trong lõi sét chống thấm là gì? Nứt gãy thủy lực xuất hiện khi áp lực nước lỗ rỗng hoặc áp lực thủy tĩnh tại một điểm vượt quá ứng suất nén hiệu quả chính nhỏ nhất của khối đất tại vị trí đó. Khi cột nước dâng cao trong giai đoạn tích nước, nếu vùng lõi sét bị suy giảm ứng suất do hiệu ứng vòm treo tải, áp lực dòng thấm chỉ cần chênh lệch từ 2 đến 4 mét cột nước cũng có thể xé rách khe nứt và gây xói ngầm nghiêm trọng.

Phương pháp đường cấp phối song song giải quyết khó khăn gì trong thí nghiệm vật liệu đá đắp? Vật liệu đá đắp ngoài hiện trường thường có kích thước hạt cực đại lên tới 1200 mm, vượt quá dung tích buồng nén 3 trục tiêu chuẩn trong phòng thí nghiệm. Phương pháp đường cấp phối song song tạo ra các mẫu thí nghiệm có cỡ hạt khống chế dưới 50 mm với đường cong thành phần hạt đồng dạng, giúp xác định chuẩn xác góc ma sát trong và các mô đun biến dạng tương đương của khối đá đắp.

Tại sao phải chia quá trình đắp đập thành hơn 8 bước tính toán trong phần mềm phần tử hữu hạn? Chia quá trình thi công thành nhiều bước giúp mô phỏng chính xác sự tích lũy biến dạng và quá trình hình thành trường ứng suất tự nhiên theo từng lớp đắp. Nếu tính toán toàn bộ khối đập được đắp tức thời trong một bước, chuyển vị lớn nhất sẽ bị gán sai lên đỉnh đập, trong khi thực tế chuyển vị lún cực đại luôn phân bố ở khoảng giữa thân đập từ 1/2 đến 2/3 chiều cao công trình.

Kết luận

Luận văn đã giải quyết toàn diện bài toán phân tích trạng thái ứng suất - biến dạng của đập vật liệu địa phương theo mô hình cơ học đất phi tuyến hiện đại với 5 kết quả nổi bật:

  • Chứng minh tính ưu việt vượt trội của mô hình đàn dẻo tăng cứng Hardening Soil so với các mô hình tuyến tính cổ điển trong việc mô phỏng ứng xử phi tuyến của đập đất đá.
  • Xác định chính xác quy luật phân bố trường ứng suất hiệu quả và biến dạng tích lũy theo từng giai đoạn thi công đắp lớp cho đập thủy điện Khao Mang Thượng cao 68.0 mét.
  • Khẳng định hiệu quả kỹ thuật của khối gia tải hạ lưu giúp cắt giảm tới 66.2% chuyển vị ngang tại các mặt cắt địa hình xung yếu, duy trì độ ổn định tổng thể cho tuyến đập.
  • Đánh giá an toàn điều kiện chống nứt gãy thủy lực trong lõi sét khi hồ chứa tích nước đến cao trình bình thường 915.00 mét, đảm bảo công trình vận hành an toàn tuyệt đối.
  • Đề xuất quy trình đồng bộ từ công tác đầm nén hiện trường, lộ trình dâng nước 4 giai đoạn đến thiết lập mạng lưới quan trắc địa kỹ thuật phục vụ phân tích ngược.

Đóng góp lớn nhất của nghiên cứu là xây dựng thành công quy trình tính toán số trị chuẩn xác, có khả năng chuyển giao và ứng dụng trực tiếp cho các dự án đập vật liệu địa phương có điều kiện địa hình phức tạp tại Việt Nam. Bước tiếp theo cần triển khai là tích hợp các số liệu quan trắc thực tế sau 12 đến 24 tháng vận hành để hoàn thiện bộ thông số cơ lý chuẩn hóa cho các dạng đập đá đổ lõi sét. Các đơn vị tư vấn và ban quản lý dự án thủy điện có thể áp dụng ngay khung phương pháp luận này để tối ưu hóa thiết kế kỹ thuật, tiết kiệm chi phí xây dựng và nâng cao độ an toàn bền vững cho công trình.