Tổng quan nghiên cứu

Sự phát triển nhanh chóng của các công trình đập bê tông khối lớn đóng vai trò then chốt trong khai thác tài nguyên nước và phát triển năng lượng bền vững. Theo số liệu từ Ủy ban Đập cao Thế giới (ICOLD), tính đến năm 2000, toàn cầu đã xây dựng khoảng 45.000 đập lớn, trong đó Trung Quốc dẫn đầu với 22.000 đập (chiếm 48%), Hoa Kỳ sở hữu 6.575 đập và Việt Nam xếp thứ 16 với 460 đập. Đối với các công trình có chiều cao trên 100m, đập bê tông trọng lực chiếm ưu thế vượt trội với tỷ lệ 38%. Tại Việt Nam, hàng loạt công trình thủy điện trọng điểm như Sơn La, Bản Vẽ, Tuyên Quang và Hủa Na đã minh chứng cho bước tiến vượt bậc về quy mô kỹ thuật xây dựng công trình thủy lợi.

Tuy nhiên, thách thức lớn nhất trong thi công bê tông khối lớn, đặc biệt là bê tông đầm lăn (RCC) và bê tông thường (CVC), là hiện tượng phát sinh ứng suất nhiệt do quá trình thủy hóa xi măng. Khi ứng suất kéo vượt quá cường độ chịu kéo của bê tông, các vết nứt bề mặt và đặc biệt là nứt xuyên khối sẽ hình thành, đe dọa nghiêm trọng đến tính ổn định và khả năng chống thấm. Thực tế tại đập thủy điện Sơn La từng ghi nhận vết nứt rộng 7cm tại khối C3 và L1, hay đập Upper Stillwater (Hoa Kỳ) cao 90m xuất hiện 13 vết nứt nhiệt dài trung bình 49m và rộng 10mm gây rò rỉ nước hạ lưu.

Nghiên cứu này tập trung phân tích trường nhiệt độ và trường ứng suất trong thân đập bê tông CVC thuộc công trình thủy điện Hủa Na (tỉnh Nghệ An) có chiều cao 94,5m. Mục tiêu cụ thể là xác định mối tương quan giữa tốc độ thi công và nhiệt độ ban đầu của hỗn hợp bê tông, từ đó đề xuất giải pháp kỹ thuật tối ưu nhằm hạ nhiệt độ đỉnh từ 49,65°C xuống dưới 45,0°C và triệt tiêu nguy cơ nứt xuyên thân đập.

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 sự tích hợp giữa lý thuyết truyền nhiệt phi ổn định và lý thuyết đàn hồi cơ học vật rắn biến dạng trong môi trường liên tục. Quá trình trao đổi nhiệt tổng thể trong thân đập được mô hình hóa thông qua ba cơ chế: dẫn nhiệt nội bộ, tỏa nhiệt đối lưu với không khí và bức xạ nhiệt mặt trời tác động lên các bề mặt hở.

Trường nhiệt độ không ổn định hai chiều được xác định thông qua phương trình vi phân truyền nhiệt tổng quát kết hợp toán tử Laplace và thông số nguồn nhiệt thủy hóa nội sinh:

$$\frac{\partial T}{\partial \tau} = a \left( \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} \right) + \frac{q_v}{C \cdot \rho}$$

Trong đó, các đại lượng vật lý đặc trưng gồm hệ số dẫn nhiệt $\lambda = 1,2 - 2,8\text{ W/m}^\circ\text{C}$, nhiệt dung riêng $C = 0,22 - 0,26\text{ kcal/kg}^\circ\text{C}$, khối lượng riêng $\rho = 2400\text{ kg/m}^3$ và nhiệt lượng tỏa ra từ các đơn khoáng xi măng (khoáng $C_3S$ tỏa 110 cal/g và $C_3A$ tỏa 188 cal/g sau 7 ngày). Điều kiện biên loại 3 của Newton-Richman được thiết lập để mô tả sự trao đổi nhiệt qua bề mặt ván khuôn và bề mặt bê tông tiếp xúc môi trường.

Về mặt cơ học, bài toán biến dạng phẳng của lý thuyết đàn hồi được áp dụng theo định luật Hooke mở rộng có xét đến biến dạng nở nhiệt ($\alpha \cdot \Delta T$). Ứng suất nhiệt bề mặt và ứng suất nhiệt đáy khối đập chịu sự khống chế của hệ số tự biến ($K$) và hệ số kiềm chế nền ($R$), quyết định trạng thái ứng suất kéo gây nứt kết cấu.

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

