祁海燕

(遼寧西北供水有限責(zé)任公司，遼寧 沈陽 110003)

1 工程簡介

恒仁水庫位于遼寧省恒仁縣境內(nèi)的渾江中游河段。恒仁水電站為渾江梯級開發(fā)中的主力電站，同時恒仁水庫也是大伙房輸水工程的重要水源地水庫。恒仁水庫始建于1958年，1967年水庫蓄水，1975年全部竣工。恒仁水庫的正常蓄水位為300m，設(shè)計洪水位302m，死水位為280m。水庫控制流域面積10400km2，多年平均徑流量約為44.80億m3，多年平均流量為142m3/s，多年平均輸沙量為125.6萬m3[1]。

樞紐工程主要由大壩、壩后式廠房、溢洪道、泄洪隧洞等建筑物組成。水庫大壩為混凝土單支墩大頭壩，壩頂高程為304.2m，最大壩高為78.5m，壩頂寬32.7m，壩頂長593m[2]。壩體設(shè)有表孔、中孔、底孔，分別擔(dān)負(fù)著泄洪、水庫正常放水以及沖砂任務(wù)。

2 監(jiān)測資料分析

2.1 環(huán)境量分析

大壩水平位移分析的建模時段為2010年1月～2014年12月，預(yù)測時段為2015年1～3月。在數(shù)據(jù)分析時，只對相應(yīng)時段的數(shù)據(jù)進(jìn)行分析和利用。

2.1.1 上下游水位分析

恒仁水庫為不完全年調(diào)節(jié)水庫，從上游水位實測值來看，一般7～9月水位較高，2～4月水位較低，呈較明顯的周期性變化特征。其中最高水位302.67m，超過設(shè)計洪水位0.68m，最低水位為286.04m，多年平均水位為295.54m。

水庫下游水位變化比較平穩(wěn)，除個別泄洪時段外，基本在247.5～249.5m之間變動。最高水位為252.51m，最低水位為247.61m，多年平均水位為248.8m。

2.1.2 溫度分析

壩址區(qū)氣溫年周期變化十分明顯，7～8月氣溫較高，12月～次年3月氣溫較低。極端最高氣溫為32.1℃，最低氣溫為-42.5℃，氣溫變幅較大，多年平均氣溫為7.1℃。水庫大壩301.00m高程以上受氣溫影響明顯，與氣溫變化有明顯的同步性，隨著水深增加，水溫受氣溫周期變化的影響逐步減小。

2.2 變形監(jiān)測分析

由于恒仁水庫大壩的視準(zhǔn)線觀測頻度為兩個月一次，不具有顯著的代表性，因此只對引張線進(jìn)行分析。引張線布置于壩頂軸線下方0+10.0m處，共設(shè)置10個觀測點，編號分別為EX1-EX10，本文僅對其中的EX5觀測點進(jìn)行分析。從2008年10月21日以來的引張線自動化監(jiān)測結(jié)果來看，壩頂順河向水平位移主要受氣溫周期性變化影響，溫度降低向下游位移，并于升高則相反。特別是7～8月水位最高時，大壩壩頂向下游的位移量最小，因此大壩水平位移受水位因素的影響較小。

3 大壩水平位移監(jiān)測混合模型

3.1 混合模型的優(yōu)勢

目前應(yīng)用于大壩位移監(jiān)測領(lǐng)域主要是統(tǒng)計模型和確定性模型[3]。其中，統(tǒng)計模型的優(yōu)勢在于模型構(gòu)建比較方便，同時還可獲得較高的精度和擬合度。統(tǒng)計模型的不足之處在于不能從力學(xué)角度解釋各分量的對效應(yīng)量的影響，同時各分量之間也不是隨機(jī)變量關(guān)系，一旦遇到突發(fā)情況，統(tǒng)計關(guān)系就不適用了。確定性模型是基于工程力學(xué)，通過計算環(huán)境量與效應(yīng)量之間的關(guān)系構(gòu)建的具有因果關(guān)聯(lián)性的模型，可以很好彌補(bǔ)統(tǒng)計模型的不足，但是確定型模型的應(yīng)用需要完整的工程資料，并需要大量的數(shù)學(xué)計算。

為了彌補(bǔ)上述兩種模型的不足，盡量發(fā)揮其優(yōu)勢，以上述兩種模型為基礎(chǔ)構(gòu)建混合模型，對溫度和時效分量采用統(tǒng)計模型擬合，而對應(yīng)力分量則采用確定性模型進(jìn)行分析計算，無疑更具有實用價值。

3.2 混合模型的構(gòu)建

根據(jù)現(xiàn)有的理論，水工構(gòu)筑物運(yùn)行過長中產(chǎn)生的位移主要與壓力載荷、溫度載荷和時效變化直接相關(guān)[4]。在建模過程中，假定大壩上某點的位移為δ，這個位移可以由空間坐標(biāo)和誘發(fā)的物理因素兩種不同方法進(jìn)行分解，其表達(dá)式為：

δ(δx/δy/δz)=δH+δT+δθ

(1)

式中，δ—大壩上某點的位移量；δx—水平位移；δy—側(cè)向位移；δz—豎直位移；δH—水壓周期性分量；δT—溫度周期性分量；δθ—時效性分量。

3.3 水壓分量及有限元計算

在混合模型構(gòu)建中，水壓分量采用基于確定性模型的有限元計算方法。同時，水壓分量由庫水位分量δH1，揚(yáng)壓力分量δH2和淤沙壓力分量δH3三部分組成。

3.3.1 庫水位分量

由上述可知，水庫大壩上某點的水平位移δH1可以表示為：

