趙激光

（中鐵上海設(shè)計(jì)院集團(tuán)有限公司，上海 200070）

0 引言

上海軟土地區(qū)修建地鐵車(chē)站，基坑止水帷幕采用地下連續(xù)墻。地下連續(xù)墻能夠阻擋潛水含水層對(duì)地鐵車(chē)站基坑穩(wěn)定的影響。因此，地下水對(duì)基坑的影響主要考慮承壓含水層。承壓水視為無(wú)補(bǔ)給來(lái)源的無(wú)限含水層，短期抽水試驗(yàn)，承壓水水位不斷下降，試驗(yàn)是非穩(wěn)定流抽水。此外，承壓含水層上下存在弱透水層，抽水過(guò)程中會(huì)出現(xiàn)越流補(bǔ)給。地下連續(xù)墻阻斷地下水的水平向流動(dòng)，地下水流動(dòng)以垂直向?yàn)橹鲗?dǎo)，抽水井附近及止水帷幕底端局部存在不適于達(dá)西定律的情況。因此，對(duì)于承壓水層進(jìn)行抽水試驗(yàn)研究是一個(gè)復(fù)雜的問(wèn)題，采用解析法計(jì)算的抽水試驗(yàn)參數(shù)比較繁瑣且誤差較大。

隨著計(jì)算機(jī)技術(shù)的發(fā)展，采用計(jì)算機(jī)對(duì)地下水滲流進(jìn)行數(shù)值模擬，使得復(fù)雜含水層的抽水試驗(yàn)的問(wèn)題得到解決。地下水?dāng)?shù)值模擬模型MODFLOW是由美國(guó)地質(zhì)調(diào)查局的McDonald和Harbaugh于20世紀(jì)80年代開(kāi)發(fā)的用于孔隙介質(zhì)的三維有限差分地下水流數(shù)值模擬模型[1]。Visual Modflow 是由加拿大Waterloo 水文地質(zhì)公司在Modflow的基礎(chǔ)上開(kāi)發(fā)研制的，是目前國(guó)際上流行的三維地下水流和溶質(zhì)運(yùn)移模擬評(píng)價(jià)的標(biāo)準(zhǔn)可視化專(zhuān)業(yè)軟件系統(tǒng)[2]。Visual Modflow可以仿真模擬深基坑降水過(guò)程中地下水位的變化情況[3]。Visual Modflow 相較其他模型而言，可視化功能強(qiáng)大，求解方法簡(jiǎn)單，適用范圍廣泛，數(shù)值模擬能力出色且操作簡(jiǎn)單[4]。

基于Visual Modflow軟件進(jìn)行地下水滲流數(shù)值模擬，反演出合理的計(jì)算模型水力學(xué)參數(shù)，逐漸成為降水試驗(yàn)的一種常規(guī)手段，廣泛應(yīng)用于軟土地區(qū)地鐵基坑水文地質(zhì)勘察中。本文就該軟件在群孔抽水試驗(yàn)中的運(yùn)用進(jìn)行介紹和分析。

1 工程概況

擬建地鐵車(chē)站工程全長(zhǎng)536m，站型為地下三層，島式車(chē)站，主體結(jié)構(gòu)為單柱雙跨現(xiàn)澆鋼筋混凝土結(jié)構(gòu)，圍護(hù)結(jié)構(gòu)為1200mm厚地下連續(xù)墻?；硬捎妹魍陧樧鞣ㄊ┕?，基坑開(kāi)挖深度約25.87～27.82m。

2 水文地質(zhì)條件

本工程位于長(zhǎng)江三角洲沖積平原區(qū)，根據(jù)勘察結(jié)果，場(chǎng)地揭露地層為第四紀(jì)全新世（Q4）濱?！涌?、晚更新世（Q3）濱?！珳\海相沉積地層，主要由黏性土、粉性土及粉砂組成。

