馬曉菲，杜明亮，劉小煜，朱豫魯

(1.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，新疆 烏魯木齊 830052；2.新疆水利工程安全與水災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室，新疆 烏魯木齊 830052)

地下水資源是干旱-半干旱地區(qū)十分關(guān)鍵的生態(tài)環(huán)境因子，其影響著當(dāng)?shù)厝嗣裆鐣?huì)和經(jīng)濟(jì)的發(fā)展。而坎兒井是吐魯番市地下水資源的最主要來(lái)源之一[1]。隨著吐魯番市灌溉面積不斷擴(kuò)張，該區(qū)域需水量也日益增長(zhǎng)，加之該區(qū)域地表水資源分布不均，導(dǎo)致該區(qū)域過(guò)度開(kāi)發(fā)利用地下水，最終致使該區(qū)域地下水位持續(xù)下降，灌區(qū)與荒漠區(qū)之間生態(tài)植被退化，荒漠化問(wèn)題日益突出[2-7]。同時(shí)由于地下水位的快速下降，使得坎兒井不斷干涸。據(jù)吐魯番市水利局水科所最新一次的調(diào)查統(tǒng)計(jì)，吐魯番市的坎兒由1957年的1 237條，減少到只有214條，平均每年減少坎兒井約18條。坎兒井衰減速率過(guò)快已經(jīng)引起當(dāng)?shù)叵嚓P(guān)部門(mén)的重視，并開(kāi)始采取相應(yīng)的保護(hù)和治理措施[8-13]。

現(xiàn)以吐魯番市為研究區(qū)，在已有地質(zhì)、水文地質(zhì)資料和長(zhǎng)期觀測(cè)資料的基礎(chǔ)上，采用數(shù)值模擬方法建立適用于吐魯番市的地下水流數(shù)值模擬模型，利用該模型對(duì)影響坎兒井出水量各因素進(jìn)行計(jì)算分析，并進(jìn)行了定量評(píng)價(jià)，該結(jié)果可為坎兒井相關(guān)研究提供理論依據(jù)，同時(shí)也可為坎兒井的保護(hù)和恢復(fù)提供實(shí)際參考價(jià)值。

1 研究區(qū)概況

吐魯番市位于烏魯木齊以東、哈密市以西的吐魯番盆地內(nèi)。地理坐標(biāo)為北緯41°12′～43°40′，東經(jīng)87°16′～91°55′。吐魯番市北部為博格達(dá)山，海拔3 500～4 000 m，南部為覺(jué)羅塔格山，海拔600～1 500 m。盆地中部為火焰山和鹽山，二者將吐魯番盆地分為北盆地和南盆地。北盆地為山前沖洪積平原，地面以較大的傾斜度向南傾斜，傾斜斜率約為0.045。南盆地為沖積細(xì)土平原，地勢(shì)低洼，坡降相交北盆地較小，海拔在10～110 m之間。

由于地質(zhì)構(gòu)造的作用，火焰山及鹽山在盆地內(nèi)部隆起，相應(yīng)的北盆地和南盆地水文地質(zhì)單元也有較大的不同。根據(jù)前人的相關(guān)水文物探工作成果表明，南部和北部盆地都分布著第四系沉積物，厚度為200～700 m，地下水資源豐富，水文地質(zhì)剖面如圖1所示，其中：

圖1 吐魯番市水文地質(zhì)剖面圖

北盆地從北向南，地下水埋藏深度由大到小變化。在600 m高程附近，地下水的埋藏深度約為100 m。往南2～4 km，地下水埋藏深度約為50 m。在火焰山地區(qū)附近，地下水的埋藏深度迅速減小，在連木沁溝、勝金溝、葡萄溝等火焰山缺口帶處，地下水以泉水的形式從地下露頭。北盆地地下水類(lèi)型主要是潛水，承壓水僅分布在火焰山附近。

