盧昱含，鮮 陽，靳孟貴，白宏偉，劉亞磊，劉延鋒

(1.中國地質(zhì)大學(xué)(武漢)盆地水文過程與濕地生態(tài)恢復(fù)學(xué)術(shù)創(chuàng)新基地，湖北 武漢 430074；2.中國地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院，湖北 武漢 430074；3.中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心,湖北 武漢 430205)

地表水與地下水聯(lián)系緊密，兩者之間的相互作用與轉(zhuǎn)化關(guān)系復(fù)雜，密切影響著潛流帶水分、物質(zhì)與能量的交換過程[1-3]。因此，研究地表水與地下水的相互作用對于水資源評價、水資源開發(fā)管理、水污染防治和生態(tài)系統(tǒng)保護均具有重要的意義[4-5]。

干旱地區(qū)人類活動引起的土地利用和渠系滲漏的變化對地表水與地下水相互作用的影響強烈，而引水灌溉的渠道入滲是干旱地區(qū)地表水與地下水相互作用的主要形式[6]。近年來，干旱地區(qū)地表水與地下水的相互作用研究受到國內(nèi)外研究機構(gòu)的廣泛重視[7-8]。干旱地區(qū)水資源開發(fā)利用與調(diào)度配置的關(guān)鍵和基礎(chǔ)是判定河流和渠道與地下水的轉(zhuǎn)化關(guān)系，量化兩者間的轉(zhuǎn)化量及其時空變化規(guī)律[9]。

干旱地區(qū)淺層地下水水溫相對恒定，而河流等地表水水溫受氣溫的影響，一般具有較為明顯的季節(jié)性變化和日變化規(guī)律，地表水與地下水發(fā)生相互作用的過程中，潛流帶溫度場會發(fā)生明顯的變化[10-12]。近年來，隨著溫度場與滲流場耦合研究的深入、溫度數(shù)據(jù)監(jiān)測成本的降低[13]以及多孔介質(zhì)水熱耦合運移模擬軟件的大量開發(fā)，運用水熱運移數(shù)值模擬定量評價地表水與地下水相互轉(zhuǎn)化關(guān)系的研究有了較大的發(fā)展。如Su等[14]和Constantz[15]利用VS2DHI軟件模擬了河岸剖面二維溫度場，并結(jié)合觀測井中監(jiān)測的季節(jié)性地下水溫度估算了含水層的水力傳導(dǎo)系數(shù)和地表水與地下水的相互轉(zhuǎn)化速率；任杰等[16]和潘維艷等[9]利用水熱運移模擬分別研究了大壩與農(nóng)渠的滲漏問題；Ju等[17]結(jié)合沙箱試驗與數(shù)值模擬，分別采用實測溫度、水頭數(shù)據(jù)分析評價了數(shù)值模擬結(jié)果，結(jié)果表明水頭與溫度數(shù)據(jù)結(jié)合可提高滲流模擬結(jié)果的可信度和精度。但目前針對干旱渠道水文動態(tài)劇烈變化條件下渠道地表水與地下水相互轉(zhuǎn)化關(guān)系和轉(zhuǎn)化水量的時空變化規(guī)律方面的研究成果較少，且轉(zhuǎn)化量評價存在很大的不確定性。

本文以新疆瑪納斯河流域炮臺鎮(zhèn)附近頭道溝渠道的排水干渠作為研究對象，通過野外原位試驗實時監(jiān)測渠道橫斷面地下水溫度、地下水水位、渠道地表水水位和不同埋深處河床沉積物的溫度變化，分析地表水與地下水相互轉(zhuǎn)化關(guān)系變化前后溫度的規(guī)律特征，并建立二維均質(zhì)多孔介質(zhì)變飽和土壤水熱運移耦合數(shù)值模型，模擬了渠道附近地下水的水熱運移，反演土壤飽和水力參數(shù)，定量計算渠道地表水與地下水的轉(zhuǎn)換量，為渠道滲漏定量估算和地表水與地下水聯(lián)合調(diào)度提供參考。