Nguồn dữ liệu đầu vào phục vụ mô hình bao gồm chuỗi số liệu quan trắc khí tượng thực tế tại khu vực tuyến công trình Hủa Na (nhiệt độ không khí, vận tốc gió trung bình năm, độ ẩm tại Bái Thượng và trực xạ mặt trời) kết hợp kết quả thí nghiệm cơ lý mẫu bê tông CVC (cường độ chịu nén, cường độ chịu kéo và mô đun đàn hồi từ 20 đến 32 GPa).

Phương pháp phần tử hữu hạn (PTHH) 2D được lựa chọn làm công cụ tính toán chủ đạo. Toàn bộ mặt cắt ngang lớn nhất của đập Hủa Na được chia lưới thành hàng nghìn phần tử tam giác 3 nút chuyển vị bậc nhất. Phương trình cân bằng nhiệt và ma trận độ cứng kết cấu được giải rời rạc hóa theo từng bước thời gian ($\Delta\tau$) bằng thuật toán sai phân nội suy tuyến tính.

Cỡ mẫu mô phỏng gồm 5 phương án thi công toàn diện (PA1 đến PA5), phản ánh các kịch bản tốc độ nâng đập khác nhau (tổng thời gian 380 ngày và 501 ngày) kết hợp việc kiểm soát hoặc không kiểm soát nhiệt độ hỗn hợp ban đầu. Phương pháp PTHH được lựa chọn vì đây là giải pháp duy nhất mô phỏng chính xác quá trình thi công nâng đập gián đoạn theo từng lớp đổ, sự thay đổi tính chất vật liệu theo tuổi bê tông và điều kiện biên thay đổi liên tục theo thời gian thực.

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

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

Mô phỏng số trường nhiệt độ và trường ứng suất trên đập Hủa Na qua 5 kịch bản thi công đã mang lại các kết quả định lượng cụ thể:

Thứ nhất, vùng tích tụ nhiệt lớn nhất ($T_{max}$) luôn tập trung ở cao độ 194m trong mọi kịch bản. Ở phương án đối chứng không xử lý nhiệt độ hỗn hợp ban đầu (PA1), nhiệt độ cực đại đạt mức 49,65°C sau 5.199 giờ (khoảng 216 ngày) thi công.

Thứ hai, việc khống chế nhiệt độ hỗn hợp bê tông ban đầu có tác dụng làm giảm đáng kể đỉnh nhiệt. Tại phương án PA2 với cùng tiến độ 501 ngày nhưng áp dụng làm mát hỗn hợp, $T_{max}$ giảm xuống còn 44,4°C sau 5.202 giờ (giảm 5,25°C, tương đương giảm 10,6% so với PA1). Tại phương án PA3, nhiệt độ đỉnh được duy trì ở mức 46,0°C.

Thứ ba, việc tăng tốc độ thi công lên 380 ngày (rút ngắn 24,1% thời gian) trong các phương án PA4 và PA5 làm mật độ phát nhiệt gia tăng nhanh chóng. Đỉnh nhiệt xuất hiện sớm hơn rất nhiều, đạt 46,5°C ở PA4 và 45,1°C ở PA5 chỉ sau khoảng 3.808 - 3.810 giờ (159 ngày).

Thứ tư, thời điểm phát sinh ứng suất kéo nguy hiểm nhất tại chân đập có sự dịch chuyển rõ rệt. Trong các phương án thi công nhanh (PA4 và PA5), ứng suất kéo cực đại xuất hiện sau khi đã đổ bê tông xong tới đỉnh đập cao độ 244m (PA4 đạt đỉnh ứng suất sau 10.302 giờ, tức 49 ngày sau khi hoàn thành đập; PA5 đạt đỉnh sau 11.447 giờ, tức 97 ngày sau khi thi công xong).

Thảo luận kết quả

Diễn biến ứng suất nhiệt trong thân đập Hủa Na phản ánh đúng bản chất vật lý của sự tương tác nhiệt - biến dạng. Trong giai đoạn tăng nhiệt ban đầu, phần lõi bê tông giãn nở tạo ra ứng suất nén ở bên trong và ứng suất kéo ở bề mặt ngoài do chênh lệch nhiệt độ với không khí. Tuy nhiên, giai đoạn nguy hiểm thực sự xảy ra trong quá trình hạ nhiệt kéo dài hàng trăm ngày sau đó.

Khi khối bê tông nguội dần và co lại, sự kiềm chế cứng từ nền đá móng và các lớp bê tông cũ đã đông cứng bên dưới (hệ số kiềm chế $R \approx 1,0$) ngăn cản biến dạng co tự do. Điều này làm chuyển hóa toàn bộ biến dạng co thành ứng suất kéo tập trung cao độ tại vùng tiếp giáp đáy.

