晏國順 陜亮 李家正

摘要：大古水電站位于西藏高寒高海拔地區(qū)，壩址區(qū)日照輻射強烈，顯著影響了環(huán)境溫差的變化，加劇了施工期乃至運行期大壩結構溫度場時空分布的不均勻性，增加開裂風險。選取大古水電站14號廠房壩段為研究對象，采用三維有限元法對大壩施工期溫度場和溫度應力進行數值模擬仿真計算，分析了日照強熱輻射對大壩的影響。結果表明：日照強熱輻射在壩體表面產生了約0.6 MPa的拉應力，拉應力區(qū)影響深度在0.6 m范圍內;采取混凝土表面保溫防護措施后，拉應力減少至約0.1 MPa，有效降低了混凝土開裂風險。

關鍵詞：日照強輻射;碾壓混凝土重力壩;溫度裂縫;三維有限元;大古水電站;高寒高海拔地區(qū);西藏

中圖法分類號：TV315 文獻標志碼：A DOI：10.15974/j.cnki.slsdkb.2021.10.009

文章編號：1006 - 0081（2021）10 - 0047 - 05

0 引 言

中國剩余水能資源開發(fā)主要分布在西部、西南部，集中于雅魯藏布江、金沙江、瀾滄江、怒江流域[1]，這些地區(qū)氣候變化劇烈，極端天氣頻繁，水電站建設難度大。碾壓混凝土筑壩技術具有澆筑面大、連續(xù)澆筑施工效率高、不設縱縫少設橫縫的優(yōu)點，因此碾壓混凝土壩成為高寒高海拔地區(qū)復雜環(huán)境下最有利的壩型[2-4]。在自然環(huán)境中的混凝土結構由于熱傳導性能較差，氣候條件對混凝土的質量影響極大，混凝土內外大溫差是導致壩體發(fā)生表面裂縫的主要原因 [5-7]。大氣的溫度變化并不是直接來源于太陽的熱量，而是來源于大地上空的空氣。太陽輻射主要以短波為主，其中大氣吸收16%，地表和海洋則吸收了51%，這51%的短波變成了熱量，又被地表以長波的形式往外輻射，這些長波被大氣吸收，影響了大氣的溫度。大古水電站壩址海拔為3 500 m，位于青藏高原高海拔地區(qū)，是中國太陽輻射最強的地區(qū)，空氣稀薄，大氣對太陽輻射的削弱作用弱，到達地面的太陽輻射強。因此，本文以大古水電站為研究對象，采用三維有限元法計算分析了日照強熱輻射對大壩的影響，旨在為類似復雜環(huán)境下的工程提供參考。

1 工程概況

大古水電站位于西藏自治區(qū)山南地區(qū)桑日縣境內，是雅魯藏布江中游桑日縣至加查縣峽谷段的第2級電站，上游距規(guī)劃的巴玉水電站約8 km，下游距規(guī)劃的街需水電站和已建成的藏木水電站分別為7 km和18 km。電站壩址控制流域面積15.74萬km2。攔河壩為碾壓混凝土重力壩，壩高126 m，電站裝機容量66萬kW，采用壩后式廠房布置。廠房壩段結構復雜，材料分組多，本文以14號廠房壩段為典型壩段進行研究分析。

2 數值模擬

2.1 結構模型

14號壩段結構模型見圖1，材料分區(qū)見圖2。坐標軸X方向為順水流方向，指向下游為正;坐標軸Y方向為壩軸線方向（橫河向），指向左岸為正;坐標軸Z方向為垂直向，向上為正，與結構實際高程一致。

本文采用長江科學院自主研發(fā)的大型水工數值分析軟件包（Ckysts1.0），進行結構溫度場和溫度應力三維有限元線彈性仿真計算，網格劃分共計335 414個節(jié)點，309 790個單元。

2.2 計算條件

依據大古水電站相關建設單位技術報告[8-11]，結合混凝土澆筑進度，參考相關類似工程經驗以及行業(yè)規(guī)程規(guī)范[12]確定計算條件。

