，，

(1.河南靈捷水利勘測(cè)設(shè)計(jì)研究有限公司,河南 南陽(yáng) 473000;2.南陽(yáng)市水利建筑勘測(cè)設(shè)計(jì)院,河南 南陽(yáng) 473000)

復(fù)雜介質(zhì)的土壩滲流計(jì)算分析

杜躍亭1，賈鵬生2，張麗麗1

(1.河南靈捷水利勘測(cè)設(shè)計(jì)研究有限公司,河南 南陽(yáng) 473000;2.南陽(yáng)市水利建筑勘測(cè)設(shè)計(jì)院,河南 南陽(yáng) 473000)

在早期修建的土石壩水利工程中，由于地質(zhì)勘查相對(duì)落后，不能準(zhǔn)確查清壩基地質(zhì)條件，加上施工過(guò)程中取材有差異，導(dǎo)致壩體、壩基局部相鄰介質(zhì)的滲透性相差較大。常規(guī)的滲流計(jì)算方法得到的滲流計(jì)算結(jié)果不合理，通過(guò)增設(shè)滲流過(guò)渡區(qū)域的手段，降低相鄰介質(zhì)滲透系數(shù)偏差，較好地解決了復(fù)雜介質(zhì)的滲流計(jì)算問題。

復(fù)雜介質(zhì)；土石壩;滲流分析

早期的土石壩，大多為“三邊”工程，加之當(dāng)時(shí)的地質(zhì)勘察技術(shù)相對(duì)落后，不能很清楚地探明壩址區(qū)域內(nèi)的地層情況，導(dǎo)致很多壩基未清理至發(fā)育較好的巖層、土層即開始筑壩，或者直接對(duì)山體進(jìn)行削坡后即作為壩體的一部分。另外，由于早期施工過(guò)程中，取材有差異，導(dǎo)致壩體材料非常復(fù)雜，相鄰的兩層介質(zhì)滲透性可能相差1～4個(gè)量級(jí)。由于早期水庫(kù)運(yùn)行水位較低，滲流問題不是很明顯，然而隨著經(jīng)濟(jì)社會(huì)的發(fā)展，用水需求日益增大，這類水庫(kù)的滲流問題逐漸凸顯。

針對(duì)這一問題的滲流場(chǎng)數(shù)值分析，目前較多采用各向同性滲流理論或沿坐標(biāo)軸的各向異性理論進(jìn)行分析[1，2]。常規(guī)的分析結(jié)果表明：相鄰兩層介質(zhì)的滲透系數(shù)比值小于103數(shù)量級(jí)時(shí)，可以直接進(jìn)行計(jì)算，否則將出現(xiàn)錯(cuò)誤的結(jié)果或迭代不收斂問題[3]。

為了避免出現(xiàn)相鄰介質(zhì)層滲透系數(shù)相差較大的情況，可通過(guò)設(shè)置過(guò)渡區(qū)域的手段加以解決。基于滲流有限元分析的基本原理, 利用Geo-studio中的SEEP /W 模塊[4],模擬壩體滲流情況。

1 計(jì)算方法

1.1 有限元計(jì)算原理

對(duì)土壩的滲透分析采用基于有限元法的SEEP/W進(jìn)行，通過(guò)賦予材料特性和邊界條件，求得土體內(nèi)部的滲流場(chǎng)和滲透力分布情況，分析最可能發(fā)生破壞的位置。

有限元滲流分析時(shí)水體在土體內(nèi)部的流動(dòng)方程符合達(dá)西規(guī)律，根據(jù)是否考慮時(shí)間效應(yīng)，將計(jì)算分析劃分為穩(wěn)態(tài)和瞬態(tài)兩種情況。穩(wěn)態(tài)為理想狀態(tài)，不考慮時(shí)間和初始內(nèi)部體積含水量變化的影響，流體的運(yùn)動(dòng)方程見式(1)，實(shí)際上這種理想狀態(tài)在入滲過(guò)程中很難達(dá)到。瞬態(tài)分析考慮了時(shí)間作用，用以分析庫(kù)水位變化時(shí)入滲過(guò)程中水位及浸潤(rùn)線等隨時(shí)間的變化過(guò)程，流體運(yùn)動(dòng)方程見式(2)，當(dāng)計(jì)算時(shí)間足夠長(zhǎng)時(shí)，瞬態(tài)分析結(jié)果接近穩(wěn)態(tài)分析結(jié)果。

穩(wěn)態(tài)方程：

(1)

瞬態(tài)方程：

(2)

上二式中H——總水頭；

kx、ky——x、y方向的滲透系數(shù)；

q——邊界流量；

θ——體積含水量；

t——時(shí)間。

