Tổng quan nghiên cứu

Thủy điện đóng vai trò then chốt trong cơ cấu năng lượng quốc gia khi cung cấp gần 40% sản lượng điện năng và chiếm khoảng 50% tổng công suất toàn hệ thống điện Việt Nam với quy mô lắp máy đạt xấp xỉ 27.000 MW. Mạng lưới sông ngòi dày đặc gồm 2.360 con sông suối có chiều dài từ 10 km trở lên mang lại nguồn trữ năng lý thuyết dồi dào ước tính khoảng 310 tỷ kWh mỗi năm. Mặc dù tiềm năng rất lớn, các dự án thủy điện vừa và nhỏ có công suất dưới 30 MW mới chỉ khai thác được khoảng 20% công suất hiện hữu, để lại dư địa phát triển lên tới 80% cần được quy hoạch và đầu tư hiệu quả.

Vấn đề cốt lõi trong việc khai thác thủy điện tại các vùng núi dốc địa hình hiểm trở thuộc khu vực Tây Bắc là lựa chọn tuyến dẫn nước đảm bảo an toàn cơ học, tối ưu thủy lực và giảm thiểu chi phí đầu tư. Việc xây dựng tuyến kênh hở thường gặp nguy cơ sạt trượt sườn dốc và bồi lắng lòng dẫn trong mùa mưa lũ. Do đó, giải pháp đường hầm dẫn nước có áp qua lòng núi trở thành lựa chọn kỹ thuật tối ưu. Đề tài tập trung nghiên cứu kết cấu hợp lý đường hầm dẫn nước vào nhà máy thủy điện Mường Kim II thuộc địa bàn huyện Than Uyên, tỉnh Lai Châu. Dự án có công suất thiết kế 10,5 MW, điện lượng bình quân hàng năm đạt 41,87 triệu kWh, lưu lượng thiết kế 27,11 m3/s, cột nước lớn nhất vượt trên 50 m và tổng chiều dài tuyến hầm khoảng 1.600 m.

Mục tiêu cụ thể của nghiên cứu là phân tích điều kiện địa chất công trình, đánh giá trạng thái ứng suất - biến dạng của hệ tương tác vỏ hầm và khối đá núi xung quanh, từ đó đề xuất hình thức mặt cắt ngang tối ưu. Ý nghĩa thực tiễn của công trình thể hiện qua việc giảm 15% đến 20% tổn thất cột nước cục bộ, tối ưu hóa khối lượng bê tông cốt thép lót vỏ hầm và đảm bảo an toàn tuyệt đối cho toàn bộ tuyến năng lượng trong suốt vòng đời vận hành trên 50 năm.

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 bộ hai khung lý thuyết cơ bản trong kỹ thuật công trình ngầm thủy công: lý thuyết cân bằng giới hạn vòm sụt lở tự nhiên và lý thuyết tương tác đàn hồi - dẻo giữa vỏ hầm với môi trường đất đá.

Lý thuyết vòm cân bằng tự nhiên của nhà khoa học Protodyakonov được sử dụng để xác định áp lực địa tầng tác dụng lên vỏ hầm. Theo nguyên lý này, khối đá xung quanh hầm được coi là một thể rời quy ước đặc trưng bởi hệ số kiên cố fk. Kích thước nhịp vòm cân bằng tự nhiên L và chiều cao vòm sụp lở h được xác định dựa trên khẩu độ đào B0, chiều cao công trình H0 và góc ma sát trong của tầng đá. Bên cạnh đó, nghiên cứu tích hợp lý thuyết biến dạng đàn hồi môi trường liên tục của B.G. Galerkin và mô hình nền đàn hồi Winkler để tính toán lực kháng đàn tính của khối đá bao quanh.

Khung phân loại địa chất ứng dụng hệ thống đánh giá chất lượng khối đá RMR (Rock Mass Rating) kết hợp chỉ số cường độ địa chất GSI theo tiêu chuẩn Hoek-Brown. Ba khái niệm chuyên ngành then chốt chi phối quá trình tính toán bao gồm: áp lực địa tầng tự nhiên thẳng đứng và nằm ngang; lực kháng đàn tính bị động của đá núi (đặc trưng qua hệ số phản lực đơn vị K0 dao động từ 200 đến 1.600 kG/cm2); và trạng thái ứng suất phẳng trong bài toán hình trụ thành dày chịu tải trọng thủy tĩnh lẫn địa kỹ thuật.

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