1 研究區(qū)域概況

瑪納斯河流域位于新疆天山北麓中段、準(zhǔn)噶爾盆地南緣，地處北緯43°27′～45°21′、東經(jīng)85°01′～86°32′，地勢自東南向西北緩緩傾斜。流域由東向西分別由塔西河、瑪納斯河、寧家河、金溝河、巴音溝河以及相連的五片沖積平原組成。從南向北地貌具有明顯的分帶性，可以劃分為5個地貌帶：山地帶、山前褶皺丘陵帶、串珠狀沖積扇帶、平原曲流帶、瑪納斯河尾閭湖泊和沙漠。

流域?qū)俚湫偷母珊祪?nèi)陸盆地氣候，冬冷夏熱，年(日)溫差大，光照充足，蒸發(fā)強烈[18]。該流域年平均氣溫為9.1℃；年平均降水量為110～200 mm，降水分布不均,垂直分帶特征明顯，主要集中在5～8月份，約占年降雨量的67%；年均蒸發(fā)量為1 500～2 000 mm[19]。流域地表水補給主要來自于降水及冰川融化，地下水補給主要來自于通過河流的垂直滲透、渠道滲漏、灌溉回灌、山前側(cè)向地下水補給[20-21]。該地區(qū)平原區(qū)地表水與地下水的相互作用同時受到自然條件和人類活動的影響，兩者之間的關(guān)系極為復(fù)雜。

頭道溝渠道的排水干渠位于炮臺鎮(zhèn)附近(N44°45′～44°50′，E85°30′～85°35′)，地處平原灌溉區(qū)，見圖1。

圖1 研究區(qū)地理位置Fig.1 Location of the research area

根據(jù)研究區(qū)左右岸2m深鉆孔的土壤巖性數(shù)據(jù)，對比白宏偉等[22]對瑪納斯河流域二維巖性結(jié)構(gòu)剖面的研究成果，認為研究區(qū)土壤質(zhì)地分布較為均一，大部分為粉質(zhì)壤土，受渠道及氣候的影響，入滲水水溫較低，低溫水入滲形成的溫度變化信號較為明顯。頭道溝渠道兼顧排水與灌溉功能，每年4月初上游開始向渠道放水，9月底停止供水，渠道地表水水位呈現(xiàn)規(guī)律性的升降變化。

2 研究方法

2.1 溫度和水位監(jiān)測

頭道溝渠道屬未襯砌的典型農(nóng)渠，本研究選定垂直頭道溝渠道的斷面，測定渠道斷面幾何形態(tài)，并布置監(jiān)測井、水位和溫度傳感器。將監(jiān)測井G1底部以下3.94 m處水平線定為基準(zhǔn)線，頭道溝渠道監(jiān)測斷面和監(jiān)測井布置見圖2。

圖2 頭道溝渠道監(jiān)測斷面和監(jiān)測井布置Fig.2 Layout of the monitoring cross sections and wells of Toudaogou canal

監(jiān)測井管為4 m長的PVC管，濾水管長為0.5 m。河岸監(jiān)測井(G1、G2、G3、G4)的濾水管布置在監(jiān)測井底部以上1 m處，每個監(jiān)測井內(nèi)均布置一個Solinst Solinst Levelogger 2參數(shù)探頭(壓力、溫度)自計水位和溫度。水位測量精度為±0.05 m，分辨率為0.001 m；溫度測量精度為±0.2℃，分辨率為0.01℃。

河床監(jiān)測井(G5、G6)的濾水管位于河床下1 m處，其中監(jiān)測井G5內(nèi)布置一個Solinst Levelogger 2參數(shù)探頭自計水位和溫度；監(jiān)測井G6內(nèi)布置一個HOBO溫度記錄儀記錄30 cm、50 cm、70 cm、100 cm埋深處河床沉積物的溫度。水位測量精度為±0.05 m，分辨率為0.001；溫度測量精度為±0.2℃，分辨率為0.02℃。

