新疆孔雀河沿岸地下水位變化的影響半徑模擬分析

李玉朋 吾買爾江·吾布力 何 宇 李衛紅

(1.中國科學院新疆生態與地理研究所荒漠與綠洲生態國家重點實驗室，新疆 烏魯木齊 830011； 2.新疆塔里木河流域管理局，新疆 庫爾勒 841000)

新疆孔雀河沿岸地下水位變化的影響半徑模擬分析

李玉朋1吾買爾江·吾布力2何 宇2李衛紅1

(1.中國科學院新疆生態與地理研究所荒漠與綠洲生態國家重點實驗室，新疆 烏魯木齊 830011； 2.新疆塔里木河流域管理局，新疆 庫爾勒 841000)

本文采用裘布依公式，以孔雀河中上游為實例，模擬了孔雀河第三分水樞紐至普惠段、普惠水庫至經濟牧場段、尉犁縣的經濟牧場至66分水閘段在不同地下水位、不同抽水降深情況下的影響半徑。結果表明：普惠水庫至經濟牧場段的影響半徑在1米、3米、5米的降深影響都是最大，且是其他兩地影響半徑的2倍甚至更大。第三分水樞紐至普惠段的單口井的抽水量較大，在棉花生長季，單口井抽水對周邊地下水影響的范圍達776米。目前，孔雀河兩岸水位嚴重下降，胡楊大片死亡，應控制耕地擴張，加快實施退耕還水計劃，堅持以水定地、以水定結構、以水定發展的原則，取締非法機電井，才能使水位逐漸恢復，利于經濟與社會的發展。

地下水埋深；水井；影響半徑；孔雀河

1 引 言

“影響半徑”R是水文地質中一個重要參數，它受到抽水量，土壤質地，釋水系數等因素的影響。影響半徑的正確估算對計算滲透系數K、開采井群間距的設計、地下水資源的評價、防護帶的圈定具有重要的意義[1-3]。確定水井影響半徑的方法很多，經驗公式是最常采用的方法有庫薩金公式、吉爾特公式、舒爾采公式。此外，還有利用觀測孔直接觀測，或根據單位涌水量、單位水位降深、含水層巖石特性等方法也可以求得。其中，抽水試驗理論方法比較成熟，對非承壓水計算國外運用最多的為Jacob修正模型[4]、Neuman模型[5]和Boulton模型[6]。其中，裘布依公式可以直接利用單孔抽水試驗的孔內涌水量來求影響半徑，由于此方法不用額外打觀察孔而廣受歡迎，應用廣泛。

2 穩定流抽水影響半徑R的計算

對任一均質含水層，滲透系數可以看作常量，它不隨流量、時間、降深而變化。但R是抽水量、降深、土壤系數、釋水系數含水層厚度以及各種補給量的函數，這些參數控制R的變化。1863年裘布衣在達西定律的基礎上，建立了穩定流單井抽水公式，即著名的裘布衣公式。奠定了穩定流的理論基礎，裘布衣公式在生產實踐中使用已有100多年的歷史，并且還將為人們繼續使用。裘布衣公式是以一定假設條件為基礎的，把抽水現場完全理想化。裘布依推導公式[7]時的假定條件是：①含水層是均質、各向同性、等厚、水平的；②地下水為層流，符合達西定律，地下水運動處于穩定狀態；③靜水位是水平的，抽水井具有圓柱形定水頭補給邊界；④對于承壓水，頂底板是完全隔水的，對于潛水，井邊水力坡度不大于1/4，底板完全隔水。對承壓井(1)與潛水井(2)計算如下：

(1)

式中，Q為抽水井流量，m3/d；K為含水層的滲透系數，m/d；M為承壓含水層厚度，m；S為抽水井中水位降深，m；R為影響半徑，m；r為抽水井的半徑，m。

(2)

式中，Q為抽水井流量，m3/d；K為含水層的滲透系數，m/d；H為潛水含水層厚度，m；S為抽水井中水位降深，m；R為影響半徑，m；r為抽水井的半徑，m。

3 案例分析

孔雀河流域地處塔里木盆地東部，近10余年來，隨著耕地面積的不斷擴大和農業灌溉用水的增加，流域地下水超采現象嚴重，地下水位大幅下降，導致綠洲外圍依靠地下水維系生存的天然植被大面積衰敗和死亡。據統計調查，孔雀河流域上、中游耕地面積擴張迅速，1990-2000年期間，耕地面積平均每年增加6.92平方公里，2000-2014年間，年平均增加34.2平方公里(新疆統計局)。地下水的嚴重超采，導致地下水位大幅下降，其中，第三分水樞紐至普惠段的地下水埋深達到-10.59m，普惠水庫至經濟牧場為-26.67m，尉犁縣的經濟牧場至66分水閘段達到-14.38m。以抽水井為中心，四周地下水在水力梯度的作用下向中心匯聚，從而形成地下漏斗，機理如圖1所示。目前，孔雀河第三樞紐至66分水閘段的地下水水位已經下降到胡楊生存的臨界地下水位[8-9]，孔雀河沿岸大片胡楊林死亡。

