李再興，李　靜，宋鵬飛，楊　芳

(河南省地質礦產(chǎn)勘查開發(fā)局測繪地理信息院，河南 鄭州 450006)

?

三維數(shù)值模擬技術(Visualmodflow)在基坑降水中的應用

李再興，李靜，宋鵬飛，楊芳

(河南省地質礦產(chǎn)勘查開發(fā)局測繪地理信息院，河南 鄭州 450006)

[摘要]基坑降水設計常采用大井法計算基坑涌水量，但是當研究區(qū)水文地質條件復雜，含水層非均值性，防滲圍堰的存在且滲透系數(shù)的不確定性時，利用“大井”解析法求解涌水量時，概化的條件較多，往往降低計算結果的可信度；而數(shù)值法不受其限制。因此本文采用了數(shù)值模擬軟件VisualModflow對基坑涌水量進行了模擬。在正確分析水文地質條件的基礎上，建立合適的水文地質概念模型并確定數(shù)學模型，通過非穩(wěn)定流計算求取基坑涌水量，并對防滲墻的防滲效果進行三種方案模擬計算，結合工程降水實際，推求防滲墻的滲透系數(shù)。

[關鍵詞]基坑降水；數(shù)值模擬；防滲帷幕

計算基坑涌水量的方法有很多[1]，如解析法、數(shù)值法和目標函數(shù)法等等,一般在基坑降水設計時，大多采用大井法，其原理易懂，方法易行。大井法一般需要對場地水文地質條件進行概化，概化成均值、等厚的含水層，再利用裘布依井流推導出來的計算公式進行求解。由于大井法概化的條件較多，概化誤差會降低參數(shù)和涌水量求解的可信度[2]。

數(shù)值模擬的優(yōu)點在于能夠刻畫均值、不等厚以及復雜的水文地質條件[3],近年來，數(shù)值模擬技術在實際應用中發(fā)揮了很大的優(yōu)越性。本文以長江流域某水電站基坑降水為例，采用三維數(shù)值模擬技術計算基坑降水涌水量，并結合工程降水實際，推求了防滲墻的滲透系數(shù)。

1研究區(qū)背景

1.1工程概述

長江流域某水電站左岸主體及導流工程分兩期施工。第一期工程先進行左岸圍堰，即在左岸灘地上修筑一期土石(防滲)圍堰，其縱向段緊靠河谷枯水期的河水位邊，上下游橫向段與左岸邊坡相連接，總體形態(tài)呈凸向右岸的弧型(見圖1)。在一期圍堰完成之后，在圍堰內(nèi)進行二期混凝土縱向圍堰上游段的施工。由于二期縱向大壩上游段地基覆蓋層平均深達45m，二期施工不能影響一期堰基穩(wěn)定和安全。通過多種處理方案的比較和分析，最后選定沉井處理方案作為二期縱向圍堰大壩上游段地基覆蓋層處理方案。

圖1　一期防滲圍堰與二期沉井群平面位置示意圖

沉井處理方案布置了若干個鄰接沉井，前期作為擋土墻及縱向圍堰堰基進行二期基坑開挖，然后作為二期圍堰的一部分，解決堰基覆蓋層的處理與二期工程施工的矛盾。

二期圍堰設計的沉井群由若干個23 m×17 m的矩形沉井組成，沉井底部進入下伏巖層，平均下沉深度45 m。沉井群依次沿一期土石圍堰成“L”型錯開布置(圖1)，采取排水下沉，沉井群下沉區(qū)可視為一個大型基坑。

設計一期防滲圍堰有多層防滲體系，包括覆蓋層的塑性砼防滲墻、下伏基巖帷幕灌漿，防滲體系的滲透系數(shù)在10-5～10-4m/d之間，防滲圍堰可視為隔水邊界。

1.2研究區(qū)地質概況

研究區(qū)兩岸地形整齊，基巖露頭良好，巖層連續(xù)分布，兩岸對應地層沒有明顯錯位。出露地層主要為三迭系上統(tǒng)(T3xj)和第四系(Qh)，第四系不整合于三疊系基巖之上。此外，三迭系中統(tǒng)雷口坡組(T21)灰?guī)r地層埋深于水電站壩基以下，埋深大于260 m，此層在地表未出露。