除圖示監(jiān)測井外，另布置一個Solinst Levelogger 2參數(shù)探頭監(jiān)測渠道地表水水位和溫度，并使用數(shù)據(jù)采集儀記錄各測點的溫度和水位，采集時間間隔為15 min，自2015年7月27日22∶00至9月2日22∶00連續(xù)監(jiān)測。

2.2 水文地質(zhì)概念模型

由于模型水力參數(shù)受河床沉積物巖性結(jié)構(gòu)和飽和度變化的影響可能會呈現(xiàn)幾個數(shù)量級的變化[10]，單純依靠水位校正數(shù)值模型會使得模型水力參數(shù)率定存在較大的不確定性，而熱傳導(dǎo)系數(shù)的變化范圍較小。因此，本文通過建立水熱運移耦合模型，利用水頭和溫度數(shù)據(jù)約束和校正模型水力參數(shù)，以降低模型的不確定性。

地表水與地下水相互作用時，土壤中水分運移與熱量傳輸會產(chǎn)生相互影響，土壤中水分運移會引起土壤的體積熱容量和熱傳導(dǎo)系數(shù)的變化，改變熱量傳輸與交換，進而改變土壤中溫度的分布。由于兩者之間的耦合作用過程非常復(fù)雜，為了研究渠道地表水與地下水的水熱運移規(guī)律，可忽略溫度變化引起水的性質(zhì)變化，將監(jiān)測斷面概化為二維均質(zhì)多孔介質(zhì)變飽和水熱運移耦合模型。

以頭道溝渠道斷面形態(tài)為上邊界，以監(jiān)測井G1埋深10.06 m處水平線為底邊界，模型深度約為10 m；以監(jiān)測井G1、G4為模型左、右邊界，模型水平長度為31.56 m。建立的頭道溝渠道研究斷面水文地質(zhì)概念模型，見圖3。

圖3 頭道溝渠道研究斷面水文地質(zhì)概念模型Fig.3 Conceptual hydrogeological model of the research cross section of Toudaogou canal

滲流場邊界條件：地表邊界B1、B3為通量邊界 (Neumann邊界)，監(jiān)測期間幾乎無降雨，邊界賦值為結(jié)合經(jīng)驗公式和實測水面蒸發(fā)E0數(shù)據(jù)計算出的邊界蒸發(fā)通量，由付秋萍等[23]研究可知，清華大學(xué)雷氏公式在新疆地區(qū)具有廣泛的適用性，參考該研究中粉壤土參數(shù)擬合成果，以函數(shù)形式賦值輸入；渠道邊界B2為已知變水頭邊界，輸入實測渠中水位值(Dirichlet邊界)；左右邊界B4、B5為已知變水頭邊界，分別輸入井孔G1、G4的實測水頭值(Dirichlet邊界)；底邊界B6為隔水邊界(Neumann邊界)。初始條件為實測壓力水頭數(shù)據(jù)的線性插值。

溫度場的邊界條件：地表邊界B1、B3為變溫度邊界，實測日平均氣溫數(shù)據(jù)(Dirichlet邊界)；渠道邊界B2為變溫度邊界，輸入實測渠道水溫數(shù)據(jù)(Dirichlet邊界)；底邊界B6為定溫度邊界，T=10℃(Dirichlet邊界)。初始條件為實測溫度數(shù)據(jù)的線性插值。

2.3 數(shù)學(xué)模型

2.3.1 土壤水分運動方程

假定研究斷面土壤為均質(zhì)各向同性，設(shè)z軸垂直向上，x軸水平向右，二維均質(zhì)各向同性土壤水分運動的偏微分方程及其定解條件可表示如下：

(1)

