孟陳銳，汪　瑩，谷帥軍，劉紅杏，周彥兆

(中國地質大學(武漢)環(huán)境學院，湖北 武漢 430074)

?

無越流含水層中完整井的井流試驗確定隔水邊界位置

孟陳銳，汪瑩，谷帥軍，劉紅杏，周彥兆

(中國地質大學(武漢)環(huán)境學院，湖北 武漢 430074)

討論了在未知無越流含水層隔水邊界的方向和位置時，通過完整井的井流試驗，利用觀測井降深-時間的數(shù)據(jù)以及Cooper-Jacob簡化公式的直線圖解法來確定隔水邊界的問題。運用數(shù)學建模的方法，分析討論了當觀測井個數(shù)為兩個時的隔水邊界位置。結果表明：使用兩個觀測孔降深-時間的數(shù)據(jù)可找出隔水邊界的兩個可能解，但無法確定哪個是真實解。因此筆者在此基礎之上提出了確定隔水邊界的一種方法：使用三個不在一條直線上的觀測井降深-時間的數(shù)據(jù)，通過直線圖解法和作圖，利用三個圓的交點近似確定隔水邊界的位置。

井流試驗；隔水邊界位置；反映法；虛井位置

在承壓含水層中進行抽水試驗時，由于巖性、地質構造等原因，經常會遇到含水層有隔水邊界的情況。確定隔水邊界的方向與位置對水資源的開發(fā)與礦坑疏干的預測有著重要意義。對于此類問題，運用反映法進行求解較為簡便。所謂反映法是對某一邊界或某些邊界來映射實際的抽水井，使得在邊界的另一邊有一個對稱的虛構的抽水井或注水井，其作用與實井等效[1]。由此邊界問題轉化為無界問題，繼而可用無限含水層中的Theis井流模型求解含水層參數(shù)。徐建國等[2]通過水文地質邊界附近非穩(wěn)定流抽水試驗資料的分析，按邊界透水性將水文地質邊界劃分為5類,并結合實例說明了邊界透水性的判別方法。這對于確定地下水系統(tǒng)邊界條件具有較為重要的意義。郭建青[3,4]等對比分析了標準曲線對比法、基于Jacob近似的直線圖解法、梯度搜索法，三種直線隔水邊界含水層抽水試驗數(shù)據(jù)的分析方法，詳細地描述了不同方法的適用性。余克林[5]等詳細描述了特定條件直線圖解法在確定隔水邊界中的應用，并結合實例闡述了復雜邊界條件下的水文地質學求參方法。陳崇希[6]等在《地下水動力學》一書中用直線圖解法近似確定兩觀察孔到虛井的距離，利用兩圓相交的方法確定直線邊界的方向和位置，但此方法的解是不唯一的。

本文指出了《地下水動力學》一書中確定隔水邊界方法的不足，在利用觀測井降深-時間的數(shù)據(jù)以及Cooper-Jacob簡化公式的直線圖解法的基礎上，運用數(shù)學建模的方法分析討論了無越流含水層中確定隔水邊界位置的方法。

1　數(shù)學模型建立與計算

1.1承壓含水層的直線圖解法

直線圖解法是將抽水試驗的實測數(shù)據(jù)(觀測井降深s，時間t)投在單對數(shù)坐標紙上并做成s-lgt曲線，此曲線在一定的區(qū)間上將呈現(xiàn)為直線，因而可以根據(jù)直線的斜率和截距來確定含水層的兩個參數(shù)(導水系數(shù)T、壓力傳導系數(shù)a)。

在s-lgt曲線中，分別延長第一、第二直線段，兩直線交于點c，此點坐標(tc，sc)，第一直線段與lgt軸的交點為t0，據(jù)此可求出虛井至觀測井的距離ρ：

(1)

式中：ρ為虛井至觀測井的距離(m)；r為觀測井至抽水井的距離(m)；tc為第一、第二直線段交點的橫坐標所對應的時間(min)；t0為第一直線段與lgt軸交點的橫坐標所對應的時間(min)。

1.2提出的疑問

按照《地下水動力學》[6]中的描述,用(1)式分別求出觀測孔1和觀測孔2至虛井的距離ρ1和ρ2，然后分別以觀測孔1、觀測孔2為圓心，ρ1和ρ2為半徑作圓，兩圓的交點即為虛井的位置。實井和虛井連線的中垂線就是隔水邊界的位置。但是兩圓相交有兩個交點，而真實的虛井只有一個，能否確定這兩個虛井中哪一個是真實的虛井呢？

此處討論的隔水邊界為已經確定存在的，并且只有一條的情況。下面我們將討論在觀測孔為兩個時，能否確定這兩個虛井中哪一個是真實的虛井。

1.3在一般情況下討論