1.3研究區(qū)水文地質條件

擬建的沉井群場地位于長江左岸的灘地之上。沉井群中心距離長江約200 m。

沉井群所在的壩址區(qū)是上覆為45厚的第四系河流沖積物(Qhal)，大致可分為三層且含水層厚度不均：0～-15 m為砂卵礫石層，-15～-30 m為砂層，-10～-45 m為砂卵礫石層；砂層又可以分為上下兩段，上層為粉細質砂，厚10 m,下層為中粗砂，厚5 m；-45～-60 m為基巖，平均厚15 m,主要為三疊系的粉細砂巖，為主要含水層，微風化為主，裂隙與解理較發(fā)育。且向下延伸至-80 m左右，含水量急劇減少，可視為不透水含水層(隔水層)。

沉井區(qū)長江河槽切割至地面以下-28 m，第四系孔隙水與長江水發(fā)生直接水力聯(lián)系，相互滲補；下伏基巖裂隙水與上部第四系孔隙水不存在相對隔水巖層，因此，沉井基坑下部裂隙水與上覆第四系孔隙水也有水力聯(lián)系，且水位基本一致。

2基坑降水實踐及存在的問題

原設計方案的基本假設和概化是：把一期防滲圍堰與山體合圍的區(qū)域(或稱圍堰基坑)作為一個封閉的地質單元(圖1)。該單元的邊界條件，一期防滲圍堰可視為隔水邊界；據(jù)壓水試驗資料，山體和基坑底部基巖的透水率比較小，岸坡段巖體和基坑底部也進行了帷幕灌漿處理，因此，山體側面和基坑底部可以視為弱隔水邊界，即概化為零流量補給邊界。原方案將圍堰基坑概化為獨立的水文地質單元，只接受上部來自大氣降水的補給。

根據(jù)以上假設，圍堰基坑的涌水量，主要來自基坑地下水的儲存量和來自大氣降水的補給。再取安全系數(shù)(富余系數(shù))1.2，經(jīng)計算圍堰基坑設計涌水量為247×104m3。

按沉井設計工期12個月計算，原方案設計月抽水量約為21×104m3，日抽水量約0.7×104m3，平均每月降深3.8 m。

此方案經(jīng)過二期施工實踐證實，設計的涌水量遠遠不能滿足工程施工的要求，設計計算值遠小于實際抽出的水量。從實際情況來看，月排水量達到原設計方案的6～8倍時，基坑水位才能夠逐步下降，達到月平均降深3.6 m左右。

從上可以看出，原設計的圍堰基坑儲量和圍堰基坑面上的降水入滲量只占基坑涌水量的少部分，圍堰基坑水量來源和計算參數(shù)的求取還存在不足。

3基坑涌水量數(shù)值模型建立

基坑涌水量計算常采用“大井法”，但遇到水文地質條件復雜等情況下，不適宜簡單、直接采用大井法，原因如下：(1)當研究區(qū)存在防滲墻時，勘察階段抽水試驗所求的影響半徑必然變化，其不能簡單代入裘布依井流公式進行求解；(2)本文研究的基坑形狀很不規(guī)則(“L”型)，同時含水層非均值，要進行條件概化，若概化不當，所求得的解可信度降低。(3)考慮到防滲墻在實際的施工中存在一些天然的或人為的不確定的因素，造成防滲墻的滲透性也不確定，在基坑降水中不能簡單的視為隔水邊界，必須對其的滲透性有充分的估計，而大井法很難推求其滲透系數(shù)的大小。

本文采用數(shù)值模擬技術可解決上述問題。

3.1模型建立技術的要點