為了查明單口井抽水對周邊地下水位變化的影響強度，我們結合孔雀河流域打井鉆孔資料和水文地質條件，采用穩定流抽水影響半徑R的計算方法，對孔雀河流域地下水位變化較大的普惠-經濟牧場地段抽水引起的地下水位變化及其影響半徑進行了計算分析，旨在為合理布局地下水井、確定地下水開采閾值提供依據。

圖1 抽水引起的周邊地下水水位下降對植被生長的影響。上面的虛線表示地下水能夠滿足植被的生長需求，下面的虛線表示地下水位下降到對植被生長產生水分脅迫

3.1 流域地質條件

土壤質地是影響土壤半徑的重要因素之一，我們在孔雀河上、中游的普惠-經濟牧場段的兩個斷面打井獲得了詳細的地質資料。每個斷面5口井，距離河道的垂直距離分別為50m，100m，200m，400m，800m。第一個斷面在普惠水庫上游，地質資料分析可見(圖2)，在距離河道200m的范圍內，從表層到深層土壤質地分別為黏土、粗沙、黏土、細沙，各層厚度略有差異，距離河道200m以外出現變化，從表層到深層土壤質地以黏土、細沙、粗沙、黏土變化，且各層出現不斷上移的趨勢。第二個斷面位于普惠水庫下游的團結荒地附件，由地質資料分析結果可知，在距離河道400米的范圍內，從表層到深層土壤質地以黏土、細沙、黏土、粗沙、細沙變化，且各層也出現不斷上移。400m以外變化更為復雜，無明顯規律，但細沙層厚度明顯變大(詳見圖2)。

(a)及不同深度的土壤質地；(b)普惠斷面；(c)團結荒地斷面圖2 鉆井斷面的位置

3.2 地下水位影響半徑計算

在調查中得知，孔雀河流域的普惠及經濟牧場一帶主要以種植棉花為主，其中，2008年以后開墾的耕地大多靠抽取地下水灌溉。結合棉花生長季400方/畝的灌溉定額(滴管)，全生長季需要灌溉16次。一次棉花生長季灌溉會導致水位下降3～5m不等。由于不同地下水埋深的差異，抽水過程中對周邊地下水的影響半徑上也是不同的。我們采用裘布依公式，對第三分水樞紐至普惠段、普惠水庫至經濟牧場段以及經濟牧場至66分水閘三個不同區段水井地下水位下降1m、3m、5m三個不同情境對周邊地下水的影響半徑進行了模擬，結果顯示，同一地下水埋深下，抽水降深越大，影響半徑越大。同一降深下，地下水埋深越大，影響半徑也越大。其中，普惠水庫至經濟牧場段由于地下水埋深較大，在灌溉期抽水對地下水的影響半徑也最大。在降深5m的情況下達到956.8m，是同情況的其他兩地影響半徑的2倍之多。

調查結果顯示，孔雀河流域第三樞紐至普惠水庫段，井灌面積約1.2×104畝，目前尚在使用的井約100眼，每年抽取地下水480×104方；普惠水庫至經濟牧場段，井灌面積約9.8×104畝，目前尚在使用的井約150眼，每年開采地下水3920×104方；經濟牧場至66分水閘段，井灌面積約5.0×104畝，地下水開采井約1000眼，每年開采地下水2000×104方。根據上述3個地段種植棉花的井灌面積，可計算出第三樞紐至普惠水庫段、普惠水庫至經濟牧場段、經濟牧場至66分水閘段三個區段單口井每年抽水量分別為4.8×104方、26×104方、20000方。根據抽水降深與當地的土壤質地可以計算出，灌溉期單口井抽水對該三個區段的影響半徑分別為776m，520m，136m。經濟牧場至66分水閘段雖然單口井的影響半徑較小，但是地下水井的分布密度較大，因此，總體對地下水位的埋深影響也是較大的。

表1 抽水降深對不同地下水埋深的影響半徑

*第三分水樞紐至普惠段的地下水位埋深，**普惠水庫至經濟牧場段的地下水位埋深，***尉犁縣的經濟牧場至66分水閘段的地下水位埋深。