式中：t為時間變量(s)；x、z均為空間坐標(biāo)(m)；h(x,z,t)為壓力水頭(m)；c(h)為比水容量,c(h)=?/?h；K(h)為非飽和土壤水力傳導(dǎo)系數(shù)(m/d)，是h的函數(shù)；h0(x,z)為滲流場初始地下水水位(m)；h1(x,z)為第一類邊界Γ1水位值(m)；ε(x,z,t)為第二類邊界Γ2通量(地表蒸發(fā)量)(m/d)；Ω為滲流區(qū)域。

2.3.2 土壤熱運移方程

采用熱傳導(dǎo)和對流方程來表征土壤熱運移過程，忽略氣態(tài)水?dāng)U散，僅考慮液態(tài)水的運動對土壤熱量傳輸?shù)挠绊?，二維均質(zhì)多孔介質(zhì)變飽和土壤的熱運移偏微分方程及其定解條件可表示如下：

(2)

2.3.3 土壤水力參數(shù)

土壤含水率θ、土壤有效飽和度Se、非飽和土壤水力傳導(dǎo)系數(shù)K這些參數(shù)均隨著壓力水頭h的變化而發(fā)生改變，需要通過數(shù)學(xué)模型來表征它們的本構(gòu)關(guān)系。本文選用Mualem-van Genuchten公式[24]進行描述，具體公式可表示如下：

(3)

(4)

(5)

式中：θr為土壤殘余體積含水率(m3/m3)；n、m和α為經(jīng)驗擬合參數(shù)，其中m=1-1/n，α為受土壤物理性質(zhì)影響的參數(shù)(1/m)；l為經(jīng)驗擬合參數(shù)，通常取平均值0.5；Se為土壤有效飽和度；Ks為土壤飽和水力傳導(dǎo)系數(shù)(m/d)。

2.4 數(shù)值模型建立與參數(shù)識別

采用地下水?dāng)?shù)值模擬軟件FEFLOW建立二維均質(zhì)多孔介質(zhì)變飽和土壤水熱運移耦合數(shù)值模型。使用advancing front剖分網(wǎng)格，得到總網(wǎng)格數(shù)為3 178、節(jié)點數(shù)為1 683。設(shè)定邊界條件，并設(shè)置土壤水力參數(shù)和土壤介質(zhì)熱特性參數(shù)初值，采用自動時間步長劃分離散時間，使用線性插值結(jié)果作為變飽和土壤非穩(wěn)定滲流場和溫度場的初始值。

數(shù)值模擬采用的土壤水力學(xué)參數(shù)：在距離監(jiān)測井G2和G3正北1 m處，鉆孔采集0 cm、-30 cm、-50 cm、-70 cm、-100 cm、-130 cm、-150 cm、-170 cm、-200 cm深度處的土壤樣品，采用中國地質(zhì)大學(xué)(武漢)材料與化學(xué)分析測試中心的Mastersizer 3000E激光衍射粒度分析儀測定土壤的粒度組成，結(jié)果表明土壤的粒度組成主要為粉粒(51.8%)和砂粒(43.84%)，并參考美國農(nóng)業(yè)部土壤質(zhì)地分類方法，將其定名為粉質(zhì)壤土；環(huán)刀采集左右岸0 cm、-50 cm、-100 cm深度處土壤樣品，測定土壤的干容重為1 650 kg/m3，利用Hydrus1D中內(nèi)嵌的Rosetta Lite V1.1模塊，使用土壤粒度組成和干容重數(shù)據(jù)確定初始土壤水力參數(shù)。

本次選用為期12 d(2015年7月27日22∶00至8月8日22∶00)的實測數(shù)據(jù)率定模型參數(shù)，使用水頭和溫度數(shù)據(jù)反演土壤水力參數(shù)，校正的模型參數(shù)見表1。

表1 反演的土壤水力參數(shù)Table 1 Inversed soil hydraulic parameters

