• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    灌排系統(tǒng)水轉(zhuǎn)化模擬模型研究

    2018-02-01 15:15:10李文暉邵東國徐保利胡能杰楊霞
    南水北調(diào)與水利科技 2018年6期
    關(guān)鍵詞:敏感性分析模型

    李文暉 邵東國 徐保利 胡能杰 楊霞

    摘要:受自然與人為因素綜合影響,灌區(qū)水轉(zhuǎn)化過程復(fù)雜多變。以湖北漳河水庫灌區(qū)為例,分別從水平和垂直方向描述典型灌排系統(tǒng)的水轉(zhuǎn)化過程,以水量平衡原理為基礎(chǔ),研發(fā)灌排系統(tǒng)邊界處置、稻作區(qū)田間多層土壤、塘壩及排水溝道之間多個水轉(zhuǎn)化模擬子模塊,由此構(gòu)建一種客觀描述灌排系統(tǒng)邊界半封閉結(jié)構(gòu)及其水平輸送、垂直運動與渠-田-溝-塘水轉(zhuǎn)化過程的耦合模擬模型,解決了應(yīng)用流域分布式水文模型求解灌區(qū)水轉(zhuǎn)化過程中面臨的灌區(qū)邊界不確定、分層土壤簡化處理等問題。利用漳河灌區(qū)譚店灌排系統(tǒng)2014-2015年水平衡測試觀測數(shù)據(jù)對該模型進行檢驗,采用拉丁超立方抽樣法對模型參數(shù)進行分層抽樣,并利用偏相關(guān)法分析參數(shù)敏感性。通過模型計算可得,支溝率定期和驗證期模擬結(jié)果的復(fù)相關(guān)系數(shù)分別為0.83、0.70,Ens系數(shù)分別為0.81、0.68,斗溝率定期和驗證期模擬結(jié)果的復(fù)相關(guān)系數(shù)分別為0.79、0.68,Ens系數(shù)分別為0.73、0.62,結(jié)果表明該模型適用于灌排系統(tǒng)水轉(zhuǎn)化過程模擬。與SWAT模型對比,對于同一研究區(qū)域,復(fù)相關(guān)系數(shù)由0.11提高到0.73,Ens系數(shù)由-0.71提高到0.70,進一步驗證了該模型的先進性,為灌排系統(tǒng)尺度水轉(zhuǎn)化模擬和機制研究提供了新方法。

    關(guān)鍵詞:灌排系統(tǒng);水轉(zhuǎn)化;模型;拉丁超立方;偏相關(guān);敏感性分析;SWAT模型

    中圖分類號:S27文獻標(biāo)志碼:A

    開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):李文暉

    Study of water transformation simulation model of irrigation and drainage system

    LI Wenhui,SHAO Dongguo,XU Baoli,HU Nengjie,YANG Xia

    (State Key Laboratory of Water Resources and Hydropower Engineering Science,Wuhan University,Wuhan 430072,China)

    Abstract:The process of water transformation in irrigation areas is complicated due to the influence of natural and artificial factors.Taking the irrigation area of Zhanghe Reservoir in Hubei Province as a case study,we described the water transformation process of a typical irrigation and drainage system from both vertical and horizontal perspectives.Based on the principle of water balance,we established a simulation model that was composed of multiple modules for boundary treatment and simulation of the water transformation between multiple soil layers of the rice field,ponds and ditches.The model was calibrated and verified based on the data of Tandian irrigation and drainage system in Zhanghe district from 2014 to 2015.After sampling the model parameters using Latin hypercube,we conducted sensitivity analysis of the parameters using partial correlation analysis.The model calculation showed that the multiple correlation coefficients of the simulation results for the branch ditches in calibration and validation periods were 0.83 and 0.70 respectively,and the Ens coefficients were 0.81 and 0.68 respectively;the multiple correlation coefficients of the simulation results for the tertiary ditches in calibration and validation periods were 0.79 and 0.68,and the Ens coefficients were 0.73 and 0.62,indicating that the model was suitable for simulating the water transformation process of irrigation and drainage systems.As compared with the SWAT model,for the same study area,the proposed model improved the multiple correlation coefficient from 0.11 to 0.73 and the Ens coefficient from -0.71 to 0.70,proving the proposed model was superior to the SWAT model.This study provided valuable information for the study of water transformation in irrigation and drainage systems.

    Key words:irrigation and drainage system;water transformation;model;Latin hypercube;partial correlation;sensitivity analysis;SWAT model

    受地形地貌、水文地質(zhì)、土壤作物、農(nóng)業(yè)水利工程等因素綜合影響,灌排系統(tǒng)結(jié)構(gòu)復(fù)雜多樣,既存在水平灌溉渠系、排水溝系及其與農(nóng)田、塘壩等之間的水量輸送,也有灌溉水量在田間作物、土壤乃至地下水等之間的垂直運動。為揭示農(nóng)田到灌區(qū)水轉(zhuǎn)化機制,國內(nèi)外開展大量的灌排試驗研究,提出了許多灌區(qū)水文過程模擬模型與水平衡機制[1-6],這些都對揭示灌區(qū)水轉(zhuǎn)化過程,促進農(nóng)業(yè)水資源高效利用,保障糧食安全具有重要意義[7-9]。

    現(xiàn)有描述灌區(qū)水文過程的模型大多依賴于國外流域分布式水文模型,根據(jù)灌區(qū)實際條件進行適當(dāng)?shù)母倪M,如代俊峰等[10-11]結(jié)合中國南方丘陵水稻灌區(qū)水文特點,改進SWAT模型中的稻田水分循環(huán)模塊,改變其陸地水文過程計算結(jié)構(gòu),增加地下水補給作用和塘壩灌溉模塊等,構(gòu)建了適應(yīng)于水稻灌區(qū)水分循環(huán)模擬的灌區(qū)分布式水文模型。也有部分學(xué)者考慮到灌溉渠道、排水溝和人工河道等人為干擾,在溝渠河網(wǎng)的提取方法、子流域劃分等方面進行了改進[12-16],但受灌排系統(tǒng)分布等影響,灌區(qū)并不像流域一樣具有明顯的分水嶺,邊界的刻畫一直是灌區(qū)水文模型研究中的一大難題,至今沒有解決。

    此外,多層土壤間的水分運動也是灌區(qū)水文過程的重要組成部分。為揭示不同土層間的水分運動規(guī)律,國內(nèi)外做了大量試驗及模擬研究[17],如甘永德等[18]為探明降雨條件下分層土壤入滲特性,通過室內(nèi)模擬降雨入滲試驗,研究了不同分層組合方式對土壤入滲特性的影響。Yangwen JIA等[19]為揭示非穩(wěn)定降雨條件下多層土壤間的入滲規(guī)律,將Green-Ampt模型推廣到多層土壤的入滲模擬,并與基于Richards方程構(gòu)建的數(shù)值模型模擬結(jié)果進行對比,兩者模擬結(jié)果較為一致,進一步驗證了所構(gòu)建的廣義Green-Ampt模型的適用性。但如何考慮多層土壤間的水分運動及其對灌區(qū)水轉(zhuǎn)化過程的影響,尚有待進一步探究。

    考慮到灌區(qū)是由多個結(jié)構(gòu)基本相似的灌排系統(tǒng)組成的復(fù)雜系統(tǒng),在自然地理條件基本相似灌區(qū),不同灌排系統(tǒng)之間也存在水轉(zhuǎn)化過程的相似性。本文以灌排系統(tǒng)為研究尺度,重點探討灌排系統(tǒng)復(fù)雜邊界條件及多層土壤對灌區(qū)水轉(zhuǎn)化過程影響問題,研發(fā)灌排系統(tǒng)邊界處置、稻作區(qū)田間多層土壤、塘壩及排水溝道之間多個水轉(zhuǎn)化模塊,構(gòu)建了一種新的灌排系統(tǒng)水轉(zhuǎn)化(Water Transformation of Irrigation and Drainage System,以下簡稱WTIDS)模型,為客觀、精細(xì)化模擬灌區(qū)水轉(zhuǎn)化過程提供了新方法。

    1WTIDS模型構(gòu)建

    1.1灌排系統(tǒng)水轉(zhuǎn)化模型概化

    灌排系統(tǒng)一般由水源、水源工程、灌溉渠系、排水溝系、各類田間工程、渠系建筑物以及容泄區(qū)組成,其平面概化示意圖見圖1。A、B為灌溉渠系,C為排水溝系。水經(jīng)各級渠道進入田間,多余水量排入溝道,部分塘壩與渠道溝道相連,發(fā)揮調(diào)蓄作用。

    垂向上,灌排系統(tǒng)水轉(zhuǎn)化過程分層從上至下可分為大氣層、作物層、土壤層以及地下水層,具體如圖2所示。水通過降雨至田間,同時蒸發(fā)散失至大氣層。水分在各層土壤間運動,部分流出底層土壤剖面補給地下水,地下水也會以毛管上升水形式補充土壤水。

    灌排系統(tǒng)區(qū)別于自然流域十分重要的一點,即自然流域擁有明確且閉合的分水嶺脊線作為流域邊界,邊界上無水量交換過程,而灌排系統(tǒng)的邊界通常為渠系溝道、堤防或是公路,由于邊界上分布有渠道或是溝道岔口,灌排系統(tǒng)邊界上存在水量輸入點。灌排系統(tǒng)邊界輸入有三種形式:(1)典型灌溉渠系按“干—支—斗—農(nóng)—毛”五級渠道劃分,從水源引水進入干渠,經(jīng)過各級渠道最終通過毛渠將水輸入田間,本模型將斗渠進水口作為控制點,編號為i(i=1,2,3,…,n1),記流量為qAi;(2)實際生活中灌溉渠系和排水溝系的布置會根據(jù)區(qū)域的地形地貌進行調(diào)整,存在直接從干渠引水進行田間灌溉的岔口,編號為j(j=1,2,3,…,n2),記流量為qj;(3)一般灌排系統(tǒng)為一套灌溉渠系和排水溝系相互配合布置,但存在兩套灌溉渠系共用一套排水溝系的情況,如圖1所示,B灌溉渠系斗渠進水口控制點編號為k(k[WTBX]=1,2,3,…,n3),記流量為qBk[WTBX]。則整個灌排系統(tǒng)的灌溉水量為

    W=∑[DD(]n1[]i=1[DD)]qAi+∑[DD(]n2[]j=1[DD)]qj+∑[DD(]n3[]k[WTBX]=1[DD)]qBk[WTBX](1)

    田間排水通過農(nóng)溝匯入斗溝,進而匯入支溝、干溝排入容泄區(qū),以斗溝與支溝的交匯口為控制點,編號為m(m=1,2,3,…,n),記流量為Qm,并在兩控制點間采用水量平衡原理,可計算出支溝排水口的徑流量。

    則灌排系統(tǒng)水平衡方程為

    SWt=SWt-1+Pt-Dt-Et-DPt+ΔWpond,t(2)

    式中:SWt、SWt-1分別表示第t天和第t-1天的土壤儲水量(mm);t表示時間(d);Pt表示第t天的降雨量(mm);Wt表示第t天的灌溉水量(mm);Dt表示第t天的排水量(mm);Et表示第t天的蒸散發(fā)量(mm);DPt表示第t天從土壤剖面底部滲出的水量(mm);ΔWpond,t表示第t天灌排系統(tǒng)內(nèi)塘壩蓄水量變化量(mm)。部分地表徑流進入排水溝道流出灌排系統(tǒng),另一部分地表徑流則進入附近塘壩,積蓄于塘壩內(nèi)。

    灌排系統(tǒng)水文過程受道路、堤防或渠溝分布等人為因素影響巨大,在劃定灌排系統(tǒng)邊界時須考慮這些人為因素的影響。灌排系統(tǒng)中,襯砌的渠道起到了隔絕水流傳輸?shù)淖饔?,公路、鐵路以及堤防在設(shè)計及建設(shè)過程中也會將路基壓實墊高,這與自然的分水嶺功能類似,因此,可以將公路、鐵路、堤防或是襯砌后的渠系作為定義灌排系統(tǒng)邊界范圍的依據(jù)。子流域劃分同樣以此為原則,以襯砌的規(guī)模較小的渠系,如支渠或斗渠,結(jié)合公路、鐵路及堤防的分布作為子流域邊界劃分的依據(jù),并將子流域作為模型水文計算的最小單元。

    為了滿足實際灌溉需求,干渠或分干渠沿線往往設(shè)置諸多分水岔口,直接引水灌溉農(nóng)田,本模型考慮邊界上水量點輸入對灌排系統(tǒng)水轉(zhuǎn)化過程的影響,在定義子流域的基礎(chǔ)上,結(jié)合實地觀測數(shù)據(jù),將某個子流域邊界上的若干輸入點的流量進行累加,使之等效為一個輸入點,流量記為Wi(i為子流域編號)。

    模型分別模擬估算各子流域的徑流量,然后通過排水溝道演算過程得到整個灌排系統(tǒng)的總徑流量。

    1.2灌排系統(tǒng)水轉(zhuǎn)化過程描述

    1.2.1田間水轉(zhuǎn)化

    本模型定義田間儲水深度,即土壤耕作層儲水量與田間淹沒水層的水量之和,并以田間儲水深度作為狀態(tài)變量,根據(jù)水量平衡方程,建立其與各輸入輸出變量間的田間水平衡方程,方程表示為:

    St=St-1+Pt+IRt+CRt-ETplant,1,t-ETzj,t-DPt(3)

    式中:St-1、St分別為第t-1天和第t天的田間儲水深度;Pt為第t天降雨量;IRt為第t天灌溉量;CRt為第t天補充耕作層的毛管水上升量;ETplant,1,t為第t天耕作層作物根系吸水量;ETzj,t為第t天株間蒸發(fā)量;DPt為第t天滲漏至犁底層的水量。以上各變量單位均為mm。

    (1)田間排水量。

    DRt=[JB({]St-θsat,1·hsoil,1[KG-*6]-[KG-*6]Hp,k[WTBX][KG*2]St-θsat,1·hsoil,1≥Hp,k[WTBX]

    0St-θsat,1·hsoil,1<Hp,k[WTBX][JB)](4)

    式中:St為第t天田間儲水深度(mm);θsat,1為耕作層土壤體積飽和含水率;hsoil,1為耕作層厚度(mm);Hp,k為水稻第k個生育階段對應(yīng)的田間最大蓄水深度(mm)。對于水田,僅當(dāng)田間淹沒水層超過最大蓄水深度時,才會通過排水孔進行排水,假設(shè)水田排水全部流入排水溝。

    (2)植株蒸騰量與株間蒸發(fā)量。

    植株蒸騰量與株間蒸發(fā)量合稱作物需水量,本模型采用國內(nèi)常用的作物系數(shù)法計算作物實際需水量[20],計算方法為

    ETt=kc·ET0t=ETplant,t+ETzj,t(5)

    式中:ETt為第t天的作物實際需水量(mm);kc為作物系數(shù),無量綱;ET0t為第t天參考作物需水量(mm);ETplant,t為第t天水稻植株蒸騰水量(mm);ETzj,t為第t天水稻株間蒸發(fā)水量。

    植株蒸騰的水分來源于作物根系吸水,由于水稻根系層深度一般大于耕作層厚度(20 cm),因此水稻蒸騰水量須在不同土層間進行分配,根據(jù)相關(guān)研究[20],水稻根系吸水總量ETplant可以表示為

    ETplant=∫.Lr0ETplant([SX(]1.8[]Lr[SX)]-[SX(]1.6[]Lr.2[SX)]z)dz(6)

    式中:z為根系深度(m);Lr為水稻有效根層的深度(m)。水稻為淺根系作物,一般90%以上的根系分布于地表以下30 cm以內(nèi)的土層中,取Lr=0.3 m。則根據(jù)式(6),可計算得到耕作層根系吸水量ETplant,1為38/45ETplant,犁底層根系吸水量ETplant,2為7/45ETplant。

    一般水稻生育初期植株小,地面裸露大,水分消耗以株間蒸發(fā)為主;隨著植株生長,葉面覆蓋率增大,植株蒸騰逐漸大于株間蒸發(fā)。本模型以式(7)來模擬水稻每日的株間蒸發(fā)量。

    ETzj,t=rzj·ET0t(7)

    式中:rzj為比例系數(shù),根據(jù)相關(guān)文獻[21],分蘗期取0.35,孕穗期、抽穗期取0.15,乳熟期取0.32,黃熟期取0.45,返青期由于植株較小,將株間蒸發(fā)量與當(dāng)日水面蒸發(fā)量近似相等處理。

    (3)灌溉需水量。

    DIt=[JB(]hmaxt-hmintht≤hmint

    0ht>hmint[JB)](8)

    式中:ht為田間淹沒水層深度(mm);hmint為田間適宜水層下限(mm);hmaxt為田間適宜水層上限(mm)。

    (4)毛管上升水量。

    在地下水與土壤底面的交界處,因毛管力作用,地下水將以毛管上升水的形式對土壤水進行補充,本模型在前人經(jīng)驗公式基礎(chǔ)上加以改進,以此來模擬毛管上升水量。具體如下式所示:

    CRt=ETt·exp(-b·d)+su(9)

    式中:b為反映土壤輸水能力的常數(shù),經(jīng)驗系數(shù),對于砂土、壤土和黏土可分別取2.1、2.0和1.9[22];d為地下水埋深(m);su為常數(shù),可根據(jù)排水口實測徑流過程反演得到。

    (5)非水田徑流模擬。

    對于旱地、城鎮(zhèn)用地以及林草地等非水田土地利用類型,本模型采用SCS曲線數(shù)法[23,24]進行徑流量模擬計算,不同土地利用類型的曲線數(shù)CN根據(jù)相關(guān)文獻選定初始值[10,24],見表1。

    1.2.2多層土壤間的水分運動

    各層土壤的物理性質(zhì),如飽和體積含水率、田間持水量、飽和水力傳導(dǎo)度等不同,導(dǎo)致各個土層內(nèi)的水平衡過程及不同土層間的水轉(zhuǎn)化關(guān)系存在差異。本模型根據(jù)湖北荊門試驗研究[25],將土壤分為三層,包括耕作層(地表以下0~20 cm)、犁底層(地表以下20~35 cm)及底土層(地表下35 cm及以下),分別模擬計算各層土壤的水平衡過程,并假設(shè)水分在土層中均勻分布。各層土壤的物理參數(shù)均由實地試驗測得,以保證模型計算結(jié)果的準(zhǔn)確性。則各層土壤含水量為

    θi=[SX(]Ssoil,i[]hsoil,i[SX)](10)

    式中:θi為某天第i層土壤的含水量,%;Ssoil,i為當(dāng)天第i層土壤的儲水深度(mm);hsoil,i為第i層土壤厚度(mm)。令i = 1表示耕作層,i=2表示犁底層,i=3表示底土層。

    (1)耕作層水平衡。

    本模型將耕作層土壤與田間淹沒水層綜合考慮,引入田間儲水量對其水平衡過程進行描述,具體如式(3)所示。

    (2)犁底層水平衡。

    犁底層接受耕作層滲漏水量補充,同時以作物根系吸水、側(cè)滲和向底土層滲漏形式損失水量,建立其水量平衡方程,具體表示為

    Ssoil,2,t=Ssoil,2,t-1+DP1,t-ETplant,2,t-DP2,t-DPhf,2,t(11)

    式中:Ssoil,2,t、Ssoil,2,t-1分別為第t天、第t-1天犁底層土壤的儲水量;DP1,t為第t天耕作層土壤垂直入滲量;ETplant,2,t為第t天犁底層作物根系吸水量;DP2,t為第t天犁底層垂直入滲至底土層的水量;DPhf,2,t為第t天犁底層側(cè)滲水量。以上變量的單位均為mm。

    犁底層土質(zhì)黏重,特別是中部和下部土層對水分下滲有較強的抑制作用,減少深層下滲,起到保水作用。

    (3)底土層水平衡。

    由于犁底層的保水作用,底土層接受犁底層滲漏補給的水量很少,根據(jù)前期研究成果表明[25],底土層土壤含水率基本保持不變,本模型默認(rèn)底土層土壤含水率為該層田間持水量。然后考慮底土層的側(cè)滲和向淺層地下水層的滲漏,可建立底土層水平衡方程,具體為

    DP2,t=DP3,t+DPhf,3,t(12)

    式中:DP3,t為第t天底土層垂直入滲至地下水層的水量(mm);DPhf,3,t為第t天底土層側(cè)滲水量(mm)。

    (4)滲漏過程。

    在計算各土壤層的滲漏量時,如果某層的含水量超過其田間持水量,且下層未飽和時,水分將會滲漏。假設(shè)當(dāng)土壤飽和時,滲漏量為定值,當(dāng)土壤含水率在田間持水量和飽和含水率之間時,滲漏量與土壤含水率成線性關(guān)系[26]。具體表示為

    DPable,i=[JB({]DP0,iθi-θsat,i

    [SX(]θi-θfc,i[]θsat,i-θfc,i[SX)]·[KG-*2]DP0,i[KG*2]θfc,i<θi<θsat,i[JB)](13)

    式中:DP0,i為第i層土壤飽和時的日滲漏量,通過試驗確定(mm);θi為某天該層土壤體積含水率,θfc,i為該層土壤田間持水量。

    下層土壤可接受的垂直入滲量為下層土壤達到飽和時的儲水量與實際儲水量之差,具體表示為

    DPacceptable,i=(θsat,i-θi)·hsoil,i(14)

    式中:DPacceptable,i為第i層土壤某天可接受上層土壤的垂直滲漏量(mm)。

    則某天第i層土壤的垂直滲漏量為

    DPi=MIN(DOable,i,DPacceptable,i+1)(15)

    上層土壤滲漏至下層土壤的水量如果超過了下層土壤可接受量,則超過部分以側(cè)滲形式流出該層土壤,則某天第i層土壤側(cè)滲水量為

    DPhf,i,t=r·Ssoil,i,t+DPable,i,t-DPi,t(16)

    式中:r為側(cè)滲系數(shù),本模型取0.01[27]。

    1.2.3淺層地下水轉(zhuǎn)化

    淺層地下水層的水量平衡方程為

    aqsh,t=aqsh,t-1+w[WTBX]supply,sh,t-qgw[WTBX],t-CRt(17)

    式中:aqsh,t為第t天淺層地下水層的儲水量;aqsh,t-1為第t-1天淺層地下水層的儲水量;w[WTBX]supply,sh,t為第t天淺層地下水層的補給量;qgw[WTBX],t為第t天匯入排水溝的地下水徑流量;CRt為第t天因土壤水分不足而進入土壤帶的水分。以上變量的單位均為mm。

    (1)補給。

    水分滲漏通過土壤最底層進入地下水層。水分流出土壤剖面與地下水層補給過程中的時間延遲,本模型應(yīng)用指數(shù)衰減權(quán)重函數(shù)來考慮[24]。第t天地下水層的總補給量計算如下:

    w[WTBX]supply,t=[JB([]1-exp[JB((]-[SX(]1[]δgw[WTBX][SX)][JB))][JB)]]DP3,t+exp[JB((]-[SX(]1[]δgw[WTBX][SX)][JB))]wsupply,t-1(18)

    式中:wsupply,t為第t天地下水層的總補給量(mm);δgw[WTBX]為延遲時間(d);DP3,t為第t天底土層垂直入滲至地下水層的水量(mm);wsupply,t-1為第t-1天地下水層的總補給量(mm)。

    地下水層的總補給量進入淺層地下水層,而后其中一部分會進入深層地下水層,則某天淺層地下水層的補給量為

    w[WTBX]Lt=(1-αdeep)w[WTBX]supply,t(19)

    式中:αdeep為進入深層地下水層的水量比例系數(shù)。

    進入到深層地下含水層的水量不參與水量收支計算,可以看作損失部分。

    (2)地下水徑流。

    淺層地下水可以補給子流域內(nèi)的排水溝道,本模型中假設(shè)僅當(dāng)淺層地下水層的儲水量超過設(shè)定的水位閾值aqshq時,進行補給,補給量為當(dāng)天淺層地下水層儲水量與水位閾值間差值。即

    qgw[WTBX],t=aqsh,t-aqshq(20)

    式中:qgw[WTBX],t為第t天淺層地下水層補給排水溝的水量(mm);aqsh,t為第t天淺層地下水層的儲水量(mm);aqshq為設(shè)定的淺層地下水層水位閾值(mm)。

    1.2.4塘壩水轉(zhuǎn)化

    塘壩具有積蓄雨水、攔截地表徑流和灌溉回歸水的功能,發(fā)揮著重要的調(diào)蓄作用,灌區(qū)塘壩數(shù)量多、分布廣,除了從渠道引水灌溉農(nóng)田,當(dāng)?shù)卮迕襁€抽取塘壩水灌溉旱地,部分塘壩直接分布在排水溝系上,排水對塘壩進行庫容補充。

    本模型以塘壩蓄水量作為狀態(tài)變量,建立其與各輸入輸出變量間的水平衡方程,沿排水溝水流方向,對流經(jīng)的各子流域依次進行計算,塘壩水平衡方程表示為

    Vt=Vt-1+Wint-WELt-Wirrt(21)

    式中:Vt、Vt-1分別為第t天、第t-1天塘壩蓄水量;Wint為第t天塘壩來水量,包括塘面集雨量、非水田地表徑流量以及上游排水溝充庫水量;WELt為第t天塘壩損失水量,包括塘面蒸發(fā)損失量和塘壩滲漏損失量;Wirrt為第t天塘壩供水量。各變量單位均以m.3表示。

    (1)塘壩入流量。

    Wint=∑[DD(]n[]i=1[DD)]Dit·Afi/1000+Pt·Ap/1000+Wck[WTBX],t(22)

    式中:Dit為第t天第i種非水田土地利用類型地表徑流量(mm);Afi為第i種非水田土地利用類型面積(m.2);Pt為第t天降雨量(mm);Ap為塘壩面積(m.2);Wckt為第t天上游排水溝對塘壩的充庫量(m.3)。

    (2)塘壩損失水量。

    WELt=E0t·Ap/1000+Wlosst((23)

    式中:E0t為第t天水面蒸發(fā)量(mm);Wlosst為第t天塘壩滲漏量(m.3),由實際觀測數(shù)據(jù)求得。

    (3)塘壩供水量。

    根據(jù)實地調(diào)查,研究區(qū)域內(nèi)旱地灌溉用水由塘壩提供,則塘壩供水量為

    Wirrt=DIdt·Ad/1000(24)

    式中:DIdt為第t天旱地灌溉需水量(mm);Ad為旱地面積(m.2)。

    (4)塘壩泄水量。

    當(dāng)Vt-1+Wint-WETt-DIdt>Vm(Vm為塘壩最大蓄水容量)時, 塘壩泄水量Xt=Vt-1+Wint-WETt-DIdt-Vm(Xt單位為m.3)。

    1.2.5排水溝道演算過程

    水田排水以及塘壩泄水排入溝道,向下游運動,溝道中的水分因溝床滲漏而損失,通過直接降落在溝道水面的雨水進行補充。本模型通過建立排水溝道水量平衡方程,模擬計算排水口徑流過程,水量平衡方程表示為

    Qout=QPt+Xt+QDRt-QLt(25)

    式中:Qout為第t天排水口徑流量;QPt為第t天降雨補充溝道水量;Xt為第t天塘壩泄水量;QDRt為第t天水田田間排入溝道的水量;QLt為第t天溝道輸水損失水量。單位均為m.3。

    (1)降雨補充溝道水量。

    QPt=Pt·A/1000(26)

    式中:A為溝道接收降雨面積(m.2)。

    (2)溝道輸水損失水量。

    溝道輸水損失水量用經(jīng)驗公式來估算,常用經(jīng)驗公式為

    qLt=σ·L·qint(27)

    式中:qLt為第t天溝道輸水損失流量(m.3/s);L為溝道長度(km);σ為每公里溝道輸水損失系數(shù);qint為第t天溝道入流量(m.3/s)。

    每公里溝道輸水損失系數(shù)σ亦采用經(jīng)驗公式估算,常用經(jīng)驗公式為

    σ=[SX(]A[]100q.mint[SX)](28)

    式中:A為溝床土壤透水系數(shù),取A=1.9;m為溝床土壤透水指數(shù),取m=0.4[20]。

    用公式計算出來的輸水損失量是在不受地下水頂托影響條件下的損失水量,若研究區(qū)域地下水位較高,渠道滲漏受地下水壅阻影響,實際滲漏水量比計算結(jié)果小,故須給計算結(jié)果乘上一個修正系數(shù)加以修正,即

    q′Lt=γ·q′Lt(29)

    式中:γ為地下水頂托修正系數(shù),取γ =0.82。

    則溝道輸水損失水量為

    QLt=86400q′Lt(30)

    式中:q′Lt為考慮地下水頂托影響后的溝道輸水損失流量(m.3/s);86 400為時間系數(shù)(s)。

    2實例應(yīng)用

    2.1研究區(qū)概況

    譚店灌排系統(tǒng)位于湖北省荊門市掇刀區(qū),112.16°E、30.91°N,屬漳河灌區(qū),面積2.66 km.2,北依漳河灌區(qū)總干渠,東鄰鳳凰水庫,距漳河水庫13.6 km,灌溉水源充足,經(jīng)實地調(diào)研,當(dāng)?shù)卮迕駜?yōu)先引渠道水進行水田灌溉,旱地則引塘壩水進行灌溉。本區(qū)屬長江中下游亞熱帶季風(fēng)氣候類型,區(qū)內(nèi)多年平均氣溫15.8 ℃,多年平均降雨量903.3 mm,多年平均蒸發(fā)量1 413.9 mm。本區(qū)大部分耕地質(zhì)地粘重,透水性較差,保水、保肥、抗旱能力較強,易于種植水稻。

    2.2模型輸入

    2.2.1空間數(shù)據(jù)

    (1)DEM(數(shù)字高程模型)與子流域劃分。

    所用數(shù)字高程數(shù)據(jù)從地理空間數(shù)據(jù)云網(wǎng)站上獲得,為GDEMV2 30M分辨率數(shù)字高程數(shù)據(jù),投影坐標(biāo)系為Beijing_1954_3_Degree_GK_CM_111E,區(qū)域DEM見圖3。根據(jù)研究區(qū)實際情況,人為定義子流域邊界,生成反映真實渠系溝系的分布圖,本灌排系統(tǒng)共劃分為10個子流域,見圖4。

    (2)土地利用類型數(shù)據(jù)。

    研究區(qū)土地利用類型根據(jù)實地調(diào)研成果,借由Arcgis平臺繪制而成,為簡化計算,若某種土地利用類型面積所占研究區(qū)總面積比例小于3%,則忽視該土地利用類型,將其整合到相近土地利用類型,整合后的研究區(qū)土地利用類型分為水稻田、旱地、塘壩、城鎮(zhèn)用地以及林地等5類,詳見表2。

    (3)土壤類型。

    譚店灌排系統(tǒng)土壤耕作層、犁底層以及底土層各項物理參數(shù),由試驗測定,具體見表3。

    2.2.2水文氣象數(shù)據(jù)

    譚店灌排系統(tǒng)有一條支溝和斗溝由北向南穿越全境,還有一條農(nóng)溝自灌排系統(tǒng)中部延伸至南部邊界。支溝和斗溝的進水流量由實測數(shù)據(jù)確定,農(nóng)溝由于無接收上游溝道水量,故模型中進水流量為0。灌排系統(tǒng)邊界上各水量輸入點的流量數(shù)據(jù)根據(jù)2014-2015年暑期在當(dāng)?shù)乇O(jiān)測得到。

    本文所用的氣象資料來源于漳河團林試驗站2014-2015年監(jiān)測的逐日降雨量、最高氣溫、最低氣溫、相對濕度、大氣壓力、相對濕度、風(fēng)速、日照時數(shù)等數(shù)據(jù)。

    2.2.3作物系數(shù)

    譚店灌排系統(tǒng)以水稻作為主要農(nóng)作物,其作物實際需水量采用作物系數(shù)法計算得到,下表為湖北省中稻作物需水系數(shù)kc值。

    2.2.4水稻各生育期淹灌水層深度

    本模型中需要輸入水稻各生育期天數(shù)以及淹灌水層深度,用于計算水文循環(huán)過程各組分量。本文輸入常規(guī)淹灌方式下的水層深度,具體見表5。[HJ1][HJ1.4mm]

    2.3模型結(jié)果分析

    2.3.1參數(shù)敏感性分析

    在進行參數(shù)敏感度分析時,輸入?yún)?shù)較多且關(guān)系復(fù)雜,輸入?yún)?shù)和輸出結(jié)果之間的相關(guān)性往往受到其他輸入?yún)?shù)的影響,因此只是簡單的分析某輸入?yún)?shù)和輸出結(jié)果的相關(guān)性是不客觀的。準(zhǔn)確反映輸入?yún)?shù)和輸出結(jié)果的相關(guān)關(guān)系,須控制其他輸入?yún)?shù)的影響。偏相關(guān)分析就是控制其他參數(shù)的影響,求得輸入變量和輸出結(jié)果相關(guān)性的全局敏感性分析方法。本文采用偏相關(guān)法對模型參數(shù)進行敏感性分析,探求各參數(shù)對排水口徑流過程的影響程度。排除若干由實測或依據(jù)經(jīng)驗得到的參數(shù),以及對排水口徑流過程影響較小的參數(shù)后,選取剩下6個參數(shù)進行分析,根據(jù)研究區(qū)域?qū)嶋H情況確定參數(shù)的取值范圍以及初始值,具體見表6。

    進行敏感性分析之前,需要根據(jù)統(tǒng)計學(xué)原理對參數(shù)進行抽樣,因拉丁超立方抽樣方法能夠使抽樣的點均勻分布在取值范圍內(nèi)[28],故本文采用拉丁超立方抽樣法對選出的6個模型參數(shù)進行分層抽樣,通過Matlab編程隨機生成500組參數(shù)。將500組抽樣參數(shù)替換本模型輸入文件中對應(yīng)的參數(shù),并分別導(dǎo)入到本模型中,得到500組輸出文件。選取在模擬中有重要意義的排水口峰值流量、峰值時間以及總徑流量這三個變量,采用偏相關(guān)分析法分析參數(shù)對輸出結(jié)果的敏感性。

    借助SPSS數(shù)據(jù)分析軟件,采用偏相關(guān)分析法進行非線性分析,得出各參數(shù)與輸出變量間的偏相關(guān)系數(shù),其中偏相關(guān)系絕對值越大表明參數(shù)敏感性越大,絕對值越小則表明敏感性越小;偏相關(guān)系數(shù)為正,表明參數(shù)與輸出變量為正相關(guān)關(guān)系,反之,偏相關(guān)系數(shù)為負(fù)則為負(fù)相關(guān)。分析結(jié)果見表7,直方圖見圖5。

    由表7和圖5可知,對峰值流量而言,參數(shù)su和DP0,1相對較為敏感,偏相關(guān)系數(shù)分別為0.407和-0.178,顯著性概率值為0,結(jié)果可靠。su與峰值流量之間呈正相關(guān),su增大將使峰值流量增大,這是因為毛管上升水量增加,使得田間儲水量增大,排水量加大,從而增大了峰值流量。參數(shù)DP0,1與峰值流量之間呈負(fù)相關(guān),DP0,1增大將使峰值流量減小,這是因為耕作層土壤滲漏量增大,使得田間儲水量減小,排水量減小,從而減小了峰值流量。

    對峰值時間而言,6個參數(shù)與其的偏相關(guān)系數(shù)都較小,敏感性不高,說明以上6個參數(shù)對峰值時間的影響不大。

    對總徑流量而言,參數(shù)su最為敏感,偏相關(guān)系數(shù)為0.453,DP0,1相對較為敏感,偏相關(guān)系數(shù)為0.174,顯著性概率值均為0,結(jié)果可靠。su與總徑流量之間呈正相關(guān)關(guān)系,su增大將使總徑流量增大,這是因為毛管上升水補充耕作層土壤的水量增加,田間儲水量增大,排水量加大,從而增大了總徑流量。參數(shù)DP0,1與總徑流量之間呈負(fù)相關(guān),即DP0,1增大將使總徑流量減小,這是因為耕作層土壤滲漏量增大,使得田間儲水量減小,排水量減小,從而減小了總徑流量。

    2.3.2模型率定與驗證

    在參數(shù)敏感性分析結(jié)果基礎(chǔ)上,利用譚店灌排系統(tǒng)2014年支溝、斗溝、農(nóng)溝排水口的徑流量實際監(jiān)測數(shù)據(jù)對模型進行參數(shù)率定,率定后的主要參數(shù)取值見表8,支溝、斗溝率定期和驗證期排水口徑流過程實測值與模擬值對比效果圖分別見圖6、圖7。率定后,農(nóng)溝排水口模擬徑流過程與實際差異較大,原因可能是農(nóng)溝涉及的區(qū)域面積較小,村民加高田埂高度,抽水灌溉等田間管理措施對區(qū)域水轉(zhuǎn)化過程影響較大,而實際調(diào)查中未能及時發(fā)現(xiàn)記錄,導(dǎo)致模擬效果不佳。

    選用線性回歸系數(shù)R.2和Nash-Suttclife系數(shù)(Ens)來評估模型在率定及驗證過程中的模擬效果,要求排水口徑流量的R.2>0.6且Ens>0.5。率定期和驗證期模型模擬效果評價指標(biāo)結(jié)果見表9,表明模型在研究區(qū)域取得了較好的模擬效果,該模型適用于灌排系統(tǒng)水轉(zhuǎn)化過程模擬。

    3與SWAT模型對比

    SWAT(Soil and Water Assessment Tool)模型是一種分布式水文模型,廣泛應(yīng)用于自然流域或灌區(qū)水文過程模擬研究[29]。為進一步探究證明本文所構(gòu)建的WTIDS模型的先進性及有效性,針對同一研究區(qū)域,[HJ1.89mm]建立基于SWAT的水轉(zhuǎn)化模型,以此模擬研究區(qū)域的水轉(zhuǎn)化過程,并對比分析兩個模型的模擬效果。

    研究區(qū)域為譚店灌排系統(tǒng),SWAT模型采用傳統(tǒng)D8方法,基于數(shù)字高程數(shù)據(jù)劃定區(qū)域邊界,定義水系分布及流向,在平原及丘陵地區(qū),其邊界劃分和水系定義結(jié)果與實際情況往往有較大出入。WTIDS模型則充分考慮道路、堤防或渠溝分布等人為因素對灌排系統(tǒng)水文過程的影響,并將其作為定義灌排系統(tǒng)邊界范圍的依據(jù)。對于譚店灌排系統(tǒng),SWAT模型共劃分子流域23個,WTIDS模型共劃分子流域7個。具體劃分情況見圖8。

    分別利用兩個模型模擬計算出排水口徑流過程,并采用R.2和Ens來對比模型模擬效果,具體模擬結(jié)果及效果見圖9。由圖9可知,SWAT模型模擬評價指標(biāo)R.2為0.11,Ens為負(fù)數(shù)。而WTIDS模型模擬結(jié)果的R.2和Ens分別為0.73和0.70,達到了R.2>0.6且Ens>0.5的要求,說明相對于SWAT模型,WTIDS模型能更好刻畫灌排系統(tǒng)子流域劃分和水系分布情況,能更為真實地反映灌排系統(tǒng)的水轉(zhuǎn)化過程,適用性更加優(yōu)越,具有一定的先進性。

    4結(jié)論

    (1)本文分析了灌排系統(tǒng)水轉(zhuǎn)化過程,以水量平衡為原理,研發(fā)了灌排系統(tǒng)邊界[HJ1.3mm]處置、稻作區(qū)田間多層土壤、塘壩及排水溝道之間多個水轉(zhuǎn)化模擬子模塊,構(gòu)建了一種客觀描述灌排系統(tǒng)邊界半封閉結(jié)構(gòu)及其水平輸送、垂直運動與渠—田—溝—塘水轉(zhuǎn)化過程的耦合模擬模型,并利用研究區(qū)域支溝、斗溝徑流實測數(shù)據(jù)進行模型的率定及驗證,發(fā)現(xiàn)支溝率定期和驗證期模擬結(jié)果的復(fù)相關(guān)系數(shù)分別為0.83、0.70,Ens系數(shù)分別為0.81、0.68,斗溝率定期和驗證期模擬結(jié)果的復(fù)相關(guān)系數(shù)分別為0.79、0.68,Ens系數(shù)分別為0.73、0.62,表明該模型適用于灌排系統(tǒng)水轉(zhuǎn)化過程模擬。

    (2)與SWAT模型模擬結(jié)果進行對比,發(fā)現(xiàn)對于同一研究區(qū)域,模擬結(jié)果的復(fù)相關(guān)系數(shù)由0.11提高到0.73,Ens系數(shù)由-0.71提高到0.70,進一步驗證了該模型的先進性,為灌排系統(tǒng)水轉(zhuǎn)化過程模擬與機制研究提供了新方法。

    (3)應(yīng)用構(gòu)建的WTIDS模型模擬時,所需率定的參數(shù)較少;其他與排水口徑流過程直接或間接相關(guān)的參數(shù),如土壤體積飽和含水率、溝床土壤透水系數(shù)等都由實測或依據(jù)經(jīng)驗得到,易于獲取。經(jīng)過實例應(yīng)用證明,WTIDS模型在南方水稻灌排系統(tǒng)中具有良好的適用性,如何研究分析研究尺度擴大對水轉(zhuǎn)化關(guān)系的影響以及該模型在不同研究區(qū)域的適用性,將是下一步研究的重點。

    參考文獻(References):

    [1]XUE J,REN L.Evaluation of crop water productivity under sprinkler irrigation regime using a distributed agro-hydrological model in an irrigation district of China[J].Agricultural Water Management,2016,176:350-365.

    [2]DE MARIA A S C,RIENZNER M,F(xiàn)ACCHI A,et al.Water balance implications of switching from continuous submergence to flush irrigation in a rice-growing district[J].Agricultural Water Management,2016,171:108-119.

    [3]李穎,王康,周祖昊.基于SWAT模型的東北水稻灌區(qū)水文及面源污染過程模擬[J].農(nóng)業(yè)工程學(xué)報,2014,30(7):42-53.(LI Y,WANG K,ZHOU Z H.Simulation of drainage and agricultural non-point source pollutions transport processes in paddy irrigation district in North-East China using SWAT[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2014,30(7):42-53.(in Chinese))

    [4]馮紹元,蔣靜,霍再林,等.基于SWAP模型的春小麥咸水非充分灌溉制度優(yōu)化[J].農(nóng)業(yè)工程學(xué)報,2014,30(9):66-75.(FENG S Y,JIANG J,HUO Z L,et al.Optimization of irrigation scheduling under deficit irrigation with saline water for spring wheat based on SWAP model[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2014,30(9):66-75.(in Chinese))

    [5]NOORY H,VAN DER ZEE S E A T,LIAGHAT A M,et al.Distributed agro-hydrological modeling with SWAP to improve water and salt management of the Voshmgir Irrigation and Drainage Network in Northern Iran[J].Agricultural Water Management,2011,98:1062-1070.

    [6]張銀輝,羅毅.基于分布式水文學(xué)模型的內(nèi)蒙古河套灌區(qū)水循環(huán)特征研究[J].資源科學(xué),2009,31(5):763-771.(ZHANG Y H,LUO Y.Water cycle modeling in Hetao irrigation district of inner Mongolia based on the distributed DEHYDROS model[J].Resources Science,2009,31(5):763-771.(in Chinese))

    [7]吳洪偉.中國重點型灌區(qū)節(jié)水配套改造發(fā)展戰(zhàn)略研究[D].北京:中國農(nóng)業(yè)科學(xué)院,2012.(WU H W.Development strategy on systematic design of water saving projects in key medium-sized irrigation areas in China [D].Beijing:Chinese Academy of Agricultural Sciences,2012.(in Chinese))

    [8]康紹忠.水安全與糧食安全[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2014,22(8):880-885.(KANG S Z.Towards water and food security in China [J].Chinese Journal of Eco-Agriculture,2014,22(8):880-885.(in Chinese))

    [9]馮保清.我國不同尺度灌溉用水效率評價與管理研究[D].北京:中國水利水電科學(xué)研究院,2013.(FENG B Q.Study on the evaluation and management of irrigation water use efficiency for different scales in countrywide [D].Beijing:China Institute of Water Resources & Hydropower Research(IWHR),2013.(in Chinese))

    [10]代俊峰,崔遠來.基于SWAT的灌區(qū)分布式水文模型——Ⅰ.模型構(gòu)建的原理與方法[J].水利學(xué)報,2009,40(2):145-152.(DAI J F,CUI Y L.Distributed hydrological model for irrigation area based on SWATⅠ.Principle and method [J].Journal of Hydraulic Engineering,2009,40(2):145-152.(in Chinese))

    [11]代俊峰,崔遠來.基于SWAT的灌區(qū)分布式水文模型——Ⅱ.模型應(yīng)用[J].水利學(xué)報,2009,40(3):311-318.(DAI J F,CUI Y L.Distributed hydrological model for irrigation area based on SWAT Ⅱ.Model application [J].Journal of Hydraulic Engineering,2009,40(3):311-318.(in Chinese))

    [12]歐陽威,黃浩波,張璇,等.基于SWAT模型的平原灌區(qū)水量平衡模擬研究[J].灌溉排水學(xué)報,2015,34(1):17-22.(OU Y W,HUANG H B,ZHANG X,et al.Simulation of water balance in Hetao irrigation district using SWAT model [J].Journal of Irrigation and Drainage,2015,34(1):17-22.(in Chinese))

    [13]馬放,姜曉峰,王立,等.基于SWAT模型的亞流域劃分方法研究[J].中國給水排水,2015,31(7):53-57.(MA F,JIANG X F,WANG L,et al.Study on subwatershed delineation based on SWAT Model[J].China Water & Wastewater,2015,31(7):53-57.(in Chinese))

    [14]鄭捷,李光永,韓振中,等.改進的SWAT模型在平原灌區(qū)的應(yīng)用[J].水利學(xué)報,2011,42(1):88-97.(ZHENG J,LI G Y,HAN Z Z,et al.Application of modified SWAT model in plain irrigation district [J].Journal of Hydraulic Engineering,2011,42(1):88-97.(in Chinese))

    [15]羅運祥,蘇保林,李卉,等.平原河網(wǎng)地區(qū)SWAT模型子流域概化方法[J].北京師范大學(xué)學(xué)報(自然科學(xué)版),2011,47(4):415-418.(LUO Y X,SU B L,LI H,et al.Sub-basin delineation methods for SWAT model in plain river network region[J].Journal of Beijing Normal University(Natural Science),2011,47(4):415-418.(in Chinese))

    [16]DUKE G D,KIENZLE S W,JOHNSON D L,et al.Improving overland flow routing by incorporating ancillary road data into digital elevation models[J].Journal of Spatial Hydrology,2003,3(2):1-27.

    [17]王全九,汪志榮,張建豐,等.層狀土入滲機制與數(shù)學(xué)模型[J].水利學(xué)報,1998(S1):77-80.(WANG Q J,WANG Z R,ZHANG J F,et al.Infiltration mechanism of layered soil and its simulation model[J].Journal of Hydraulic Engineering,1998( S1):77-80.(in Chinese))

    [18]甘永德,賈仰文,仇亞琴,等.降雨條件下分層土壤入滲特性[J].水土保持學(xué)報,2012,26(5):217-219.(GAN Y D,JIA Y W,CHOU Y Q,et al.Stratified soil infiltration characteristics during rainfall[J].Journal of Soil and Water Conservation,2012,26(5):217-219.(in Chinese))

    [19]JIA Y,TAMAI N.Modeling infiltration into a multi-layered soil during an unsteady rain [J].Annual Journal of Hydraulic Engineering,JSCE,1997(41):31-36.

    [20]郭元裕.農(nóng)田水利學(xué)[M].第三版.北京:中國水利水電出版社,1997.(GUO Y Y.Irrigation and Drainage Engineering[M].The third edition.Beijing:China Water & Power Press,1997.(in Chinese))

    [21]包莉.羅江灌區(qū)雜交稻騰發(fā)量變化及控水條件對產(chǎn)量構(gòu)成的影響[D].雅安:四川農(nóng)業(yè)大學(xué),2016.(BAO L.The change of Luojiang irrigation hybrid rice evapotranspiration and the influence of water control conditions on the yield of rice[D].Yaan:Sichuan Agricultural University,2016.(in Chinese))

    [22]茆智,李遠華,李會昌.逐日作物需水量預(yù)測數(shù)學(xué)模型研究[J].武漢水利電力大學(xué)學(xué)報,1995,28(3):253-259.(MAO Z,LI Y H,LI H C.Study of mathematic model for forecasting daily crop evapotranspiration [J].Journal of Wuhan University of Hydropower & Electric Engineering,1995,28(3):253-259.(in Chinese))

    [23][JP+2]ARNOLD J G,KINIRY J R,SRINIVASAN R,等.ArcSWAT2009輸入輸出文件手冊[M].鄭州:黃河水利出版社,2012.(ARNOLD J G,KINIRY J R.,SRINIVASAN R,et al.Soil and Water Assessment Tool Input/ Output File Documentation Version 2009 [M].Zhengzhou:Yellow River water conservancy press,2012.(in Chinese))

    [24]NEITSCH S L,ARNOLD J G,KINIRY J R,等.SWAT2009理論基礎(chǔ)[M].鄭州:黃河水利出版社,2012.(NEITSCH S L,ARNOLD J G,KINIRY J R,et al.Soil and Water Assessment Tool Theoretical Documentation Version 2009[M].Zhengzhou:Yellow River water conservancy press,2012.(in Chinese))

    [25]邵東國,楊霞,徐保利,等.南方丘陵區(qū)多層土壤結(jié)構(gòu)水平和垂向滲流特征試驗與模擬研究[J].水利學(xué)報,2017,48(7):799-807.(SHAO D G,YANG X,XU B L,et al.Field experiment and simulation of lateral and vertical water flow in multi-layer soil in hilly region of southern China [J].Journal of Hydraulic Engineering,2017,48(7):799-807.(in Chinese))

    [26]胡能杰,邵東國,陳述,等.基于系統(tǒng)動力學(xué)的稻田塘堰系統(tǒng)水轉(zhuǎn)化模擬及驗證[J].農(nóng)業(yè)工程學(xué)報,2017,33(12):130-135.(HU N J,SHAO D G,CHEN S,et al.Simulation and verification of water transformation of rice paddy and pond system based on system dynamics[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2017,33(12):130-135.(in Chinese))

    [27]鄭明昇.桃園灌區(qū)之區(qū)域回歸水分析研究[D].國立中央大學(xué)水文科學(xué)研究所,1996.(ZHENG M S.Analysis of regional return water in Taoyuan irrigation district[D].Institute of Hydrology of National Central University,1996.(in Chinese))

    [28]段明印,李傳奇,熊劍智.基于偏相關(guān)法的SWMM模型參數(shù)全局敏感性分析[J].人民長江,2017,48(17):53-57.(DUAN M Y,LI C Q,XIONG J Z.Parameters sensitivity analysis of SWMM model based on partial correlation method [J].Yangtze River,2017,48(17):53-57.(in Chinese))

    [29]DESSU S B,MELESSE A M.Modelling the rainfall-runoff process of the Mara River basin using the Soil and Water Assessment Tool[J].Hydrological Processes,2012,26(26):4038-4049.

    [30]TURCOTTE R,F(xiàn)ORTIN J P,ROUSSEAU A N,et al.Determination of the drainage structure of a watershed using a digital elevation model and a digital river and lake network[J].Journal of Hydrology,2001,240:225-242.

    猜你喜歡
    敏感性分析模型
    一半模型
    p150Glued在帕金森病模型中的表達及分布
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    區(qū)域性系統(tǒng)性金融風(fēng)險影響因素研究
    區(qū)域性系統(tǒng)性金融風(fēng)險影響因素研究
    城鎮(zhèn)化下基本養(yǎng)老保險制度運行可持續(xù)性的敏感性分析
    商(2016年35期)2016-11-24 10:42:45
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    興和縣萬鵬養(yǎng)殖合作社經(jīng)濟效益評價
    三级经典国产精品| 你懂的网址亚洲精品在线观看 | 变态另类成人亚洲欧美熟女| eeuss影院久久| 我要搜黄色片| 一级av片app| 天堂网av新在线| 人人妻人人澡人人爽人人夜夜 | 日韩成人伦理影院| 国产乱人视频| 久久精品综合一区二区三区| 国产美女午夜福利| 久久久久久久久中文| 精品一区二区三区人妻视频| 国产高清视频在线观看网站| 亚洲欧美中文字幕日韩二区| 日韩视频在线欧美| 亚洲第一区二区三区不卡| 久久国产乱子免费精品| 黄色视频,在线免费观看| 久久久久网色| 久久久久久久亚洲中文字幕| 国产精品久久久久久久电影| 亚洲av男天堂| 99国产精品一区二区蜜桃av| 别揉我奶头 嗯啊视频| 亚洲国产日韩欧美精品在线观看| 波多野结衣高清无吗| 国产黄a三级三级三级人| 亚洲最大成人手机在线| 天堂网av新在线| 午夜福利在线观看免费完整高清在 | 国产乱人视频| 国产精品.久久久| 日本在线视频免费播放| 亚洲人与动物交配视频| 久久中文看片网| 岛国毛片在线播放| 禁无遮挡网站| 久久久久久伊人网av| 三级经典国产精品| 亚洲自拍偷在线| 亚洲激情五月婷婷啪啪| 日韩一本色道免费dvd| 亚洲精品亚洲一区二区| 国产黄a三级三级三级人| 久久午夜亚洲精品久久| 青春草视频在线免费观看| 国产老妇女一区| 国产蜜桃级精品一区二区三区| 国产精品久久久久久av不卡| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 久久这里只有精品中国| 久久久午夜欧美精品| 中文字幕免费在线视频6| 精品久久国产蜜桃| 在线a可以看的网站| 亚洲图色成人| 欧美性猛交黑人性爽| 在线观看av片永久免费下载| 美女脱内裤让男人舔精品视频 | 精品国产三级普通话版| 一边亲一边摸免费视频| 日韩欧美精品免费久久| 高清在线视频一区二区三区 | 日韩一本色道免费dvd| av又黄又爽大尺度在线免费看 | 日本黄色片子视频| 两性午夜刺激爽爽歪歪视频在线观看| 少妇被粗大猛烈的视频| 色播亚洲综合网| 国产久久久一区二区三区| 久久精品国产亚洲av涩爱 | 少妇的逼好多水| 国产精品永久免费网站| 色哟哟·www| 人人妻人人看人人澡| 亚洲精品久久国产高清桃花| 国产高清有码在线观看视频| 国产精品一及| 日韩精品有码人妻一区| 精品一区二区免费观看| 午夜a级毛片| 国产成人91sexporn| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 亚洲精华国产精华液的使用体验 | 欧美成人免费av一区二区三区| 少妇熟女aⅴ在线视频| av国产免费在线观看| 久久韩国三级中文字幕| 亚洲精品乱码久久久v下载方式| 看免费成人av毛片| 国产蜜桃级精品一区二区三区| 性色avwww在线观看| 亚洲在线观看片| 国产爱豆传媒在线观看| or卡值多少钱| 高清毛片免费看| 99久久九九国产精品国产免费| 国产一区二区在线观看日韩| 国产高清有码在线观看视频| 久久久久久久亚洲中文字幕| 国产精品综合久久久久久久免费| 久久久精品欧美日韩精品| 女人被狂操c到高潮| 99久久久亚洲精品蜜臀av| 亚洲美女视频黄频| 久久久成人免费电影| 欧美一区二区国产精品久久精品| 中国国产av一级| 欧美性猛交黑人性爽| 亚洲aⅴ乱码一区二区在线播放| 午夜福利视频1000在线观看| 国产91av在线免费观看| 国产乱人偷精品视频| 国内精品久久久久精免费| 国产白丝娇喘喷水9色精品| 久久久久久久久久久免费av| kizo精华| 九九在线视频观看精品| 99久久中文字幕三级久久日本| 午夜福利在线观看吧| 人妻久久中文字幕网| 99久国产av精品| 国产一级毛片在线| 午夜精品国产一区二区电影 | 成人综合一区亚洲| 女同久久另类99精品国产91| 欧美日韩精品成人综合77777| 国产精品人妻久久久影院| 亚洲精品久久国产高清桃花| 精品久久久久久久久久免费视频| 午夜福利在线在线| 亚洲高清免费不卡视频| 少妇猛男粗大的猛烈进出视频 | 亚洲av成人av| 人妻久久中文字幕网| 国产亚洲精品久久久久久毛片| 国产成人福利小说| 精品日产1卡2卡| 日本黄色视频三级网站网址| 我要搜黄色片| 级片在线观看| 国产三级在线视频| 一本久久中文字幕| 中国美白少妇内射xxxbb| 99国产极品粉嫩在线观看| 99国产精品一区二区蜜桃av| a级一级毛片免费在线观看| 亚洲乱码一区二区免费版| 我要搜黄色片| 三级毛片av免费| 国产精品美女特级片免费视频播放器| 亚洲成人精品中文字幕电影| 少妇人妻一区二区三区视频| 免费电影在线观看免费观看| 国产精品女同一区二区软件| 免费大片18禁| 最近手机中文字幕大全| 51国产日韩欧美| 久久久久久久亚洲中文字幕| 亚洲国产日韩欧美精品在线观看| 看免费成人av毛片| 哪个播放器可以免费观看大片| 亚洲欧美成人综合另类久久久 | 悠悠久久av| 欧美+亚洲+日韩+国产| 国产精品av视频在线免费观看| 亚洲av免费在线观看| 亚洲国产欧美人成| 好男人视频免费观看在线| 一个人观看的视频www高清免费观看| 久久精品国产亚洲网站| 一级av片app| 国产在线精品亚洲第一网站| 色5月婷婷丁香| 看十八女毛片水多多多| 国产精品国产高清国产av| 日韩av不卡免费在线播放| 成人国产麻豆网| 亚洲最大成人av| 一级黄片播放器| 白带黄色成豆腐渣| 嫩草影院精品99| 亚洲第一电影网av| 久久99蜜桃精品久久| 丝袜美腿在线中文| 日韩欧美在线乱码| 97超视频在线观看视频| 亚洲精品456在线播放app| 亚洲经典国产精华液单| 国产精品三级大全| 亚洲第一电影网av| 国产一区二区三区av在线 | 国产爱豆传媒在线观看| av视频在线观看入口| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线| 99热全是精品| 久久久久久久久久久免费av| 男女边吃奶边做爰视频| 亚洲久久久久久中文字幕| 一进一出抽搐gif免费好疼| 99国产精品一区二区蜜桃av| 亚洲人成网站在线观看播放| 少妇人妻一区二区三区视频| 热99re8久久精品国产| 国产精品1区2区在线观看.| 男女边吃奶边做爰视频| 日本在线视频免费播放| 亚洲精品自拍成人| 国产私拍福利视频在线观看| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 赤兔流量卡办理| 日本成人三级电影网站| 亚洲精品成人久久久久久| 69人妻影院| 国产单亲对白刺激| 国产欧美日韩精品一区二区| 99久久久亚洲精品蜜臀av| 久久人人爽人人爽人人片va| 一夜夜www| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 高清日韩中文字幕在线| 亚洲精品456在线播放app| 欧美不卡视频在线免费观看| 久久久精品欧美日韩精品| 亚洲天堂国产精品一区在线| 看非洲黑人一级黄片| 一本久久精品| 国产亚洲精品av在线| 国产精品一及| 国产精品久久电影中文字幕| 日本与韩国留学比较| 99久国产av精品国产电影| 亚洲精华国产精华液的使用体验 | 亚洲人成网站高清观看| 国产av在哪里看| 中国美女看黄片| 18+在线观看网站| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 成年版毛片免费区| 一夜夜www| 美女被艹到高潮喷水动态| www.色视频.com| 精品无人区乱码1区二区| 熟女电影av网| 久久午夜福利片| 午夜福利成人在线免费观看| 久久人人精品亚洲av| 99久国产av精品国产电影| 麻豆一二三区av精品| 亚洲av.av天堂| 乱系列少妇在线播放| 大香蕉久久网| 淫秽高清视频在线观看| 亚洲三级黄色毛片| 草草在线视频免费看| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 乱人视频在线观看| 日韩一区二区视频免费看| 亚洲成人中文字幕在线播放| 最近的中文字幕免费完整| 亚洲,欧美,日韩| 日韩一本色道免费dvd| 91在线精品国自产拍蜜月| 又黄又爽又刺激的免费视频.| 熟女人妻精品中文字幕| 人体艺术视频欧美日本| 美女xxoo啪啪120秒动态图| 级片在线观看| 日韩制服骚丝袜av| 麻豆乱淫一区二区| 国产精品一区二区性色av| 一夜夜www| 久久这里有精品视频免费| 18禁黄网站禁片免费观看直播| 禁无遮挡网站| av天堂中文字幕网| 美女内射精品一级片tv| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站| 一本精品99久久精品77| 国产蜜桃级精品一区二区三区| 久久久久性生活片| ponron亚洲| 一级毛片我不卡| 国产成人一区二区在线| 久久精品影院6| 国产av一区在线观看免费| 欧美日韩在线观看h| 国产亚洲精品久久久com| h日本视频在线播放| 偷拍熟女少妇极品色| 看片在线看免费视频| 国产精品免费一区二区三区在线| 又爽又黄a免费视频| 国产三级中文精品| 美女 人体艺术 gogo| 永久网站在线| 成人无遮挡网站| 久久鲁丝午夜福利片| 色综合站精品国产| a级一级毛片免费在线观看| 亚洲自偷自拍三级| 一本精品99久久精品77| 久久亚洲精品不卡| 亚洲图色成人| 久久国产乱子免费精品| 国产黄色视频一区二区在线观看 | 亚州av有码| 中文字幕免费在线视频6| 神马国产精品三级电影在线观看| 午夜激情欧美在线| 亚洲av中文av极速乱| 久久九九热精品免费| 小说图片视频综合网站| 午夜福利在线观看吧| 婷婷亚洲欧美| 99久久无色码亚洲精品果冻| 亚洲经典国产精华液单| 亚洲av熟女| 欧美+日韩+精品| 国产精品精品国产色婷婷| 免费一级毛片在线播放高清视频| 中文精品一卡2卡3卡4更新| 国产精品永久免费网站| 免费av毛片视频| 国产久久久一区二区三区| 国产不卡一卡二| 成年免费大片在线观看| 国产三级在线视频| 久久久国产成人精品二区| 99国产极品粉嫩在线观看| 51国产日韩欧美| 国产爱豆传媒在线观看| 久久精品夜色国产| 中文字幕av成人在线电影| 老司机影院成人| av天堂中文字幕网| 午夜免费激情av| 人体艺术视频欧美日本| 中国美女看黄片| 亚洲av二区三区四区| 久久6这里有精品| 成人高潮视频无遮挡免费网站| 日本免费a在线| 国产伦一二天堂av在线观看| 观看美女的网站| 成人永久免费在线观看视频| 搡老妇女老女人老熟妇| www.av在线官网国产| 乱人视频在线观看| 特级一级黄色大片| 欧美在线一区亚洲| 搞女人的毛片| 97超视频在线观看视频| 欧美丝袜亚洲另类| 91久久精品国产一区二区三区| 人妻系列 视频| 高清午夜精品一区二区三区 | 国产色爽女视频免费观看| 成人三级黄色视频| 精品人妻熟女av久视频| 精品不卡国产一区二区三区| 男女边吃奶边做爰视频| 久久草成人影院| 色哟哟哟哟哟哟| 欧美在线一区亚洲| 日韩成人伦理影院| 哪里可以看免费的av片| 在线播放国产精品三级| 男人舔女人下体高潮全视频| 夜夜夜夜夜久久久久| 精品熟女少妇av免费看| 听说在线观看完整版免费高清| 亚洲国产精品国产精品| 又粗又爽又猛毛片免费看| 精品免费久久久久久久清纯| 蜜桃久久精品国产亚洲av| 成年女人看的毛片在线观看| 午夜亚洲福利在线播放| 一级毛片aaaaaa免费看小| 在线观看av片永久免费下载| 亚洲高清免费不卡视频| 免费观看精品视频网站| 最新中文字幕久久久久| 日韩一区二区三区影片| 国产高清不卡午夜福利| 亚州av有码| 国产视频内射| 国产一级毛片七仙女欲春2| 亚洲av一区综合| 特级一级黄色大片| 天天躁日日操中文字幕| 成年女人看的毛片在线观看| 精品人妻熟女av久视频| 精品一区二区三区人妻视频| 热99re8久久精品国产| 只有这里有精品99| 亚洲一区高清亚洲精品| 久久精品久久久久久噜噜老黄 | 色综合色国产| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 97热精品久久久久久| 变态另类成人亚洲欧美熟女| 久久久色成人| 午夜福利在线在线| 老司机福利观看| 久久99热这里只有精品18| 国产综合懂色| 国产精品日韩av在线免费观看| 在线播放无遮挡| 一夜夜www| 国语自产精品视频在线第100页| 一区二区三区免费毛片| 成人毛片60女人毛片免费| av女优亚洲男人天堂| 日韩精品青青久久久久久| 亚洲成a人片在线一区二区| 91aial.com中文字幕在线观看| 麻豆av噜噜一区二区三区| 最好的美女福利视频网| 成年av动漫网址| 婷婷六月久久综合丁香| 成人综合一区亚洲| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区 | 在线播放国产精品三级| АⅤ资源中文在线天堂| 麻豆成人av视频| 日本爱情动作片www.在线观看| 亚洲欧美日韩东京热| 两个人视频免费观看高清| 夫妻性生交免费视频一级片| 最近手机中文字幕大全| 国产高潮美女av| 午夜激情欧美在线| 色噜噜av男人的天堂激情| 中文在线观看免费www的网站| АⅤ资源中文在线天堂| 久久久久久久久久成人| 久久韩国三级中文字幕| 三级男女做爰猛烈吃奶摸视频| 亚洲av免费高清在线观看| 国产三级中文精品| 一卡2卡三卡四卡精品乱码亚洲| 神马国产精品三级电影在线观看| 免费看日本二区| 亚洲美女视频黄频| 欧美最黄视频在线播放免费| 亚洲av一区综合| 欧美最新免费一区二区三区| 18禁裸乳无遮挡免费网站照片| 欧美精品一区二区大全| 精品无人区乱码1区二区| 五月伊人婷婷丁香| 天堂影院成人在线观看| 一级毛片电影观看 | 插逼视频在线观看| 国产视频首页在线观看| 亚洲av成人精品一区久久| 国产一区二区亚洲精品在线观看| 伦精品一区二区三区| 欧美激情在线99| av.在线天堂| 午夜爱爱视频在线播放| 日本与韩国留学比较| 国产久久久一区二区三区| 久久精品影院6| 麻豆成人av视频| 免费人成视频x8x8入口观看| 麻豆久久精品国产亚洲av| 麻豆精品久久久久久蜜桃| 亚洲精品日韩av片在线观看| 国产久久久一区二区三区| 久久人人爽人人爽人人片va| 亚洲精华国产精华液的使用体验 | 12—13女人毛片做爰片一| 亚洲一区高清亚洲精品| 免费电影在线观看免费观看| 97超碰精品成人国产| 亚洲七黄色美女视频| 日本黄色视频三级网站网址| 桃色一区二区三区在线观看| 亚洲av.av天堂| 99热精品在线国产| 日日干狠狠操夜夜爽| 成年女人看的毛片在线观看| 久久精品久久久久久久性| 少妇人妻一区二区三区视频| 最近最新中文字幕大全电影3| 国产精品日韩av在线免费观看| 国内精品一区二区在线观看| 亚洲欧美日韩高清专用| 久久国产乱子免费精品| 午夜久久久久精精品| 亚洲精品国产成人久久av| 欧美+日韩+精品| 久久久久九九精品影院| 午夜福利在线在线| 级片在线观看| 欧美性猛交╳xxx乱大交人| 日本成人三级电影网站| 久久精品国产亚洲网站| 热99在线观看视频| 波多野结衣巨乳人妻| 免费人成视频x8x8入口观看| 久久热精品热| 女的被弄到高潮叫床怎么办| 99热只有精品国产| 国产成人一区二区在线| 在线观看午夜福利视频| 亚洲真实伦在线观看| 精品人妻一区二区三区麻豆| 国产精品三级大全| h日本视频在线播放| 一进一出抽搐动态| 九九久久精品国产亚洲av麻豆| 波多野结衣巨乳人妻| 毛片一级片免费看久久久久| 国产高潮美女av| 亚洲av电影不卡..在线观看| 色哟哟·www| 欧美+日韩+精品| 国产精华一区二区三区| 91aial.com中文字幕在线观看| 高清在线视频一区二区三区 | 一区二区三区免费毛片| 欧美又色又爽又黄视频| 老师上课跳d突然被开到最大视频| 有码 亚洲区| 久久久久久久久久久免费av| 午夜爱爱视频在线播放| av国产免费在线观看| 青春草视频在线免费观看| 观看美女的网站| 亚洲一区二区三区色噜噜| 如何舔出高潮| 亚洲av免费在线观看| 伊人久久精品亚洲午夜| 亚洲国产色片| 国产精品av视频在线免费观看| 午夜免费男女啪啪视频观看| 少妇猛男粗大的猛烈进出视频 | 亚洲美女搞黄在线观看| 国产黄a三级三级三级人| 国产精品.久久久| av女优亚洲男人天堂| 成人国产麻豆网| 成人欧美大片| 久久6这里有精品| 日韩人妻高清精品专区| 亚洲国产精品国产精品| 国产精品爽爽va在线观看网站| 在线观看免费视频日本深夜| 久久人人爽人人片av| 九九久久精品国产亚洲av麻豆| 亚洲精华国产精华液的使用体验 | 只有这里有精品99| 一本精品99久久精品77| 一个人免费在线观看电影| 国产精品一区二区性色av| 一级av片app| 五月玫瑰六月丁香| 色综合色国产| 午夜爱爱视频在线播放| 国产不卡一卡二| 看黄色毛片网站| 午夜老司机福利剧场| 中文欧美无线码| 日本一本二区三区精品| 久久6这里有精品| 此物有八面人人有两片| 十八禁国产超污无遮挡网站| av在线天堂中文字幕| 午夜福利在线观看免费完整高清在 | 在线天堂最新版资源| 婷婷色av中文字幕| 精品午夜福利在线看| 精品少妇黑人巨大在线播放 | 免费人成在线观看视频色| 夜夜爽天天搞| 亚洲欧美精品专区久久| 国产综合懂色| 亚洲av熟女| 国产精品久久久久久亚洲av鲁大| 十八禁国产超污无遮挡网站| 天美传媒精品一区二区| 午夜福利高清视频| 国产乱人偷精品视频| 听说在线观看完整版免费高清| 亚洲一区二区三区色噜噜| 高清毛片免费看| 噜噜噜噜噜久久久久久91| 午夜福利在线观看免费完整高清在 | 欧美+亚洲+日韩+国产| 久久人妻av系列| 丝袜美腿在线中文| 国产亚洲5aaaaa淫片| av在线老鸭窝| 日本爱情动作片www.在线观看| 国内精品美女久久久久久| 99久久人妻综合| 欧美性猛交╳xxx乱大交人| 12—13女人毛片做爰片一| 亚洲天堂国产精品一区在线|