4 結 論

影響半徑R是水文地質中的重要參數，在計算滲透系數K、開采井群間距的設計、地下水資源的評價、防護帶的圈定上具有重要的作用。裘布依公式根據抽水量，水位降深可以很方便求出在承壓井與潛水井的影響半徑R。

本文以孔雀河中上游為實例，模擬了第三分水樞紐至普惠段、普惠水庫至經濟牧場段、尉犁縣的經濟牧場至66分水閘段在不同地下水位埋深和不同抽水降深的影響半徑。結果表明：普惠水庫至經濟牧場段的影響半徑在1米、3米、5米的降深影響都是最大，且是第三分水樞紐至普惠段、普惠水庫至經濟牧場段影響半徑的2倍，甚至更大。在棉花全生長季，第三分水樞紐至普惠段的單口井的抽水量最大，對周邊地下水影響的范圍可達776m。

孔雀河流域的農業開墾面積已超水資源承載力，因此，應控制耕地擴張，加快實施退耕還水計劃，堅持以水定地、以水定結構、以水定發展的原則，取締非法機電井，實施壓采平衡。在生態輸水背景下，爭取盡快恢復地下水位，以維系沿河兩岸胡楊林生態系統的健康。

[1]張茹，劉水，馬璐瑤.淺談單井抽水的影響半徑計算[J].科技風，2015，10：126.

[2]胡月，劉國東，暴偉越，李俊.地下水環境影響評價中的“影響半徑”及其確定[J].山西科技，2015，06：59-62.

[3]陳雨孫，姚道聖，花仁榮，朱益春，周長瑚.試論“影響半徑”[J].勘察技術資料，1976，06：1-29.

[4]Jacob C E.Correlation of ground-water levels and precipitation on Long Island，New York，Part Theory[J].Am.Geophys.Union，1944(24)：564-573.

[5]Carrera J，Neuman S.P.Estimation of aquifer parameters under transient and steady state conditions：2.Uniqueness，stability，and solution algorithms[J].Water Resour.Res，1986，22 (2)：211-227.

[6]Boulton N S.Analysis of data from non-equilibrium pumping tests allowing for delayed yield from storage[C].Ice Virtual Library，1963：469-482.

[7]薛禹群.地下水動力學[M].北京：地質出版社，2009.

[8]陳亞寧，陳亞鵬，李衛紅，張宏鋒.塔里木河下游胡楊脯氨酸累積對地下水位變化的響應[J].科學通報，2003，09：958-961.

[9]劉加珍，陳亞寧，李衛紅，陳永金.荒漠河岸植被的受損過程與受損機理分析[J].地理學報，2006，09：946-956.

Simulation Analysis on the Radius of Influence of Ground Water Level Change Along the Kongque River，Xinjiang

LI Yupeng1Wumairjiang Wubli2HE Yu2LI Weihong1

(1.State Key Laboratory of Desert and Oasis Ecology，Xinjiang Institute of Ecology and Geography， Chinese Academy of Sciences，Xinjiang 830011 China； 2. Xinjiang Tarim River Basin Management Bureau，Korla 841000，China)

This paper through the Dupuit formula in the upper reaches of Kongque River as an example，simulation of the third water diversion project to the Pratt &；Whitney section，Pratt &；Whitney reservoir to economic pasture section，Yuli County Economy pasture to 66 sluice in different stages of different underground water level drawdown radius of Influence. The results show that：inclusive reservoir to economic pasture segment radius of influence in 1 meters，3 meters，the impact of the depth of 5 meters is the largest，and is the other two affected radius of more than 2 times or even greater. In the whole cotton growing season，third water diversion project to the pumping irrigation amount per mu inclusive maximum range，the impact on the surrounding groundwater is 776 meters. At present，the Peacock River water level decline，Populus death should be controlled，and cultivated land，accelerate the implementation of returning farmland to water plan，adhere to the water，with water，with water development，banning illegal electromechanical wells，in order to make the water level gradually restored，conducive to the development of economy and society.

Groundwater；Well；radius of Influence；Kongque River

國家科技支撐計劃項目(2014BAC15B02)和中國科學院“率先行動”計劃課題(TSS-2015-014-FW-2-3)

李玉朋，在讀博士生，從事氣候變化與水文過程研究

李衛紅，研究員，從事氣候變化與水文過程研究

X21

A

1673-288X(2017)02-0145-03

引用文獻格式：李玉朋 等.新疆孔雀河沿岸地下水位變化的影響半徑模擬分析[J].環境與可持續發展，2017，42(2)：145-147.