Dữ liệu tính toán có thể được minh họa trực quan qua biểu đồ đường cong lịch sử nhiệt độ - ứng suất theo thời gian tại các cao độ 194m, 227m, 239m và 244m, kết hợp sơ đồ phân bố đẳng ứng suất dạng bản đồ nhiệt trên mặt cắt ngang thân đập. So sánh với các công trình quốc tế như đập Liễu Khê (Trung Quốc) hay đập Upper Stillwater (Hoa Kỳ), việc đẩy nhanh tốc độ nâng đập mà không đi kèm biện pháp hạ nhiệt hỗn hợp đổ sẽ tích lũy nhiệt lượng khổng lồ, khiến gradient nhiệt độ vượt ngưỡng an toàn ($> 20^\circ\text{C}$), từ đó kích hoạt các vết nứt xuyên thân đập trong giai đoạn bắt đầu tích nước vận hành.

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

Nhằm đảm bảo an toàn tuyệt đối cho công trình đập CVC Hủa Na và các dự án đập bê tông trọng lực tương tự, 4 nhóm giải pháp kỹ thuật cần được triển khai đồng bộ:

  • Khống chế nhiệt độ hỗn hợp bê tông ban đầu: Chủ đầu tư và nhà thầu trộn bê tông cần thiết lập trạm làm lạnh cốt liệu bằng nước lạnh và sử dụng nước đá vảy khi trộn trong mùa nắng nóng, đặt mục tiêu duy trì nhiệt độ bê tông trước khi đổ không vượt quá 20 - 22°C nhằm giới hạn đỉnh nhiệt thân đập dưới 45,0°C. Timeline thực hiện: xuyên suốt từ giai đoạn sản xuất bê tông đến khi kết thúc đợt đổ.
  • Tối ưu hóa quy trình phân lớp và nhịp độ đổ: Đội ngũ kỹ sư thi công hiện trường cần khống chế chiều dày mỗi lớp đổ từ 1,5m đến 2,0m và duy trì khoảng cách thời gian giữa hai lớp đổ liên tiếp từ 5 đến 7 ngày nhằm tạo điều kiện cho bề mặt tỏa nhiệt tự nhiên, hạn chế chênh lệch nhiệt độ trong - ngoài khối đập dưới 15°C. Timeline thực hiện: duy trì liên tục trong suốt tiến độ 380 - 501 ngày.
  • Sử dụng xi măng ít tỏa nhiệt và phụ gia khoáng hoạt tính: Ban chỉ đạo vật liệu cần thay thế 15% đến 30% lượng xi măng Poóc-lăng bằng phụ gia tro bay hoặc Puzolan hoạt tính để giảm nhiệt thủy hóa 7 ngày xuống dưới 50 - 60 cal/g, đồng thời bổ sung phụ gia hóa dẻo chậm đông kết. Timeline thực hiện: kiểm định cấp phối 60 ngày trước khi triển khai thi công đại trà.
  • Tăng cường bảo dưỡng giữ nhiệt và làm mát bề mặt: Đơn vị thi công đập cần duy trì chế độ phun nước bảo dưỡng liên tục và che phủ bạt cách nhiệt trên bề mặt bê tông mới đổ trong tối thiểu 28 ngày đầu để tránh hiện tượng sốc nhiệt khi thời tiết chuyển lạnh đột ngột. Timeline thực hiện: thực hiện ngay sau khi tháo dỡ ván khuôn cho đến 90 ngày tuổi của từng khối đổ.

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

Nội dung và kết quả nghiên cứu mang lại giá trị thực tiễn cao cho nhiều nhóm chuyên môn trong ngành xây dựng công trình thủy lợi - thủy điện:

  • Kỹ sư thiết kế công trình thủy lợi - thủy điện: Tiếp cận khung phương pháp tính toán mô phỏng nhiệt - ứng suất tin cậy, làm cơ sở xây dựng chỉ dẫn kỹ thuật thi công và tiêu chuẩn nghiệm thu nhiệt độ bê tông khối lớn.
  • Chỉ huy trưởng và kỹ sư quản lý thi công hiện trường: Sử dụng các thông số định lượng về tốc độ đổ và nhiệt độ ban đầu để lập tiến độ thi công chi tiết, điều hành trạm trộn và bố trí mặt bằng công nghệ hợp lý.
  • Chuyên gia tư vấn giám sát và quản lý chất lượng dự án: Tham khảo quy trình theo dõi diễn biến trường nhiệt và các ngưỡng ứng suất giới hạn để thiết lập hệ thống quan trắc nhiệt độ thân đập bằng cảm biến chuyên dụng.
  • Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Xây dựng Thủy điện/Địa kỹ thuật: Sử dụng làm tài liệu tham khảo học thuật giá trị về ứng dụng phương pháp phần tử hữu hạn trong giải quyết bài toán cơ - nhiệt phi tuyến phức tạp.

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

