孫寶亮，梁俊紅，崔衛(wèi)利，吳義松，馬廣智

(1.東北大學(xué)資源與土木工程學(xué)院，沈陽 110004; 2.中南冶金地質(zhì)研究所，湖北 宜昌 310014;3.中冶北方工程技術(shù)有限公司，遼寧 大連 116600)

?

孫寶亮1，梁俊紅1，崔衛(wèi)利1，吳義松2，馬廣智3

(1.東北大學(xué)資源與土木工程學(xué)院，沈陽 110004; 2.中南冶金地質(zhì)研究所，湖北 宜昌 310014;3.中冶北方工程技術(shù)有限公司，遼寧 大連 116600)

文章提出并采用水力坡度分析方法，對潛水回水初始階段進(jìn)行了定性-半定量分析，為了解潛水回水過程提供了新的視角。通過分析潛水回水初始時刻的水力坡度，河渠側(cè)的界面被分為2段：河間潛水位之上的鉛直段和之下的水平段。水平段上點(diǎn)的水力坡度最大且相等，方向?yàn)樗较蛴?；鉛直段上點(diǎn)的水力坡度方向趨于鉛直向下，越向上水力坡度越小。經(jīng)過微小時間段Δt后，水力坡度均相應(yīng)變小，鉛直段相鄰點(diǎn)間地下水運(yùn)動發(fā)生沖突，水力坡度方向由趨于鉛直向下向右偏轉(zhuǎn)。將水力坡度分析方法用于分析潛水井流的過程和特征，水躍的成因可解釋為：在相同水力坡度和過水?dāng)嗝婷娣e下，井內(nèi)大空間水流速度遠(yuǎn)大于孔隙中的潛水流速度，井內(nèi)水位急速下降，從而形成水躍；水躍越大，進(jìn)入井內(nèi)的流量就越大。

潛水；回水；水力坡度；定性分析；井；水躍(滲出面)

0 引言

現(xiàn)有大多數(shù)地下水動力學(xué)教材中，潛水回水非穩(wěn)定運(yùn)動數(shù)學(xué)模型基于布西涅斯克方程建立，之后線性化布西涅斯克方程，再用有限傅里葉正弦變換得到該模型的級數(shù)解析解[1-5]。問題在于，布西涅斯克方程應(yīng)用的是裘布依假設(shè)，而實(shí)際上在河附近的垂直流動分量大且不可忽視，裘布依假設(shè)并不成立[6-7]。河水位在大幅度變動的情況下，潛水厚度(M)不能代表真實(shí)水位高度，導(dǎo)水系數(shù)(T=KM)與實(shí)際不符，因而線性化的布西涅斯克方程解的誤差較大[7-8]。鑒于此，本文提出利用水頭和水力坡度分析方法，對潛水回水初始階段非穩(wěn)定運(yùn)動進(jìn)行定性-半定量分析；由于潛水井流與潛水回水初始階段有相似之處，可將水力坡度分析方法用于分析描述潛水井流，有助于深入了解潛水井流的過程和特征，獲得較為準(zhǔn)確的流網(wǎng)，從而解釋水躍產(chǎn)生的原因。水頭和水力坡度分析方法為了解潛水回水過程及理解不同層位潛水運(yùn)動的狀態(tài)和特征提供了新的視角。

1 潛水回水的物理模型

河水面的抬高、河渠地表水的側(cè)滲作用來補(bǔ)充地下水，會引起潛水位相應(yīng)抬高，稱為潛水回水[5]。其物理模型為：含水層均質(zhì)各向同性，隔水底板水平；初始狀態(tài)左河渠與河間潛水有相同水位h2；左河渠水位汛時上升為定值h1，發(fā)生瞬時回水，河渠間潛水非穩(wěn)定運(yùn)動；不考慮地表水垂向補(bǔ)給和非飽和水的運(yùn)動。進(jìn)而將此模型虛擬為：受可移動的鉛直無厚度隔水板PP′的阻擋，左河渠不能瞬時回水，當(dāng)迅速抽出隔水板PP′后，才會發(fā)生瞬時回水(圖1)。

