李 磊 陳 干 唐 沛 姚德華 黨峰榮

(1.中鐵第五勘察設(shè)計(jì)院集團(tuán)有限公司，102627，北京;2.南京大學(xué)地球科學(xué)與工程學(xué)院，210023，南京∥第一作者，高級(jí)工程師)

滲透系數(shù)是基坑疏降水、地下水滲流控制、地質(zhì)環(huán)境影響評(píng)價(jià)等水文地質(zhì)計(jì)算中的關(guān)鍵參數(shù)［1-2］之一。獲取水文地質(zhì)參數(shù)的方法很多，如:直接采用《水文地質(zhì)手冊(cè)》［3］及GB 50287—1999《水利水電工程地質(zhì)勘察規(guī)范》［4］中羅列的經(jīng)驗(yàn)值;通過(guò)室內(nèi)試驗(yàn)進(jìn)行測(cè)定［3］;通過(guò)現(xiàn)場(chǎng)抽水試驗(yàn)或注水試驗(yàn)，根據(jù)抽、注水量和監(jiān)測(cè)水位的動(dòng)態(tài)關(guān)系，通過(guò)解析公式或數(shù)值法進(jìn)行反演計(jì)算［5］。文獻(xiàn)［3-4］中給定的不同巖性的經(jīng)驗(yàn)參數(shù)往往是一個(gè)取值范圍;室內(nèi)試驗(yàn)測(cè)定的參數(shù)應(yīng)用到場(chǎng)地工程存在尺度效應(yīng)問(wèn)題;而解析法引入了大量的假設(shè)條件，如含水層幾何形狀規(guī)則、均質(zhì)各向同性、邊界條件單一等，其使用條件苛刻，難以滿足實(shí)際復(fù)雜水文地質(zhì)條件下參數(shù)反演的要求。采用上述方法求解得到的水文地質(zhì)參數(shù)結(jié)果存在一定偏差，難以滿足工程計(jì)算的要求［6］。因此，有必要采用解析法與數(shù)值法相結(jié)合的方法進(jìn)行含水層水文地質(zhì)參數(shù)的反演［7-9］。通過(guò)對(duì)多種方法反演的水文地質(zhì)參數(shù)進(jìn)行對(duì)比，相互驗(yàn)證水文地質(zhì)參數(shù)的合理性，以降低單種方法求解水文地質(zhì)參數(shù)的不確定性。這對(duì)于水文地質(zhì)條件的正確判別，以及地下水控制方案的合理制定具有重要的工程應(yīng)用意義。

哈爾濱地鐵4號(hào)線一期第1標(biāo)段全段位于松花江以北、阿什河以南的漫灘區(qū)，漫灘區(qū)含水層厚度較大。以往諸多地下水控制工程，如基坑降水、地鐵抗浮等工程，更多地憑借工程人員的經(jīng)驗(yàn)進(jìn)行設(shè)計(jì)，造成了大量的工程地質(zhì)問(wèn)題［10-12］。因此，本文分別采用解析法和數(shù)值法［8］求解研究區(qū)段主要含水層的水文地質(zhì)參數(shù)，如滲透系數(shù)K、導(dǎo)數(shù)系數(shù)T等，對(duì)哈爾濱漫灘區(qū)地鐵工程地下水滲控具有借鑒意義。

1 科技五街站抽水試驗(yàn)工程概況

1.1 抽水試驗(yàn)場(chǎng)地條件

抽水試驗(yàn)場(chǎng)地位于哈爾濱市松北區(qū)地鐵4號(hào)線科技五街站附近。抽水試驗(yàn)井采用“5抽+3觀測(cè)”的組合。其中，5口抽水井沿東西向直線布置，間距均為10 m;3口觀測(cè)井與抽水井連線垂直布置，觀測(cè)井1(觀1)與抽水井3(主井3)間距為10 m，與觀測(cè)井2(觀2)間距為20 m，與觀測(cè)井3(觀3)間距為30 m。抽水試驗(yàn)井位置詳見(jiàn)圖1。單井抽水試驗(yàn)設(shè)計(jì)參數(shù)見(jiàn)表1。

圖1 科技五街站抽水試驗(yàn)井布置圖Fig.1 Location of pumping test wells around Kejiwujie Station

表1 抽水試驗(yàn)場(chǎng)地單井抽水試驗(yàn)設(shè)計(jì)參數(shù)表Tab.1 Design parameters table of single well pumping test Station

