郭張軍,徐建光,韓建新,馮超,彭鵬

（1.陜西電力科學(xué)研究院,陜西西安710054；2.安康水力發(fā)電廠，陜西安康725000；3.河海大學(xué),江蘇南京210098）

0 引言

消力池是水工建筑物經(jīng)常采用的消能設(shè)施之一，其長(zhǎng)期穩(wěn)定和安全運(yùn)行對(duì)確保大壩正常發(fā)電和汛期行洪安全至關(guān)重要，特別是消力池底板安全穩(wěn)定尤為關(guān)鍵。大量工程實(shí)踐表明，長(zhǎng)期處于水下的消力池底板混凝土在滲水壓力的作用下，材料中的水泥水化物溶出和分解并不斷被滲水帶走，導(dǎo)致混凝土微觀結(jié)構(gòu)改變和孔隙率增加，滲透性增大，化學(xué)侵蝕逐步加重，混凝土因而損失膠凝性，強(qiáng)度和耐久性能逐漸下降。因此，消力池底板混凝土耐久性的減弱是滲流場(chǎng)、化學(xué)場(chǎng)等多場(chǎng)耦合作用下的結(jié)果[1，2]。

目前，國(guó)內(nèi)很多學(xué)者已對(duì)壩址區(qū)的多場(chǎng)耦合作用做了很多研究[3，6]，但是考慮有關(guān)水化學(xué)侵蝕的還比較少。水化學(xué)侵蝕是消力池底板混凝土強(qiáng)度及耐久性能劣化的重要原因，必須予以重視。本文通過建立化學(xué)損傷的混凝土-水-化學(xué)耦合模型，數(shù)值模擬底板混凝土化學(xué)侵蝕的過程，分析化學(xué)侵蝕對(duì)消力池底板混凝土強(qiáng)度和耐久性的影響，預(yù)測(cè)化學(xué)侵蝕對(duì)已破損消力池底板今后運(yùn)行的影響，并為今后消力池底板修復(fù)中如何提高混凝土抗化學(xué)侵蝕提供依據(jù)，工程實(shí)例表明，該耦合模型有較高的精度。

1 多場(chǎng)耦合模型

水文地球化學(xué)的角度來講，消力池底板混凝土的化學(xué)侵蝕是在滲水壓力作用下，由混凝土-水系統(tǒng)狀態(tài)變化引起的化學(xué)組分在混凝土固相和液相（滲水溶液）兩相之間重新分布的結(jié)果，是一種滲入水溶液在混凝土中流動(dòng)、遷移并與周圍混凝土土體不斷發(fā)生一系列化學(xué)作用的混凝土-水化學(xué)作用過程。消力池底板混凝土-水-化學(xué)耦合模型就是在適當(dāng)?shù)臈l件下，通過滲流模型、溶質(zhì)遷移模型和地球化學(xué)模型相互耦合而建立的。

1.1 基本假設(shè)

研究區(qū)為等效連續(xù)介質(zhì)，水流為飽和流，其運(yùn)動(dòng)服從達(dá)西定律；滲流和組分遷移之間的耦合是從滲流到遷移的單向耦合，即滲流狀態(tài)的變化對(duì)水溶液中組分濃度的變化有顯著影響，而組分濃度變化引起的水密度、粘性等物理特性改變對(duì)滲流的影響可以忽略不計(jì)；混凝土材料孔隙度、滲透系數(shù)等物理特性只是空間位置的函數(shù)，不隨時(shí)間而變化。

1.2 多場(chǎng)耦合作用模型

根據(jù)上述假定，認(rèn)為運(yùn)行工況下消力池底板混凝土-水-化學(xué)相互作用的多場(chǎng)耦合作用模型由多個(gè)子模塊構(gòu)成。包括滲流模塊、化學(xué)反應(yīng)模塊[7]、溶質(zhì)遷移模塊等，見式（1）。