地下水流數(shù)值模擬的核心是擬合，進行擬合的過程就是不斷修改參數(shù)，不斷修正數(shù)學模型，使數(shù)學模型與水流原型達到逼真的相似，即達到與原型等效[4]。對于地下水流問題，要使模型和原型擬合很好，必須對原型有足夠的認識，即對含水層的結構、邊界條件，地下水的補給、排泄、徑流以及運動規(guī)律等有正確的認識并進行一系列合理的概化。主要包括：模擬計算區(qū)域的確定，含水層的概化，邊界條件確定和概化，主要水文地質參數(shù)的確定等等?；佑克康臄?shù)值模擬方法是建立在研究區(qū)水文地質條件認識清楚、水文地質模型概化合理、數(shù)學模型建立和運用正確、所有參數(shù)具有代表性這些基坑之上，劃分水均衡區(qū)，通過水均衡的計算來預測基坑涌水量的全過程。

3.2水文地質條件的概化

研究區(qū)的主要含水、透水介質即T32巖組，其相對隔水底板是T31巖組，相對隔水頂板是T33巖組。研究區(qū)主要含水透水巖層是河床基巖(T32)和第四系覆蓋層(Qhal)。

3.2.1巖性概化和參數(shù)選取

1)第四系覆蓋層基巖概化(Qhal)

研究區(qū)第四系覆蓋層自地表往下主要可分為三大層有：砂卵礫石層8～13 m，下界面高程249～252 m；砂層厚14～24 m，下界面高程237～253 m;含崩塊石的砂卵礫石層厚8～22 m?；鶐r(T32)可分為兩層：強風化基巖層和弱風化基巖層。

2)水文地質參數(shù)分區(qū)

由上面的分析含水層組可概化為5大層，為了模型的精確性，第四系地層要進一步細分，分成8層。三疊系T32地層是含水透水地層，對基坑涌水有一定的影響，分成兩層。即模型細分為十層。分別給這十層含水層進行滲透系數(shù)分區(qū)巖土的滲透系數(shù)取決于巖土的成因和物質組成；基巖的滲透系數(shù)還取決于其風化程度，裂隙發(fā)育程度及連通性等。在對滲透系數(shù)的采用上，根據(jù)勘察資料水文鉆孔的抽水試驗和取土樣進行室內(nèi)滲透試驗以及利用的資料成果，結合工程經(jīng)驗，推薦滲透系數(shù)(K)值如表1 所示。

表1　模擬區(qū)滲透系數(shù)取值

3.2.2邊界條件

沉井群段基坑位于一期土石圍堰內(nèi)，其南是長江，其距基坑中心距離約100 m，認為此處地下水水位不會隨著基坑開挖而變化，視為定水頭邊界,取值270 m。距離基坑中心約240 m處的北面是山體，含有基巖裂隙水，視為定水頭邊界，取值271 m。東西方向距離基坑中心分別約為460，200，此兩處水位受長江影響較大，取水位值271～270 m。因此模型四周是定水頭邊界。模型上邊界接受大氣降水的補給，模型的下邊界是零流量邊界，位于河床基巖中等偏弱透水巖體之下，視為隔水底板。

3.3數(shù)學模型

模型采用潛水三維非穩(wěn)定流模擬[5,6]：

式中：ω為源匯相，u為給水度，K為滲透系數(shù)，Z為含水層標高；H為水頭；B1為第1類邊界；H1(x,y,z)為第1類邊界定水頭值。B2為模型下邊界零流量邊界。

3.4模型的建立

模擬區(qū)域為一矩形，面積1 000 m×400 m,模擬最大深度為176 m,頂部最大標高為280.26 m。 模型分為10層，每層100列，80行，共80 000個單元。

由于計算模型是三維的，事先知道每一層的頂面高程和底面高程。本文通過鉆孔勘察資料，記錄每個鉆孔點的坐標和同一層的高程，利用surfer 軟件克里格插值 ，得每一層的頂面高程和底面高程。在模擬時，賦給visualmodflow建立的模型。

3.5初始水頭的確定

模擬區(qū)的初始水頭有3種方法來確定：一是根據(jù)鉆孔中地下水的埋深資料，利用插值方法或趨勢面分析得到；二是先利用類比方法確定該地區(qū)的水力坡度，利用該水力坡度和現(xiàn)今的河流水位推算出各單元的地下水位埋深；三是對模型進行穩(wěn)定流模擬，得到后綴名.HDS水頭文件的水頭值，把此結果作為非穩(wěn)定流模型的初始水位。