2.2.1 環(huán)境溫度

（1）年氣溫作用。年氣溫作用是指氣溫年變化對混凝土結構的影響，其變化相對緩慢。根據大古氣象站氣象要素月平均氣溫統計，采用下式擬合：

[Tc（t）=10.2+7.9sin2π365（t-100）]? ? ? ? ? ?（1）

式中：[Tc（t）]為隨時間變化的氣溫;t為計算時刻至1月1日的時間。

（2）日照溫度作用。日照溫度變化主要是太陽輻射作用引起，大壩混凝土表面受到日照直射輻射、散射輻射的作用，表面溫度要高于該時刻的大氣溫度。日照影響相當于氣溫增高值[ΔTα]，由下式計算：

[ΔTα=R/β]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（2）

式中：R為物體表面吸收的熱輻射，材料吸收系數取0.65;β為物體表面放熱系數。太陽輻射統計和受日照影響的氣溫增高換算值列于大古氣象站氣象要素統計表1中。

由表1看出，大古水電站受日照影響相當于年平均氣溫增加10.5℃。5月份日照最強，相當于氣溫增高12.5℃;12月份日照最弱，相當于氣溫增高7.6℃;年變幅A=（12.5-7.6）/2=2.4℃。因此，在計算分析中按年平均氣溫提高10.5℃，氣溫年變幅增加2.4℃的情況進行大古水電站日照溫度作用影響的敏感性分析，擬合式如下：

[Tc（t）=20.7+10.3sin2π365（t-100）]? ? ? ? （3）

2.2.2 材料力學和熱學性能

（1）混凝土及基巖彈性模量?；炷粮鼾g期彈性模量用式（4）擬合，擬合參數見表2。

[E（t）=E0（1-e-AtB）]? ? ? ? ? ? ? ? ? ? ? ? ? （4）

式中：[E0]為混凝土最終彈性模量，GPa;[t]為混凝土齡期，d;A，B為擬合參數。

基巖變形模量取值為13.5 GPa，泊松比取0.2，不計自重?；炷撩芏龋耗雺夯炷罜9015為2.42×103

kg/m3，其他常態(tài)混凝土取為2.37×103 kg/m3，泊松比0.167。

（2）混凝土徐變度。溫度應力計算中，混凝土徐變度采用式（5）擬合：

[C（t， τ）=C1（τ）[1-e-k1（t-τ）]+C2（τ）[1-e-k2（t-τ）]] （5）

式中： C（t，τ）為混凝土徐變度，×10-6/MPa;t為混凝土齡期，d;τ為混凝土加荷齡期，d;[C1（τ）=A1B1τm] ，[C2（τ）=A2B2τn]，A1，A2，B1，B2，n，m，k1，k2為擬合系數，見表3。

（3）熱學性能參數?；炷?、基巖熱學性能參數見表4?；炷两^熱溫升通過膠凝材料水化熱值以及混凝土配合比的膠凝材料用量計算得出，采用[θ（t）=θ0tA+t]擬合，式中：[θ（t）]為混凝土絕熱溫升，℃;t為混凝土齡期，d;[θ0]為混凝土最終絕熱溫升，℃。

2.2.3 邊界條件與初始條件

（1）邊界條件。溫度場計算中，基礎各側面、底面、壩段對稱面取絕熱邊界，壩體暴露在空氣中的面取氣溫邊界。壩體與氣溫接觸的邊界，按第三類邊界條件處理，放熱系數β取10.73 W/（m2·℃），等效熱交換系數為927 kJ/（m2·d·℃）。壩體廊道孔洞內取氣溫邊界，考慮到空氣流通不暢，孔洞內表面與大氣接觸， 放熱系數β取5.0 W/（m2·℃），等效熱交換系數為452 kJ/（m2·d·℃）。

應力計算中，基礎左右兩側面、基礎下游面取法向約束，基礎底面取三向約束;基礎上游面取為自由面?？紤]到壩體側面與相鄰壩段之間有橫縫，取為自由面。

