李　飛，王科奇

(合肥工業(yè)大學土木與水利工程學院，安徽 合肥 230009)

?

基于arcgis的淮北平原區(qū)域潛水蒸發(fā)量計算

李飛，王科奇

(合肥工業(yè)大學土木與水利工程學院，安徽 合肥 230009)

潛水蒸發(fā)是陸面水量損耗的重要部分，對農作物生長有重要影響。區(qū)域潛水蒸發(fā)量的計算對區(qū)域水資源的合理配置及水均衡研究有重要意義。本文利用安徽省淮北平原151個淺層孔隙水水位監(jiān)測點及周邊分布的8個水面蒸發(fā)站點某月監(jiān)測數(shù)據(jù)，通過arcgis并插值生成柵格文件。用阿維揚諾夫公式計算淮北平原該月區(qū)域潛水蒸發(fā)量為1.70E8 m3；利用Visual Modflow建立淮北平原區(qū)域地下水數(shù)值模型，計算數(shù)值模型區(qū)域潛水蒸發(fā)量。兩者進行比對分析，驗證利用arcgis進行不同時間尺度和空間尺度上的區(qū)域潛水蒸發(fā)量這一方法的可行性。

arcgis；區(qū)域潛水蒸發(fā)；阿維揚諾夫公式

潛水含水層及潛水面以上的非飽和帶處于的近地表環(huán)境，與土壤、水、空氣以及人類等生物聯(lián)系最為密切，這些要素相互作用的三類過程包括生態(tài)過程、生物地球化學過程和水文過程[1]。其中水文過程通過潛水的水位變動及土壤水汽運移使得物質和能量在時間和空間上得意重新分布，潛水蒸發(fā)就是該過程的重要環(huán)節(jié)。

潛水蒸發(fā)是指潛水通過毛管作用向非飽和帶輸送水分，并通過土壤蒸發(fā)進入大氣的過程。這一過程伴隨著能量和物質傳遞，蒸發(fā)的水分可以補給非飽和帶農作物根系[2]，同時是區(qū)域水均衡中重要的排泄項之一[3]，這對農作物生長及區(qū)域潛水含水層地下水資源量的評價、增加區(qū)域水均衡中其他各項的計算的可靠性有重要作用。但是目前潛水蒸發(fā)計算的研究多是在特定時間和地域條件下，在不同時間尺度和空間尺度上的潛水蒸發(fā)量計算還需要更多的研究。

生產實踐中常用的潛水蒸發(fā)經驗公式有以下幾種[4]：

(1)

葉水庭公式：E=E0e-aΔ

(2)

(3)

式中：E為潛水蒸發(fā)強度，m/d；E0為水面蒸發(fā)強度，m/d；Δ為埋深，m；Δ0為潛水蒸發(fā)極限埋深，m；n為與包氣帶土質、氣候有關的蒸發(fā)指數(shù)，一般取1～3；k為與土質有關的潛水蒸發(fā)經驗系數(shù)；a是經驗常數(shù)，與土質和地下水埋深有關；r、b是擬合系數(shù)。其中，阿維里揚諾夫公式、葉水庭公式被《中國水資源調查和評價工作技術細則》推薦為計算潛水蒸發(fā)的經驗公式。

本文選用阿維揚諾夫公式[5]，根據(jù)研究區(qū)監(jiān)測站點的數(shù)據(jù)，利用arcgis將研究區(qū)域的監(jiān)測點的水面蒸發(fā)強度和潛水埋深進行插值生成柵格文件，通過柵格轉點，得到單位柵格面積的平均水面蒸發(fā)強度和潛水埋深，利用屬性表進行字段計算，最終求得區(qū)域潛水蒸發(fā)量。

1　研究區(qū)基本概況