δH1=δ1H+δ2H+δ3H

(2)

由于恒仁水庫大壩為混凝土重力壩，根據(jù)相關(guān)文獻(xiàn)和有限元計算結(jié)果可得：δH1與上游水深H以及H2、H3的多項式之間是正比關(guān)系，其表達(dá)式如下：

(3)

式中，a—待定回歸系數(shù)；H—上游水深。

3.3.2 揚(yáng)壓力分量

揚(yáng)壓力分量包括壩基揚(yáng)壓力和壩身揚(yáng)壓力[5]。根據(jù)相關(guān)文獻(xiàn)和工程力學(xué)計算，壩基揚(yáng)壓力的計算公式為：

(4)

式中，δH2f—壩基揚(yáng)壓力；hd—壩高；Ec—壩體砼的彈性模量；m—下游壩面坡度；d—觀測點距壩頂?shù)呢Q直距離；H—上游水深。

壩身揚(yáng)壓力對大壩水平位移的影響主要與上下游的水頭差線性相關(guān)，其表達(dá)式為：

(5)

式中,變量的含義同上。

3.3.3 淤沙壓力分量

由于泥沙淤積對大壩水平位移的作用比較復(fù)雜，同時這種作用還會隨著泥沙的不斷淤積而變化，因此這種影響作用并不能通過簡單的數(shù)學(xué)參數(shù)予以表達(dá)[6]。考慮到其影響不能忽略不計，因此在建模過程中將其視為時效因子考慮。

3.3.4 有限元計算

根據(jù)上述計算公式，對有限元計算中的壩基和壩體相關(guān)參數(shù)假定見表1。

表1 大壩不同部位力學(xué)參數(shù)假定數(shù)值

計算過程中的下游水位區(qū)平均值248.8m。上游水位選擇包括最高水位、最低水位、設(shè)計和校核水位的原則，選取304m、301m、299m、297m、295m、293m、291m、289m、287m，共9組水位進(jìn)行計算。

選擇EX5所在的壩體段為整體建立有限元模型。其研究邊界為上游2倍壩高，下游1.5倍壩高，豎直1.5倍壩高。壩體兩側(cè)為自由邊界，基巖兩側(cè)限制軸向運(yùn)動，基巖上下游面施加水流方向約束，基巖底部三向固定約束。整個模型劃分為23716個節(jié)點，22171個單元。

利用有限元模型計算不同上游水位條件下壩頂?shù)乃轿灰屏恳姳?。

表2 不同上游水位條件下EX5段壩頂位移計算結(jié)果

將上述表格中的數(shù)據(jù)進(jìn)行擬合，并引入壩體混凝土平均彈性模量的的修正系數(shù)X后，獲得EX5水壓分量表達(dá)式：

δH=X(28.71-4.265ΔH2+25.103ΔH3)

(6)

3.4 混合模型表達(dá)式

3.4.1 溫度分量

由于恒仁水庫大壩監(jiān)測只設(shè)置了六個溫度監(jiān)測點，且其中三個監(jiān)測點由于溫度計不能正常工作而沒有完整的監(jiān)測資料，因此溫度資料極為缺乏，不便于采用有限元模型計算，只能采用統(tǒng)計模型進(jìn)行擬合。鑒于多溫度模型的擬合效果明顯優(yōu)于但溫度模型，故采用多溫度模型進(jìn)行擬合[7]，并得到如下成果：

(7)

式中，b—待定回歸系數(shù)；T—觀測日大壩溫度；t—距初始觀測日的時間，d。

3.4.2 時效分量

時效分量利用統(tǒng)計模型擬合，結(jié)果為：

(8)

式中，C1，C2—待定系數(shù)；t—距離基準(zhǔn)日觀測天數(shù)。

綜合上述成果，可得大壩水平位移的混合模型表達(dá)式：

(9)

式中，符號所代表的含義同上。

3.5 模型參數(shù)分析

運(yùn)用SPASS軟件進(jìn)行回歸計算[8]，并剔除異常點，獲得的各項系數(shù)見表3。

由表3的結(jié)果可知，上述混合模型的回歸方程為：

(10)

表4 EX5實測值與預(yù)測值對比表

注：正值表示向下游偏移，負(fù)值表示向上游偏移。

3.6 模型效果分析

通過對2015年2月至2016年1月的實測值和預(yù)測值的對比，其中部分?jǐn)?shù)值見表4，可以發(fā)現(xiàn)兩者具有較好的擬合度，模型的剩余標(biāo)準(zhǔn)差為0.375，能夠基本描述發(fā)展趨勢，具有理想的預(yù)測效果。因此，混合模型兼具統(tǒng)計模型和確定性模型的優(yōu)勢，具有一定的實際用用價值。

4 結(jié)語

通過人為的、定期的、科學(xué)的大壩監(jiān)測，不僅可以實現(xiàn)對大壩現(xiàn)實安全性實施監(jiān)控，還有利于對安全隱患未爆發(fā)時采取必要的工程措施，具有重要的實際應(yīng)用價值。我國的許多大壩建于20世紀(jì)50，60年代，運(yùn)行時間已經(jīng)達(dá)到50年以上，運(yùn)營管理單位急需對大壩進(jìn)行各方面的監(jiān)測，持續(xù)關(guān)注其運(yùn)行狀態(tài)。本文基于對恒仁水庫大壩的監(jiān)測研究，提出了具有一定精確度，便于操作的大壩水平位移混合模型，對大壩運(yùn)行管理者及時判斷大壩的運(yùn)行狀況具有重要幫助。但是，研究模型仍然存在待定系數(shù)多、計算量大的問題，今后需要在此方面進(jìn)行深入研究。