圖1 潛水瞬時回水的物理模型(據(jù)文獻(xiàn)[5]修改)Fig.1 The conceptual model of a stream recharging into unconfined aquifer

2 潛水回水初始階段的定性-半定量分析

用水力坡度法定性-半定量分析潛水回水初始階段非穩(wěn)定運(yùn)動的過程。①比較隔水板PP′左側(cè)點(diǎn)A，B，C，D，E，F(xiàn)和右側(cè)水平對應(yīng)點(diǎn)A′，B′，C′，D′，E′，F(xiàn)′共12個點(diǎn)測壓水頭的大小(圖1，A，D′分別在隔水板PP′左、右側(cè)的液面位置，F(xiàn)在隔水底板上，B，C在AD間且B在C之上，E在DF間)；②在迅速抽出隔水板PP′的瞬間(初始時刻)，A，B，C，D，E，F(xiàn)點(diǎn)水力坡度的大小和方向如何？③微小時間段Δt后，A，B，C，D，E，F(xiàn)點(diǎn)的水力坡度的大小和方向如何變化？

2.1 測壓水頭大小分析

靜止水體內(nèi)各點(diǎn)的測壓水頭相等。A，B，C，D，E，F(xiàn)各點(diǎn)位于靜止左河內(nèi)，其測壓水頭相等，等于左河的水位h1。D′，E′，F(xiàn)′處于靜止河間潛水內(nèi)，測壓水頭相等，A，D′點(diǎn)處于水面或潛水位上，A′，B′，C′在潛水位之上，這6點(diǎn)只有位置水頭沒有壓強(qiáng)水頭。A，A′處于同一高度，B′，C′，D′高程依次降低，所以h1=HA=HA′>HB′>HC′>HD′。比較12個點(diǎn)的測壓水頭的大小，結(jié)果為：

h1=HA=HB=HC=HD=HE=HF=HA′>HB′>HC′>HD′=HE′=HF′=h2

在隔水板PP′左側(cè)點(diǎn)的測壓水頭相等，高于右側(cè)除A′點(diǎn)外的其他點(diǎn)，二者間存在水頭差，可形成大致向右下方的水力坡度。

2.2 初始時刻水力坡度的大小和方向

分析水力坡度的大小和方向可判斷地下水的運(yùn)動趨勢，地下水總是選擇最大水力坡度的方向運(yùn)動。當(dāng)選擇左側(cè)界面AD上的一點(diǎn)，同一方向上選擇不同的滲流路徑長度，其水力坡度大小不等，長度越小，水力坡度越大；當(dāng)終點(diǎn)M在右側(cè)液面運(yùn)動時，越靠近D′點(diǎn)水力坡度越大。比較AD上不同點(diǎn)的水力坡度大小時，如果以選定點(diǎn)至D′點(diǎn)作為滲流路徑長度，則越在下面的點(diǎn)其水力坡度越大；如果選擇任意的滲流路徑長度，則不同點(diǎn)的水力坡度大小不定(圖2)。

圖2 初始時刻左河側(cè)與河間地塊界面上點(diǎn)的水力坡度分析Fig.2 The initial hydraulic gradient analysis at point on the interface between the left stream and the inter-stream arear被放大。右上部分：微小時間段Δt后，鉛直段點(diǎn)的滲流發(fā)生沖突

以r為定長滲流路徑長度，對所有點(diǎn)以r為半徑在右側(cè)畫半圓或1/4圓(圖2)，規(guī)定水平向右方向?yàn)榱愣?，順時針旋轉(zhuǎn)為正方向。這種方法適合尋找最大水力坡度方向，也可比較不同點(diǎn)的水力坡度。在迅速抽出隔水板PP′的瞬間(初始時刻)，依據(jù)河間地塊潛水位h2，將左河側(cè)與河間地塊的界面AF分為兩段：高于h2的AD段，低于等于h2的DF段。界面AF上的A，B，C，D，E，F(xiàn)點(diǎn)具有代表性，只需分析這6個點(diǎn)的水力坡度，就可了解界面AF上所有點(diǎn)的情況。