科技五街站抽水試驗(yàn)井揭示，研究區(qū)段的地層主要為厚45～55 m的第四系松散堆積物，為現(xiàn)代河漫灘相堆積，分布較為廣泛。抽水試驗(yàn)場(chǎng)地自上而下主要為:地表上覆2.0～6.0 m的亞黏土或亞砂土，低洼處為淤泥質(zhì)黏土;全新統(tǒng)下段細(xì)砂、中粗砂;下更新統(tǒng)為中粗砂、砂卵礫石層;下覆白堊系嫩江組泥巖、砂巖及泥質(zhì)砂巖，構(gòu)成了松散堆積物底板。

研究區(qū)段含水層為典型的非均質(zhì)微承壓含水層，地下水主要補(bǔ)給來(lái)源于大氣降水的滲入補(bǔ)給、豐水期松花江江水側(cè)向滲入補(bǔ)給，以及地下水的側(cè)向徑流補(bǔ)給;地下水的排泄主要以蒸發(fā)為主，其次為地下水側(cè)向徑流排泄，枯水期排泄方式為向松花江排泄和人工開采。地下水位埋深僅為2.0～4.5 m，年變幅一般為1.0～3.0 m。

1.2 抽水試驗(yàn)數(shù)據(jù)分析

本次抽水試驗(yàn)分單井抽水試驗(yàn)和群井抽水試驗(yàn)兩種類型。采用單井抽水試驗(yàn)分別進(jìn)行3次不同水位降深S的抽水。其中，最大水位降深發(fā)生在采用最大抽降能力抽水時(shí)，其余2次水位降深分別為最大水位降深值的1/3和2/3。3次抽水落程采用的抽水流量Q分別為221.1 m3/h、154.0 m3/h和78.5 m3/h。

群井抽水試驗(yàn)中，5口抽水井的抽水流量全部采用200 m3/h，并對(duì)其進(jìn)行一次水位降深的抽水，同時(shí)讀取各觀測(cè)井水位。現(xiàn)場(chǎng)均采用穩(wěn)定流進(jìn)行抽水，為了滿足非穩(wěn)定流的計(jì)算要求，結(jié)合非穩(wěn)定流要求進(jìn)行數(shù)據(jù)觀測(cè)。單井和群井抽水試驗(yàn)結(jié)果如圖2～3所示。

圖2 單井和群井抽水試驗(yàn)Q-t關(guān)系圖Fig.2 Q-t relationship between single well and group well pumping tests

2 采用解析法求解抽水試驗(yàn)場(chǎng)地水文地質(zhì)參數(shù)

2.1 計(jì)算方法

本次非穩(wěn)定流水文地質(zhì)參數(shù)的計(jì)算，分別選取穩(wěn)定流Dupuit計(jì)算公式、非穩(wěn)定流Theis計(jì)算公式和水位恢復(fù)法反演水文地質(zhì)參數(shù)［3］。

2.1.1 穩(wěn)定流Dupuit計(jì)算公式

抽水試驗(yàn)采用完整井進(jìn)行抽水。場(chǎng)地含水層近似水平、均質(zhì)、各向同性，等水頭面是共軸的圓柱面，和過(guò)水?dāng)嗝嬉恢?。抽水井附近達(dá)到穩(wěn)定或近似穩(wěn)定后，可采用Dupuit公式進(jìn)行計(jì)算。

圖3 單井和群井抽水試驗(yàn)S-t關(guān)系圖Fig.3 S-t relationship between single well and group well pumping tests

式中:

K——滲透系數(shù)，m/d;

H——微承壓含水層厚度，m;

Q——抽水井出水量，m3/d;

S1——1號(hào)觀測(cè)井的水位降深，m;

S2——2號(hào)觀測(cè)井的水位降深，m;

r1——1號(hào)觀測(cè)井與抽水井的距離，m;

r2——2號(hào)觀測(cè)井與抽水井的距離，m。

2.1.2 非穩(wěn)定流Theis計(jì)算公式和水位恢復(fù)法

利用群孔(完整井)在抽水試驗(yàn)場(chǎng)地進(jìn)行定流量非穩(wěn)定流抽水試驗(yàn)。試驗(yàn)中，微承壓含水層近似均質(zhì)、等厚、側(cè)向無(wú)限延伸，含水層中水頭損失符合達(dá)西定律。采用Theis計(jì)算公式進(jìn)行抽水時(shí)段水文地質(zhì)參數(shù)估算;或利用抽水恢復(fù)階段的觀測(cè)資料，采用水位恢復(fù)法進(jìn)行水文地質(zhì)參數(shù)計(jì)算。