Nguồn dữ liệu đầu vào của luận văn được thu thập từ tài liệu khảo sát địa chất công trình thực tế dọc theo tuyến hầm dài 1.600 m của thủy điện Mường Kim II. Dữ liệu khảo sát bao gồm 12 lỗ khoan thăm dò địa chất sâu từ 35 m đến 120 m, kết hợp các thí nghiệm cơ lý mẫu đá tại hiện trường và trong phòng thí nghiệm. Phương pháp lấy mẫu tuân thủ nguyên tắc phân tầng địa chất theo hai đới thạch học chủ yếu: đá phun trào trachit thuộc phức hệ Nậm Kim và đá trầm tích gắn kết thuộc hệ tầng Mường Trai.

Nghiên cứu áp dụng phương pháp phân tích cơ học kết cấu giải tích cổ điển (mô hình vòm cao, vòm thấp, vòm khép kín trên nền đàn hồi) kết hợp phương pháp số phần tử hữu hạn (FEM - Finite Element Method). Lý do lựa chọn kết hợp phương pháp số và giải tích nhằm giải quyết bài toán tiếp xúc phi tuyến phức tạp, đồng thời kiểm chứng chéo độ tin cậy của các kết quả nội lực. Mô hình tính toán kiểm tra độ bền và độ ổn định theo hai trạng thái giới hạn với hệ số lệch tải nghiêm ngặt: hệ số lệch tải trọng lượng đá trong vòm cân bằng n = 1,5; áp lực ngang n = 1,2; áp lực nước trong và hiện tượng nước va n = 1,0; áp lực nước ngầm n = 1,1; và áp lực phụt vữa n = 1,2. Toàn bộ quá trình nghiên cứu, xử lý số liệu địa chất và mô phỏng số được thực hiện liên tục trong chu kỳ 24 tháng.

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

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

Quá trình tính toán giải tích và mô phỏng số đã mang lại bốn phát hiện khoa học và kỹ thuật quan trọng:

Thứ nhất, điều kiện địa chất dọc tuyến hầm có sự phân hóa rõ rệt về khả năng chịu lực. Đới đá trachit tươi và phong hóa nhẹ có hệ số kiên cố fk đạt từ 6 đến 10, mô đun kẽ nứt thấp dưới 1,5 kẽ nứt/mét dài, cho phép đá tự tạo vòm cân bằng tốt và áp lực địa tầng tác dụng lên đỉnh hầm chỉ dao động trong khoảng 0,05 đến 0,10 MPa. Ngược lại, đới đá trầm tích nứt nẻ mạnh của hệ tầng Mường Trai có hệ số kiên cố fk giảm xuống dưới 4, làm chiều cao vòm sụp lở tăng lên và tạo áp lực hông đáng kể lên vách hầm từ 0,20 đến 0,30 MPa.

Thứ hai, kết quả so sánh giữa các dạng hình học mặt cắt cho thấy: đối với đường hầm có áp chịu cột nước làm việc trên 50 m, mặt cắt hình tròn thể hiện ưu thế vượt trội về mặt chịu lực so với mặt cắt dạng móng ngựa hay vòm đáy bằng. Mặt cắt tròn giúp triệt tiêu hoàn toàn các điểm tập trung ứng suất mômen uốn cục bộ tại các góc chân tường, phân phối đều ứng suất tiếp tuyến dọc theo toàn bộ chu vi vỏ hầm, giúp giảm giá trị mômen uốn cực đại từ 42 kNm/m ở mặt cắt vòm xuống dưới 12 kNm/m ở mặt cắt tròn.

Thứ ba, sự tham gia của lực kháng đàn tính từ khối đá đóng vai trò quyết định trong việc giảm ứng suất kéo của bê tông vỏ hầm. Khi xét đến hệ số phản lực nền K0 của đá trachit (đạt từ 800 đến 1.400 kG/cm2), ứng suất kéo lớn nhất ở thớ trong của vỏ hầm giảm đi khoảng 35% đến 40% so với trường hợp tính toán kết cấu vỏ hầm làm việc độc lập không có nền đàn hồi.