(1)AD段(高于h2)，即鉛直段。

對于A點(diǎn)，當(dāng)θ→90°時，JA→1，達(dá)到極大值，方向趨于鉛直向下，所以潛水面上點(diǎn)的水力坡度最大只能為1。D之上附近點(diǎn)的滲流路徑長度小于半徑r，水力坡度更大，方向趨于鉛直向下。

(2)DF段(低于等于h2)，即水平段。

D點(diǎn)及其以下附近點(diǎn)，其右側(cè)半圓周與河間潛水位相交，高于潛水位點(diǎn)的水頭大，不取；F點(diǎn)及其附近點(diǎn)，其右側(cè)半圓周與隔水底板相交；因此，D，F(xiàn)及其附近點(diǎn)僅取在河間潛水內(nèi)的半圓周部分。

從數(shù)學(xué)角度看，滲流路徑r可取無窮小，除A點(diǎn)外的點(diǎn)的水力坡度為無窮大。實(shí)際上，達(dá)西或非線性滲透定律中的水力坡度都是從宏觀滲流的角度獲得的，r取無窮小不在滲流研究范圍。

2.3 微小時間段Δt后水力坡度的大小、方向及滲流速度

(2)鉛直段AD。初始時刻，鉛直段AD上點(diǎn)的最大水力坡度方向趨于鉛直向下。經(jīng)過微小時間段Δt后，任選上面B點(diǎn)滲流沖撞下面C點(diǎn)處的滲流，而受下面點(diǎn)的影響，B點(diǎn)最大水力坡度調(diào)整，方向向右偏轉(zhuǎn)和變小。

如果考慮Δt時間長一些，D點(diǎn)之上附近點(diǎn)的滲流接觸右側(cè)河間潛水面后，轉(zhuǎn)折為水平向右運(yùn)動，可以認(rèn)為，D點(diǎn)之上附近點(diǎn)有右傾直線接近平行右側(cè)液面的水力坡度，這樣的水力坡度對D點(diǎn)之上的點(diǎn)均有影響。

經(jīng)過微小時間段Δt后，所有點(diǎn)的水力坡度都減小，鉛直段的水力坡度方向向右偏轉(zhuǎn)。

3 對水躍的正面解釋

3.1 現(xiàn)有教材對水躍的解釋所存在的問題

現(xiàn)有大多數(shù)地下水動力學(xué)教材中對水躍的解釋，沿用1955年前蘇聯(lián)教材中的以反證法解釋潛水井水躍(滲出面)的形成[2-3,5,9]。我們認(rèn)為該解釋因所依據(jù)的流網(wǎng)不準(zhǔn)確，未能從根本上解釋水躍形成的原因。

井內(nèi)水面下是大空間中的水，井壁應(yīng)看作是一個等水頭面[10]。井壁等水頭面與上游的等水頭面ab相交于a，就說明流網(wǎng)是錯誤的，因?yàn)榱骶W(wǎng)中的流面或等勢面不能相交(圖3)。井壁也可以不是等水頭面，而與不同值的等水頭面相交，ab等水頭面與井壁間的區(qū)域(圖3中陰影部分)不是等勢體。這一點(diǎn)可從半無限厚承壓含水層與隔水頂板相連的井壁進(jìn)水不完整井獲得支持，Muskat得出結(jié)論：等水頭面是對稱于z軸的半旋轉(zhuǎn)橢球面，井的過濾器不是等水頭面[5]。這種情況下，地下水可通過ab等水頭面流向井，證明有側(cè)向補(bǔ)給的潛水井流下部的水力坡度大于上部，而不是只有a點(diǎn)過水。

圖3 潛水井水躍(滲出面)示意圖Fig.3 Sketch of the presence of a seepage face in a pumped well

退一步說，如果流量不通過ab等水頭面向井，曲面ab轉(zhuǎn)為流線。產(chǎn)生水躍后，如果井內(nèi)水位hW還是有類似ab等水頭面圍成的等勢體，而不能進(jìn)水，潛水只能從上面的滲出面入井，這也與實(shí)際情況不符(圖3)。