有限元軟件SEEP/W材料本構(gòu)中引入非飽和土特性，程序中反映非飽和材料特性的參數(shù)有：土體的滲透系數(shù)、基質(zhì)吸力、土體體積含水量等；其中非飽和土體的滲透系數(shù)是土體體積含水量和基質(zhì)吸力的函數(shù)[5]，土體體積含水量的變化受到土體應(yīng)力狀態(tài)和土性的影響，土體應(yīng)力狀態(tài)主要由(σ-ua)和(ua-uw)控制(其中：σ1土體總應(yīng)力，在滲透分析中保持不變；ua為孔隙氣壓力)[6]，假定孔隙氣壓力在滲流分析中也保持恒定，則土體體積含水量只與ua-uw有關(guān)，體積含水量與時(shí)間的關(guān)系可由式(3)表示：

(3)

式中mw——由負(fù)孔隙水壓力引起的體積變化系數(shù)，為土水特征曲線斜率。

將式(3)代入式(2)得到非飽和土體的流體運(yùn)動(dòng)方程：

(4)

1.2 過(guò)渡區(qū)域滲透系數(shù)及厚度的確定

由于過(guò)渡區(qū)域的滲透性和兩種介質(zhì)都有關(guān)系，同時(shí)也為防止因兩種介質(zhì)滲透性差異過(guò)大導(dǎo)致模型求解不收斂問題，采用下式確定過(guò)渡區(qū)域的滲透系數(shù)[3]：

(5)

式中k——過(guò)渡區(qū)域的滲透系數(shù)；

k1、k2——相鄰介質(zhì)1、2的滲透系數(shù)。

過(guò)渡帶的厚度根據(jù)單元尺寸取0.5m。

2 工程概況

南陽(yáng)西部某水庫(kù)于1968年10月開工興建，水庫(kù)控制流域面積58km2,總庫(kù)容2447萬(wàn)m3，工程等別為Ⅲ等，主要建筑物級(jí)別為3級(jí)。洪水標(biāo)準(zhǔn)為1000年一遇校核，校核洪水位258.65m(黃海高程系,下同)，相應(yīng)庫(kù)容2182萬(wàn)m3；50年一遇設(shè)計(jì)，設(shè)計(jì)洪水位257.12m，相應(yīng)庫(kù)容1973萬(wàn)m3；興利水位253.20m，興利庫(kù)容1250萬(wàn)m3。主壩為黏土均質(zhì)壩，壩頂高程261.30m，壩長(zhǎng)630m，最大壩高36.3m。目前，該水庫(kù)是一座以防洪、灌溉為主，兼顧發(fā)電、養(yǎng)魚等綜合利用的水利樞紐工程。

該水庫(kù)大壩左臺(tái)地段原為一條很厚的土丘嶺，修建水庫(kù)時(shí)在左臺(tái)地上、下游取土，并按照設(shè)計(jì)壩坡坡度將天然土丘嶺進(jìn)行削坡，最終形成類似壩體形式的壩基(根據(jù)工程實(shí)際，為了便于模型建立和滲流分析計(jì)算，簡(jiǎn)化為壩體進(jìn)行分析)。壩體材料自上而下依次為：壩體填土、泥質(zhì)中砂、粉質(zhì)壤土、泥質(zhì)砂礫石(詳見圖1)。

壩體填土厚15.8～18.3m,黃褐、棕褐色，填土性質(zhì)較復(fù)雜，以含礫重粉質(zhì)壤土、含礫粉質(zhì)黏土為主，含風(fēng)化石渣、礫石及砂粒，分布不均勻。泥質(zhì)中砂以透鏡體形式分布于241.90～244.50m高程，層厚0.8～4.2m, 滲透系數(shù)4.65×10-5m/s，強(qiáng)透水層，局部貫通形成滲漏通道。粉質(zhì)壤土層厚9.7～14.0m，含少量礫石、砂粒，砂粒分布不均勻。泥質(zhì)砂礫石層頂高程229.1～234.1m,層厚2.0～7.3m，滲透系數(shù)9.2×10-4m/s，強(qiáng)透水層，泥質(zhì)含量21.4%，砂粒含量50.4%，卵、礫石含量28.2%，卵石粒徑5～10cm，成分以灰白色砂礫巖為主。

由于水庫(kù)除險(xiǎn)加固前一直為低水位運(yùn)行，未形成滲流破壞。近年來(lái)隨著水資源日益緊缺，工程除險(xiǎn)加固后，水庫(kù)開始長(zhǎng)期高水位運(yùn)行，強(qiáng)透水層滲流問題逐漸顯現(xiàn)，水位在254.60m時(shí)，散浸面積達(dá)780m2，水位在250.00m時(shí)，散浸面積也有150m2。為此，設(shè)計(jì)采取塑性混凝土防滲墻工程措施來(lái)進(jìn)行除險(xiǎn)加固。

圖1 最大壩高斷面地層分布

3 計(jì)算模型及參數(shù)

結(jié)合土石壩壩址區(qū)地形地質(zhì)條件,以壩段最大壩高剖面為對(duì)象,建立平面有限元計(jì)算模型，采用 3節(jié)點(diǎn)和4節(jié)點(diǎn)平面等參單元對(duì)壩體及壩基劃分單元,單元總數(shù)8951個(gè), 節(jié)點(diǎn)總數(shù)9165個(gè),二維有限元模型見圖2。各分區(qū)材料的滲透系數(shù)見下表。

