徐芬芳，王芊予，段亦多，降亞楠，2，潘嘉培

（1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院，陜西楊凌712100；2.西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室，陜西楊凌712100）

0 引言

我國(guó)是一個(gè)農(nóng)業(yè)大國(guó)，灌區(qū)在我國(guó)國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展中發(fā)揮重要作用，灌區(qū)用占全國(guó)40%的耕地生產(chǎn)了75%的糧食和90%以上的經(jīng)濟(jì)作物，為我國(guó)社會(huì)經(jīng)濟(jì)發(fā)展和糧食安全提供了重要的基礎(chǔ)保障[1]。隨著人口增長(zhǎng)、居民生活水平提高和氣候變化影響的加劇，人類社會(huì)對(duì)水資源的需求急劇增加，水資源供需矛盾日益突出，導(dǎo)致部分區(qū)域農(nóng)業(yè)用水供需矛盾加劇。尤其是干旱半干旱的西北地區(qū)，由于缺乏足夠的地表徑流來(lái)滿足灌溉用水，發(fā)展了大規(guī)模的井灌區(qū)與原有的渠灌形成了渠井結(jié)合灌溉模式，以井灌代替排水，形成了渠井雙灌模式[2]。

地表水和地下水聯(lián)合利用提高了灌溉用水保證率，增加了灌區(qū)應(yīng)對(duì)氣候變化的韌性，對(duì)保障灌區(qū)生態(tài)環(huán)境的良性循環(huán)和經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展發(fā)揮了重要作用。受氣候變化和人類活動(dòng)的雙重影響，我國(guó)部分渠井結(jié)合灌區(qū)的水資源供需矛盾日益突出。由于地下水資源的開發(fā)利用缺乏有效監(jiān)管，致使一些灌區(qū)地下水超采嚴(yán)重，渠井灌溉用水比例失調(diào)、地下水采補(bǔ)失衡，進(jìn)而引發(fā)了一系列農(nóng)業(yè)生態(tài)環(huán)境問(wèn)題，直接影響灌區(qū)水資源的高效安全利用和灌溉農(nóng)業(yè)的可持續(xù)發(fā)展[3]。灌區(qū)是氣候變化的敏感區(qū)和人類活動(dòng)的密集區(qū)，渠井結(jié)合灌區(qū)的水循環(huán)受諸多人類活動(dòng)的影響[4]，因此過(guò)程復(fù)雜且難以量化。傳統(tǒng)的研究方法（如統(tǒng)計(jì)模型）[5]可以得到水量平衡關(guān)系，但難以定量計(jì)算人類活動(dòng)因素對(duì)水循環(huán)的影響，難以做到對(duì)地下水循環(huán)過(guò)程的準(zhǔn)確描述[6]。

當(dāng)前部分井渠結(jié)合管理中存在諸多問(wèn)題，如運(yùn)行管理費(fèi)缺乏有效保障，灌溉水水價(jià)偏低，農(nóng)戶對(duì)灌區(qū)地下水涵養(yǎng)認(rèn)知不足，渠井分屬不同的管理部門等導(dǎo)致許多灌區(qū)沒(méi)有真正實(shí)現(xiàn)“渠井雙灌”，而是走向了“渠井分灌”之路[7]，使得地表水和地下水聯(lián)合調(diào)度優(yōu)勢(shì)無(wú)法發(fā)揮，灌區(qū)水資源高效利用和農(nóng)業(yè)可持續(xù)發(fā)展受到嚴(yán)峻挑戰(zhàn)。

隨著計(jì)算機(jī)技術(shù)的發(fā)展和灌區(qū)水資源管理水平的提高，模型在水循環(huán)機(jī)理研究和水資源管理利用方面發(fā)揮著越來(lái)越重要的作用[8]。當(dāng)前灌區(qū)水資源管理方式的改進(jìn)需要深入了解利益相關(guān)者用水行為的驅(qū)動(dòng)因素，即農(nóng)戶的灌溉行為應(yīng)該作為模型開發(fā)中重要的一環(huán)，國(guó)際水文科學(xué)協(xié)會(huì)也提倡應(yīng)多關(guān)注水文學(xué)和社會(huì)之間的相互作用和反饋以加深對(duì)水循環(huán)的理解[9]。在傳統(tǒng)灌區(qū)水文模型中，往往將人類的水資源管理活動(dòng)視作水循環(huán)動(dòng)力系統(tǒng)的外部因子，雖然模擬了復(fù)雜的物理過(guò)程，但往往忽略了水文系統(tǒng)與人類活動(dòng)之間的互饋機(jī)制。