南盆地從北向南，地下水的埋藏深度也是由大到小變化。地下水埋藏深度大于50 m的區(qū)域主要是在火焰山南部一帶；地下水埋藏深度大于10～30 m的區(qū)域主要分布在南盆地的中部，例如鄯善縣魯克沁鎮(zhèn)、吐峪溝鄉(xiāng)、二堡鄉(xiāng)、三堡鄉(xiāng)等；地下水埋藏深度小于10 m的區(qū)域主要分布在南盆地南部靠近艾丁湖的部分區(qū)域。

2 數(shù)值模型

根據(jù)吐魯番市水文地質(zhì)資料，建立該區(qū)域的水文地質(zhì)概念模型。本次模型共建立三層：上部為潛水含水層，中部為潛水、承壓水混合含水層，下部為承壓水含水層。對(duì)模型進(jìn)行500×500 m網(wǎng)格剖分，南北向劃分為212行，東西向劃分為619列，共131 228個(gè)網(wǎng)格區(qū)域，根據(jù)邊界條件將部分區(qū)域設(shè)置為無(wú)效單元，計(jì)算范圍為有效單元，有效區(qū)單元數(shù)為48 813個(gè)。

吐魯番市的主要補(bǔ)給來(lái)源為山前側(cè)向和地表補(bǔ)給量，其中地表補(bǔ)給量主要為降雨入滲量、河道和渠系滲漏量、田間入滲量、庫(kù)塘入滲量以及井灌回歸補(bǔ)給量。用通用水頭邊界對(duì)山前側(cè)向補(bǔ)給量進(jìn)行賦值模擬；用River模塊對(duì)河道滲漏補(bǔ)給進(jìn)行模擬；用Recharge模塊模擬面狀補(bǔ)給量。

吐魯番市主要排泄量為潛水蒸發(fā)蒸騰量、人工開(kāi)采量、坎兒井排泄及泉水溢出。采用Well模塊模擬人工開(kāi)采量；采用Evapotranspiration模塊模擬潛水蒸發(fā)蒸騰量；使用Drain模塊模擬坎兒井排泄量及泉水溢出量，見(jiàn)圖2所示。

圖2 數(shù)值模型概化圖

3 結(jié)果分析與討論

3.1 坎兒井結(jié)構(gòu)設(shè)計(jì)與出水量關(guān)系

3.1.1 坎兒井斜率與出水量關(guān)系

坎兒井的斜率會(huì)影響坎兒井與地下水位相交的位置，進(jìn)而影響坎兒井的出水量。根據(jù)計(jì)算結(jié)果坎兒井出水量隨著其斜率的減小而增大，尤其是在坎兒井斜率減小到0.003后，其相對(duì)出水量變化驟然加大，由2.37×104m3增加到7.03×104m3(見(jiàn)表1和圖3)，但考慮坎兒井施工水平和地下水自流因素，一般坎兒井斜率不會(huì)比0.003小太多。

表1 坎兒井斜率與出水量關(guān)系表

圖3 坎兒井斜率與出水量變化關(guān)系曲線

3.1.2 坎兒井井徑與出水量關(guān)系

增加坎兒井井徑意味著可以增加坎兒井的過(guò)水?dāng)嗝婷娣e，相應(yīng)也會(huì)增加坎兒井的出水量。計(jì)算結(jié)果表明隨著坎兒井的井徑的不斷增大，從0.2 m增加到1 m，其出水量略有增加，由21.51×104m3至21.62×104m3(表2和圖4)，但總體變化不大，可認(rèn)為井徑的變化對(duì)坎兒井出水量影響不大。

圖4 坎兒井井徑與出水量變化關(guān)系曲線

表2 坎兒井間距與出水量關(guān)系表

3.1.3 坎兒井間距與出水量關(guān)系