模型中，μs為儲(chǔ)水率；H為水頭；K為滲透系數(shù)張量；t為時(shí)間；f0、f1、f2均為已知初始函數(shù)和邊界的函數(shù)；Γ1、Γ2分別為已知水頭和已知流量邊界；nΓ為Γ2邊界外法線方向；C為地下水中含有目標(biāo)離子的總濃度；Ci為可溶礦物溶解目標(biāo)離子的濃度；XEq是由化學(xué)熱力學(xué)引起的目標(biāo)濃度隨時(shí)間的變化量；Nc為水溶液中組分總數(shù)；ai為第i個(gè)離子或絡(luò)合物的活度；Ki為第i個(gè)離子或絡(luò)合物組分的化學(xué)計(jì)量數(shù)；νij為第i個(gè)反應(yīng)中第j個(gè)組分的計(jì)量數(shù)；aj為地下水溶液中第j個(gè)組分的活度；Tj為第j個(gè)組分的總摩爾濃度；?Kin為由動(dòng)力學(xué)引起的目標(biāo)濃度隨時(shí)間的變化量；kp為第P個(gè)反應(yīng)的速率常數(shù)；Sp為第P個(gè)礦物物種反應(yīng)表面積；Kp為第P個(gè)反應(yīng)的平衡常數(shù)；下標(biāo)Eq和Kin表示相關(guān)量分別對(duì)應(yīng)于熱力學(xué)反應(yīng)和動(dòng)力學(xué)反應(yīng)。u為地下水實(shí)際流速；DL為垂直方向水動(dòng)力彌散系數(shù)；DT為水平方向水動(dòng)力彌散系數(shù)；ux為地下水水平方向?qū)嶋H流速；uy為地下水垂直方向?qū)嶋H流速；W為源匯項(xiàng)，是由于化學(xué)作用引起的目標(biāo)組分濃度變化；B1、B2分別為已知邊界條件；（n1，n2）為B2邊界外法線方向。

1.3 模型的求解

采用弱耦合方式進(jìn)行求解，弱耦合又稱順序耦合或者間接耦合[8]，它是多場(chǎng)按照一定的順序排列后，再逐個(gè)進(jìn)行計(jì)算分析，即將前一個(gè)場(chǎng)計(jì)算分析的結(jié)果作為已知條件或者邊界條件，帶入到第二個(gè)場(chǎng)分析中進(jìn)行耦合，依次類推，這樣可以簡(jiǎn)化計(jì)算步驟，減小計(jì)算量。以鈣離子為研究對(duì)象[9]，其主要步驟如下：

（1）首先在給定庫(kù)水位的情況下，計(jì)算得到Δt時(shí)刻消力池底板滲流場(chǎng)各點(diǎn)水頭值（滲流量）的分布。

（2）根據(jù)滲流場(chǎng)各點(diǎn)水頭值（滲流量）的分布計(jì)算各個(gè)部位由化學(xué)反應(yīng)引起的溶解（或沉淀）量，可以得到Δt時(shí)刻分別由化學(xué)熱力學(xué)和化學(xué)動(dòng)力學(xué)貢獻(xiàn)鈣離子的反應(yīng)量。

（3）將化學(xué)反應(yīng)計(jì)算得到的鈣離子反應(yīng)量與水中原有的鈣離子含量相疊加，并代入溶質(zhì)遷移模型，利用有限元的求解方法得到Δt時(shí)刻鈣離子含量在滲流場(chǎng)中的分布。

（4）根據(jù)鈣離子的遷移量，計(jì)算Δt時(shí)刻消力池底板鈣離子的含量。

（5）根據(jù)模擬時(shí)間要求，反復(fù)進(jìn)行步驟（2）～（4）直到滿足一定的時(shí)間精度。

2 實(shí)例應(yīng)用

2.1 工程概況

