吐魯番盆地典型剖面的地下水徑流特征

張明江，蘇明磊，胡伏生，楚華兵(.新疆維吾爾自治區地質礦產勘查開發局第一水文工程地質大隊，新疆 烏魯木齊 83009；.中國地質大學，北京 00083)

吐魯番盆地典型剖面的地下水徑流特征

張明江1，蘇明磊1，胡伏生2，楚華兵1
(1.新疆維吾爾自治區地質礦產勘查開發局第一水文工程地質大隊，新疆 烏魯木齊 830091；2.中國地質大學，北京 100083)

為研究吐魯番盆地南北斷面方向上地下水的流動規律，沿吐魯番盆地中部南北剖面建立了二維地下水流動模型，采用二維地下水流數值模擬方法對剖面模型進行識別，在此基礎上分析了吐魯番盆地第四系地下水流動系統在剖面上的循環特征、流速分布特征、剖面上地下水的平均滯留時間和徑流強度。模擬結果表明，吐魯番盆地地下水流動系統存在局部流動系統、中間流動系統和區域流動系統，局部流動系統循環深度北盆地為600 m，南盆地為300 m，中間流動系統循環深度達1 000 m，區域流動系統循環深度達1 200 m。

吐魯番盆地；典型剖面；地下水徑流特征

1　研究區概況

1.1水文地質條件

據吐魯番盆地中部NS向地質剖面上3個典型勘探孔可知：北盆地TW-NC勘探孔揭露的第四系巖性自上而下依次為卵礫石、卵石、砂卵石；剖面中部TW-SE勘探孔揭露的第四系巖性自上而下依次為粉土、細砂、粉土、粉質粘土、粉土、粉質粘土；剖面南部TW-SS勘探孔揭露的第四系巖性自上而下依次為粉質粘土、粉土、粉質粘土、粉土、粉質粘土、粉土、粉質粘土、粉土、粉質粘土（圖1）。據地層層序和巖相分析，將剖面上第四系自上而下分為A層、B層、C層?日本國際航業株式會社.新疆吐魯番盆地水資源可持續利用項目最后報告書，2003。南盆地地層巖性和巖相具較好的連續性，采用上述地層層序劃分方法將剖面上的水文地質單元從上至下劃分為A層、B層、C層（圖2）。從剖面上看，地下水類型有潛水和承壓水，由北向南，含水層富水性從水量豐富到水量貧乏各等級均有分布?新疆地礦局第一水文地質大隊.吐魯番盆地煤炭基地地下水勘查，2010-2013。北盆地為單一結構的砂礫石潛水含水層分布區，大面積分布水量豐富的潛水含水層，水量中等和水量貧乏的含水層主要分布在盆地東部的局部地區。火焰山-鹽山構造缺口及其以南地區為多層結構的砂巖類含水層分布區，向南延伸至艾丁湖，含水層富水性由水量豐富逐漸過渡到水量貧乏，水量豐富和中等的含水層分布較廣。

圖1　吐魯番盆地典型剖面地層層序劃分圖Fig.1 Division of stratigraphic sequence of typical profile in Turpan basin

地下水從北部和中部獲得補給后向南徑流，沿

項目資助:2010年自治區地質勘查基金項目“新疆東疆地區煤炭基地地下水勘查吐魯番-哈密盆地地下水勘查”(2010005)資助

圖2　吐魯番盆地典型剖面水文地質單元劃分Fig.2 Division of hydrogeological units of typical profile in Turpan basin

1.2研究的必要性

Tóth地下水流動系統理論將區域地下水流動系統劃分為3個級次：局部流動系統、中間流動系統和區域流動系統[1]。吐魯番盆地降水稀少，區域水資源主要依靠出山口河流供給。從剖面上看，吐魯番地下水盆地缺少支持Tóth地下水流動系統理論的大氣降水的入滲補給，其地下水循環模式不能完全套用該理論，需根據發展完善后的區域地下水流理論進行研究[2]。對于地下水局部流動系統，地下水的補給、徑流、排泄的周期短，更新速度快；對于地下水區域流動系統，地下水補給、徑流、排泄的周期長，更新速度最慢；中間流動系統介于上述二者之間。研究吐魯番盆地地下水在斷面上的徑流特征，劃分地下水流動系統，對地下水開采方案的制定和地下水資源的保護具有重要的指導意義。

2　地下水位的變幅