將抽水井用Pw表示，觀測井用k1、k2表示，虛井用P1、P2表示。再以觀測井k1、k2為圓心半徑為ρ1、ρ2作圓，連接他們的圓心k1k2和交點P1P2(如圖1所示)，則k1k2是P1P2的中垂線。

圖1　線段k1k2、P1P2垂直示意圖

以k1k2為X軸，P1P2為Y軸建立坐標系(如圖2所示)。設P1(0，a)、P2(0，-a)為可在Y軸上移動的點(P1、P2關于X軸對稱)，k1(b，0)、k2(c，0)為可在X軸上移動的點，抽水井Pw的坐標為(m，n)。

圖2　抽水井Pw取值范圍示意圖

Pw應滿足：(m-b)2+n2

(2)

連接PwP1、PwP2，并作他們的中垂線(如圖3所示)，中垂線即為隔水邊界。因為抽水井、觀測井應在隔水邊界的一側，而虛井應在隔水邊界的另一側，所以隔水邊界不能與線段k1k2相交。

若PwP1、PwP2的中垂線與線段k1k2無交點，則圓的兩個交點都為虛井，隔水邊界有兩種可能：若PwP1、PwP2的中垂線與線段k1k2有1個交點，則圓的兩個交點中只有1個交點為虛井，隔水邊界只有1種可能；若PwP1、PwP2的中垂線與線段k1k2有兩個交點，則沒有隔水邊界。

圖3　k1、k2位于X軸異側時隔水邊界位置討論圖

經計算，我們發(fā)現(xiàn)PwP1、PwP2這兩條中垂線在m0時交于X軸上同一點，并且這一點不在線段k1k2上，在m=0平行X軸，與線段k1k2也沒有交點。證明如下：

2016年，山東省有6個果蔬區(qū)域公用品牌成功入選為第一批著名農產品區(qū)域公用品牌名單，膠州大白菜作為其中一個。2017年，11個區(qū)域公用品牌，如膠州大白菜被選為最受歡迎的農產品區(qū)域公用品牌。膠州大白菜協(xié)會繼續(xù)研發(fā)白菜品種、栽培技術和白菜深加工產品的創(chuàng)新和開發(fā)形成較強的市場競爭力。

直線PwP1、PwP2中垂線的方程分別為：

(3)

(4)

在上述方程中分別令y=0，求得其兩個解x1、x2為：

(5)

又由(2)式可得：

m2+n2-a2<2mb且a2+b2-a2<2mc

(6)

當m=0時，這兩條中垂線平行于線段k1k2。

補充：1.因為以上證明沒有對k1(b，0)、k2(c，0)的值進行限定，他們是在X軸上任意移動的點，所以當k1k2位于X軸的同側時證明依然成立，m>0時(如圖4)，m<0時(如圖5)，PwP1、PwP2的中垂線與線段k1k2沒有交點。

圖4　k1、k2位于X軸同側且m>0時隔水邊界位置示意圖

圖5　k1、k2位于X軸同側且m<0時隔水邊界位置示意圖

由此可推出，這兩條隔水邊界都可使抽水井、觀測井在隔水邊界的一側，虛井在隔水邊界的另一側，兩條隔水邊界都有可能是真實解，需要其他方法來作進一步確定。

1.4特殊情況

該問題存在一種特殊情況，當k1、k2和虛井剛好在一條直線上并且虛井不能在線段k1k2上(否則不滿足抽水井、觀測井在隔水邊界的一側，而虛井在隔水邊界的另一側)。所以虛井不在線段k1k2上，但滿足三點k1、k2和虛井在一條直線上，這時作圖會得到兩個內切的圓(如圖6所示)。

圖6　特殊情況下隔水邊界位置示意圖

此時兩圓的切點即為虛井的位置。但只有當觀測井k1、觀測井k2和虛井在一條直線上時才會出現(xiàn)此情況，也等同于觀測井k1、k2兩口井的連線剛好垂直于隔水邊界。這種情況在我們運用此方法求解隔水邊界位置時出現(xiàn)的概率很小。

1.5結論驗證

為了進一步證明我們的結論，我們用動畫進行了演示。我們在不改變抽水井Pw、觀測井k1的位置的情況下，僅改變觀測井k2位置，觀察兩條隔水邊界與線段k1k2的交點。在動畫前期，兩條隔水邊界與線段k1k2有兩個相同的交點，所以這兩條隔水邊界都不符合，但是實際問題中，應存在一個隔水邊界，所以這種情況不符合實際。當觀測井k2運動到臨界位置(k2剛好位于兩條隔水邊界的交點時)，觀測井位于隔水邊界上也不符合實際。當觀測井k2繼續(xù)運動時，兩條隔水邊界與線段k1k2沒有交點，說明兩條隔水邊界都有可能是真的，符合我們之前的結論。

2　隔水邊界的確定

