梁 越，陳建生，陳 亮

隨著社會的發(fā)展，工農(nóng)業(yè)及居民生活淡水用量日益增大，地下水的開采量也隨之增加。在沿海地區(qū)，如果地下水抽水量超過了它的補給量，使海岸帶附近地下水水頭下降，當?shù)w水頭低于附近海水楔形體的水頭時，咸淡水界面就要向內(nèi)陸推進，直至形成新的平衡，這樣就造成了海水的入侵［1］。目前，全世界已經(jīng)有幾十個國家和地區(qū)的幾百個地方發(fā)現(xiàn)了海水入侵問題，并給各國沿海地區(qū)帶來了嚴重危害，造成了巨大的經(jīng)濟損失，嚴重阻礙經(jīng)濟社會的持續(xù)發(fā)展。全球范圍的海水入侵問題引起了各國的關注，有關國家積極開展了海水入侵問題的研究與治理［2］。

造成海水入侵的根本原因是地下水超采等人為因素引起的地下水位降低，從而使得海水回倒灌入地下含水層，因此，升高地下水的水位形成水力帷幕可以阻止海水的入侵。升高地下水水位的方法有很多，許多研究者對此做了研究［3－5］；本文討論注水井水力帷幕這一阻止海水入侵形式，即通過注水井注水升高地下水水位，通過調(diào)節(jié)注水井的注水量及注水井的間距形成水力帷幕，達到阻止海水入侵的目的。通過計算分析得出了這一方法的可行性與有效性，并提出了這種帷幕形式的初步優(yōu)化設計方法。

1 注水井水力帷幕阻止海水入侵原理

1．1 研究海水入侵的常用模型

由于海水與淡水是可混溶的，因此兩者間的交界面實際上是一個過渡帶，但當過渡帶范圍較窄，或為了計算方便，過渡帶可以近似看作為突變界面，所以，研究海水入侵常用的概化模型可以分為突變界面模型與過渡帶模型兩種。突變界面模型是一種理想化的模型，由它所得的結果是一種近似解；而過渡帶模型是滲流過程與溶質運移過程的耦合，由一個偏微分方程組來描述，因此很難得到解析解，只能利用數(shù)值方法得到解答。早期，研究人員采用突變界面模型對海水入侵的機理與發(fā)展進行了大量的研究［6，7］。本文提出了注水井水力帷幕的簡單化模型，研究它阻止海水入侵的可行性、有效性及設計方法，并通過突變界面模型給出模型的近似解答。

1．2 基于突變界面模型的注水井水力帷幕阻止海水入侵原理

利用注水井向地層中注水會引起地下水水位的升高，每隔一定距離設置注水井，調(diào)節(jié)井間距及注水量，在注水井沿線可以形成一個完整的水力帷幕?；谕蛔兘缑婺Ｐ?，注水井水力帷幕阻止海水入侵的原理可以由圖1來描述。

圖1中，L為注水井間距；a為注水井到海岸等效邊界的距離；h為淡水面在海平面之上的高度；Z為突變界面在海平面之下的高度。根據(jù)突變界面模型，只要h升高到一定的程度，就可以阻止海水的繼續(xù)入侵，如果繼續(xù)使h升高，那么海水將會在淡水的驅使下向海洋后退，因此，海水的前鋒剛好位于注水井沿線是阻止海水入侵的臨界狀態(tài)。為了研究這一臨界狀態(tài)，考慮一種比較理想的狀態(tài)，假設一均質地層，不透水邊界為水平面，每個注水井到海岸等效邊界的距離a都是相等的，注水井的間距及單井的注水量也是相等的，那么注水井沿線的水位變化可以由圖2來描述。

圖1 水力帷幕阻止海水入侵原理圖Fig．1 Principle of the hydraulic curtain stopping seawater intrusion

圖2 等間距注水井水力帷幕形成前后的地下水水位變化Fig．2 Ground water level change infected by the equal space water injection hydraulic curtain