安康水電站位于漢江上游安康城西南18 km處，是以發(fā)電為主、兼顧航運(yùn)、防洪、旅游及養(yǎng)殖等綜合效益的水利樞紐工程。電站樞紐由攔河壩、右岸壩后廠房、左岸泄流及垂直升船機(jī)等建筑物組成。壩體為折線型混凝土重力壩，水庫(kù)控制流域面積35600 km2,平均年徑流量192億m3,水庫(kù)正常蓄水位330 m，正常發(fā)電尾水位325 m左右，汛限水位325 m，設(shè)計(jì)洪水35700 m3/s,泄洪建筑物包括5個(gè)表孔、5個(gè)中孔和4個(gè)底孔。表孔布置在河床中部11～15壩段，采用寬尾墩與消力池結(jié)合的消能方式，表孔消力池處于左右底孔之間。5個(gè)表孔中，左右兩側(cè)表孔（即1號(hào)和5號(hào)）采用非對(duì)稱寬尾墩，中間為對(duì)稱寬尾墩。右底孔為窄縫式異型鼻坎，表孔和右底孔泄流均納入表孔消力池內(nèi)消能。安康水電站不但承擔(dān)著下游城市和長(zhǎng)江下游段的防洪安全，而且是陜西電網(wǎng)重要的調(diào)峰電源，其安全穩(wěn)定至關(guān)重要。安康水電站的表孔消力池底板雖經(jīng)4次修補(bǔ)，仍屬“帶病”運(yùn)行，隱患沒有得到根治，嚴(yán)重危脅著廠房、大壩的安全。根據(jù)歷次表孔消力池的外觀檢查以及強(qiáng)度、耐久性檢測(cè)成果來看，底板表面破損較嚴(yán)重，基礎(chǔ)廊道滲水量較大，混凝土強(qiáng)度及耐久性能產(chǎn)生不同程度的降低。從滲流水質(zhì)成分化驗(yàn)、析出物成分檢驗(yàn)結(jié)果來看，消力池底板滲水具有較強(qiáng)的軟水侵蝕性作用，且堿性較強(qiáng)，鈣離子濃度較高，個(gè)別廊道內(nèi)析出物鈣含量較多。

2.2 模擬結(jié)果及分析

2.2.1 滲流場(chǎng)模擬結(jié)果及分析

滲流計(jì)算所獲的消力池底板內(nèi)各單元滲流速度是反應(yīng)-遷移計(jì)算的重要組成部分。安康水電站表孔消力池滲流模型以正常水頭差約20 m作為上、下表面的定水頭邊界條件進(jìn)行穩(wěn)定流計(jì)算，其結(jié)果見圖1和圖2。

圖1 表孔消力池底板滲流場(chǎng)等表面圖（單位：m）Fig.1 The surface map of seepage at plate of stilling basin

圖2 表孔消力池底板高程228 m（上、下兩層接觸面）處滲流場(chǎng)剖面圖（單位：m）Fig.2 The surface map of seepage at plate of stilling basin at the elevation 228 m(contact interface)

從滲流場(chǎng)計(jì)算結(jié)果可以看出：中間局部水頭值變化較快，表明這些部位徑流條件較好，滲透率較大，滲流速度較快，為消力池底板混凝土的化學(xué)侵蝕提供了有利條件。從圖2表孔消力池底板高程228 m（混凝土上下兩層間接觸面）處三維滲流場(chǎng)剖面圖來看，消力池底板抗沖層上表層的水頭值變化較快，滲流場(chǎng)相對(duì)較為活躍，尤其是原型檢測(cè)成果中發(fā)現(xiàn)裂縫較多的部位，這可能與其上表面分布的裂隙較多有關(guān)，滲透性較大，這為消力池底板混凝土的化學(xué)侵蝕提供了有效的水力條件。從鉛直方向來看，消力池底板混凝土抗沖層滲流場(chǎng)相對(duì)較為活躍，基礎(chǔ)混凝土層滲流條件相對(duì)滯緩，表明目前上層混凝土損傷相對(duì)較為嚴(yán)重，導(dǎo)致混凝土強(qiáng)度和耐久性的降低，而下層混凝土滲流條件相對(duì)滯緩，水頭變化較慢，反映出混凝土內(nèi)部完整性和抗?jié)B性較好。

2.2.2 耦合模型模擬結(jié)果及分析