2.1.2.1 非穩(wěn)定流Theis計(jì)算公式

非穩(wěn)定流Theis計(jì)算公式如下:

式中:

W(u)——泰斯井函數(shù)，u=r2/(4 at);

S (r，t)——距離抽水井r處，t時(shí)刻的水位降深值，m;

a——壓力系數(shù)，m2/d;

T——導(dǎo)水系數(shù)，m2/d;

t——抽水延續(xù)時(shí)間，s。

2.1.2.2 水位恢復(fù)法

水位恢復(fù)法的計(jì)算公式如下:

令:

則式(3)變?yōu)?

式中:

h——井內(nèi)恢復(fù)水位高度，m;

ha——抽水穩(wěn)定時(shí)井內(nèi)水位高度，m;

R——影響半徑，m;

r——抽水井半徑，m。

停抽穩(wěn)定后，與t1、t2相應(yīng)的h1、h2分別為:

求解式(5)—式(6)，得到:

2.2 計(jì)算結(jié)果分析

基于不同水位降深的單井抽水試驗(yàn)和群孔非穩(wěn)定流抽水試驗(yàn)數(shù)據(jù)，通過(guò)Aquifer test軟件，采用穩(wěn)定流Dupuit計(jì)算公式、非穩(wěn)定流Theis計(jì)算公式和水位恢復(fù)法分別反演水文地質(zhì)參數(shù)，計(jì)算結(jié)果如表2所示。非穩(wěn)定流Theis函數(shù)配線見(jiàn)圖4。通過(guò)不同解析法求得的滲透系數(shù)為35.9～71.3 m/d，導(dǎo)水系數(shù)為1.61×103～3.21×103m2/d。

表2 各解析法下抽水試驗(yàn)場(chǎng)地水文地質(zhì)參數(shù)統(tǒng)計(jì)表Tab.2 Statistical table of hydrogeological parameters of pumping test sites

圖4 非穩(wěn)定流Theis函數(shù)配線圖Fig.4 Distribution diagram of unsteady flow Theis function wirings

3 采用數(shù)值法求解抽水試驗(yàn)場(chǎng)地水文地質(zhì)參數(shù)

為了更加準(zhǔn)確地刻畫抽水試驗(yàn)場(chǎng)地的水文地質(zhì)特征，本文利用FEFLOW軟件進(jìn)行三維地下水流建模，并采用抽水試驗(yàn)動(dòng)態(tài)數(shù)據(jù)對(duì)模型進(jìn)行識(shí)別，通過(guò)反演得到該場(chǎng)地最佳的水文地質(zhì)參數(shù)。

3.1 抽水試驗(yàn)場(chǎng)地地下水流概念模型

根據(jù)抽水試驗(yàn)場(chǎng)地的水文地質(zhì)條件，對(duì)含水層結(jié)構(gòu)進(jìn)行概化。地下水流概念模型的地層由上到下劃分為3層:

1)第四系覆蓋層(素填土、粉土、粉質(zhì)黏土層)厚3.0 m;

2)第四系孔隙微承壓層厚39.0 m;

3)基巖隔水底層厚6.8 m。

由于研究區(qū)段范圍較小，且模擬區(qū)主要含水層的巖性變化較小，因此將含水層概化為均質(zhì)、各向異性的介質(zhì)。

為減少邊界選取的隨意性給計(jì)算結(jié)果帶來(lái)誤差，通過(guò)解析法試算抽水試驗(yàn)的影響半徑。以抽水井(主井1—主井5)為中心，向四周長(zhǎng)、寬方向各延伸1 400 m的區(qū)域?yàn)槟M區(qū)。排除抽水井對(duì)邊界的影響，模擬區(qū)范圍為2 853 m×2 853 m。

該模型共分為3層。其中，上邊界為粉質(zhì)黏土層，為自由面，由于不考慮大氣降水入滲、蒸發(fā)排泄等因素，故上邊界不作處理，默認(rèn)為隔水邊界;下邊界為基巖底板，將其設(shè)置為隔水邊界;模擬區(qū)周界未受到抽水影響;將側(cè)向模擬邊界設(shè)為定水頭補(bǔ)給邊界;給定水頭值為現(xiàn)場(chǎng)實(shí)測(cè)的含水層初始水位111.22 m。

3.2 抽水試驗(yàn)場(chǎng)地地下水流數(shù)學(xué)模型

根據(jù)上述建立的概念模型，在不考慮水的密度變化的前提下，將該模型概化為三維微承壓含水層中的非穩(wěn)定流數(shù)學(xué)模型:

式中:

Ω——滲流區(qū)域;

h——含水層的水位標(biāo)高;

Kx、Ky、Kz——分別為x、y、z方向的滲透系數(shù);

μ——含水層的給水度或者貯水系數(shù);

ε——含水層的源匯項(xiàng);

h0——含水層的初始水位分布;

?！獫B流區(qū)域的給定流量邊界，包括滲流區(qū)域的側(cè)向流量或隔水邊界和含水層底部的隔水邊界;

Kn——邊界面法向方向的滲透系數(shù);

q(x，y，z，t)——定義為給定流量邊界的單位面積流量，本次為隔水邊界，設(shè)為0。

3.3 抽水試驗(yàn)場(chǎng)地地下水流數(shù)值模型

圖5所示為抽水試驗(yàn)場(chǎng)地三維數(shù)值模型圖。該模型在豎向分為3層，在橫向采用三角形網(wǎng)格剖分(抽水井和觀測(cè)井附近網(wǎng)格需加密)，最小的網(wǎng)格面積為0.04 km2，每層共剖分為11 438個(gè)網(wǎng)格，形成5 852個(gè)節(jié)點(diǎn)。

圖5 抽水試驗(yàn)場(chǎng)地三維地下水流數(shù)值模型圖Fig.5 Schematic diagram of 3D ground water flow numerical model at pumping test spots

非主要含水層初始滲透系數(shù)和貯水系數(shù)等參數(shù)的取值主要根據(jù)含水層巖性，并結(jié)合《水文地質(zhì)手冊(cè)》進(jìn)行查取，而主要含水層的初始參數(shù)則通過(guò)解析法進(jìn)行求取(見(jiàn)表2)。本次模擬區(qū)巖性分布均勻，故在橫向上不再對(duì)其進(jìn)行參數(shù)分區(qū)。

3.4 水文地質(zhì)參數(shù)模擬反演結(jié)果分析

選用群孔抽水試驗(yàn)時(shí)段為模擬期，對(duì)孔隙含水層水文地質(zhì)參數(shù)進(jìn)行反演，其中，時(shí)間步長(zhǎng)通過(guò)誤差迭代自行控制。在參數(shù)反演中，依據(jù)觀測(cè)井水文觀測(cè)值和計(jì)算水位的差值來(lái)調(diào)整滲透系數(shù)。通過(guò)反復(fù)調(diào)試，得到該參數(shù)的最佳擬合效果，如圖6所示。通過(guò)數(shù)值法反演所得的含水層滲透系數(shù)為44.5 m/d。

圖6 基于數(shù)值法的地下水水位擬合曲線和流場(chǎng)圖Fig.6 Groundwater level fitting curves and flow field diagram based on numerical method

通過(guò)構(gòu)建數(shù)值模型，可模擬在不同抽水時(shí)刻試驗(yàn)區(qū)地下水流動(dòng)態(tài)。如在抽水開始后1 800 min，各抽水主孔產(chǎn)生的漏斗已相互疊加，在觀測(cè)孔處模擬水位與觀測(cè)水位誤差為0.1～0.2 m，精確地反映了群井抽水過(guò)程中地下水位的動(dòng)態(tài)變化。

本次采用解析法求解得到的滲透系數(shù)為35.9～71.3 m/d，平均值為54.1 m/d;采用數(shù)值法反演得到的滲透系數(shù)為44.5 m/d?？紤]在實(shí)際工程中，解析法對(duì)工程條件進(jìn)行了假設(shè)和簡(jiǎn)化，且擬合的隨意性較大，而數(shù)值法能更精確地描述含水層的實(shí)際結(jié)構(gòu)和地下水的真實(shí)運(yùn)動(dòng)狀態(tài)。

4 結(jié)語(yǔ)

1)利用解析法計(jì)算得到孔隙微承壓含水層滲透系數(shù)為35.9～71.3 m/d，平均值為54.1 m/d。

2)數(shù)值法模擬了群孔非穩(wěn)定流抽水試驗(yàn)過(guò)程中地下水流動(dòng)態(tài)過(guò)程，通過(guò)擬合得到主要含水層的滲透系數(shù)為44.5 m/d。

3)解析法和數(shù)值法計(jì)算結(jié)果的相互驗(yàn)證，大大增加了水文地質(zhì)參數(shù)反演結(jié)果的可靠度，可為后期地鐵工程地下水滲控方案的制定和實(shí)施提供參數(shù)。