圖2 中H0為水力帷幕設置前的地下水水位；hf為兩注水井間距中點位置在水力帷幕設置后的水位升高。從圖2中可以看出，帷幕中注水井處的水位最高，兩注水井中點處水位最低，為保證帷幕阻止海水入侵的有效性，就需要兩注水井間距中點位置處具有阻止海水入侵的能力，并以此處阻止海水入侵的臨界狀態(tài)作為整個帷幕有效的臨界狀態(tài)。根據(jù)突變界面模型［6］，H0與hf之間的關系可以表示為

式中：γf為淡水密度；γs為海水密度。

如果取 γs＝1．025 g／cm3，γf＝1．000 g／cm3，那么H0，也就是說，只要保證兩注水井間中點處的水位升高大于H0，就能保證帷幕的有效性，達到阻止海水入侵的目的。

2 注水井水力帷幕的設計

對于上述的帷幕形式，只要保證兩注水井間中點處的水位升高大于，就能達到阻止海水入侵的目的，因此，設計一條阻止海水入侵的注水井水力帷幕，就是合理選擇帷幕的各參數(shù)，在保證帷幕有效的前提下，達到經(jīng)濟適用的目的。

2．1 注水井水力帷幕參數(shù)間的關系

對于某一特定的地層，在確定了水力帷幕設置的位置之后，確定帷幕形式的參數(shù)只有2個：單井注水量以及注水井間距。注水量越小越經(jīng)濟，而擴大注水井間距可以節(jié)省注水井個數(shù)以節(jié)約成本，注水井個數(shù)同樣影響著總的注水量的大小，2個參數(shù)之間是相關的，因此，首要的任務是理清單井注水量與注水井間距之間的關系。

注水井作用下的地下水滲流場可以用水頭勢函數(shù)來描述。對于單個注水井影響下的滲流場勢函數(shù)可以利用鏡像法及疊加原理得到［8］（圖3）：

圖3 單注水井作用下滲流場分布Fig．3 Seepage field under the single water injection well

滲流場中任一點P的水頭勢函數(shù)可以表示為

式中：Q為單井注水量；r1為點P到注水井的距離；r2為點P到虛擬抽水井的距離；C為常數(shù)，其值可以根據(jù)邊界條件確定。

對于由無窮多個注水井組成的水力帷幕影響下的滲流場，可以看作是單注水井作用下的疊加（圖4）。

滲流場中某點A（x，y）的勢函數(shù)可以表示為

圖4 鏡像法計算注水井水力帷幕滲流場勢函數(shù)示意圖Fig．4 Potential function of the seepage field based on the image method and superposition principle

式中：rf0為（0，a）處注水井到點A的距離，rf0＝）處虛擬抽水井到點A的距離為編號±n的注水井到點A的距離為編號±n的虛擬抽水井到點A的距離

將各參數(shù)代入式（3），得

根據(jù)復函數(shù)性質x2＋y2＝（x＋yi）（x－yi）可以將式（4）化簡為

利用式（5）減式（6），得到：

考查兩注水井間距中點處的最低水頭勢，即在x＝nL／2，y＝a處：

對于在潛水層中的注水井，

式中K為地層的滲透系數(shù)。

如前所述，要保證帷幕有效，只要保證在兩注水井中點處hf＝H0，代入式（10），得在兩注水井的水頭勢為

將式（11）代入式（9），得到保證帷幕有效性時的單井注水量：

式（12）是在無限長水力帷幕情況下得到的單井注水量，對于足夠長的水力帷幕，也可以用式（12）計算注水量。設水力帷幕的長度為D（D足夠長，D＞＞L），那么整個帷幕所需注水井的個數(shù)為

整個水力帷幕所需的注水量為

2．2 注水井水力帷幕的優(yōu)化設計