通過之前的研究我們發(fā)現(xiàn)圓的兩個交點都符合條件，無法判斷哪個是真實的虛井，也無法判斷哪個是真實的隔水邊界，我們猜想能否通過增加一個觀測井的方式來確定隔水邊界的真實位置。

通過直線圖解法可得出觀測井k3到虛井的距離ρ3，以k3為圓心ρ3為半徑作圓。

圖7　k1、k2、k3在一條直線上時隔水邊界位置示意圖

圖8　k1、k2、k3不在一條直線上時隔水邊界位置示意圖

若三個觀測井k1、k2、k3在同一條直線上，以k3為圓心ρ3為半徑作圓，則此圓會交于虛井p1、p2，同樣無法從兩個交點中判斷出哪個是真實的虛井(如圖7)；若三個觀測井不在同一直線上時，因為觀測井k3不在虛井p1p2的中垂線上，所以觀測井k3到兩虛井的距離不相等，若以k3為圓心ρ3為半徑作圓，則此圓必與其中的一個虛井相交(如圖8)，這三個圓的交點就是虛井的位置，實井與虛井連線的中垂線就是隔水邊界的位置。第三個觀測井的取值應在兩隔水邊界的內部。

以此類推，當使用三個以上觀測井的降深數(shù)據(jù)時，只要這些觀測井不在同一條直線上就可確定隔水邊界的位置。

3　結語

在用完整井的井流試驗方法確定無越流含水層中的隔水邊界時，若未知隔水邊界的方向和位置，使用兩個觀測孔的降深觀測數(shù)據(jù)可找出隔水邊界的兩個可能解，但無法確定哪個是正解(只有當兩觀測井和虛井剛好在同一條直線上時才可以直接確定虛井的位置)；使用三個及以上觀測井的降深數(shù)據(jù)時，只要這些觀測井不在同一條直線上就可確定隔水邊界的位置。

[1]趙運德,安進強.幾種不同邊界附近水井的滲流計算模式[J].地下水.2001,4:198-201.

[2]徐建國.利用水文地質邊界附近抽水試驗資料判定邊界透水性[J].勘察科學技術.1998,6:12-14.

[3]郭建青,周宏飛,李彥,等.隨機搜索算法在確定含水層參數(shù)中的應用[J].中國農村水利水電.2010,12:48-51.

[4]郭建青,鄭力,雷玉平,等.分析直線隔水邊界附近抽水試驗數(shù)據(jù)的新方法[J].工程勘察.2006,1:26-29.

[5]余克林.特定條件直線圖解法在水文地質參數(shù)確定中的應用[J].現(xiàn)代礦業(yè).2013,2:58-59.

[6]陳崇希，林敏.地下水動力學[M].武漢：中國地質大學出版社.1999:135-137.

[7]葛曉云.水文地質邊界性質s～lgt直線斜率關系的探討[J].地下水.1992,4:194-196.

[8]張力春,肖長來,梁秀娟.穩(wěn)定井流直線圖解法確定含水層參數(shù)[J].世界地質.2006,1:67-70.

[9]劉英龍.關于水文地質虛擬邊界的討論及應用[J].今日科苑.2009,2:31.

[10]葛孝椿,陳光臨,張立爭,等.直線補給邊界深井降水滲流分析的鏡像反射方法[J].巖土工程學報.2005,7:812-818.

[11]李同貴,王德文,劉鐵平,等.正交直線供水邊界完整(井)干擾井穩(wěn)定流計算公式推導及應用[J].黑龍江水專學報.2009,2:60-63.

[12]楊科元,趙建宇.映射法在隔水邊界附近對地下水水位預報的應用——以馬鬃山西段營毛沱水源地為例[J].甘肅地質.2010,3:60-65.

Determination of impervious boundary with a pumping test in a confined aquifer

MENG chen-rui，Wangying，GU shuai-jun，LIU hong-xing，ZHOU yan-zhao

(School of Environmental Studies,China University of Geosciences,Wuhan,Hubei,430074,China)

With the unknown of the direction and position of impervious boundary,the developed method can be used to determine the position of impervious boundary through the pumping test data associated with the Cooper Jacob method.Using the method of mathematical modeling,the position of the impervious boundary was analyzed and discussed for the condition of with two observation wells.The results indicated that the drawdown-time data of two observation wells can be used to find out the two possible solutions of the impervious boundary`s position,but it is not sure which is the real solution.Based on this,we proposed a method to determine the position of impervious boundary:using drawdown-time data in three observation wells that do not have a straight line and through the Cooper Jacob method and matching,using the intersection of the three circles to determine the position of impervious boundary.

well flow experiment；impervious boundary′s position；reflection method；imagined well

2016-04-26

孟陳銳(1995-)，男，湖北武漢人，主攻方向：環(huán)境工程。

P641.2

B

1004-1184(2016)05-0019-04