數(shù)值模擬采用的土壤介質(zhì)熱特性參數(shù)有：干土的比熱容Cs、干土的熱導(dǎo)率λs、水的比熱容Cw、水的熱導(dǎo)率λw、土壤縱向熱擴散度DL、土壤橫向熱擴散度DT。利用Lu等[25]研究模型中的常溫土壤熱導(dǎo)率公式計算干土熱導(dǎo)率；參考Bridger等[26]和龍燾等[27]的研究成果設(shè)定干土比熱容的數(shù)值；參考張盼盼[28]的研究數(shù)據(jù)設(shè)定土壤熱運移的縱向擴散度和橫向擴散度，土壤介質(zhì)的熱特性參數(shù)見表2。

表2 土壤介質(zhì)的熱特性參數(shù)Table 2 Soil thermal parameters

3 結(jié)果與討論

保持率定后的模型參數(shù)不變，使用為期37 d(2015年7月27日22∶00至9月2日22∶00)的野外監(jiān)測數(shù)據(jù)驗證模型的可靠性。本文利用均方根誤差RMSE和決定系數(shù)R2來評價模擬值與實測值的吻合程度，其計算公式如下：

(6)

(7)

不同埋深處河床沉積物溫度模擬值與實測值的對比結(jié)果，見圖4。

由圖4可見，30 cm、50 cm、70 cm、100 cm埋深處河床沉積物溫度模擬值與實測值的差異較小，總體趨勢一致，且隨著埋深的增加，溫度波動的振幅減?。粶囟饶M值與實測值的總均方根誤差RMSE為0.220，決定系數(shù)R2均大于0.95，認為河床沉積物溫度梯度的模擬結(jié)果合理可靠。

圖4 不同埋深處河床沉積物溫度模擬值與實測值的對比Fig.4 Comparisons between simulated and observed bed sediment temperature at different depths

河岸監(jiān)測井G2、G3水頭模擬值與實測值的對比結(jié)果，見圖5。

圖5 河岸監(jiān)測井G2、G3水頭模擬值與實測值的對比Fig.5 Comparisons between simulated and observed hydraulic head for monitoring wells G2 and G3

由圖5可見，盡管河岸監(jiān)測井G2、G3處水頭的模擬值與實測值在某些觀測點上存在一定的誤差，但總體趨勢基本一致；水頭模擬值與實測值的總均方根誤差RMSE為0.033，決定系數(shù)R2均大于0.98，認為滲流場的模擬結(jié)果合理可靠，部分誤差可能是由實際渠道斷面土壤的空間變異性造成的。

河岸監(jiān)測井G2、G3溫度模擬值與實測值的對比結(jié)果，見圖6。

由圖6可見，河岸監(jiān)測井G2、G3溫度模擬值與實測值的差異較小，總體趨勢一致；總均方根誤差RMSE為0.115，決定系數(shù)R2均大于0.98，認為溫度場的模擬結(jié)果合理可靠。

綜上分析可見，研究區(qū)地溫和水溫的波動大，有明顯的晝夜變化。渠道地表水與地下水相互作用時，渠道地表水的溫度信號影響著河床沉積物的溫度分布。本文針對不同埋深處河床沉積物的溫度數(shù)據(jù)，提取其日周期波動的正弦函數(shù)擬合振幅值見圖7。

圖6 河岸監(jiān)測井G2、G3溫度模擬值與實測值的對比Fig.6 Comparison of simulated and monitored temperatures for monitoring wells G2 and G3

圖7 不同埋深處河床沉積物溫度曲線正弦函數(shù)擬合振幅Fig.7 Sinusoidal fitting amplitudes of bed sediment temperature curves at different depths