安徽省淮北平原面積37 411 km2。地處黃淮海平原南部，包括淮河以北和霍邱、壽縣北部。為北亞熱帶和暖溫帶半濕潤季風氣候區(qū)。冬季干旱少雨，夏季炎熱多雨，四季分明，氣候溫和。該區(qū)多年平均氣溫14℃～15℃,由南向北逐漸降低。年際變化不大，年內最冷1月份為0℃左右,最熱7月份在27℃左右；光照充足,多年平均日照時數(shù)為2 085 h；無霜期較長,多年平均212 d。由于該區(qū)處在南北分界線上。屬冷暖氣團經常交鋒過渡地帶，年均降水量700～900 mm，年均蒸發(fā)量1 000～1 300 mm。汛期6～9月降水量占年總量的63%左右,多年蒸發(fā)量為1 181.3 mm。地下水資源豐富，以松散巖類孔隙含水巖組分布最廣。一般厚200～600 m，東部小于100 m，西部可達800 m。地表水屬淮河水系。淺部地下水資源豐富，埋藏淺。平原區(qū)松散巖類淺層孔隙潛水主要埋藏于全新世古河道砂層中，枯水期自西北往東南由3～4 m減至1～2 m，河間1～3 m，濱河2～4 m，黃河古道和山麓地帶3～8 m，豐水期可普遍上升1～2 m。目前，淺層孔隙水仍然是淮北平原農業(yè)灌溉和農村人畜飲水的重要水源。

2　利用arcgis計算潛水蒸發(fā)量

2.1潛水水位變化及分布

根據(jù)淮北平原分布的151個淺層孔隙水水位監(jiān)測點的監(jiān)測數(shù)據(jù)，某月的月平均水位利用arcgis插值成柵格文件，繪成潛水水位等值線(圖1a)，利用地表高程與潛水位求得淮北平原的潛水埋深分布，生成埋深等值線及其柵格文件，(圖2a)根據(jù)淮北平原及周邊分布的8個蒸發(fā)站點監(jiān)測數(shù)據(jù)，繪成蒸發(fā)量等值線(圖1b)，利用arcgis插值成柵格文件(圖2b)。

(a)潛水埋深等值線圖

(b)水面蒸發(fā)等值線圖

2.2字段計算

每個柵格文件的面積均為992.47× 992.47 m2。將得到潛水埋深等值線和蒸發(fā)量等值線的柵格文件轉為點文件，每個柵格中心的一個點的蒸發(fā)量和潛水埋深代表該柵格所覆蓋位置的單位蒸發(fā)量和潛水埋深?；幢逼皆恋闹饕愋蜑樯敖谕梁忘S潮土，其中砂姜黑土約占54%，經過大量的潛水蒸發(fā)實驗，淮北平原的潛水發(fā)埋深的經驗值為6～7 m，此處設定為7 m。阿維揚諾夫公式中理論潛水蒸發(fā)極限埋深為4 m。蒸發(fā)指數(shù)定為6.81。各柵格文件分別求的單位柵格的潛水蒸發(fā)量，然后求和計算得淮北平原潛水蒸發(fā)總量為1.70E8 m3。

(a)潛水埋深柵格圖

(b)水面蒸發(fā)柵格圖

3　VisualModflow模型驗證

本文采用Visual MODFLOW 建立數(shù)值模擬模型。Visual Modflow中潛水線性蒸發(fā)的方程表述為：

(4)

式中：E為蒸發(fā)強度，mm/d；h為地下水水位標高，m；hs為地面標高，m；d為蒸發(fā)極限埋深，m。

3.1水文地質概念模型

模型以安徽境內淮河以北區(qū)域為模擬區(qū)，模擬區(qū)域面積約為37 411 km2，全區(qū)幾乎為松散巖類孔隙含水巖組所覆蓋，以埋深40 m且分布穩(wěn)定之粘性土為界，大致可分為淺層和深層兩個部分。淺層大部由上更新統(tǒng)亞粘土、亞砂土、粉砂和細砂組成，僅山丘坡麓為上更新統(tǒng)下部之粘性土，基巖侵蝕基準面上覆有下更新統(tǒng)礫砂，北部黃泛區(qū)及河谷地帶為全新統(tǒng)砂性土和粘性土。根據(jù)研究區(qū)的水文地質條件，將模擬區(qū)南側淮河定為水頭邊界。

3.2數(shù)學模型

根據(jù)水文地質概念模型及含水層水力性質，將模擬區(qū)地下水流概化成非均質各向異性非穩(wěn)定三維地下水流系統(tǒng)，并建立相應的數(shù)學模型。

(5)