Thứ tư, phân tích các tổ hợp tải trọng xác định được trường hợp bất lợi nhất đối với vỏ hầm là trạng thái xả cạn kiệt hầm để phục vụ công tác kiểm tra, duy tu định kỳ. Lúc này, áp lực nước bên trong bằng 0 trong khi áp lực đá núi ngoài đạt giá trị tối đa với hệ số lệch tải n = 1,5 kết hợp áp lực thủy tĩnh của nước ngầm n = 1,1, tạo ra trạng thái nén lệch tâm uốn lớn tại khu vực hông và đáy hầm.

Thảo luận kết quả

Nguyên nhân cơ bản tạo nên sự khác biệt về ứng suất là do cơ chế tương tác đàn tính giữa kết cấu vỏ mỏng và môi trường vô hạn. Khi hầm vận hành đầy nước dưới áp lực cao, vỏ hầm có xu hướng biến dạng giãn nở hướng tâm ra phía ngoài, ép chặt vào khối đá nguyên khối. Khối đá xung quanh phản ứng lại bằng lực kháng đàn tính thụ động, giữ vai trò như một vòng đai gia cường tự nhiên chịu phần lớn lực kéo thay cho bê tông.

Khi so sánh với các công trình thủy điện lớn đã xây dựng như Nậm Chiến (hầm dài 10 km, chịu cột nước cao) hay Huội Quảng (7 tuyến hầm áp lực dài hơn 4 km), kết quả của nghiên cứu này hoàn toàn tương thích về mặt phân bố trường ứng suất tiếp tuyến và ứng suất pháp. Sự tương đồng này khẳng định việc áp dụng lý thuyết nền đàn hồi Winkler và lời giải Galerkin cho đường hầm thủy điện Mường Kim II là hoàn toàn chính xác và tin cậy.

Toàn bộ tập hợp dữ liệu tính toán được trực quan hóa thông qua biểu đồ quan hệ phi tuyến giữa hệ số kiên cố fk và lực kháng đàn tính đơn vị K0, kết hợp biểu đồ bao nội lực mômen uốn và lực dọc tại 24 tiết diện kiểm tra trên chu vi hầm. Đồng thời, bảng tổng hợp phân bố ứng suất chính σ1 và σ3 được thiết lập chi tiết cho từng trường hợp tải trọng, giúp kỹ sư thiết kế xác định chuẩn xác các vị trí xuất hiện ứng suất kéo vượt quá cường độ chịu kéo tính toán của bê tông để bố trí mật độ cốt thép gia cường phù hợp.

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

Dựa trên các kết quả phân tích cơ học và thủy lực, bốn nhóm giải pháp cụ thể được đề xuất nhằm tối ưu hóa công trình:

Thứ nhất, tối ưu hóa thiết kế hình học mặt cắt ngang. Đơn vị tư vấn thiết kế cần áp dụng mặt cắt hình tròn đường kính lọt lòng 2,8 m đến 3,2 m cho toàn bộ các đoạn hầm có áp chịu cột nước trên 35 m đi qua đới địa chất phong hóa trung bình và nứt nẻ. Đối với các đoạn hầm đi qua vùng đá trachit nguyên khối kiên cố cao fk > 8 ở đầu tuyến năng lượng, có thể chuyển tiếp sang mặt cắt vòm đáy bằng để thuận lợi cho việc vận hành cơ giới hóa khi đào hầm, giúp tiết kiệm 15% khối lượng bê tông lót trong giai đoạn thiết kế bản vẽ thi công.

Thứ hai, gia cố chủ động khối đá bao quanh bằng công nghệ phụt vữa và neo đá. Nhà thầu thi công cần triển khai phụt vữa xi măng - bentonite áp lực cao từ 0,5 đến 1,0 MPa tại vùng đỉnh vòm và hai bên tường hầm trên các đoạn qua hệ tầng Mường Trai để điền đầy các kẽ nứt nẻ, nâng chỉ số RMR của khối đá lên thêm 20 đến 25 điểm trước khi đổ bê tông vỏ hầm. Thời gian thực hiện phải hoàn thành ngay sau chu kỳ đào từng gương hầm.