灌區(qū)大規(guī)模引、抽、灌、排等人類活動(dòng)顯著改變了灌區(qū)水循環(huán)各環(huán)節(jié)的水分通量[10]，因此研究必須要定量考慮人類活動(dòng)的影響。Srinivasan 等[11]建議將用水戶的社會(huì)行為納入水循環(huán)的建模范圍，在此基礎(chǔ)上進(jìn)一步發(fā)展水文模型，考慮人類行為與水文過(guò)程的密切作用，構(gòu)建社會(huì)水文學(xué)模型來(lái)解決水資源管理問(wèn)題[12,13]。Pankaj 等[14]探討了如何將人與水之間的關(guān)系聯(lián)系起來(lái)，以管理當(dāng)?shù)氐挠盟枨螅瑫r(shí)減輕對(duì)水文循環(huán)的不利影響。Nazemi等[15]發(fā)現(xiàn)目前大量模型在代表人類用水需求的能力相當(dāng)有限，強(qiáng)調(diào)了傳統(tǒng)水資源管理模型的局限性。英國(guó)帝國(guó)理工學(xué)院的O′Keeffe 等[16]基于印度灌區(qū)的野外調(diào)查提出了一個(gè)包括農(nóng)民灌溉用水行為的社會(huì)水文學(xué)框架，提供了一個(gè)“自下而上”構(gòu)建灌區(qū)水文模型的方法，通過(guò)模型概化將農(nóng)戶的行為和決策與水文模型相結(jié)合，定量分析不同灌溉情景對(duì)灌區(qū)地下水位的影響，以幫助決策者進(jìn)行科學(xué)的決策。當(dāng)前盡管有部分考慮人類活動(dòng)的社會(huì)水文學(xué)模型，但針對(duì)渠井結(jié)合灌區(qū)考慮農(nóng)戶灌溉用水的社會(huì)水文學(xué)模型較少。

O′Keeffe 等的研究結(jié)果表明在模型開發(fā)中直接納入利益相關(guān)者的見(jiàn)解和經(jīng)驗(yàn)可以增加模型的適用性，但是文獻(xiàn)中沒(méi)有提供模型的源代碼，且印度與中國(guó)在灌區(qū)管理模式、基本政策和農(nóng)民行為存在許多差異，而且其模型設(shè)置的灌溉方式僅是渠灌或井灌的單一模式，并不符合渠井結(jié)合灌區(qū)的實(shí)際情況。

因此，本研究以典型渠井結(jié)合灌區(qū)——陜西省寶雞峽灌區(qū)為例，進(jìn)行了大量的田間調(diào)查，收集了作物種植類型、產(chǎn)量和種植成本，灌溉水源選擇、灌溉用水量等基礎(chǔ)數(shù)據(jù)，并將觀察到的農(nóng)戶灌溉行為概念化和程序化，基于Python 語(yǔ)言開發(fā)了一個(gè)農(nóng)戶灌溉用水模型。該模型將農(nóng)戶的灌溉用水行為與水文物理過(guò)程耦合，包含水文、作物產(chǎn)量、農(nóng)民生計(jì)三大模塊，在塬上、塬下分別設(shè)置了多種情景，模擬了不同渠井灌溉用水比例對(duì)地下水位動(dòng)態(tài)、作物產(chǎn)量和農(nóng)戶生計(jì)的影響，分區(qū)確定地下水取水用量和水位控制指標(biāo)，為準(zhǔn)確量化渠井結(jié)合灌區(qū)的地下水取用量提供依據(jù)，為灌區(qū)地表水地下水聯(lián)合高效安全利用提供強(qiáng)有力的技術(shù)支撐。