式中：Kx、Ky、Kz為滲透系數(shù)在x、y、z方向的分量(m/d)，(假定滲透系數(shù)主軸方向與坐標軸的方向一致)；H為地下水水頭(m)；M為含水層厚度(m)；W為單位體積流量，用以代表流進源或流出匯的水量；μs為含水層的儲水率(1/m)；H0為初始水頭(m)； t為時間(d)；D為模擬區(qū)范圍；Γ為二類邊界；q為邊界流量。

3.3模型識別與驗證

通過模型參數(shù)的調整，六個潛水位觀測孔的實測值和計算值的基本擬合，模型識別曲線見圖3。

圖3　模型擬合曲線

可以利用模型水均衡模塊來驗證基于arcgis的區(qū)域潛水蒸發(fā)值。根據(jù)模型運行結果，該月淮北平原潛水蒸發(fā)量為1.47E8 m3。

4　結果分析

根據(jù)計算，已知基于arcgis利用阿維揚諾夫公式求得淮北平原2010年12月區(qū)域潛水蒸發(fā)量為1.7E8 m3，利用Visual Modflow模型求得的區(qū)域同期潛水蒸發(fā)量為1.47×108m3，兩者相差13.5%。Visual Modflow蒸發(fā)模塊中的潛水蒸發(fā)計算公式均是關于水面蒸發(fā)強度、水位埋深、蒸發(fā)極限埋深的經驗公式，兩者產生差異的原因主要是阿維揚諾夫公式中引入潛水蒸發(fā)指數(shù)，這一指數(shù)讓該公式適用范圍更廣。因此，基于arcgis利用阿維揚諾夫公式計算區(qū)域潛水蒸發(fā)量準確性有所增加。

本文以淮北平原為例，基于arcgis利用阿維揚諾夫公式計算了淮北平原某一月份的區(qū)域潛水蒸發(fā)量，為計算不同時間尺度和空間尺度上的區(qū)域潛水蒸發(fā)量提供了一種方法。經其與識別過的數(shù)值模型計算的蒸發(fā)量值比較，具有一定若可獲得計算區(qū)域內土地類型及土壤分類資料，根據(jù)實測資料分別為潛水蒸發(fā)經驗公式的參數(shù)賦值，利用arcgis進行字段分類計算，可使區(qū)域潛水蒸發(fā)量的計算更加精確。

[1]楊建鋒,張翠光.地球關鍵帶:地質環(huán)境研究的新框架[J].水文地質工程地質.2014,41.

[2]雷志棟,胡和平,楊詩秀,et al.塔里木盆地綠洲耗水分析[J].水利學報.2006,37(12):1470-1475.

[3]岳衛(wèi)峰,楊金忠,高鴻永,et al.內蒙河套灌區(qū)義長灌域水均衡分析[J].灌溉排水學報.2004,23(6):25-28.

[4]束龍倉,荊艷東,黃修東,et al.改進的無作物潛水蒸發(fā)經驗公式[J].吉林大學學報 (地球科學版).2012,42(6):1859-1865.

[5]劉淼,郝振純,王加虎,et al.淮北平原裸地潛水蒸發(fā)計算公式比較[J].吉林大學學報：地球科學版.2008,38(4):649-653.

Calculation on Huaibei Plain Regional Groundwater Evapotranspiration Based on Arcgis

LiFei，Wang Ke-qi

(School of Civil Engineering,Hefei University of Technology，Hefei 230009，Anhui)

The phreatic water evaporation is an important part of the land surface loss of water,it have a major effect on crop growth.It is important that Regional phreatic water evaporation calculation to the rational allocation of regional water resources.In this paper,we use Huaibei Plain 151 phreatic water level monitoring stations and 8 water surface evaporationmonitoring stationsdata a month and interpolated to generate raster files by arcgis.Then we get theHuaibei Plain region phreatic water evaporation the month is 1.7E8m3calculatedwith Aviriyanover’s empirical formula; Compared the value with the result calculated by regional groundwater numerical modelestablished with Visual Modflow,we can verify the feasibility of the phreatic water evaporationcalculate method.

Arcgis；region phreatic evaporation and Aviriyanover’s empirical formula

2016-05-03

李飛(1991-)，男，內蒙古鄂爾多斯人，在讀碩士研究生，主攻方向：地下水資源評價與保護。

P641.2

A

1004-1184(2016)05-0010-02