根據(jù)鉆孔觀測水位資料分析，研究區(qū)的水位變化不大，基本保持在262 m左右，因此本次模擬取的初始水位262 m。

3.6對防滲墻圍堰的處理

考慮防滲墻圍堰的滲透性存在不確定性，在基坑降水設計中不能簡單將它視為隔水邊界。

在即坑降水三維數(shù)值模擬時，把它作為擋墻邊界(Wall)來處理，處理時，需要輸入兩個參數(shù)：擋墻的厚度和滲透系數(shù)。擋墻的厚度是一定的取0.8 m，滲透系數(shù)在一定范圍內(nèi)取值，來考慮對基坑涌水量的影響。

3.7補給邊界

MODFLOW用補給程序包(RCH)來處理地下水補給。補給程序包用來模擬地下水來自大氣的補給。大多數(shù)情況下，大氣補給表現(xiàn)為降水向地下水系統(tǒng)的入滲結果。本文模型的上邊界就是一個補給邊界，把大氣降水量按照模擬的時間逐月分配到模型中。根據(jù)基坑開挖的時間，模擬的時間定為一年，從第一年5月18日5時40分到第二年年5月18日。根據(jù)多年平均降水量得下表2。

表2　模擬時間內(nèi)每月平均降水量

表3　應力期的分布及其初末水頭值

3.8應力期

在利用visualmodflow進行非穩(wěn)定流模擬時，模擬時間分為多個時間段，在每個時間段內(nèi)系統(tǒng)上所有的應力(邊界條件、抽水速率等)保持不變的時間期，這個時間期叫應力期。

考慮到基坑開挖的進度，模擬期模擬12個時間段(共一年時間)，即12個應力期。每個應力期的初末水頭如表3所示。

圖2　均衡區(qū)的形狀大小示意圖

3.9模擬時對基坑的處理

把基坑底部的水頭設為變化的定水頭，具體輸入的數(shù)值按表3進行。另外把基坑所在的空間滲透系數(shù)設置很大，視為空的區(qū)域，地下水流進基坑時是無障礙流動。

4基坑涌水量的預測

模型是通過水均衡的計算來求基坑內(nèi)涌水量。在模擬過程中，模型將設置兩個均衡區(qū)，一個從基坑的周邊進水，設置均衡區(qū)1；另一個從基坑的底部進水設置均衡區(qū)2。均衡區(qū)的形狀等同于沉井群段基坑的形狀(見圖2)。

4.1預測方案的選取

考慮到防滲墻在實際施工中存在一些天然的或人為的不確定的因素，造成防滲墻的滲透性的不確定，因此模型在模擬時考慮防滲墻滲透性，制定了多個預測方案，見表4。

表4　預測方案

4.2方案的模擬結果

三種方案的模擬結果分別見表5、表6和表7。

表5　方案一涌水量計算結果記錄表

表6　方案二涌水量計算結果統(tǒng)計表

表7　方案三涌水量計算結果統(tǒng)計表

4.3模擬結果的分析

從表5，表6和表7的三種方案的模擬計算結果來看，如果防滲墻的防滲效果好，防滲墻的滲透系數(shù)(0.08 m/d)與基巖的滲透系數(shù)相當(如方案三)，從基坑四周流入量和底部流入量相差不大，且基坑總涌水量小。如果防滲圍堰效果不好，那么基坑四周進入的水量增幅很大，占總量的60%以上。

當防滲墻滲透系數(shù)k為8 m/d時，模擬計算出降深10 m時日排水量達到了約8×104m3；

當防滲墻滲透系數(shù)k為0.8 m/d時，模擬計算出降深10 m時日排水量達到了約5×104m3；

當防滲墻滲透系數(shù)k為0.08 m/d時，模擬計算出降深10 m時日排水量達到了約1.4×104m3。

而實際排水降到10 m時，大約需要100 d，則模擬的日排水量換算成月排水量分別為約240×104m3、152×104m3、42×104m3；而降水實際表明此時月排水在120～160×104m3之間時，沉井段基坑水位能夠逐步下降。則兩種結果對比較表明防滲墻的防滲效果達到了2～0.8 m/d。方案二與實際情況比較一致。