在一定范圍內(nèi)若存在多條坎兒井，它們會(huì)不均等地分配地下入滲補(bǔ)給水以及地下水徑流量，以至于出現(xiàn)不同的坎兒井出水量不一的現(xiàn)象，因此進(jìn)行計(jì)算分析。分析結(jié)果表明：坎兒井的出水量隨著井間距的減小而減小。當(dāng)井間距介于200～500 m之間，坎兒井出水量的變化比率最大；當(dāng)井間距大于1 km時(shí)，其出水量變化比率值很小(表3和圖5)，其對(duì)坎兒井的影響可以忽略不計(jì)。

表3 坎兒井間距與出水量關(guān)系表

圖5 坎兒井間距與出水量關(guān)系曲線

3.2 水文地質(zhì)參數(shù)與出水量的關(guān)系

3.2.1 滲透系數(shù)與出水量關(guān)系

滲透系數(shù)反映含水層透水性能，會(huì)影響地下水量和水位的相應(yīng)變化。根據(jù)分析結(jié)果，坎兒井的出水量隨著水平滲透系數(shù)的不斷增大，出現(xiàn)先增大后減小的現(xiàn)象，并在水平滲透系數(shù)達(dá)到1 m/d時(shí)，坎兒井的出水量達(dá)到最大，最大值為56.34×104m3(表4和圖6)；滲透系數(shù)增大時(shí)，坎兒井出水量減小，主要由于滲透系數(shù)增大，從而引起地下水位下降。

表4 滲透系數(shù)與出水量關(guān)系表

圖6 水平滲透系數(shù)與出水量關(guān)系曲線

3.2.2 給水度與出水量關(guān)系

給水度反映含水層的給水能力，該參數(shù)會(huì)直接影響坎兒井的出水量。根據(jù)分析結(jié)果，當(dāng)給水度從0.01變化到0.2時(shí)，坎兒井出水量從16.17×104m3增加到36.98×104m3，二者呈正相關(guān)關(guān)系(表5和圖7)。

表5 給水度與出水量關(guān)系表

圖7 給水度與出水量關(guān)系曲線

3.3 地表水體入滲與出水量關(guān)系

吐魯番市地表水體入滲補(bǔ)給量主要包括降雨入滲量、渠系入滲量、田間入滲量、地表河流入滲量以及井灌回歸補(bǔ)給量，吐魯番市降雨量稀少，該項(xiàng)可以忽略不計(jì)。分析結(jié)果表明，地表水體入滲量越多，坎兒井的出水量越大；但灌溉入滲量的變化值遠(yuǎn)比坎兒井出水量的變化值大得多，當(dāng)灌溉入滲量由9.12×104m3增加到456.25×104m3，坎兒井出水量?jī)H僅變化了10.62×104m3(表6和圖8)，說(shuō)明灌溉入滲量對(duì)坎兒井出水量雖有影響，但是影響相對(duì)并不大。

表6 灌溉入滲量與出水量關(guān)系表

圖8 地表水體入滲量與出水量關(guān)系曲線

3.4 地下水變化與出水量關(guān)系

3.4.1 側(cè)向補(bǔ)給量與出水量關(guān)系

側(cè)向補(bǔ)給量作為山區(qū)進(jìn)入平原區(qū)的不重復(fù)補(bǔ)給量，在吐魯番市地下水補(bǔ)給量中占有一定的比例，而側(cè)向補(bǔ)給量隨著山區(qū)降雨入滲量隨季節(jié)變化，因此也會(huì)對(duì)坎兒井出水量產(chǎn)生一定的影響。分析結(jié)果表明，側(cè)向補(bǔ)給量與坎兒井出水量呈正相關(guān)關(guān)系(圖9)，但是由于側(cè)向補(bǔ)給量占總補(bǔ)給量的比例不大，所以從總體上看對(duì)坎兒井出水量的影響并不大，表7表明側(cè)向補(bǔ)給量從16.31×104m3增加到90.61×104m3，坎兒井出水量只變化了2.3×104m3，因此側(cè)向補(bǔ)給量對(duì)坎兒井出水量的影響較小。