根據2011年9月、2012年2月、2012年7月3次地下水位統測數據分析地下水水位的變化幅度。

從吐魯番盆地地下水2012年2月統測水位減去2011年9月統測水位所得差值來看，盆地周邊2月份的水位低于9月份水位，盆地中心則相反。盆地周邊是補給區，冬春季節地下水補給強度弱，在天然水力坡度作用下，地下水向盆地中心匯集，水位下降。盆地中心是地下水的天然排泄區和開采區，冬春季節地下水的開采量和蒸發量小，在盆地周邊地下水的補給作用下，盆地中心水位上升。在火焰山與艾丁湖之間，地下水位上升幅度大，一般大于10 m。從2012年7月統測水位減去2012年2月統測水位所得差值來看，盆地周邊夏季水位普遍高于冬季水位，盆地中心則相反。地表河流不發育的鄯善和艾丁湖南區地下水位上升點多，說明存在較強的暴雨洪流的補給作用。

對比上述水位變幅發現，火焰山南部的鄉鎮存在夏季地下水位降幅大于冬季水位升幅的現象，出現了局部超量開采問題。

3　典型剖面水文地質條件的概化

3.1典型剖面的選取

模擬要求所選取剖面的地下水流動特征要對吐魯番盆地的地下水流動系統具代表意義，剖面既能反映基本地形地質特征又必須經過區域地下水主要補給區和排泄區，此外，剖面水流流向需體現區域地下水的水流方向，即選取的剖面應盡量與等水位線垂直。據鉆孔和水文地質資料，選取的剖面如圖3。

3.2剖面水文地質條件概化

從區域水文地質條件和水文地質剖面可知，雖然剖面含水層巖性滲透系數不同，富水性差異較大，但無明顯承壓含水層存在，因此，將剖面巖層概化為潛水含水層和隔水基底兩層。典型剖面的巖性具明顯的水平分布特征，在垂向上巖性有分層性，將剖面概化為非均質各向異性。模擬僅研究剖面的水流流動特征，不進行水位動態的表現，模型概化為均衡穩定模型。

3.3剖面邊界條件概化

據模擬軟件本身特點，剖面模型的表層，在主要補給區（點）設置recharge和well項進行補給量的調控，在排泄區（點）設置general head項進行排泄量的控制，還有一些特殊邊界依據地勢條件和地表狀況設置為special head提供進行補排量，其他表層邊界僅設置為活動單元格。

為此，在吐魯番盆地典型剖面的表層邊界，山前補給帶設為活動單元格，以recharge和well的形式輸入到模型中，在火焰山溢出帶附近，設置為general head，艾丁湖附近設為special head，其他邊界區僅設為活動單元格。據巖性分布，將剖面南北部側邊界及下邊界定為隔水邊界。

圖3　吐魯番典型剖面水文地質剖面圖Fig.3 Hydrogeological section of typical profile in Turpan basin

4　建立模型

4.1數學模型

依據上述水文地質概念模型及邊界條件，建立非均質各向異性二維穩定流數學模型來模擬剖面含水層水流流動。

式中：Ω——剖面滲流區域；

h——含水層水位標高（m）；

Kx——水平方向的滲透系數（m/d）；

Kz——垂向上的滲透系數（m/d）；

w——源匯項（m3/d）;

h1——水頭邊界上水頭分布（m）；

T——含水層導水系數（m2/d）;

B1——水頭邊界；

B2——流量邊界；

q(x,z)——流量邊界單寬流量（m3/d），流入為正，流出為負，隔水邊界為0。

4.2研究區剖分

考慮到剖面橫縱差及單位距離內水流方向垂向上比橫向上變化大，剖分時網格縱向應比橫向小得多，網格大小設定為水平方向1 km，垂向20 m。Y軸厚度為50 m，斷面模型的總網格數為85列×90層＝7 650個，參與模擬的網格共3 164個。

4.3滲透系數分區

在參考水文地質剖面圖、鉆孔柱狀圖及抽水試驗的基礎上進行水平滲透系數分區，將特殊巖性互層段進行合并，水平滲透系數取相應的等效值，最后形成10個滲透系數分區。

含水層的各向異性通過水平滲透系數與垂向滲透系數的比值Kh/Kv表示，通過控制補給區向下滲流的深度來影響整個水流系統的循環深度。據不同巖性進行分區試驗確定最終的Kh/Kv值：在北部山區基巖和第四系用50、第三系用1 000、南部裸露基巖用200作為模型的垂向滲透設置。