擬建工程場(chǎng)地屬于太湖流域黃浦江水系。地下水根據(jù)成因類(lèi)型可分為孔隙潛水和（微）承壓水。潛水主要賦存于淺部填土、黏性土、粉性土中。微承壓水主要分布在⑤2層黏質(zhì)粉土中，⑤2層與⑦2層水力聯(lián)系弱，基坑止水帷幕埋深在⑤2層層底，基坑范圍采用疏干降水處理。承壓水主要分布在⑦2層粉砂、⑧2層粉質(zhì)黏土夾粉砂和⑨層粉砂中。⑦2層與⑧2層兩層承壓水之間⑧1層作為隔水層。⑦2層與⑧2層水力聯(lián)系弱。本次群孔抽水試驗(yàn)對(duì)相鄰含水層越流補(bǔ)給影響的⑧2層水文地質(zhì)參數(shù)進(jìn)行確定。

3 抽水試驗(yàn)

抽水試驗(yàn)設(shè)置82Y1～82Y3三個(gè)抽水試驗(yàn)孔，82G1～82G4四個(gè)觀(guān)測(cè)孔，孔深為65～68m?？孜徊贾靡?jiàn)圖1。

圖1 抽水孔及觀(guān)測(cè)孔布置示意圖

群孔抽水期間，82Y1～82Y3作為抽水孔，剩余為觀(guān)測(cè)孔。抽水歷時(shí)約8d，各觀(guān)測(cè)孔水位穩(wěn)定，滿(mǎn)足群孔抽水試驗(yàn)停抽要求。

4 數(shù)值模擬反演

4.1 含水層的結(jié)構(gòu)特征

根據(jù)抽水試驗(yàn)數(shù)據(jù)建立數(shù)值模型，考慮基坑面積及抽水影響半徑，確定模型尺寸800×800（m）。行列間距按照20m間距設(shè)置，核心區(qū)按照2m進(jìn)行加密，模型層厚按照100m設(shè)置。網(wǎng)絡(luò)最外邊緣按定水頭考慮。網(wǎng)格立體剖分圖見(jiàn)圖2。

圖2 離散模型網(wǎng)絡(luò)三維圖

4.2 模型調(diào)整優(yōu)化

根據(jù)抽水試驗(yàn)所得參數(shù)，對(duì)模型進(jìn)行初步賦值。模型的建立往往會(huì)忽略一些次要因素，因此模型不能準(zhǔn)確代表實(shí)際系統(tǒng)，此時(shí)需要通過(guò)模擬結(jié)果與實(shí)測(cè)結(jié)果進(jìn)行反演，實(shí)際模型參數(shù)最終通過(guò)不斷地調(diào)整及優(yōu)化確定。

4.3 水力特征

天然孔隙含水層中的地下水流的Reynolds數(shù)和裂隙中地下水流的水力坡度，遠(yuǎn)小于臨界Reynolds數(shù)和臨界水力坡度，天然地下水多處于層流狀態(tài)[5]。地下含水層隨地層呈層狀分布，不同的含水層滲透系數(shù)不同。因此，一般建立模型時(shí)，可以根據(jù)場(chǎng)地實(shí)際地層確定含水層。地層相同，地質(zhì)參數(shù)相近，可視為均質(zhì)含水層，在水平方向上視為各向同性介質(zhì)。垂直方向上地層不同，滲透系數(shù)各異，例如黏性土滲透系數(shù)很小，視為隔水層，阻礙地下水的垂向運(yùn)動(dòng)，因此垂向視為各向異性介質(zhì)。

地下水的流動(dòng)在含水層中以水平向?yàn)橹?，垂向?yàn)檩o。相鄰含水層為透水層或弱透水層，在抽水的情況下，會(huì)產(chǎn)生明顯的層間越流，形成三維地下水滲流場(chǎng)。地下水流是連續(xù)的，取地下水微單元體進(jìn)行分析，根據(jù)地下水質(zhì)量守恒關(guān)系，進(jìn)入離散的單元體水的質(zhì)量與流出單元體質(zhì)量差等于單元體內(nèi)液體質(zhì)量的變化。根據(jù)地下水流連續(xù)方程及邊界條件，建立地下水三維非穩(wěn)定流模型。

4.4 模擬期及應(yīng)力期確定