Thứ ba, chuẩn hóa quy trình kiểm soát vật liệu và chiều dày vỏ hầm. Đơn vị giám sát và thi công phải sử dụng bê tông mác M250 đến M300 có phụ gia chống thấm cấp B8, chiều dày vỏ hầm khống chế từ 25 cm đến 35 cm tùy theo từng phân đoạn địa chất, đồng thời duy trì hệ số lệch tải nghiêm ngặt từ 1,2 đến 1,5 cho toàn bộ các cấu kiện chịu lực trong thời gian thi công 18 tháng.

Thứ tư, thiết lập hệ thống quan trắc ứng suất - biến dạng tự động. Ban quản lý dự án cần đầu tư lắp đặt chuỗi cảm biến quan trắc bao gồm 8 trạm đo áp lực tiếp xúc giữa vỏ hầm và đá, 6 thiết bị đo áp lực nước ngầm kẽ rỗng và hệ thống đo độ mở rộng khe nứt dọc theo 1.600 m tuyến hầm. Dữ liệu quan trắc cần được truyền liên tục 24/7 về trung tâm điều hành nhà máy trong suốt quá trình tích nước thử nghiệm và vận hành thương mại.

Đố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 lại giá trị thực tiễn cho bốn nhóm đối tượng chuyên môn:

Nhóm thứ nhất là các kỹ sư tư vấn thiết kế công trình thủy lợi - thủy điện. Luận văn cung cấp quy trình mẫu về tính toán kết cấu hầm áp lực từ việc thiết lập sơ đồ tính, xác định tải trọng bất lợi theo quy chuẩn hiện hành, đến việc áp dụng mô hình giải tích và phương pháp phần tử hữu hạn để tối ưu kích thước kết cấu vỏ hầm.

Nhóm thứ hai là chủ đầu tư và các ban quản lý dự án năng lượng tái tạo. Báo cáo cung cấp căn cứ kỹ thuật rõ ràng để so chọn phương án tuyến dẫn nước, quản lý kiểm soát khối lượng đào đắp - bê tông, giảm thiểu rủi ro phát sinh kinh phí và dự báo các biến động địa chất công trình phức tạp tại vùng núi phía Bắc.

Nhóm thứ ba là các nhà thầu thi công công trình ngầm. Doanh nghiệp xây lắp có thể tham khảo trực tiếp các thông số kỹ thuật về áp lực phụt vữa chèn, sơ đồ cắm neo gia cố và quy trình xử lý chống thấm cho hầm dẫn nước đường kính nhỏ và vừa trong điều kiện địa chất đá nứt nẻ.

Nhóm thứ tư là giảng viên, học viên cao học và sinh viên chuyên ngành Xây dựng Công trình thủy, Thủy điện và Địa kỹ thuật công trình. Tài liệu là nguồn tham khảo học thuật giá trị, kết nối mạch lạc giữa lý thuyết cơ học môi trường liên tục kinh điển với bài toán thiết kế kỹ thuật thực tiễn.

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

Tại sao dự án thủy điện Mường Kim II lựa chọn đường hầm dẫn nước có áp thay vì đào kênh hở?

Khu vực huyện Than Uyên có địa hình đồi núi dốc đứng và bị chia cắt sâu. Nếu xây dựng kênh hở dài 1.600 m, khối lượng đào đắp sườn đồi sẽ rất lớn và dễ gây sạt lở vào mùa mưa lũ. Tuyến hầm ngầm giúp rút ngắn chiều dài dẫn nước, triệt tiêu nguy cơ đất đá bồi lấp như từng xảy ra ở bậc thang thủy điện Mường Kim phía trên năm 2012, đồng thời bảo đảm lưu lượng ổn định 27,11 m3/s.

Phương pháp xác định áp lực đá núi của Protodyakonov được áp dụng cụ thể như thế nào?

Phương pháp Protodyakonov áp dụng trực tiếp cho các đới đá nứt nẻ có hệ số kiên cố fk nhỏ hơn 4. Nghiên cứu xác định kích thước vòm sụp lở tự nhiên dựa trên nhịp hầm và góc ma sát trong của đá. Với hệ số lệch tải n = 1,5 cho trọng lượng khối đá vòm, phương pháp cho phép định lượng chính xác tải trọng thẳng đứng phân bố đều tác dụng trực tiếp lên đỉnh vòm hầm.