4.4初始水頭

在穩定流計算里無需給初始水頭變量賦值，但設定一個與計算結果相近的初始值可提高計算效率，并可防止在模型運行中含水層的疏干。這里將剖分網格每一列的初始水頭設置成對應列的第一個網格中心點的高程值。

4.5模型運行與參數調整

在GMS軟件中利用MODFLOW模塊來運行模型以獲得剖面水頭分布[3]。調參過程中主要遵循：①各觀測點的計算水位與實測水位基本一致；②符合水均衡原則，即補排項和存量的變化平衡；③流場分布與實際水文地質特征相對應；④調參后的水文地質參數既要以實驗參數為依據，也要符合剖面巖性變化規律。

4.6模型識別

計算水位與實際水位對比 模型運行最終得到的典型剖面計算等水位線分布圖(圖4)。等水位線值的變化趨勢符合剖面總體地下水補徑排趨勢，即從補給區（點）到排泄區（點）水頭值減小，說明模型的模擬結果具一定可信度。

圖4　吐魯番典型剖面計算等水頭線分布Fig.4 Distribution of calculational equivalent water head of typical profile in Turpan basin

調整后的滲透系數與剖面巖性變化規律對比調整后的水平滲透系數分布情況與原來的水平滲透系數分區相比，總體上變化不大，僅在個別區域有變化。說明參數調整適當，數值模擬結果合理。

模型識別結果利用GMS中的MODFLOW模塊進行計算模擬，經多次參數調試和模擬識別，模型參數基本合理，能夠反映典型剖面流場特征。

5　剖面水流系統

5.1剖面地下水循環特征

采用MODPATH程序模塊生成地下水流線分布，經剖面地下水流系統數值模擬，吐魯番盆地剖面流場分布和地下水循環特征的分布情況見圖5。

沿用Tóth的理論，將剖面地下水流系統分為：局部循環、中間循環和區域循環3級循環系統。由圖3可見，火焰山以北和艾丁湖以南為地下水單一循環區，火焰山以南至艾丁湖之間為多級循環區。多級水流系統中的局部水流系統是灌溉水回滲與蒸發過程在空間上的表現。

圖5　吐魯番典型剖面流場和地下水循環系統Fig.5 Flow field and groundwater circulation system of typical profile in Turpan basin

5.2剖面流速分布特征

在局部流動系統范圍內，流速大；在南部基巖低山區以及剖面底部基巖區域的區域和中間流動系統范圍內，流速較小。地下水水平滲透速度較大地區分布在山前補給區及火焰山附近溢出區，最大可達8.0 m/d；其次分布在艾丁湖附近排泄區，約為5.0 m/d；在綠洲分布地區的小的局部循環范圍內，也有較大的滲透速度，約2.0 m/d。

5.3剖面地下水平均滯留時間分析

橫向上，在北部山前補給區、火焰山溢出帶以及艾丁湖地下水排泄區，地下水滯留時間短；在南部基巖低山區上部裂隙發育區域的地下水補給區，地下水滯留時間較短；在剖面底部基巖，巖層滲透系數小，地下水滯留時間長。垂向上，與循環系統分布圖對比，火焰山兩側地下水補給和排泄區、綠洲區、艾丁湖排泄區的局部循環系統內，地下水滯留時間短；南部基巖低山區上部地下水補給區局部循環系統內，地下水滯留時間較長；在剖面底部基巖的中間和區域流動系統內，地下水滯留時間長。

5.4剖面地下水徑流強度分析

潛水面水量交換 潛水面凈補排量曲線的變化情況和剖面的補排特征相對應：在北部山前，地下水接受暴雨洪流與河水的入滲補給，約占總補給量的35%；從山前到洪積扇扇沿，有河水或干渠的滲漏補給，補給強度相對較弱；地下水在火焰山北麓大量溢出，約占總排泄量的45%；在火焰山南麓，之前溢出的地下水重新轉化為地下水，約占總補給量的30%；在火焰山以南的小洪積扇扇緣也有一部分地下水溢出，約占溢出帶的20%；其它區域潛水凈補排量小，且補排交替出現。

地下徑流與深度 據剖面流場水循環特征及潛水面的主要補給、徑流、排泄點，選取了8個典型斷面，采用GMS中的BUDGET程序模塊進行地下水徑流強度分析。從斷面的徑流強度看，山前砂礫石層水流強度大；隨著山前入滲的補給，在斷面6處的流量達到最大；艾丁湖附近斷面1和斷面2的徑流強度最弱。