消力池底板混凝土中主要化學(xué)成分是鈣，其在水化產(chǎn)物中也占重要地位。鈣類水化物在具有軟水特征的滲入液的作用下最易被溶解，造成混凝土孔隙結(jié)構(gòu)的改變，使裂縫擴(kuò)展，從而導(dǎo)致混凝土強(qiáng)度和耐久性能降低。由此可見，混凝土強(qiáng)度及耐久性能的劣化主要取決于鈣類水化產(chǎn)物的溶解過程，因此鈣離子濃度是混凝土化學(xué)侵蝕一個(gè)很好的指示劑，可以把Ca2+濃度作為化學(xué)侵蝕狀態(tài)變量，則消力池底板混凝土侵蝕過程就可以用孔隙溶液中鈣離子濃度變化來模擬。本研究針對(duì)安康水電站表孔消力池化學(xué)侵蝕的特點(diǎn)，立足于Ca2+濃度變化，從時(shí)間上和空間上定量研究混凝土化學(xué)侵蝕程度。

以安康水電站表孔消力池底板混凝土上表面滲入液中Ca2+的水化學(xué)檢測(cè)濃度和消力池底板化學(xué)侵蝕模擬計(jì)算活性離子濃度分別作為反應(yīng)-遷移模型的邊界條件和初始條件，以7天作為時(shí)間步伐，以10年作為模擬總時(shí)間，利用算子分裂法進(jìn)行整個(gè)消力池底板區(qū)域內(nèi)的反應(yīng)-遷移計(jì)算，獲得了由于混凝土-水間的化學(xué)反應(yīng)（反應(yīng)溫度在計(jì)算區(qū)均取15℃）引起的消力池底板混凝土與其中的滲水溶液之間物質(zhì)交換量在空間和時(shí)間上的變化過程，其結(jié)果見圖3、圖4和表1（圖表單位均為mg/L）。

圖3 表孔消力池底板混凝土滲水溶液鈣離子等表面圖Fig.3 The distribution of calcium ion concentration at plate of stilling basin

圖4 10年后消力池底板混凝土滲水溶液鈣離子等表面圖Fig.4 The distribution of calcium ion concentration after 10 years at plate of stilling basin

表1 安康水電站表孔消力池混凝土-水-化學(xué)耦合模型水化學(xué)資料統(tǒng)計(jì)表Table 1:Chemical compositions of the water in the con?crete-water-chemical coupling model

從模擬計(jì)算結(jié)果可以看出，消力池底板混凝土中諸如Ca(OH)2一類的水化產(chǎn)物在整個(gè)模擬期間（10年）始終處于溶解狀態(tài)，并隨著時(shí)間的增加而增加。其中，消力池底板混凝土上部抗沖層溶解強(qiáng)度相對(duì)較大，部分壩段消力池底板混凝土鈣離子溶解相對(duì)較多，濃度變化較大，最大鈣離子濃度為145 mg/L，這表明該部位消力池混凝土化學(xué)侵蝕較為嚴(yán)重，混凝土強(qiáng)度和耐久性能相對(duì)較差。從橫向分布來看，消力池底板12壩段和14壩段部位混凝土中鈣離子溶解相對(duì)較大，其次為13壩段和15壩段部位，但是其化學(xué)侵蝕程度也不容忽視；11壩段部位化學(xué)侵蝕程度相對(duì)較輕，其強(qiáng)度及耐久性相對(duì)較好。從縱向分布來看，消力池底板2池和3池部位混凝土中鈣離子溶解量相對(duì)較大，最大溶解量達(dá)120 mg/L，并隨著時(shí)間發(fā)展，混凝土中鈣離子含量越來越少，化學(xué)侵蝕對(duì)混凝土強(qiáng)度及耐久性的影響由“不顯著”變?yōu)椤帮@著”。從發(fā)展趨勢(shì)來看，目前鈣離子溶解量較大的部位，隨著鈣離子溶解量的增大（即混凝土中鈣的物質(zhì)含量減少），10年后，混凝土內(nèi)部形成集中的滲漏通道，將嚴(yán)重降低混凝土強(qiáng)度和耐久性，進(jìn)而危害消力池安全運(yùn)行。

3 結(jié)論及建議

（1）安康水電站表孔消力池環(huán)境水質(zhì)具有軟水性侵蝕作用的特征，正是通過這種環(huán)境水與混凝土中固相介質(zhì)發(fā)生化學(xué)反應(yīng)，使混凝土材料中可溶性物質(zhì)被溶失，造成混凝土裂縫擴(kuò)展和孔隙率增加，使混凝土強(qiáng)度和耐久性受到損害。安康水電站表孔消力池底板混凝土中鈣類水化物最易被溶解流失，最大溶解量約139 g/d，而此類物質(zhì)的流失將使混凝土孔隙率增加，滲透性增大，進(jìn)而降低混凝土強(qiáng)度和耐久性，形成化學(xué)損傷。