如果將本帷幕的成本簡單地歸結為帷幕所需要的注水量與帷幕建設維護成本兩部分，由式（13）及式（14）可知這兩部分都可以表示成注水井間距L的函數(shù)。帷幕每年總成本可以表示為

式中：MT為每年總成本；f為每年注水量成本系數(shù)；g為注水井及注水設備每年建設維護成本系數(shù)。

要確定注水井間距L的最優(yōu)值，只需對式（15）求極小值，得到MT取得極小值時對應的注水井間距。

3 注水井水力帷幕設計實例

下面以一簡單算例說明注水井水力帷幕的優(yōu)化設計方法。假設一條長50 km的海岸線遭受到海水入侵，利用注水井水力帷幕對海水入侵進行阻止。以每年的總成本確定帷幕的優(yōu)化設計方案。每年注水量成本系數(shù)f以及注水設備每年建設維護成本系數(shù)g可以根據(jù)材料價格及施工費用得到。本例中涉及的其參數(shù)見表1。

表1 計算參數(shù)表Table 1 Parameters used for the example

將表1中的數(shù)據(jù)代入式（15），得到每年總成本的表達式：

年成本與注水井間距的關系如圖5。

圖5 年成本與注水井間距關系曲線Fig．5 Relation between the total cost of each year and the space between wells

由圖5可以看出，年成本隨著注水井間距的增加先減小再增大，在L＝26．2 m處取得極小值，極小值為萬元。也就是說，在本例所涉及的注水井止水帷幕中，注水井間距為26．2 m時所需的年成本最少，為最優(yōu)選擇。將L的值代入式（12）可以得到單井的注水量 Q為0．16 m3／d；將 L值代入式（14）得到總注水量 QT為305．3 m3／d；注水井總數(shù)為 N＝50 000／26．2≈1 908個。

4 結 論

（1）升高地下水位是治理海水入侵的有效方法，通過調(diào)節(jié)注水井的間距及單井注水量，所形成的水力帷幕可以有效地阻止海水入侵；

（2）通過推導注水量與注水井間距的關系可以將帷幕運行的年成本表示為注水井間距的函數(shù)，并通過求解年成本的極小值可以確定最優(yōu)的注水井間距，實現(xiàn)注水井水力帷幕的優(yōu)化設計。

（3）需要指出的是，實際的地層情況比算例中的情況要復雜得多，注水井間距及其到計算邊界的距離等參數(shù)可能不再是常數(shù)，同時突變界面模型也影響了計算的精度，因此，對于這些復雜的情況，解析方法將不能滿足需求，需要利用數(shù)值方法得到解答。

［1］ 薛禹群，吳吉春，謝春紅，等．萊州灣沿岸海水入侵與咸水入侵研究［J］．科學通報，1997，11（22）：2360－2367．

［2］ 郭占榮，黃奕普．海水入侵問題研究綜述［J］．水文，2003，623（3）：10－15．

［3］ 劉青勇，董廣清，耿樹德，等．淡水帷幕防治海（咸）水入侵的物理模擬試驗研究［J］．水利學報，1998，2（2）：30－33．

［4］ 張 奇．海水入侵的實驗研究［J］．水文地質與工程地質，2005，4：43－47．

［5］ 高學平，楊建華，涂向陽，等．帷幕防治海水入侵的數(shù)值模擬研究［J］．海洋環(huán)境科學，2006，11，25（4）：55－58．

［6］ BEAR J．Hydraulics of Groundwater［M］．New York：McGraw－Hill，1979．

［7］ MOOR Y H，STOESELL R K，EASLEY D H．Fresh－wa－ter／See－water Relationship within a Groundwater Flow System，Northeastern Coast of the Yucatan Peninsula［J］．Ground Water，1992，30（3）：343－350．

［8］ 陳建生，李興文，趙維炳．堤防管涌產(chǎn)生集中滲漏通道機理與探測方法研究［J］．水利學報，2000，9（9）：48－54．