1 研究區(qū)域概況和調(diào)查方法

1.1 研究區(qū)域概況

寶雞峽灌區(qū)位于陜西省關(guān)中盆地，屬大陸性季風(fēng)氣候區(qū)，總面積2 355 km2。降水年內(nèi)分配不均，冬春夏多旱，秋季多澇，7-9月降水量占全年的50%左右。寶雞峽灌區(qū)是一典型的渠井結(jié)合灌區(qū)[17]，自渭河左岸引水，灌溉范圍包括寶雞、咸陽(yáng)、楊凌、西安4 個(gè)市（區(qū)）的14 個(gè)縣（區(qū)）。灌溉面積約為19.44 萬(wàn)hm2，糧食總產(chǎn)量和商品糧供應(yīng)分別占全省的1/7 和1/4，是陜西省最大的國(guó)有灌區(qū)，也是全省糧、油、菜、果的生產(chǎn)基地。主要農(nóng)作物有小麥、玉米、棉花和油菜等，按自然地理和工程布局分為塬上和塬下兩大灌溉系統(tǒng)[18]。

目前，寶雞峽灌區(qū)的供水模式以地表水與地下水聯(lián)合供水為主，形成了農(nóng)業(yè)、工業(yè)、城鄉(xiāng)鎮(zhèn)生活用水多元化的灌區(qū)用水格局。但灌區(qū)人均水資源可用量不足500 m3，水資源供需矛盾極為突出，在社會(huì)經(jīng)濟(jì)快速發(fā)展的大背景下，一系列水資源問(wèn)題接踵而至，如灌區(qū)渠井用水比例失調(diào)，地下水超采，枯水期渭河徑流量減少，河流沖沙能力降低等，水資源的高效安全利用存在很大隱患[19,20]。

1.2 調(diào)查方法

本研究采用社會(huì)科學(xué)領(lǐng)域采用的半結(jié)構(gòu)化訪談來(lái)收集灌溉水源選擇、農(nóng)戶灌溉行為特征和其他用水驅(qū)動(dòng)因素的定性和定量信息。傳統(tǒng)的結(jié)構(gòu)化訪談以相同方式向所有受訪者提出的一系列固定問(wèn)題，而半結(jié)構(gòu)化訪談則是圍繞某一個(gè)主題對(duì)受訪者進(jìn)行面試，能夠讓受訪者更充分地表達(dá)自己對(duì)相關(guān)問(wèn)題的理解。半結(jié)構(gòu)化訪談?dòng)兄诳焖?、有效地在?shù)據(jù)匱乏的地區(qū)收集豐富和詳細(xì)的數(shù)據(jù)，能更好地揭示各種事件的驅(qū)動(dòng)因素和農(nóng)戶用水決策背后的動(dòng)機(jī)[21]。