Bê tông khối lớn trong đập trọng lực CVC được định nghĩa như thế nào?
Theo tiêu chuẩn ACI 116R-90 của Hoa Kỳ, bê tông khối lớn là kết cấu bê tông có kích thước hình học đủ lớn đòi hỏi phải áp dụng các biện pháp kiểm soát nhiệt sinh ra do quá trình thủy hóa xi măng. Mục đích chính là hạn chế sự biến đổi thể tích gây nứt vỡ trong thân đập.

Tại sao vết nứt xuyên khối lại đặc biệt nguy hiểm đối với đập dâng nước?
Khác với nứt bề mặt có thể tự khép kín, nứt xuyên khối chia cắt tính toàn khối của thân đập, làm suy giảm khả năng chịu lực trọng lực và tạo đường rò rỉ nước áp lực cao từ thượng lưu về hạ lưu như từng xảy ra tại đập Upper Stillwater (13 vết nứt dài 49m), đe dọa sự an toàn của toàn bộ cụm công trình.

Tốc độ thi công ảnh hưởng như thế nào đến ứng suất nhiệt trong đập Hủa Na?
Khi đẩy nhanh tiến độ thi công từ 501 ngày xuống 380 ngày (phương án PA4 và PA5), nhiệt lượng thủy hóa bị giữ lại nhiều hơn trong thân đập, dẫn đến việc ứng suất kéo cực đại xuất hiện muộn hơn (từ 49 đến 97 ngày sau khi đổ xong đập) tại vùng tiếp giáp chân móng bị kiềm chế.

Làm thế nào để khống chế hiệu quả nhiệt độ hỗn hợp bê tông ban đầu?
Biện pháp hiệu quả nhất là làm mát cốt liệu lớn bằng nước lạnh, sử dụng đá vảy thay thế một phần nước trộn và phối trộn 15% - 30% tro bay hoặc Puzolan. Giải pháp này giúp hạ nhiệt độ thủy hóa của xi măng sau 7 ngày xuống dưới ngưỡng 60 cal/g.

Phương pháp phần tử hữu hạn giải bài toán nhiệt - ứng suất đập bê tông bằng cách nào?
Phương pháp chia miền mặt cắt đập thành lưới phần tử tam giác, sử dụng hàm nội suy bậc nhất để giải hệ phương trình vi phân truyền nhiệt từng bước theo thời gian ($\Delta\tau$), sau đó tích hợp trường nhiệt độ vào bài toán biến dạng phẳng đàn hồi để tìm phân bố ứng suất kéo.

Kết luận

  • Luận văn đã mô phỏng thành công trường nhiệt độ và trường ứng suất nhiệt 2D cho đập bê tông CVC Hủa Na cao 94,5m trên 5 phương án thi công khác nhau.
  • Xác định cao độ 194m là vị trí tích tụ nhiệt nguy hiểm nhất với nhiệt độ đỉnh dao động từ 44,4°C đến 49,65°C tùy theo chế độ làm mát ban đầu.
  • Chứng minh rằng việc thi công nhanh trong 380 ngày làm đỉnh nhiệt xuất hiện sớm sau 159 ngày và kéo dài thời gian xuất hiện ứng suất kéo đáy đập tới 97 ngày sau khi hoàn thành.
  • Khẳng định tính ưu việt và độ chính xác của phương pháp phần tử hữu hạn trong việc mô hình hóa bài toán phi ổn định cơ - nhiệt cho công trình thủy lợi lớn.
  • Đề xuất giải pháp kiểm soát nhiệt độ hỗn hợp ban đầu dưới 22°C kết hợp phụ gia khoáng nhằm triệt tiêu hoàn toàn nguy cơ nứt nhiệt thân đập.

Đóng góp chính của công trình là cung cấp cơ sở khoa học và định lượng chuẩn xác giúp các nhà quản lý dự án cân bằng tối ưu giữa bài toán đẩy nhanh tiến độ kinh tế và chi phí xử lý nhiệt độ. Kế hoạch tiếp theo cần tập trung vào việc số hóa dữ liệu quan trắc nhiệt tự động trong suốt quá trình tích nước hồ chứa. Các kỹ sư và nhà quản lý xây dựng công trình thủy hãy ứng dụng ngay mô hình tính toán nhiệt - ứng suất PTHH vào giai đoạn thiết kế thi công để bảo đảm chất lượng công trình trường tồn cùng thời gian.