通過分析圖4和圖7可見，河床沉積物的平均溫度從上至下逐漸降低，埋深小于30 cm處河床沉積物的溫度受到渠道水溫波動的影響較大，在埋深為50 cm及以下處河床沉積物的溫度受渠道水溫波動的影響不明顯；當(dāng)渠道地表水補給地下水時，渠道水溫日周期波動信號向下傳遞，河床沉積物溫度日周期波動的振幅(正弦函數(shù)擬合)隨埋深增加而減??；當(dāng)?shù)叵滤a給渠道地表水時，低溫地下水向上進行熱量交換，緩沖了由渠道地表水引起的河床沉積物溫度波動，50 cm埋深處河床沉積物的溫度曲線信號隨之變化，其日周期正弦波動的振幅減小?？梢?，通過分析50 cm埋深以上處河床沉積物的溫度信號，即使在無水位數(shù)據(jù)的情況下，也能分析得出渠道地表水與地下水的相互作用關(guān)系。

為了研究渠道地表水與地下水的相互作用對渠道垂直剖面溫度分布規(guī)律的影響，本文利用上述模擬結(jié)果，選取2015年8月12日22∶00、8月24日22∶00、8月25日22∶00、8月29日22∶00，即模擬時長分別為t=16 d、t=28 d、t=29 d和t=33 d 這4個不同時刻的溫度場模擬結(jié)果，對比分析渠道地表水與地下水不同轉(zhuǎn)化關(guān)系時，渠道研究斷面的溫度分布規(guī)律，見圖8。

圖8 不同時刻渠道研究斷面的溫度分布Fig.8 Distribution of temperature at the research cross section of the cannal at different times

由圖8可見，渠道研究斷面的溫度分布由下至上大致可以分為低溫區(qū)(10～17℃)、中溫區(qū)(17～20℃)和高溫區(qū)(20～23℃)。其中，低溫區(qū)分布較為穩(wěn)定，無論斷面是渠道地表水補給地下水還是地下水補給渠道地表水，低溫區(qū)的面積變化不大；高溫區(qū)分布于受氣溫、水溫影響較大的淺部非飽和層，其分布范圍較小，面積變化差異大；中溫區(qū)位于兩者之間，受到渠道地表水與地下水相互作用中熱傳導(dǎo)作用的影響。

當(dāng)渠道地表水補給地下水時，受溫度日變異性大的渠道地表水入滲的影響，中溫區(qū)的溫度梯度整體隨埋深的增加而逐漸增大，但局部的溫度梯度方向可能相反；當(dāng)?shù)叵滤a給渠道地表水時，中溫區(qū)的溫度梯度方向一致，隨著埋深的增加中溫區(qū)的溫度梯度逐漸減小。這在一定程度上表征了渠道研究斷面渠道地表水與地下水的相互作用關(guān)系。

為了定量評價渠道地表水與地下水的轉(zhuǎn)化水量，本文選取渠道邊界單元，繪制出渠道研究斷面單位長度累計入滲量的變化曲線，見圖9。

圖9 渠道研究斷面單位長度累計入滲量的變化曲線Fig.9 Cumulative recharge per unit length of the research cross section of the canal

由圖9可見，0～28.6 d內(nèi)渠道研究斷面單位長度累計入滲量隨時間增大，渠道地表水補給地下水，隨著渠道地表水水位的降低，渠道研究斷面單位長度累計入滲量曲線斜率減??；28.6～37 d內(nèi)渠道研究斷面單位長度累計入滲量隨時間減小，渠道地表水水位低于地下水水位，轉(zhuǎn)化為地下水向渠道排泄，隨著地下水水位的降低，渠道單位長度累計入滲量曲線斜率的絕對值減小。0～28.6 d內(nèi)渠道單位長度累計入滲量為35.3 m2/d；28.6～37 d內(nèi)渠道研究斷面單位長度累計排泄量為9.626 m2/d；37 d內(nèi)渠道單位長度累計入滲量為25.674 m2/d。