3.2 潛水井流的特征和流網(wǎng)

有滲出面的潛水井流與潛水回水，經(jīng)微小時間段Δt后水流有相似之處。而存在的差別有2個方面：①潛水井井壁是柱面，遠(yuǎn)小于前述河渠與河間地塊界面的面積；②抽水初期的非穩(wěn)定流動階段，井壁外側(cè)界面和井內(nèi)水面都不是定水頭；到穩(wěn)定流動階段，井壁外側(cè)界面和井內(nèi)水面都是定水頭。潛水回水中左河面是定水頭，河間地塊不是定水頭。可將用水力坡度方法對潛水回水初始階段定性-半定量分析的過程和結(jié)論用于分析潛水井流，以了解潛水井流的過程和特征，獲得更為準(zhǔn)確的流網(wǎng)，解釋產(chǎn)生水躍(滲出面)的原因。

井壁外側(cè)界面以井內(nèi)水面為準(zhǔn)，可分成2段：水面之下的水平段和水面之上的鉛直段(對應(yīng)滲出面)(圖4)。水平段點(diǎn)的水力坡度近相等最大且水平向井，即等水頭面是以井軸為心的圓柱面，相應(yīng)有最大的滲流速度，相當(dāng)于承壓含水層。鉛直段對應(yīng)滲出面，水力坡度小于水平段，越向上水力坡度越小；方向略偏離鉛直向下，即等水頭面略偏離垂直井壁，對應(yīng)有越小的滲流速度。這個結(jié)論與近來用數(shù)值和實(shí)驗(yàn)方法模擬所獲的結(jié)果較吻合。

圖4 井壁外側(cè)界面滲流速度方向及分段Fig.4 Segment division and flow velocities of the soil-well interface

圖5 潛水面保持恒定高程，有滲出面抽水完整井的流網(wǎng)(據(jù)文獻(xiàn)[12]修改)Fig.5 The equipotentials and flowlines of pumped well with a seepage face in an unconfined aquifer for which the overlying water table is assumed to remain at a coant elevation

5篇文獻(xiàn)[11-15]模擬不同情況下(完整井、不完整井、定高程潛水面、非零常流量地面、考慮不飽和水運(yùn)動等)有滲出面的潛水完整井穩(wěn)定流，所獲得的流網(wǎng)和滲流速度等特征為：對應(yīng)井內(nèi)水柱部分，含水層中等勢線幾乎平行于井壁，流線水平進(jìn)入井內(nèi)。但沿著滲出面，流線逐漸由傾斜趨于水平線(圖5)。在滲出面底部滲流速度達(dá)到最大值，井水面下的滲流速度接近相等，大于滲出面上的滲流速度(圖6)。所以進(jìn)井流量按長度比例來說主要來自井水面下井壁。

圖6 有滲出面的潛水完整井穩(wěn)定流數(shù)學(xué)模型及進(jìn)井流量(據(jù)文獻(xiàn)[12]修改)Fig.6 Math model flow rate of stable flow into pumped well with a seepage face in an unconfined aquifer

3.3 對水躍(滲出面) 的正面解釋

式中，d為管子直徑；v為斷面平均流速；R為斷面的水力半徑；C為Chezy系數(shù)[16]。

舉例，多孔介質(zhì)粒徑中砂級(0.5 mm)，井直徑300 mm[17]，假設(shè)井壁內(nèi)垂直井軸每層有12個圓管，管直徑為77.6 mm，圓管緊密堆積后，這些圓管總橫截面積略小于外井壁面積，井壁內(nèi)側(cè)的水力坡度略大于井壁外側(cè)。井內(nèi)大空間的和多孔介質(zhì)的d或R差距達(dá)150倍以上。不論采用上述哪個公式，水流速度v的差距都至少在12倍以上。潛水速度還要轉(zhuǎn)化為滲流速度，按孔隙度n=0.35，井中水流速度是潛水滲流速度的35倍以上，即在相同水力坡度和過水?dāng)嗝婷娣e下，井內(nèi)水流速度遠(yuǎn)遠(yuǎn)大于潛水滲流速度，所以抽水早期，抽水來自井儲水，由于井內(nèi)水位急速下降，井壁形成滲出面[18]。即使對于承壓井，井外水頭高于井內(nèi)水位，存在井損。潛水井形成水躍大于井損。