表7 地下水變化與出水量關(guān)系

圖9 側(cè)向補(bǔ)給量與出水量關(guān)系

3.4.2 機(jī)井與出水量關(guān)系

機(jī)井對(duì)坎兒井出水量的影響主要表現(xiàn)在以下兩個(gè)方面：其一是機(jī)井的抽水量越大，在吐魯番市這樣水資源匱乏的城市，當(dāng)開(kāi)采量大于可開(kāi)采量時(shí)，會(huì)造成地下水位持續(xù)下降，進(jìn)而導(dǎo)致坎兒井出水能力大幅下降甚至干涸。根據(jù)分析結(jié)果，當(dāng)機(jī)井抽水量增加206%時(shí)，坎兒井出水量減小39%(表8)，二者呈線性相關(guān)的關(guān)系，見(jiàn)圖10。

表8 機(jī)井抽水量與坎兒井出水量關(guān)系表

圖10 機(jī)井抽水量與出水量關(guān)系曲線

其次，機(jī)井和坎兒井之間的距離也會(huì)影響坎兒井的出水量，機(jī)井離坎兒井集水段的距離越近，坎兒井集水段附近的地下水位越低，進(jìn)而坎兒井的出水量越小。根據(jù)分析結(jié)果，集水段的機(jī)井距離坎兒井越遠(yuǎn)，坎兒井的出水量越大，當(dāng)距離大于1 km時(shí)，對(duì)坎兒井的出水量沒(méi)有影響。如果機(jī)井在輸水段，無(wú)論距離遠(yuǎn)近都對(duì)坎兒井出水量影響不大，具體見(jiàn)表9和圖11。

表9 機(jī)井、坎兒井距離與坎兒井出水量關(guān)系表

圖11 機(jī)井抽水量與出水量關(guān)系曲線

4 結(jié)語(yǔ)

(1)影響坎兒井出水量的主要因素可歸納概括為坎兒井結(jié)構(gòu)、水文地質(zhì)參數(shù)、地表水體入滲以及地下水變化四個(gè)方面。對(duì)坎兒井出水量的影響由大至小依次為：坎兒井結(jié)構(gòu)、地下水變化、水文地質(zhì)參數(shù)、地表水體入滲。

(2)坎兒井結(jié)構(gòu)設(shè)計(jì)中對(duì)出水量影響較大的主要是坎兒井斜率和坎兒井間距?？矁壕甭试?.003附近時(shí)對(duì)出水量影響最為顯著。坎兒井間距在1 km以上時(shí)，對(duì)出水量幾乎沒(méi)有影響。

(3)水文地質(zhì)參數(shù)對(duì)出水量影響較大的主要是給水度和滲透系數(shù)。給水度與坎兒井出水量呈正線性相關(guān)關(guān)系。滲透系數(shù)小于1 m/d時(shí)與出水量呈正相關(guān)關(guān)系，大于1 m/d時(shí)與出水量呈負(fù)相關(guān)關(guān)系。

(4)地表水體入滲對(duì)出水量影響較大的主要是地表水體入滲補(bǔ)給量。地表水體入滲量增加，坎兒井出水量也會(huì)相應(yīng)增加，二者呈線性相關(guān)關(guān)系。但入滲量對(duì)坎兒井出水總量的影響相對(duì)并不大。

(5)地下水變化對(duì)出水量影響較大的因素依次為機(jī)井位置、機(jī)井抽水量以及側(cè)向補(bǔ)給量的變化。機(jī)井位于坎兒井集水段時(shí)對(duì)坎兒井出水量影響最大，但距離大于1 km時(shí)幾乎沒(méi)有影響。機(jī)井抽水量、側(cè)向補(bǔ)給量與坎兒井出水量呈正線性相關(guān)關(guān)系。