假設(shè)渠道研究斷面為典型斷面，可利用渠道研究斷面單位長度累計入滲量估算整個渠道地表水與地下水的轉(zhuǎn)化量。已知研究區(qū)斷面位于新疆石河子農(nóng)八師瑪納斯河灌區(qū)頭道溝渠道的排水干渠上，干渠全長為12.6 km，南北走向，垂直等高線分布,控制灌溉面積為2萬hm2。通過分析得出:2015年7月27日22∶00至2015年8月25日12∶30，頭道溝干渠渠道地表水入滲補給地下水，地表水入滲補給水量約為4.448×105m3，2015年8月25日12∶30至9月2日22∶00，頭道溝干渠排泄地下水，地下水排泄水量約為1.213×105m3。

4 結(jié) 論

本文以新疆瑪納斯河流域炮臺鎮(zhèn)附近頭道溝渠道的排水干渠為研究對象，通過野外原位試驗實時監(jiān)測渠道橫斷面地下水溫度、地下水水位、渠道地表水水位和不同埋深處河流沉積物的溫度變化，分析了渠道地表水與地下水相互轉(zhuǎn)化關(guān)系變化前后溫度的規(guī)律特征，并建立二維均質(zhì)多孔介質(zhì)變飽和土壤水熱運移數(shù)值模型，模擬分析與計算了頭道溝干渠渠道地表水與地下水的相互轉(zhuǎn)化關(guān)系和轉(zhuǎn)化量，并得到如下結(jié)論：

(1) 淺層河床沉積物(埋深為50 cm以上)的溫度曲線可示蹤表征渠道地表水與地下水的相互轉(zhuǎn)化關(guān)系。在觀測計算時段內(nèi)，2015年7月27日22∶00至2015年8月25日12∶30，渠道地表水補給地下水，2015年8月25日12∶30至9月2日22∶00，地下水向渠道排泄。河床沉積物可大致分為低、中、高三個溫區(qū)，當(dāng)渠道水補給地下水時，受渠道地表水入滲的影響，中溫區(qū)溫度梯度隨埋深的增加而逐漸增大，可能出現(xiàn)淺層沉積物溫度低于深層沉積物溫度的現(xiàn)象；當(dāng)?shù)叵滤蚯琅判箷r，中溫區(qū)溫度梯度方向一致，隨著埋深的增加溫度梯度逐漸減小。河床沉積物的溫度日周期波動振幅(正弦函數(shù)擬合)隨著埋深的增加而減小，當(dāng)渠道地表水補給地下水時，50 cm埋深處河床沉積物的溫度日周期波動振幅隨時間增大；當(dāng)?shù)叵滤蚯琅判箷r，50 cm埋深處河床沉積物的溫度日周期波動振幅隨時間減小。

(2) 使用土壤水熱運移耦合模型反演得到試驗區(qū)土壤的飽和水力傳導(dǎo)系數(shù)為7.05 m/d。使用監(jiān)測數(shù)據(jù)驗證模擬結(jié)果，溫度模擬值與實測值的RMSE為0.220，決定系數(shù)R2均大于0.95，水頭模擬值與實測值的RMSE為0.033，決定系數(shù)R2均大于0.98，模擬值和實測值之間具有較好的一致性，表明率定后的土壤水熱運移模型能較好地仿真渠道地表水與地下水相互轉(zhuǎn)化的動態(tài)變化規(guī)律。

(3) 渠道研究斷面單位長度累計入滲量為35.3 m2/d，發(fā)生在2015年8月25日12∶30時，模擬結(jié)束時渠道研究斷面單位長度累計入滲量為9.626 m2/d。2015年7月27日22∶00至2015年8月25日12∶30期間，頭道溝干渠渠道地表水入滲補給地下水，12.6 km長渠道地表水入滲補給量約為4.448×105m3；2015年8月25日12∶30至9月2日22∶00期間，頭道溝干渠排泄地下水，地下水排泄量約為1.213×105m3。