（2）鈣質(zhì)的析出是安康水電站表孔消力池底板混凝土強(qiáng)度和耐久性降低的重要原因，因此鈣離子濃度是混凝土化學(xué)侵蝕一個(gè)很好的指示劑，可以把Ca2+濃度作為化學(xué)侵蝕狀態(tài)變量。安康水電站表孔消力池底板混凝土中鈣類物質(zhì)目前一直處于流失狀態(tài)，并隨著時(shí)間發(fā)展逐漸增加，化學(xué)侵蝕速率加快，相對(duì)來說，消力池底板混凝土抗沖層遭受化學(xué)侵蝕作用程度相對(duì)較嚴(yán)重。

（3）安康水電站表孔消力池底板混凝土化學(xué)侵蝕程度的加重，將給其強(qiáng)度和耐久性帶來一系列不利影響，10年后，化學(xué)侵蝕較嚴(yán)重的部位將形成明顯的貫穿性裂縫或缺陷，影響消力池底板的整體性及穩(wěn)定性，同時(shí)也影響混凝土的強(qiáng)度和耐久性。綜合分析說明，安康水電站表孔消力池底板化學(xué)侵蝕對(duì)其強(qiáng)度及耐久性的影響是緩慢的、潛移默化的，但局部化學(xué)侵蝕作用相對(duì)較強(qiáng)，且這種作用是長(zhǎng)期存在的，對(duì)其工作性能的弱化作用具有不可逆性。

（4）建議對(duì)消力池底板抗沖層混凝土外觀缺陷進(jìn)行有效處理，以降低底板混凝土化學(xué)侵蝕速率，減小消力池底板混凝土安全運(yùn)行風(fēng)險(xiǎn)。建議在今后消力池底板修復(fù)中，通過在混凝土中添加一定量的摻合料、外加劑或改變不同材料的配合比等，降低鈣類等易溶失物質(zhì)的含量，以提高或改善混凝土的抗化學(xué)侵蝕性，盡可能提高混凝土的使用性能和延長(zhǎng)使用壽命。■

[1]宋漢周.大壩環(huán)境水文地質(zhì)研究[M].北京:中國(guó)水利水電出版社,2007.

[2]沈振中,陳小虎,徐力群.重力壩應(yīng)力-滲流相互作用的無單元耦合分析[J].巖土力學(xué),2008,29(11)：74～78.

[3]陳曉平,茜平一,梁志松，等.非均質(zhì)土壩穩(wěn)定性的滲流場(chǎng)和應(yīng)力場(chǎng)耦合分析[J].巖土力學(xué),2004,25(6)：860～864.

[4]頂秀麗,盛謙.三峽大壩左廠房3號(hào)壩段壩基滲流與應(yīng)力場(chǎng)耦合分析[J].巖石力學(xué)與工程學(xué)報(bào),2000,19(增)：1001～1005.

[5]盛金昌,速寶玉,趙堅(jiān)，等.滲流應(yīng)力耦合分析在溪洛渡電站壩址區(qū)的應(yīng)用[J].巖土力學(xué),2000,21(4)：410～415.

[6]柴軍瑞,李守義.大柳樹水利樞紐區(qū)滲流場(chǎng)與應(yīng)力場(chǎng)耦合分析[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(增1)：2322～2325.

[7]阿里木·吐爾遜.壩基老化巖－水－化學(xué)作用數(shù)值模擬研究[D].河海大學(xué)博士論文，2005.

[8]盛金昌.多孔介質(zhì)流-固-熱三場(chǎng)全耦合多孔介質(zhì)數(shù)學(xué)模型及數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報(bào),2006,25(增1)：3028～3033.

[9]Gérard,G.Pijaudier-Cabot and C.La Borderie.Coupled Dif?fusion-Damage Modeling and the Implications on Failure due to Strain Localization.International Journal of Solids and Struc?tures[J].1998,35(30):4107～4120.