在此防滲效果下，模擬預測降深-45 m時，基坑涌水量的計算結果與后期降水實踐比較一致。

5結語

本文以長江流域某水利基坑降水為例，采用Visualmodflow三維數(shù)值模擬技術對基坑涌水量進行計算，并以此指導工程實踐，為降水設計與工程實踐提供依據(jù)。

(1)在基坑降水設計時，要全面的分析研究區(qū)的水文地質條件，當含水層出現(xiàn)防滲墻時，難以引入正確的“影響半徑”時，而采用數(shù)值模擬技術不受其限制。數(shù)值模擬技術不僅可以計算基坑涌水量，模擬出整個開挖過程中非穩(wěn)定流過程，而且還可以預測基坑周邊每一點的水位降深。

(2)防滲墻在實際施工中存在一些天然的或人為的不確定的因素，造成防滲墻的滲透性的不確定性，在基坑降水設計中不能簡單將它視為隔水邊界，必須對防滲墻滲透性的不確定性有充分的估計，采用大井法計算很難處理這一問題，而數(shù)值模擬技術能夠解決此問題。

本文利用數(shù)值模擬軟件(VisualModflow)對防滲墻的防滲效果進行了三種方案的非穩(wěn)定流模擬計算，結合降水實際情況，得出該工程防滲效果為2～0.8 m/d。

(3)建議在基坑涌水量預測計算中應采用多種模型計算，并結合地區(qū)經(jīng)驗類比，不應完全依賴某種計算結果。在施工過程密切注意涌水量情況，根據(jù)涌水量變化資料修正數(shù)值模型條件和參數(shù)，動態(tài)預測涌水量，以提高預測的可信度和安全性；在高滲透性含水層和補給豐富的基坑降水工程中，尤其是特大型深基坑，必須充分考慮防滲墻的滲透性，采用安全性的評估方法，進行多方案模擬計算，以保證降水的安全施工。

參考文獻

[1]李再興.有關基坑降水方法的探討.地下水[J].2008(2),72-74.

[2]李再興，梁杏，郭付三，等.大井法在基坑涌水量計算中的應用[J].人民長江.VOL. 40(15):50-55.

[3]陳崇希，成建梅，孫紅林.深基坑地下水的排降——地面沉降模型的探討[J].勘察科學技術.1998，(2)：3-7.

[4]李艷梅.遼河油田地區(qū)地下水水質模擬與預測[D].吉林大學碩士學位論文.2003.

[5]陳崇希，林敏.地下水動力學[M].武漢：中國地質大學出版社.1996.

[6]陳崇希，唐仲華.地下水流動問題數(shù)值方法[M].武漢：中國地質大學出版社.1990.

Application of Numerical Modeling to Dewatering Design for Foundation Pit

LI Zai-xing，LIJing，SONG Peng-fei，YANGFang

(Institute Of Surveying Mapping And Geoinformation of Henan ,Zhengzhou 450006, Henan)

Abstract:The total water displacement is obtained through the generalized large well method in the dewatering design for foundation pits,but the hydrology geology condition in the research area is more complicated, heterogeneity of aquifer, exist seepage control cofferdam and filter efficiency is uncertain, making use of " generalized large " resolution method to solve stream amount of water, the condition of generalization is more, usually cause the credibility of compute result lower. So my paper have used VisualModflow to simulate the amount of water. Build on the suitable hydrology geology concept model and assurance mathematics model, got stream amount of water of foundation pit dewatering by the calculation unstabilized flow. The influence of impervious effect of penetrating wall for stream amount of water is very greatly, take three kinds of projects to calculation the impervious effect, combining with practical engineering, have inquired into the permeability coefficient of seepage proof curtain.

Key words:Foundation pit dewatering；numerical modeling and seepage proof curtain

[中圖分類號]P641.4+1

[文獻標識碼]A

[文章編號]1004-1184(2016)02-0022-04

[作者簡介]李再興(1981-)，男，安徽安慶人，工程師，主要從事工程地質、水文地質與環(huán)境地質的勘測與研究工作。

[收稿日期]2015-11-16