圖2 SEEP/W單元網(wǎng)格模型

壩體及過(guò)渡區(qū)域材料滲透系數(shù)表

4 計(jì)算成果及分析

本文僅選取正常水位工況下壩體滲流場(chǎng)的計(jì)算結(jié)果進(jìn)行分析，計(jì)算結(jié)果見圖3、圖4。

圖3 未設(shè)置滲流過(guò)渡區(qū)域滲流場(chǎng)分布 單位：m

圖4 設(shè)置滲流過(guò)渡區(qū)域滲流場(chǎng)分布 單位：m

由圖3、圖4正常水位工況下壩體壩基的滲流場(chǎng)可以發(fā)現(xiàn)，當(dāng)壩體未設(shè)置滲流過(guò)渡區(qū)域時(shí)，浸潤(rùn)線經(jīng)塑性混凝土心墻時(shí)折減較小，壩體浸潤(rùn)線呈不規(guī)則分布，背水坡中部出現(xiàn)封閉的滲流場(chǎng)，且出現(xiàn)二次逸出的現(xiàn)象，出現(xiàn)該種情況的原因是相鄰兩層介質(zhì)的滲透系數(shù)太大，導(dǎo)致有限元計(jì)算時(shí)形成的滲透系數(shù)矩陣出現(xiàn)奇異。當(dāng)壩體設(shè)置滲流過(guò)渡區(qū)域時(shí)，可以看出浸潤(rùn)線經(jīng)塑性混凝土心墻時(shí)折減較大，滲流場(chǎng)分布合理，浸潤(rùn)線呈光滑漸變，逸出點(diǎn)位置合理，說(shuō)明增設(shè)滲流過(guò)渡區(qū)域可以較好地解決相鄰介質(zhì)滲透系數(shù)相差較大的滲流計(jì)算問題。另外，通過(guò)在壩體設(shè)置的浸潤(rùn)觀測(cè)管觀測(cè)數(shù)據(jù)，其形成的浸潤(rùn)線與增設(shè)過(guò)渡區(qū)域后的滲流計(jì)算結(jié)果接近，這也進(jìn)一步說(shuō)明本文的結(jié)果是合理的。

5 結(jié) 語(yǔ)

早期修建的土石壩，地質(zhì)勘查技術(shù)相對(duì)落后，未準(zhǔn)確查清地質(zhì)情況即開始筑壩，導(dǎo)致清基不徹底，壩體比較復(fù)雜，相鄰介質(zhì)的滲透性相差較大，常規(guī)的滲流計(jì)算方法得到的滲流計(jì)算結(jié)果大多不合理甚至錯(cuò)誤。本文通過(guò)在滲透性相差較大的相鄰介質(zhì)中間增設(shè)滲流過(guò)渡區(qū)域的手段，降低相鄰介質(zhì)滲透系數(shù)偏差，較好地解決了復(fù)雜介質(zhì)的滲流計(jì)算問題。

[1]Brunone B.，Romano N.，Santini A.Numerical Simulations of one-dimensional infiltration into layered soils with Richards equation using different estimates of the interlayer conductivity [J].Vadose Zone Journal，2003(2): 193-200.

[2]Lima-Vivancos V.，Voller V.R.Two numerical methods for modeling variably saturated flow in layered Media [J].Vadose Zone Journal，2004(3): 1031-1037.

[3]劉昌軍,丁留謙,寧保輝,張順福.心墻土石壩滲流場(chǎng)的無(wú)單元法模擬[J]. 水力發(fā)電學(xué)報(bào)，2012(3).

[4]Fredlund D G， Rahardjo H.Soil Mechanics for Unsaturated Soils[M].New York:Wiley,1993.

[5]Fredlund D G， Morgenstern N R. Constitutive relations for volume change in unsaturated soils [J]. Can Geotech J.,1976(13)：261-276.

CalculationAnalysisofEarthDamSeepageinComplicatedMedium

DU Yue-ting1,JIA Peng-sheng2, ZHANG Li-li1

(1.HenanLingjieHydraulicInvestigationandDesignResearchInstituteCo.,Ltd.,Nanyang473000,China;

2.NanyangHydraulicBuildingInvestigationandDesignInstitute,Nanyang473000,China)

During early process of constructing earth and rock dam hydraulic project, geological condition on the foundation of the dam can not be found precisely due to geological investigation was relatively lagged behind, plus there was difference in using materials during the process of construction, causing greater differences of penetration in local neighbored medium in dam body and the foundation of dam. Seepage calculation result obtained from regular seepage calculation methodology is not reasonable, and seepage calculation issue of complicated medium can be better solved with means of adding seepage transition area and reducing deviation of penetration coefficient of neighbored medium.

complicated medium; earth and rock dam; seepage analysis

TV641

A

1673-8241(2014)09-0067-04