根據(jù)有效應(yīng)力原理，隨著水位的下降，水的應(yīng)力減小，土的應(yīng)力增加，總應(yīng)力保持不變。抽水試驗(yàn)抽水時(shí)間為8d，因此群孔試驗(yàn)?zāi)M期及應(yīng)力期為8d。為了準(zhǔn)確反應(yīng)模擬曲線(xiàn)變化，本次分三個(gè)周期，每個(gè)周期外部源匯項(xiàng)的強(qiáng)度保持不變。

4.5 源匯項(xiàng)

源匯項(xiàng)指承壓含水層中流入（源）或流出（匯）項(xiàng)。在Visual Modflow軟件中可以對(duì)抽水試驗(yàn)孔、觀(guān)測(cè)孔的參數(shù)進(jìn)行設(shè)置。抽水孔可以設(shè)置過(guò)濾管的位置，觀(guān)測(cè)孔可以設(shè)置水位值。研究區(qū)范圍邊界在抽水試驗(yàn)影響半徑之外，邊界條件即初始水位，為固定值。

將抽水試驗(yàn)得到的涌水量值輸入模擬系統(tǒng)，反演獲得模擬水位值，通過(guò)模擬水位值與實(shí)測(cè)值比較，判斷模擬是否有效。

5 模擬結(jié)果分析

通過(guò)三維數(shù)值反演分析，獲取的模型參數(shù)如表1所示。

表1 群孔試驗(yàn)?zāi)Ｐ蛯訁?shù)

通過(guò)群孔抽水試驗(yàn)4個(gè)觀(guān)測(cè)孔的實(shí)測(cè)降深曲線(xiàn)與三維數(shù)值模擬取得的計(jì)算降深曲線(xiàn)進(jìn)行對(duì)比分析，觀(guān)測(cè)孔水位對(duì)比見(jiàn)圖3～6（點(diǎn)：代表實(shí)測(cè)數(shù)據(jù)，線(xiàn)：代表理論計(jì)算數(shù)據(jù)）。

圖3 82G1觀(guān)測(cè)孔實(shí)測(cè)水位曲線(xiàn)與模擬水位曲線(xiàn)對(duì)比圖

圖4 82G2觀(guān)測(cè)孔實(shí)測(cè)水位曲線(xiàn)與模擬水位曲線(xiàn)對(duì)比圖

圖5 82G3觀(guān)測(cè)孔實(shí)測(cè)水位曲線(xiàn)與模擬水位曲線(xiàn)對(duì)比圖

根據(jù)以上圖像分析可知，三維數(shù)值模擬取得的計(jì)算降深曲線(xiàn)與實(shí)測(cè)水位曲線(xiàn)較接近，說(shuō)明模型建立較為合理。數(shù)值模擬水位降深值與實(shí)測(cè)水位降深值見(jiàn)表2。

圖6 82G4觀(guān)測(cè)孔實(shí)測(cè)水位曲線(xiàn)與模擬水位曲線(xiàn)對(duì)比圖

表2 數(shù)值模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)對(duì)比

理論上觀(guān)測(cè)孔距離抽水孔越遠(yuǎn)，越接近初始水位，模擬值也一定越接近實(shí)測(cè)值。本次模擬從82G1觀(guān)測(cè)孔的模擬降深值與實(shí)測(cè)水位降深的偏差10.13%，到82G4觀(guān)測(cè)孔的模擬降深值與實(shí)測(cè)水位降深的偏差0.68%，呈逐漸減小的趨勢(shì)也證實(shí)了這一點(diǎn)。實(shí)測(cè)水位降深值與模擬水位降深值平均偏差2.04%，滿(mǎn)足工程需要。

6 結(jié)束語(yǔ)

通過(guò)三維數(shù)值計(jì)算得⑧2層滲透系數(shù)Kh=0.18m∕d，Kv=0.01m∕d，貯水率Ss=3.31×10-51∕m。數(shù)值模擬水位變化與實(shí)測(cè)抽水試驗(yàn)水位變化相符，說(shuō)明模型建立較為合理，反演參數(shù)較為準(zhǔn)確，反演計(jì)算的水文地質(zhì)參數(shù)能代表場(chǎng)地水文地質(zhì)特征。