若抽水量Q非常小，則無水躍或水躍不明顯。若抽水量Q大到一定程度，就會產(chǎn)生明顯的水躍。保持井內(nèi)外大水頭差，才能保證潛水既有大的平均水力坡度，又有大的過水?dāng)嗝婷娣e，從而產(chǎn)生大進(jìn)井流量。由此可證明水躍越大，單寬流量越大。

設(shè)井外水頭位于A點(diǎn)，井內(nèi)水位等于D點(diǎn)，水躍AD=l，井外水頭AF=h0。以井軸為z軸，計算井壁AF上進(jìn)井單寬流量q(圖4)。單寬流量q等于水平段DF的流量q1與鉛直段AD流量q2之和，應(yīng)用達(dá)西定律有：

當(dāng)水躍l為變量時，l∈(0,h0)，單寬流量q是l的增函數(shù)，即水躍l越大進(jìn)入井內(nèi)的單寬流量q越大，流量就越大。井壁附近可能不適用達(dá)西定律，而適用非線性滲透定律。另外，把AD段的水力坡度考慮成h/r+1，也與實(shí)際情況略有差距，但這樣處理能簡化問題，結(jié)論也仍可作為定性分析的依據(jù)。

4 結(jié)論

(1)通過潛水回水初始時刻水力坡度分析，將左河側(cè)界面分為2段：河間潛水位之上的鉛直段和之下的水平段。水平段的點(diǎn)的水力坡度最大且相等，方向?yàn)樗较蛴?；鉛直段的點(diǎn)的水力坡度方向趨于鉛直向下，越向上水力坡度越小。

(2)相比于初始時刻，潛水回水經(jīng)過微小時間段Δt后，水力坡度都相應(yīng)變小。鉛直段相鄰點(diǎn)間地下水運(yùn)動發(fā)生沖突，水力坡度方向由趨于鉛直向下向右偏轉(zhuǎn)。可將潛水回水初始階段定性-半定量分析的過程和結(jié)論用于分析潛水井流，有助于更好地了解潛水井流的過程和特征。

(3)地下水動力學(xué)教材中對水躍的解釋多未從根本上解釋水躍形成的原因。形成水躍的根本原因是：在相同水力坡度和過水?dāng)嗝婷娣e下，井內(nèi)大空間水流速度遠(yuǎn)大于孔隙中潛水流速度。井內(nèi)水位急速下降，從而形成水躍。水躍越大進(jìn)入井內(nèi)的流量就越大。

致謝：在寫作過程中得到2008級姬祥、閆振、張光亮同學(xué)和2009級許江、張忠杰等同學(xué)給予的啟發(fā)和幫助，在此一并表示感謝！

[1] 郭東屏，宋焱勛，錢會，等. 地下水動力學(xué)[M]. 西安：陜西科學(xué)技術(shù)出版社，1994：112-115，132-133.

[2] 何俊杰，王明偉，王廷國. 地下水動力學(xué)[M]. 北京：地質(zhì)出版社，2009：34-39，56.

[3] 蔣輝，曾波，潘宏雨. 地下水動力學(xué)[M]. 北京：地質(zhì)出版社，2009：50-54，66-67.

[4] 吳吉春，薛禹群. 地下水動力學(xué)[M]. 北京：中國水利水電出版社，2009：53-57.

[5] 薛禹群，吳吉春. 地下水動力學(xué)[M]. 北京：地質(zhì)出版社，2010：62-66，79-80，167-168.

[6] Guo W X. Transient groundwater flow between reservoirs and water-table aquifers[J]. Journal of Hydrology, 1997, 195: 370- 384.

[7] Serrano S E，Workman S R，Srivastava K，et al．Models of nonlinear stream aquifer transients[J]. Journal of Hydrology, 2007, 336: 199-205.