從徑流強度與高程關系看，隨著循環深度增加，徑流強度明顯減弱。北盆地積極循環帶的深度在600 m以淺，南盆地400 m以淺；局部循環系統徑流強度較大；區域循環系統徑流強度較小。

6　結論

（1）所選典型剖面能夠代表吐魯番盆地第四系地下水的循環特征。

（2）剖面上地下水流動系統劃分為局部流動系統、中間流動系統和區域流動系統。局部流動系統位于南盆地中間流動系統內；北盆地和南盆地地下水流動系統屬中間流動系統；從盆地北部山前至艾丁湖之間的深層地下水流動系統屬區域地下水流動系統。

（3）火焰山以北和艾丁湖以南為地下水單一循環區，火焰山以南至艾丁湖之間為多級循環區。

（4）在局部流動系統范圍內，流速大；在南部基巖低山區及剖面底部基巖區域的區域和中間流動系統范圍內，流速較小。

（5）橫向上，北部山前補給區、火焰山溢出帶及艾丁湖地下水排泄區的地下水滯留時間短；垂向上，火焰山兩側地下水補給和排泄區、綠洲區、艾丁湖排泄區的局部循環系統內，地下水滯留時間短。

（6）山前砂礫石層水流強度大，艾丁湖附近的徑流強度最弱。

（7）區域流動系統的循環深度達1 200 m；中間流動系統的循環深度達1 000 m；局部流動系統循環深度北盆地為600 m，南盆地為300 m。

[1]Tóth.A theoretical analysis of groundwater flow in small drainage basins[J].Journal of Geophysical Research，1963.

[2]蔣小偉，萬力，王旭升.區域地下水流理論進展[M].北京:地質出版社，2013.

[3]McDonald M.G.，A.W.Harbaugh.A Modular Three-Dimensional finite-difference ground-water flow model,U.S.Geological Survey Techniques of Water Resources Investigations Book 6.United States Government Printing Office,Washington,D.C Chapter A1, 1988.

Groundwater Runoff Characteristics of Typical Profile in Turpan Basin

Zhang Mingjiang1,Su Minglei1,Hu Fusheng2,Chu Huabing1
(1.The 1st Hydro-engineering geological Brigade,Xinjiang Bureau of geological and mineral resources,Urumqi, Xinjiang,830091,China;2.China university of Geosciences,Beijing,100083,China)

To study the groundwater flow law which in the north and south section of Turpan basin,two-dimensional groundwater flow model is established.It is along north and south section of the central Turpan basin.Using two-dimensional groundwater flow numerical simulation method to identify the profile model.Based on the model，this paper analyzes many Quaternary groundwater flow system characteristics in Turpan basin,such as circulation characteristics in profile,velocity distribution characteristics,average residence timeand groundwater runoff intensity in profile.Simulation results show that groundwater flow system in Turpan basin have three types,including local flow system，intermediate flow system and regional flow system.Local flow system circulation depth is 600 m in the north basin and 300 m in the south basin,the circulation depth of intermediate flow system and regional flow system is 1000 m and 1200 m respectively.

Turpan basin;Typical section;Groundwater runoff characteristic

1000-8845(2016)03-418-05

P641.2

A

2016-03-11;

2016-05-16;作者E-mail:728870395@qq.com

張明江（1965-），男，重慶人，高級工程師，1987年畢業于長春地質學院水文地質專業，現從事水文地質、工程地質、環境地質勘查工作途排泄于開采和蒸發蒸騰。第四系地下水的主要補給來源有：河流、渠系、水庫、田間灌溉水的入滲、補給，河谷潛流的側向流入補給，山區基巖裂隙水的側向流入補給，暴雨形成的洪水入滲補給。地下水的排泄方式有：機電井開采，坎兒井開采，泉水排泄，地下水蒸發蒸騰，艾丁湖是地下水的排泄中心。對吐魯番盆地地下水的研究，目前僅局限于從平面上對水文地質條件進行分析，尤其是對地下水徑流的研究僅從編制地下水等水位線圖的角度了解地下水在平面上的流向和水力坡度等參數，尚未開展過斷面地下水徑流特征的模擬分析，對剖面上地下水流動系統的劃分需做進一步研究。