訪談主要在咸陽(yáng)、禮泉、寶雞、乾縣、扶風(fēng)5 個(gè)地區(qū)進(jìn)行，通過(guò)查閱文獻(xiàn)和典型農(nóng)戶問(wèn)詢?cè)O(shè)計(jì)了半結(jié)構(gòu)化訪談?wù){(diào)查問(wèn)卷，對(duì)信義、段家灣、常興、杏林、貞元等15 個(gè)管理站的200 個(gè)典型農(nóng)戶進(jìn)行了面對(duì)面談話，典型農(nóng)戶的選擇標(biāo)準(zhǔn)為：①種植作物為典型作物，作物包括種植冬小麥、夏玉米或果樹等；②灌溉經(jīng)驗(yàn)豐富；③長(zhǎng)期在家務(wù)農(nóng)，熟悉農(nóng)田灌溉情況（見(jiàn)圖1）。

調(diào)查的內(nèi)容主要分為4個(gè)方面：①農(nóng)戶作物種植情況（作物類型、種植面積、產(chǎn)量、經(jīng)營(yíng)情況）；②灌溉水源情況（渠、井灌作物類型及其用水量、灌溉方式、灌溉成本）；③土地流轉(zhuǎn)情況；④節(jié)水灌溉開展情況。

2 模型構(gòu)建和開發(fā)

2.1 模型概化

為充分了解農(nóng)戶灌溉用水行為對(duì)灌區(qū)水循環(huán)的影響，對(duì)野外調(diào)查的調(diào)查問(wèn)卷進(jìn)行分析總結(jié)，將研究區(qū)農(nóng)民灌溉用水情況總結(jié)，為后續(xù)構(gòu)建水文模型提供理論依據(jù)。野外調(diào)查的主要信息可概括為：①灌區(qū)地勢(shì)平坦，側(cè)向徑流量??；②渠井結(jié)合灌溉是灌區(qū)的主要灌溉方式，農(nóng)戶根據(jù)所處地理位置和水源條件選擇灌溉水源，在渠道設(shè)施完好時(shí)傾向用渠水灌溉糧食作物，渠系水量不足時(shí)用井水灌溉經(jīng)濟(jì)價(jià)值較高的經(jīng)濟(jì)作物；③農(nóng)戶通常遵循傳統(tǒng)灌溉實(shí)踐模式進(jìn)行灌溉，灌區(qū)管理局按照灌溉制度在相應(yīng)的時(shí)段開閘放水；④部分渠系工程年久失修，地下水滲漏量大；⑤農(nóng)戶反映灌溉時(shí)間增加時(shí)井灌區(qū)的地下水位會(huì)不斷下降，農(nóng)戶需支付更多電費(fèi)抽取相同數(shù)量的地下水，導(dǎo)致灌溉成本增加，當(dāng)渠灌比例上升后，灌區(qū)地下水位會(huì)有所回升。基于以上幾點(diǎn)，本研究構(gòu)建了農(nóng)戶灌溉用水模型，該模型包括水文、作物產(chǎn)量和農(nóng)民生計(jì)3個(gè)模塊。

2.2 水文模塊

針對(duì)作物根系土層進(jìn)行水量平衡分析，設(shè)植物根區(qū)耗散量為Dr，對(duì)作物根區(qū)土層建立日尺度的水量平衡方程：

式中：ETc為作物實(shí)際蒸散量，mm；R為作物根系對(duì)地下含水層的補(bǔ)給量，mm；P為降雨，mm；I為灌溉量，mm；t和T分別代表日和年。

假設(shè)當(dāng)作物根部區(qū)域的含水量超過(guò)田間持水量時(shí)，發(fā)生從根部區(qū)域到地下含水層的補(bǔ)給，補(bǔ)給量設(shè)為R，計(jì)算如下：

干旱缺水條件下，作物遭受水分脅迫，此時(shí)作物的實(shí)際蒸散量可以用下式進(jìn)行估算：

式中：ET0為參考作物蒸散量，用Penman-Monteith 公式估算；Kc為作物系數(shù)，反映了作物自身的生物學(xué)特性，它與作物種類、品種、生育期、作物群體葉面指數(shù)等因素有關(guān)，對(duì)于每種作物的作物系數(shù)的具體值可以參閱FAO文獻(xiàn)。

Ks是水分脅迫系數(shù)，計(jì)算如下：

其中，TAW是作物最大可利用的總有效土壤水量，不同作物及同一作物不同的生育階段對(duì)土壤含水量有不同的要求，通常將田間持水量作為土壤有效水分的上限，凋萎系數(shù)作為土壤有效水分的下限，因此TAW計(jì)算如下：

式中：θFC為田間持水量；θwp為凋萎系數(shù)；Zr為最大根深，mm。

RAW稱為容易被植物利用的土壤有效水分，是在水分脅迫發(fā)生前作物從根區(qū)內(nèi)吸收的土壤水分，用下式計(jì)算：

式中：ρ為水分脅迫發(fā)生前根區(qū)消耗的土壤水分占總有效土壤水的比例。

ρ值隨作物種類、作物蒸發(fā)蒸騰量ETc的大小而變化。FAO-56 提供了不同作物在ETc=5 mm/d 時(shí)的ρ值，在該條件下冬小麥、夏玉米的推薦ρ值為0.55。

模型運(yùn)行過(guò)程中，渠井灌溉用水比例(PCW)設(shè)為渠道灌溉水量占總灌溉量的比例。則含水層存儲(chǔ)量H的變化計(jì)算如下：

式中：H為單位面積含水層存儲(chǔ)量，mm；ε為渠道滲漏系數(shù)；Ic為渠灌用水量；Iw為井灌用水量。

則地下水埋深Hi(m)為前一天的埋深減去逐日地下水位變化量，計(jì)算如下：

2.3 作物模塊

在該模塊中，作物產(chǎn)量是農(nóng)民生計(jì)和農(nóng)業(yè)用水的重要聯(lián)系。它使用Doorenbos and Kassam (1979)開發(fā)的作物生產(chǎn)和蒸發(fā)之間的關(guān)系來(lái)計(jì)算，這可以表示為：

式中：Yx(T)為最大產(chǎn)量，kg/hm2，根據(jù)實(shí)地調(diào)查的田間報(bào)告確定；Ya(T)為實(shí)際產(chǎn)量，kg/hm2；ky為產(chǎn)量響應(yīng)因子；ETx為作物最大蒸發(fā)量，mm；ETc為作物實(shí)際蒸發(fā)量，mm；ds為播種日；dh為收獲日。

2.4 農(nóng)戶生計(jì)模塊

農(nóng)民一年的凈收入L(元)用農(nóng)民收入與支出的差額確定：

假設(shè)農(nóng)民收入全部來(lái)自于作物的售賣，則農(nóng)民收入為：

式中：nc為一個(gè)年度種植的農(nóng)作物數(shù)量；p為農(nóng)作物售賣的毛利潤(rùn)，元/kg；A為種植面積，hm2。

由于種植成本（化肥、種子，雇工等）每年基本不變，因此公式中已經(jīng)考慮了種植成本。農(nóng)民支出僅考慮灌溉成本（渠灌和井灌），則農(nóng)民每年支出計(jì)算如下：

式中：Zc(T)為渠灌支出；pc為渠道水水價(jià)，元/m3；Zw(T)為井灌支出，農(nóng)民僅支付電費(fèi)。

將1 t 水提升1 m 消耗的電量為0.002 72 kWh，η為電泵機(jī)效率，灌區(qū)抽水站實(shí)際運(yùn)行效率約為60%～70%，V為日提水量（t），f為機(jī)井用電綜合單價(jià)（元/kWh），Hn為提水凈泵揚(yáng)程（m），計(jì)算如下：

式中：μ為給水度；H0為地下水年初始埋深，m。

2.5 模型開發(fā)

農(nóng)戶灌溉用水模型采用Python（3.8.2）語(yǔ)言開發(fā)，使用的第三方庫(kù)主要為具有強(qiáng)大數(shù)據(jù)處理和分析能力的NumPy（1.18.4）、Sympy（1.0）和Pandas（1.0.3），將前期準(zhǔn)備好的數(shù)據(jù)以CSV 文件保存，以Excel 為可視化的操作界面，基于Pandas 庫(kù)的Dataframe 對(duì)數(shù)據(jù)進(jìn)行分析計(jì)算，并用Matplotlib（3.2.1）進(jìn)行圖形繪制。由于目前的模型版本（v1.0）需要通過(guò)手工調(diào)參的方式來(lái)校準(zhǔn)，故在模型開發(fā)時(shí)用字典分別為模型的3 個(gè)子模塊設(shè)置了參數(shù)集dic1、dic2、dic3，方便后續(xù)調(diào)參。

3 模型應(yīng)用

3.1 模型初始化和模型校準(zhǔn)

對(duì)灌區(qū)氣象站連續(xù)30年（1986-2016）的實(shí)測(cè)氣象資料系列進(jìn)行了水文頻率分析計(jì)算后，得到相應(yīng)的年降水頻率曲線，最終選取2003年作為豐水年，2006年作為平水年，2013年作為枯水年，其降雨量分別為771.6、445.0、360 mm。根據(jù)作物生長(zhǎng)期，水文年設(shè)定為10月初到翌年10月初。

根據(jù)調(diào)查，興平區(qū)和秦都區(qū)近十年來(lái)的灌溉模式都是渠井結(jié)合灌溉，故設(shè)置模型初始化參數(shù)時(shí)選取的兩個(gè)觀測(cè)井分別位于塬上的興平（066號(hào)井）和塬下的秦都（013號(hào)井），通過(guò)查閱《陜西省地下水埋深年鑒》和相關(guān)文獻(xiàn)資料，兩個(gè)典型灌域的模型初始參數(shù)設(shè)置如表1所示。

表1 模型初始化參數(shù)設(shè)置值表Tab.1 The model initializes the parameter setting value table

將2001-2011年秦都區(qū)013 號(hào)井地下水埋深的實(shí)測(cè)值和模擬值繪制如圖，并用納什效率系數(shù)NSE對(duì)模型進(jìn)行檢驗(yàn)，計(jì)算如下：

式中：yi為實(shí)測(cè)值；為模型預(yù)測(cè)值；NSE值越接近1，表明模型預(yù)測(cè)能力越好。

輸入典型灌區(qū)模型的相關(guān)參數(shù)和數(shù)據(jù)，輸出模型模擬結(jié)果和實(shí)測(cè)值對(duì)比如圖2所示，計(jì)算得NSE=0.58，表明估計(jì)誤差的方差較小，模型具有可靠性。由對(duì)比圖也可以看出，模擬值較好地模擬出了實(shí)測(cè)點(diǎn)的整體趨勢(shì)。

3.2 渠井用水比例變化情景設(shè)置

本文設(shè)置的情景為豐水年的地下水位在一定范圍內(nèi)上升，發(fā)揮地下水庫(kù)的調(diào)蓄作用；枯水年的地下水位在一定范圍內(nèi)下降，發(fā)揮地下水的抗旱作用；平水年的地下水位需要在一定范圍內(nèi)維持穩(wěn)定。模型設(shè)定了不同的渠井用水比例，保持灌溉總量不變，改變渠井用水比例，使得豐水年，平水年和枯水年的地下水位呈現(xiàn)不同幅度的變化，為決策者制定灌區(qū)水資源管理政策提供依據(jù)。模擬過(guò)程中設(shè)置的情景如表2所示。

表2 模型情景設(shè)置Tab.2 Model scenario settings

3.3 數(shù)據(jù)來(lái)源和數(shù)據(jù)輸入

根據(jù)寶雞峽灌區(qū)的灌溉制度給農(nóng)民分配灌溉水量，為了考慮灌溉的時(shí)空異質(zhì)性，模型在進(jìn)行編程時(shí)根據(jù)正態(tài)分布確定每日的灌溉量，正態(tài)分布的均值和標(biāo)準(zhǔn)差依據(jù)田間調(diào)查資料和灌溉制度率定，灌區(qū)小麥和玉米充分灌溉的灌溉制度如表3所示[22]。

表3 灌區(qū)小麥和玉米充分灌溉時(shí)的灌溉制度Tab.3 Irrigation system under the condition of full irrigation of wheat and corn in irrigation areas

通過(guò)查閱FAO 文獻(xiàn)[23]，計(jì)算每種作物的作物系數(shù)的平均值見(jiàn)表4。

表4 小麥、玉米的作物系數(shù)Tab.4 Crop coefficients for wheat and corn

氣象數(shù)據(jù)主要包括氣溫、降雨、土壤三方面的數(shù)據(jù)，來(lái)源于國(guó)家科學(xué)氣象數(shù)據(jù)中心，用水量數(shù)據(jù)來(lái)自《陜西省水利統(tǒng)計(jì)年鑒》，模型使用的社會(huì)經(jīng)濟(jì)數(shù)據(jù)主要來(lái)自于田間實(shí)地調(diào)查。

4 結(jié)果

4.1 分灌季的適宜渠井用水比例方案

為方便決策者管理，模型分灌季制定渠井用水比例方案。由灌溉定額，夏灌的灌溉量為1 200 m3/hm2，冬、春灌分別為825 和750 m3/hm2，因此夏灌對(duì)于地下水位的影響更大，為簡(jiǎn)化方案設(shè)置，將冬灌與春灌的PCW 設(shè)為定值，只改變夏灌的渠井用水比例，觀察地下水埋深的變化情況。根據(jù)“豐蓄枯補(bǔ)”原則，允許豐水年的地下水位在一定范圍內(nèi)上升，枯水年的地下水位在一定范圍內(nèi)下降，平水年的地下水位需要在一定范圍內(nèi)維持穩(wěn)定，因此，模擬過(guò)程中從方案1 到方案5 逐漸提高井灌所占比例，使地下水位從上升到穩(wěn)定再到下降：S1、X1 使地下水位上升了2 m；S2、X2 使地下水位上升了1 m；S3、X3 使地下水位保持穩(wěn)定；S4、X4 使地下水位下降了1 m；S5、X5 使地下水位下降了2 m（圖3 和圖4）。比例設(shè)置見(jiàn)表5和表6。

表5 塬上渠井用水比例Tab.5 The proportion of water used on the upper part of Baojixia irrigation district

表6 塬下渠井用水比例Tab.6 The proportion of water used on the lower part of Baojixiairrigation district

由于豐水年降雨量對(duì)灌區(qū)地下水動(dòng)態(tài)變化影響比枯水年更為明顯，所以豐水年即使全部井灌，地下水位也下降得也并不明顯[見(jiàn)圖3(a)]，故塬上豐水年只設(shè)置了4 種方案，塬下豐水年只設(shè)置了3 種方案。因此豐水年可以適度開采地下水，加大井灌比例避免地下水過(guò)度抬升避免土地鹽堿化。而枯水年可以加大渠灌比例涵養(yǎng)地下水。

塬上塬下3 個(gè)典型年不同方案下的地下水埋深變化如圖3和圖4所示。

4.2 不分灌季的適宜渠井用水比例方案

為更好地反映地下水埋深特征，將全年的渠井用水比例統(tǒng)一，同樣對(duì)于塬上塬下的不同典型年分別制定5種方案：渠井用水比例S1、X1 使地下水位上升了2 m；S2、X2 使地下水位上升了1 m；S3、X3 使地下水位保持穩(wěn)定；比例S4、X4 使地下水位下降了1 m；S5、X5 使地下水位下降了2 m（圖5和圖6）。

由于即使全部井灌，在豐水年時(shí)塬上塬下地下水位都下降不多，故塬上豐水年只設(shè)置了4種方案，塬下豐水年只設(shè)置了3種，比例見(jiàn)表7和表8。

表7 塬上不分灌季的渠井用水比例Tab.7 The proportion of water used on the upper part of Baojixia irrigation district during both irrigation season

表8 塬下不分灌季的渠井用水比例Tab.8 The proportion of water used on the lower part of Baojixiairrigation district during both irrigation season

塬上塬下3 個(gè)典型年不同方案下的地下水埋深變化如圖5和圖6所示。

綜上所述，根據(jù)模擬結(jié)果，在相同的地下水埋深變化下，塬上塬下渠井用水比例區(qū)別不大，但是當(dāng)埋深從上升2 m到下降2 m時(shí)，塬下井灌增加的幅度要更大，比如塬上平水年井灌比例從62%提升到82%，塬下從57%提升到90%（表7 和表8），即塬上的地下水埋深對(duì)PCW更加敏感。

這一現(xiàn)象與模型設(shè)置的初始參數(shù)和塬上塬下地理位置都有關(guān)，但主要原因是塬上的給水度小于塬下，因此輸入的灌溉總量相同時(shí)，塬上含水層存儲(chǔ)量的變化量大于塬下，因此PCW變化幅度相同時(shí)，塬上年地下水埋深變化更大。

4.3 作物產(chǎn)量與農(nóng)戶生計(jì)

由于模型模擬中，假設(shè)農(nóng)民嚴(yán)格按照灌溉制度進(jìn)行灌溉，各方案并未改變灌溉總量，只改變了渠井用水比例，因此作物產(chǎn)量固定不變，則農(nóng)戶的總收入不變，當(dāng)井灌增加時(shí)灌溉成本增加，反而導(dǎo)致農(nóng)民凈收入下降。但是作物產(chǎn)量會(huì)隨降雨量改變，各比例下的每公頃凈收入以及各水平年的作物產(chǎn)量如表9所示。

表9 作物產(chǎn)量和農(nóng)戶收入的模擬結(jié)果Tab.9 Simulation results of crop yields and farmers'income

根據(jù)模型模擬結(jié)果，隨著井灌占比增大導(dǎo)致灌溉成本增加，塬上和塬下灌區(qū)畝均凈收入也隨之下降，但塬上灌區(qū)下降幅度遠(yuǎn)大于塬下，塬下凈收入下降幅度較小。這是由于塬上灌區(qū)地下水埋深較大，地下水灌溉的費(fèi)用更大，因此井灌成本遠(yuǎn)大于塬下，而塬下地下水的埋深相對(duì)較小，農(nóng)民的井灌成本相應(yīng)也小。

5 結(jié)論

本文在大量野外調(diào)查的基礎(chǔ)上基于Python 語(yǔ)言開發(fā)了一個(gè)適用于渠井結(jié)合灌區(qū)的農(nóng)戶灌溉用水模型，該模型包括水文、作物和農(nóng)戶生計(jì)3個(gè)模塊。該模型將農(nóng)戶的灌溉用水行為與水文過(guò)程緊密耦合，模型在典型區(qū)域應(yīng)用的納什效率系數(shù)為0.58，表明模型可以用來(lái)模擬渠井結(jié)合灌區(qū)的地下水位，并能夠量化不同渠井用水比例對(duì)地下水位、作物產(chǎn)量和農(nóng)戶生計(jì)的影響，為灌區(qū)的水資源可持續(xù)管理提供了新的技術(shù)支持，同時(shí)突出了將農(nóng)戶灌溉用水納入灌區(qū)水資源管理政策制定的重要性。

根據(jù)模型模擬結(jié)果，若要維持地下水穩(wěn)定，塬上豐、平、枯水年推薦的全年渠井用水比例分別為2∶98、28∶72、30∶70，塬下豐、平、枯水年推薦的渠井用水比例分別為全井灌、28∶72、29∶71；根據(jù)豐蓄枯補(bǔ)原則，豐水年地下水位上升2 m 時(shí)，塬上的渠井比例為10∶90，塬下的渠井比例為17∶83；枯水年地下水位下降2 m時(shí)，塬上的渠井比例為21∶79，塬下的渠井比例為13∶87，總體來(lái)說(shuō)塬下用水比例要比塬上的變化幅度更大，即塬下地下水埋深對(duì)PCW 的變化更加敏感，故塬上塬下應(yīng)進(jìn)行分區(qū)管理，確定不同的地下水取水用量和水位控制指標(biāo)。

納入農(nóng)民行為可以使模型模擬結(jié)果更符合實(shí)際，本文通過(guò)向模型輸入不同的渠井結(jié)合用水比例，來(lái)分析其對(duì)地下水位動(dòng)態(tài)的影響，同時(shí)得到不同比例下的每畝凈收入，由于灌溉總量不變，不同PCW 下作物產(chǎn)量沒(méi)有改變，而畝凈收入下降。在模型的實(shí)際運(yùn)用過(guò)程中，可以根據(jù)實(shí)際情況改變灌溉總量，從而分析其對(duì)作物產(chǎn)量和農(nóng)戶收入的影響，以此為依據(jù)在小尺度范圍優(yōu)化農(nóng)民灌溉決策。

該模型也可以量化其他多種因素對(duì)灌溉用水的影響，如種植結(jié)構(gòu)變化，節(jié)水灌溉技術(shù)應(yīng)用等。比如渠井結(jié)合灌區(qū)采取滴灌等節(jié)水措施時(shí)，水源都取自于地下水，隨著節(jié)水灌溉技術(shù)的推行導(dǎo)致井灌用水增加，從而導(dǎo)致PCW 和灌溉用水總量都隨之改變，即可根據(jù)模型模擬結(jié)果分析節(jié)水灌溉技術(shù)對(duì)地下水位、作物產(chǎn)量和農(nóng)戶生計(jì)的影響，因此本模型在渠井結(jié)合灌區(qū)具有較好的應(yīng)用前景。