Lực kháng đàn tính của khối đá xung quanh giúp tiết kiệm vật liệu vỏ hầm ra sao?

Dưới tác dụng của cột nước áp lực hơn 50 m, vỏ hầm biến dạng nở hông và ép vào vách đá. Khối đá sinh ra phản lực đàn tính K0 từ 200 đến 1.600 kG/cm2 chống lại biến dạng này. Sự tương tác đó hấp thụ từ 35% đến 40% ứng suất kéo phát sinh trong kết cấu, cho phép kỹ sư giảm chiều dày lớp bê tông cốt thép lót từ 50 cm xuống còn khoảng 30 cm.

Điểm khác biệt căn bản giữa phương pháp cơ học kết cấu và phương pháp phần tử hữu hạn là gì?

Phương pháp cơ học kết cấu đơn giản hóa bài toán thành hệ thanh, vòm hoặc dầm trên nền đàn hồi một chiều với các giả thiết phân bố lực parabôn. Trong khi đó, phương pháp phần tử hữu hạn mô phỏng trường ứng suất và biến dạng liên tục hai chiều, phản ánh chính xác trạng thái tiếp xúc phi tuyến và tính dị hướng phức tạp của các đới đá trachit và trầm tích.

Tổ hợp tải trọng nào được xác định là nguy hiểm nhất đối với kết cấu đường hầm Mường Kim II?

Tổ hợp nguy hiểm nhất là giai đoạn tháo cạn nước trong hầm để kiểm tra định kỳ trong mùa khô. Lúc này không còn áp lực nước trong để cân bằng, vỏ hầm phải chịu toàn bộ áp lực thẳng đứng của đá núi (hệ số lệch tải n = 1,5), áp lực đẩy ngang của đất đá (n = 1,2) và áp lực thủy tĩnh của nước ngầm (n = 1,1), gây ra ứng suất uốn nén cục bộ lớn nhất.

Kết luận

  • Luận văn đã phân tích toàn diện đặc điểm địa hình, địa chất và thủy lực của công trình thủy điện Mường Kim II với quy mô công suất 10,5 MW và chiều dài tuyến hầm 1.600 m.
  • Xác lập cơ sở khoa học chứng minh mặt cắt hình tròn là hình thức kết cấu hợp lý nhất cho đoạn hầm có áp chịu cột nước trên 50 m, giúp triệt tiêu mômen uốn tập trung và tối ưu hóa phân bố ứng suất tiếp.
  • Lượng hóa chính xác vai trò của lực kháng đàn tính từ nền đá với hệ số K0 đạt tới 1.400 kG/cm2, giúp giảm 35% đến 40% ứng suất kéo nguy hiểm trong bê tông vỏ hầm.
  • Đề xuất các giải pháp kỹ thuật đồng bộ về xử lý phụt vữa áp lực cao từ 0,5 đến 1,0 MPa, bố trí chiều dày vỏ bê tông M250 dày 25 đến 35 cm và lắp đặt mạng lưới quan trắc tự động 24/7.
  • Đóng góp hệ thống phương pháp luận tính toán kết hợp giữa giải tích cơ học kết cấu và mô hình số phần tử hữu hạn cho các dự án thủy điện vừa và nhỏ tại Việt Nam.

Đóng góp chính của công trình là cung cấp một giải pháp thiết kế đường hầm dẫn nước có áp chuẩn mực, bảo đảm tính bền vững kết cấu và hiệu quả kinh tế cao cho các dự án năng lượng vùng Tây Bắc. Kế hoạch tiếp theo cần tập trung mở rộng quan trắc số liệu đo đạc thực địa trong chu kỳ 5 đến 10 năm khai thác để hiệu chỉnh các mô hình tương tác địa cơ học. Các cơ quan tư vấn, ban quản lý dự án và doanh nghiệp xây dựng hãy liên hệ và áp dụng ngay khung giải pháp kỹ thuật này để nâng cao chất lượng và tối ưu hóa chi phí cho các công trình ngầm thủy công hiện nay.