（2）初始條件。基巖初溫擬定為10.2℃，基巖上、下游表面在氣溫邊界條件下計算20 a至混凝土開始澆筑時所得的溫度場作為基巖的初始溫度場?；炷翝仓r的澆筑溫度作為混凝土的初始溫度。

2.2.4 澆筑溫度與通水冷卻

廠房壩段建基面呈多層臺階式，混凝土相應允許澆筑溫度見表5。混凝土下料時即開始一期通水冷卻，通水水溫12℃，通水時長30 d;高溫季節(jié)澆筑的混凝土入冬前進行必要的中期通水冷卻降溫，從當年9月初開始，通水時長60 d，通水水溫12℃，冷卻目標溫度16℃。

2.2.5 混凝土表面保護

高寒地區(qū)碾壓混凝土壩表面保溫方案多樣[13]，大古水電站14號廠房壩段混凝土表面分為保溫、不保溫2種情況計算分析，保溫效果參數見表6。

3 結果分析

在自然環(huán)境中的混凝土結構，由于熱傳導性能較差，氣候條件對混凝土的質量影響極大，日照熱輻射等作用使混凝土表面溫度迅速變化，而其結構內部溫度仍處于原來狀態(tài)，從而形成較大溫度梯度，產生溫度變形，這樣的變形在受到結構內、外約束阻礙時，會產生相當大的溫度應力，施工階段由于混凝土水化熱作用，更容易形成混凝土內外表面較大溫度差，引起裂縫的產生。

3.1 溫度場

壩體內部混凝土與壩體外表面混凝土最高溫度見表7，日照影響的混凝土外表面溫度歷程見圖3。溫度場的仿真計算表明：混疑土外表面溫度受太陽日照的影響較大，混凝土澆筑層內部區(qū)域受太陽日照的影響較小，在混凝結構中造成較大的溫度梯度。與不考慮太陽日照影響相比，大壩混凝土表面保溫條件下，考慮日照影響的外表面混凝土最高溫度增加了10.61℃。結構內部碾壓混凝土、排沙廊道常態(tài)混凝土、壩下0+74前常態(tài)混凝土最高溫度分別增加了1.80℃，0.34℃，0.47℃;大壩混凝土表面不保溫條件下，考慮日照影響的外表面混凝土最高溫度增加了12.56℃，結構內部碾壓混凝土、排沙廊道常態(tài)混凝土、壩下0+74前常態(tài)混凝土最高溫度分別升高2.62℃，1.50℃，0.71℃。

3.2 溫度應力

壩體內部與外表面混凝土橫河向拉應力對比分別見表8和表9，壩體外表面混凝土橫河向應力歷程見圖4。溫度應力的仿真計算表明：日照強輻射作用對壩體內部混凝土應力基本無影響，而對外表面混凝土拉應力影響較大，增大了外表面混凝土的開裂風險。壩體內部混凝土，采取表面保溫措施時，相對于不考慮太陽日照影響，混凝土橫河向拉應力，碾壓混凝土增加了0，排沙廊道常態(tài)混凝土減小0.03 MPa，壩下0+74前常態(tài)混凝土減小0.03 MPa;表面不保溫情況下，內部碾壓混凝土增加了0.04 MPa，排沙廊道常態(tài)混凝土減小0.23 MPa，壩下0+74前常態(tài)混凝土減小0.07 MPa。壩體表面混凝土采取表面保溫措施時，相對于不考慮日照影響，混凝土橫河向拉應力增加0.12 MPa，拉應力區(qū)影響深度在0.6 m以內;表面不保溫情況下，橫河向拉應力增加0.54 MPa，拉應力區(qū)影響深度在1.2 m以內。

4 結 語

（1）高寒高海拔地區(qū)氣候環(huán)境復雜，太陽日照對直接接受太陽輻射熱的外表面混凝土拉應力影響較大，增加了混凝土表面開裂風險。