[8] Moutsopoulos K N. The analytical solution of the Boussinesq equation for flow induced by a step change of the water table elevation revisited[J]. Transp Porous Med, 2010, 85: 919-940.

[9] 陳崇希，林敏. 地下水動力學(xué)[M]. 武漢：中國地質(zhì)大學(xué)出版社，1999：61-62.

[10] 王大純，張人權(quán)，史毅紅，等．水文地質(zhì)學(xué)基礎(chǔ)[M]．北京：地質(zhì)出版社，1986：34-35.

[11] Simpson M J, Clement T P, Gallop T A. Laboratory and numerical investigation of flow and transport near a seepage-face boundary[J]. Ground Water, 2003, 41(5): 690-700．

[12] Rushton K R. Significance of a seepage face on flows to wells in unconfined aquifers[J]. Quarterly Journal of Engineering Geology and Hydrogeology, 2006, 39: 323-331.

[13] Chenaf D, Chapuis R P. Seepage face height, water table position, and well efficiency at steady state[J]. Ground Water, 2007, 45(2): 168-177.

[14] Yakirevich A, Gish T J,l. Potential impact of a seepage face on solute transport to a pumping well[OL]. Vadose Zone Journal, 2010 (9): 686-696[2010-08-03]. http://www.VadoseZoneJournal.org

[15] Behrooz-Koohenjani S，Samani N，Kompani-Zare M．Steady flow rate to a partially penetrating well with seepage face in an unconfined aquifer[J]. Hydrogeology Journal, 2011, 19: 811-821.

[16] 薛禹群，朱學(xué)愚，吳吉春，等．地下水動力學(xué)[M]．北京：地質(zhì)出版社，1997：213-219.

[17] Fetter C W. Applied hydrogeology[M]．Englewood Cliffs, New Jersey: Prentice Hall, 2001: 123.

[18] Rushton K R. Groundwater hydrology[M]．England: John Wiley & Sons Ltd, 2003: 151.

The qualitative-semiquantitative analysis of the unconfined aquifer’s recharging from stream at the initial stage and explanation for presence of a seepage face in a pumped well

SUN Baoliang1, LIANG Junhong1, CUI Weili1, WU Yisong2, MA Guangzhi3

(1．CollegeofResourcesandCivilEngineering，NortheastChinaUniversity，Shenyang110004，China;2．CentralSouthResearchInstituteofMetallurgicalGeology,Yichang443003,Hubei，China;3．NorthernEngineeringandTechnologyCorporation，MetallurgicalCorporationofChinaLtd.，Dalian116600，Liaoning,China)

We use hydraulic gradient method to analyze the initial unconfined aquifer’s recharging from streams qualitatively-semiquantitatively thus provide a new visual angle for recharging process. Analysis of the hydraulic gradient of initial recharging, the left stream interface is divided into two segments, the vertical segment and the horizontal segment above and under the water table in the inter-stream area. The hydraulic gradient of the points in the segment is maximum and equal. It is to right in the horizontal direction. The hydraulic gradient of the points in the vertical segment tends downward in the vertical direction and it gets smaller upward. After a tiny period, the hydraulic gradients of any point become smaller, waters coming from adjacent points in the vertical segment collide, and the direction of the hydraulic gradient is deflected to the right. The qualitative-semiquantitative analytical method can also be employed to analyze unconfined radial flow to a pumped well. Presence of a seepage face can be explained in the way: under same hydraulic gradient and cross-sectional area, the flow velocity in big space of the well is much greater than in the pore. The water level in the well rapidly declines thus a seepage face is developed. The higher the seepage face is, the greater the rate of discharge into the well is.

unconfined aquifer; recharge; hydraulic gradient; qualitative analysis; well; seepage-face

2015-06-19； 改回日期： 2016-06-21； 責(zé)任編輯： 趙慶

孫寶亮(1970—)，男，講師，從事地下水、沉積學(xué)方面研究。通信地址：沈陽市東北大學(xué)265信箱；郵政編碼：110004；E-mail: sunbaoliang@mail.neu.edu.cn

10. 6053/j. issn.1001-1412. 2016. 03. 020

P641.2

A