（2）壩體混凝土表面采用保溫防護措施十分必要，可有效降低工程在施工期面臨的開裂風險。

（3）在較高溫季節(jié)日照強烈時澆筑混凝土，有必要采取遮陽板、遮陽棚等溫控措施，減小太陽輻射造成環(huán)境溫度高于壩體溫度情況而產生熱量倒灌現象的影響。

參考文獻：

[1] 范國福. 西藏水電規(guī)劃綜述[J]. 水力發(fā)電，2008，38（7）：14-20.

[2] 鈕新強.? 中國水電工程技術創(chuàng)新實踐與新挑戰(zhàn)[J].? 人民長江， 2015， 46（19）：13-17.

[3] 李響， 彭松濤， 田德智，等.? 相同設計強度下凝灰?guī)r粉對水工混凝土性能的影響[J].? 人民長江， 2019，50（增2）：216-219.

[4] 王松春.? 中國碾壓混凝土筑壩技術及其發(fā)展[J].? 水利水電科技進展， 2001， 21（2）：16-18，22.

[5] 邱冬梅， 熊濤， 蔡暢，等.? DG水電站大壩工程碾壓混凝土倉面設計[J].? 水電與新能源， 2019（12）：12-16.

[6] 龔召雄. 水工混凝土的溫控與防裂[M]. 北京：中國水利水電出版社，1999.

[7] 朱伯芳. 大體積混凝土的溫度應力與溫度控制[M]. 北京：中國電力出版社，2012.

[8] 華東勘測設計研究院有限公司. 西藏大古水電站大體積混凝土溫控防裂施工技術要求[R]. 杭州：華東勘測設計研究院有限公司，2018.

[9] 華東勘測設計研究院有限公司. 西藏大古水電站大壩混凝土熱力學性能試驗檢測報告[R]. 杭州：華東勘測設計研究院有限公司，2019.

[10] 中國水利水電第九工程局有限公司. DBTJ-進度-2019-007 2019年實施性進度計劃（已審批）[R]. 貴陽：中國水利水電第九工程局有限公司，2019.

[11] 中國水利水電第九工程局有限公司. 西藏大古水電站大壩土建及金屬結構安裝工程施工總進度計劃[R]. 貴陽：中國水利水電第九工程局有限公司，2018.

[12] NB/T 35092-2017 混凝土壩溫度控制設計規(guī)范[S]

[13] 尚層.? 高寒地區(qū)碾壓混凝土壩上游面保溫方案比選[J].? 人民長江， 2016， 47（1）：79-82.

（編輯：江 文）

Influence of sunlight strong heat radiation on RCC gravity dam of

Dagu Hydropower Station in Tibet

YAN Guoshun1， SHAN Liang2，3， LI Jiazheng2，3

（1. Huadian Tibet Hydropower Development Co.， Ltd， Lasa 856200， China;? ?2. Material and Engineering Structure Department， Changjiang Scientific Research Institute， Wuhan 430010， China;? ?3. Research Center of Water Engineer Security and Disease Prevention Under Ministry of Water Resources， Changjiang Scientific Research Institute， Wuhan 430010， China）

Abstract： Dagu hydropower station is located in the high cold and high altitude area of Tibet. The sunshine radiation in the dam site area has a strong impact on the change of environmental temperature difference， which aggravates the uneven distribution of temperature field in the construction period and even operation period， and increases the risk of cracking. The No.14 plant dam section is selected as the research object. The three-dimensional finite element method is used to simulate the temperature field and temperature stress during the construction period of the dam， and the influence of the strong heat radiation of sunshine on the dam is analyzed. The results show that the tensile stress of the dam surface is about 0.6MPa， and the influence depth of the tensile stress area is within 0.6m. After taking the thermal insulation protection measures on the concrete surface， the tensile stress decreases to about 0.1 MPa， which effectively reduces the risk of concrete cracking.

Key words： strong sunlight radiation; compacted concrete gravity dam; temperature crack; 3D finite element; Dagu hydropower station; High altitude and cold area; Tibet