譚 倩，王淑萍，張?zhí)镦?/p>

基于實(shí)證數(shù)學(xué)規(guī)劃模型的農(nóng)業(yè)水價(jià)政策效應(yīng)模擬

譚 倩，王淑萍，張?zhí)镦?/p>

（中國農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，北京 100083）

水價(jià)政策的合理制定對水資源高效利用意義重大。水價(jià)政策的決策過程往往主觀性過強(qiáng)，缺乏水價(jià)政策效應(yīng)的量化研究方法。針對多水源多灌溉方式下的農(nóng)業(yè)灌區(qū)，基于實(shí)證數(shù)學(xué)規(guī)劃方法構(gòu)建農(nóng)業(yè)水價(jià)政策效應(yīng)模擬方法。將建立的方法應(yīng)用于甘肅省民勤縣進(jìn)行實(shí)證研究，分別模擬了3種水價(jià)政策下用水效益、灌溉用水量、種植占地和節(jié)水灌溉工程面積等對水價(jià)上漲的響應(yīng)機(jī)制，分析了不同政策下的水價(jià)上漲閾值和不同水價(jià)水平下的適宜政策。結(jié)果表明：在單一水價(jià)和兩部制水價(jià)政策下，民勤縣農(nóng)業(yè)水價(jià)宜處于0.24～2.10元/m3之間；在差別水價(jià)政策下，水價(jià)宜處于0.24～1.50元/m3之間；當(dāng)計(jì)量部分水價(jià)低于2.10元/m3或高于6元/m3時(shí)，差別水價(jià)政策的綜合效益最高；當(dāng)計(jì)量部分水價(jià)在2.10～6元/m3之間時(shí)，單一水價(jià)政策最為理想；實(shí)行差別水價(jià)政策且計(jì)量部分水價(jià)為1.50元/m3時(shí)取得的單方水效益最高。該文證實(shí)了民勤縣目前實(shí)行的差別水價(jià)政策是適宜的，水價(jià)的穩(wěn)步上漲不會降低用水的綜合效益，但應(yīng)在水價(jià)上調(diào)的同時(shí)考慮當(dāng)?shù)貙?shí)際情況，且注意加大對節(jié)水灌溉工程的投入和建立配套的農(nóng)戶收入補(bǔ)貼政策。該文建立的農(nóng)業(yè)水價(jià)政策效應(yīng)模擬方法具有較廣泛的適用性。

農(nóng)業(yè)；水；模型；實(shí)證數(shù)學(xué)規(guī)劃；水價(jià)政策；用水效益

0 引 言

水資源是最珍稀的自然資源之一，關(guān)系著人類的生存與發(fā)展[1-2]。水價(jià)是水資源管理的重要手段。一方面，價(jià)格杠桿有望緩解農(nóng)業(yè)用水緊張的問題[3-4]。另一方面，水價(jià)政策的制定與農(nóng)戶收入、農(nóng)業(yè)生產(chǎn)、能源消耗和糧食安全等其他經(jīng)濟(jì)、社會與環(huán)境要素密切相關(guān)[5-6]。水價(jià)過低會導(dǎo)致農(nóng)戶節(jié)水積極性下降、灌溉用水量增高、水資源利用效率降低等諸多負(fù)面問題；水價(jià)過高會影響正常的農(nóng)業(yè)生產(chǎn)活動(dòng)，降低農(nóng)戶的收入與生活質(zhì)量，不利于農(nóng)業(yè)的可持續(xù)發(fā)展。因此，在充分理解農(nóng)業(yè)水價(jià)政策效應(yīng)的基礎(chǔ)上合理制定水資源價(jià)格政策，具有重要的科學(xué)研究價(jià)值和應(yīng)用前景。

長期以來，中國實(shí)行無償或低價(jià)農(nóng)業(yè)水價(jià)政策，導(dǎo)致水資源被過度開發(fā)和大量浪費(fèi)[7]。當(dāng)前，國家已開始在部分地區(qū)推行農(nóng)業(yè)水價(jià)綜合改革，以期得到更合理的農(nóng)業(yè)水價(jià)政策實(shí)施方案。然而，農(nóng)業(yè)灌溉系統(tǒng)包含資源、技術(shù)、經(jīng)濟(jì)和社會等諸多互相關(guān)聯(lián)的子系統(tǒng)和驅(qū)動(dòng)因素，并涉及多種水源、灌溉技術(shù)與作物類型，其水價(jià)改革對農(nóng)業(yè)行為以及經(jīng)濟(jì)收益等帶來的影響錯(cuò)綜復(fù)雜[8-9]?，F(xiàn)有研究與實(shí)踐工作對水價(jià)政策的制訂多采取經(jīng)驗(yàn)與類比方法，主觀性過強(qiáng)，且鮮有針對多水源、多灌溉方式類農(nóng)業(yè)灌區(qū)的水價(jià)政策進(jìn)行效應(yīng)模擬、預(yù)測與比較[10-11]。

近年來，國內(nèi)外已出現(xiàn)一些研究水價(jià)政策效應(yīng)的量化方法。其中，計(jì)量經(jīng)濟(jì)學(xué)方法難以預(yù)測新政策的實(shí)施效果、無法真正反映用戶面對新政策的決策機(jī)制和過程，且需要大量的數(shù)據(jù)[12-13]。一般可計(jì)算均衡模型對樣本數(shù)量要求極高且方程參數(shù)難以準(zhǔn)確確定[14-15]。除這2種常用的經(jīng)濟(jì)學(xué)模型外，國外也有學(xué)者采用數(shù)學(xué)規(guī)劃方法來模擬政策效應(yīng)[16]，該方法通過添加約束來限定決策變量的變動(dòng)范圍，雖能保證模型計(jì)算結(jié)果與實(shí)際觀測值相近，但降低了模型的靈活性。實(shí)證數(shù)學(xué)規(guī)劃（positive mathematical programming，PMP）是一種能夠自校準(zhǔn)的政策效應(yīng)模擬方法，它對樣本數(shù)量要求低，可反映決策機(jī)制和過程，能將模型的輸出結(jié)果精確校準(zhǔn)至實(shí)際觀測值并在其基礎(chǔ)上預(yù)測新政策的實(shí)施效果[17-18]。Heckelei等[19]認(rèn)為PMP模型彌補(bǔ)了經(jīng)濟(jì)學(xué)模型和傳統(tǒng)數(shù)學(xué)規(guī)劃模型的不足，在這2類模型間起著橋梁作用。

目前PMP方法在農(nóng)業(yè)經(jīng)濟(jì)、環(huán)境、資源政策分析與評估中有了較為廣泛的應(yīng)用[20]。Iglesias等[21]首次將PMP方法應(yīng)用在水價(jià)政策效應(yīng)的研究中，模擬了計(jì)量水價(jià)政策下，灌溉用水、農(nóng)戶收入、肥料使用量等指標(biāo)隨水價(jià)上升的變化情況。隨后，Gallego-Ayala等[22-26]在Iglesias的基礎(chǔ)上，分別利用 PMP方法從不同角度對水價(jià)政策效應(yīng)開展了進(jìn)一步的研究。然而，以往基于PMP的水價(jià)政策效應(yīng)研究尚存在一些不足：1）在模型結(jié)構(gòu)上僅涵蓋單一水源，且未考慮糧食、能源供需方面的限制，對采用多水源、多灌溉方式的農(nóng)業(yè)灌區(qū)來說并不適用；2）對水價(jià)政策效應(yīng)規(guī)律的探討不夠深入。在對水價(jià)政策進(jìn)行模擬時(shí)，盡管部分研究得出了某一政策下各效應(yīng)指標(biāo)隨水價(jià)上升的變化趨勢，但未能進(jìn)一步分析響應(yīng)規(guī)律并提出水價(jià)閾值。在進(jìn)行不同水價(jià)政策間的橫向比較時(shí)，現(xiàn)有研究僅主觀選擇零星的水價(jià)政策方案或水平，且往往局限于農(nóng)戶收入、灌溉用水量等少數(shù)幾個(gè)效應(yīng)指標(biāo)，未能在多個(gè)效應(yīng)指標(biāo)下對一系列水價(jià)政策方案進(jìn)行量化比較后得出不同水價(jià)水平所對應(yīng)的適用政策（即政策適用范圍）。

基于以上背景，本研究針對多水源多灌溉方式下的灌區(qū)，構(gòu)建基于PMP的農(nóng)業(yè)水價(jià)政策效應(yīng)模擬方法，在灌區(qū)農(nóng)業(yè)水價(jià)政策制定前，對構(gòu)想新政策的實(shí)施效果進(jìn)行預(yù)測。以農(nóng)業(yè)水價(jià)綜合改革試點(diǎn)區(qū)甘肅省民勤縣為例，對不同農(nóng)業(yè)水價(jià)政策的實(shí)施效果進(jìn)行模擬和比較，進(jìn)一步得出不同政策制度下的水價(jià)上漲閾值和不同水價(jià)水平下的適宜政策。本研究旨在科學(xué)地揭示水價(jià)上漲的政策效應(yīng)，以系統(tǒng)最優(yōu)化為目標(biāo)，從理論角度探索不同政策下的水價(jià)上漲閾值和不同水價(jià)水平下的適宜政策，并以此來減少水價(jià)政策實(shí)施可能造成的經(jīng)濟(jì)、資源風(fēng)險(xiǎn)，而非囿于提出具體的水價(jià)建議值。這是由于，水價(jià)的制定從來都是自然、經(jīng)濟(jì)和社會、政治、歷史等多維因素共同作用的結(jié)果。在厘清水價(jià)上漲對各效應(yīng)指標(biāo)的客觀影響機(jī)制之后，還需征求多方意見，并綜合考慮當(dāng)?shù)噩F(xiàn)有水價(jià)、農(nóng)民承受能力、政府財(cái)政補(bǔ)貼等多種因素來確定具體的實(shí)施水價(jià)。

1 PMP方法建立

1.1 PMP核心思想

PMP模型研究始于20世紀(jì)90年代初，是在Heady等早期研究基礎(chǔ)上開展的[27-29]。最常用PMP模型由Richard教授于1995 年正式提出，其核心思想在于，假定觀測到的農(nóng)戶生產(chǎn)行為是最優(yōu)化的結(jié)果，認(rèn)為農(nóng)戶生產(chǎn)決策不但取決于可觀測的成本收益，還會受到不可觀測因素的影響，如農(nóng)民知識水平、土壤養(yǎng)分狀況等，并通過基準(zhǔn)期的校準(zhǔn)將這些不可觀測因素體現(xiàn)在模型中，使得模型輸出結(jié)果與實(shí)際觀測值一致[17,30]。本研究針對多水源多灌溉方式下的灌區(qū)，構(gòu)建基于PMP的農(nóng)業(yè)水價(jià)政策效應(yīng)模擬方法。在現(xiàn)行水價(jià)政策下，利用基準(zhǔn)期實(shí)際觀測數(shù)據(jù)和PMP方法對農(nóng)業(yè)水價(jià)線性規(guī)劃模型進(jìn)行校準(zhǔn)，計(jì)算模型中的未知參數(shù)。將未知參數(shù)代入模型，先驗(yàn)證模型在基準(zhǔn)期的最優(yōu)解與實(shí)際種植面積保持一致，再用于后續(xù)不同水價(jià)政策的效應(yīng)模擬。

1.2 PMP模型構(gòu)建及求解

1.2.1 建立并求解灌區(qū)農(nóng)業(yè)水價(jià)線性規(guī)劃模型得到校正約束對應(yīng)的對偶值

以農(nóng)戶純收入最大為目標(biāo)，以不同地區(qū)不同水源不同作物不同灌溉方式的種植面積為決策變量，構(gòu)建灌區(qū)農(nóng)業(yè)水價(jià)線性規(guī)劃模型，其中，約束式（3）為保證決策變量實(shí)際觀察值和計(jì)算結(jié)果一致的校正約束，其他約束均為資源約束，具體模型如下所示。

目標(biāo)函數(shù)：

約束條件：

1）決策變量實(shí)際觀察值校正約束

2）最大最小總種植面積的約束

耕作作物的總面積應(yīng)小于可用耕地的總面積，大于保留耕種的最小面積。

式中LMIN為當(dāng)?shù)乇Ａ舾N的最小種植面積，hm2；LMAX為當(dāng)?shù)乜捎酶氐目偯娣e，hm2。

3）地表水可用水量約束

用于各作物灌溉的地表水之和，不高于農(nóng)業(yè)灌溉的地表水可用總量。

式中SW為地表水總可用水量，m3；為地表水總可用量中可用于農(nóng)業(yè)灌溉的比例，%。

4）地下水可用水量約束

用于各種作物灌溉的地下水之和，不高于農(nóng)業(yè)灌溉的地下水可用總量。

式中GW為地下水總的可用水量，m3；為地下水總可用量中可用于農(nóng)業(yè)灌溉的比例，%。

5）可用電量約束

抽取地下水和滴灌節(jié)水技術(shù)的用電總量，不高于農(nóng)業(yè)灌溉的最大可用電量。

式中X,i,k,2為區(qū)域內(nèi)以地下水為水源的作物在灌溉方式下的種植面積，hm2；GE為抽取地下水灌溉的單位耗電量，kWh/m3；X,2,k,l為區(qū)域內(nèi)以水源的作物在滴灌方式下的種植面積，hm2；DE為滴灌技術(shù)灌溉的單位耗電量，kWh/m3；EN為農(nóng)業(yè)灌溉的最大可利用電量，kW·h。

6）糧食需求約束

為確保研究區(qū)各種作物的自給自足，引入了糧食安全約束，它根據(jù)總?cè)丝诤腿司魑镄枨蟠_定。

式中D為人均作物需求，每人每年糧食占有量不得少于400 kg，其他作物的最小需求量由《中國居民膳食指南》[31]查得，kg/(人·a)；TP為研究區(qū)人口總量，人。

7）人均最小耕地面積約束

為了維持人口和農(nóng)民的收入，確定最小耕地面積約束。

式中MAL為人均最小耕地面積，hm2/人。

8）決策變量非負(fù)約束

通過求解上述線性模型，可得到對應(yīng)校正約束的對偶值，該值間接提供了與成本相關(guān)的信息。

1.2.2 利用校正約束對應(yīng)的對偶值對灌區(qū)農(nóng)業(yè)水價(jià)線性規(guī)劃模型進(jìn)行校準(zhǔn)

根據(jù)邊際收益遞減的經(jīng)濟(jì)學(xué)假設(shè)，將成本函數(shù)設(shè)為一元二次[17]形式：

為保證模型的靈活性，去掉校正約束式（3），其他約束和1.2.1節(jié)中的模型保持一致，得到含未知參數(shù)的新模型。為求解未知參數(shù)矩陣α,i,k,s和γ,i,k,s，針對1.2.1節(jié)和1.2.2節(jié)的模型，利用庫恩塔克定理，對式（1）和式（12）求偏導(dǎo)可得：

式中1c,i,k,s和2c,i,k,s分別表示校正約束和資源約束對應(yīng)的對偶值，(x,i,k,s)表示在x,i,k,s點(diǎn)處起作用的校準(zhǔn)約束向量，(x,i,k,s)為在x,i,k,s點(diǎn)處起作用的資源約束向量集。

由式（14）和式（15）可以得到

將單位面積成本觀察值（不包括灌溉水費(fèi)）和模型聯(lián)系起來，使得不同區(qū)域不同水源不同作物不同灌溉方式的平均成本與基期成本觀察值相等[13,32]，則得到

結(jié)合式（17）和式（18），可得α,i,k,s和γ,i,k,s的最終值為

將參數(shù)α,i,k,s和γ,i,k,s代入式（12）中，可得到校準(zhǔn)后的灌區(qū)農(nóng)業(yè)水價(jià)實(shí)證數(shù)學(xué)規(guī)劃模型。

基于該模型在基準(zhǔn)期的最優(yōu)解與實(shí)際種植面積，對模型進(jìn)行驗(yàn)證，驗(yàn)證后用于不同水價(jià)政策的效應(yīng)模擬。改變式（2）中水價(jià)政策變量WP,i,k,s、B的值，用校準(zhǔn)后的模型模擬不同水價(jià)政策的實(shí)施效果。

2 實(shí)證研究

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

民勤縣（102°45′～103°55′E，38°20′～39°10′N）位于甘肅省西北方位，石羊河流域下游（圖1），是沙漠中的綠洲，鄰近地區(qū)的天然生態(tài)屏障[33]。民勤縣面積1.58萬km2，2015年常住人口24.12萬，耕地面積5.86萬hm2，縣內(nèi)包含紅崖山、環(huán)河、昌寧3個(gè)灌區(qū)，主要的糧食作物有小麥和玉米，主要的經(jīng)濟(jì)作物有葵花、蔬菜、瓜類、棉花等。主要的水源為地表水、地下水。主要的節(jié)水灌溉方式是膜下滴灌和溫室滴灌，傳統(tǒng)灌溉方式包括畦灌、溝灌、管灌和大水漫灌。

民勤縣年均降雨量113 mm，蒸發(fā)量2 644 mm，是中國境內(nèi)最干旱的地區(qū)之一。縣內(nèi)長期以種植高耗水糧食作物為主，農(nóng)業(yè)灌溉用水量占總用水量的70%以上，水資源供需矛盾突出。為緩解農(nóng)業(yè)用水緊張的態(tài)勢，民勤縣于2003年開始實(shí)施水價(jià)改革，逐步提高農(nóng)業(yè)水價(jià)。2013年，農(nóng)業(yè)灌溉的地表水價(jià)從0.108元/m3調(diào)整至0.24元/m3。地下水實(shí)行兩部制水價(jià)，即基本水價(jià)與計(jì)量水價(jià)相結(jié)合。其中，基本水價(jià)是為維持水利工程單位正常運(yùn)轉(zhuǎn)而向用水戶收取的最低基本費(fèi)用，按補(bǔ)償供水工程的直接固投、管理費(fèi)用和50%的折舊費(fèi)、修理費(fèi)的原則制定，以灌溉面積收取水費(fèi)；計(jì)量水價(jià)按實(shí)際灌溉用水量收取水費(fèi)。兩部制水價(jià)制度既能利用基本水價(jià)部分使供水成本得到補(bǔ)償，保證供水工程持續(xù)、穩(wěn)定運(yùn)轉(zhuǎn)；又能利用計(jì)量水價(jià)部分防止水資源浪費(fèi)，促進(jìn)農(nóng)業(yè)節(jié)約用水，目前已在國內(nèi)外許多大、中型農(nóng)業(yè)灌區(qū)實(shí)施[34-35]。民勤現(xiàn)行的基本水價(jià)為30元/hm2，地下水計(jì)量水價(jià)為0.235元/m3。此外，2013年武威市和民勤縣人民政府分別出臺了《武威市人民政府關(guān)于深化水價(jià)改革的實(shí)施意見》（武政發(fā)2013 91 號）和《民勤縣人民政府關(guān)于印發(fā)民勤縣深化水價(jià)改革實(shí)施意見的通知》（民政發(fā) 2013 109 號），提出推進(jìn)差別水價(jià)的實(shí)行[36]。

圖1 甘肅省民勤縣示意圖

2.2 情景設(shè)置與輸入數(shù)據(jù)

本研究考慮3種水價(jià)上漲政策情景，分別為：1）單一水量計(jì)價(jià)政策情景（以下簡稱“單一水價(jià)情景”）：單純按水量計(jì)價(jià)；2）兩部制水價(jià)政策情景（以下簡稱“兩部制水價(jià)情景”）：按計(jì)量水價(jià)和基本水價(jià)相結(jié)合的方式計(jì)價(jià)；3）兩部制水價(jià)耦合差別水價(jià)政策情景（以下簡稱“差別水價(jià)情景”）：在兩部制水價(jià)的基礎(chǔ)上實(shí)行差別水價(jià)，對于設(shè)施農(nóng)業(yè)和實(shí)施滴灌的大田節(jié)水作物收取較低的水價(jià)，對采用傳統(tǒng)種植方式的糧食作物收取較高的水價(jià)。

為提高政策的公眾接受度和易操作度，且為保持效應(yīng)模擬時(shí)的一致性從而更好地反映趨勢變化，本模型研究的基本假設(shè)和簡化如下：1）為促進(jìn)水價(jià)收取的公平性，在研究區(qū)實(shí)行同樣的水價(jià)政策；2）由于現(xiàn)行水價(jià)政策下地表水和地下水計(jì)量水價(jià)接近，所以模擬中不區(qū)分地下水和地表水價(jià)，耕地的基本水價(jià)設(shè)定為30元/hm2，計(jì)量水價(jià)部分以現(xiàn)狀水價(jià)0.24元/m3為起漲點(diǎn)；3）作物灌溉方式分為滴灌和傳統(tǒng)灌溉2種，不考慮傳統(tǒng)灌溉方式間成本、灌溉用水和作物單產(chǎn)的差異，滴灌灌溉方式主要考慮田間滴灌帶的投資成本；4）結(jié)合《武威市深化水權(quán)水價(jià)改革實(shí)施意見》和民勤差別水價(jià)實(shí)施情況，按照上漲后的水價(jià)，將滴灌灌溉的水價(jià)優(yōu)惠50%，將傳統(tǒng)種植方式下糧食作物的水價(jià)上浮50%；5）模擬水價(jià)上漲過程中，基本水價(jià)以及計(jì)量部分差別水價(jià)的優(yōu)惠與上浮幅度不變，僅提升計(jì)量部分的水價(jià)數(shù)值。

本研究基于民勤縣2015年的種植面積、種植成本、灌溉用水量、作物產(chǎn)量、人均作物需求量等各項(xiàng)基礎(chǔ)數(shù)據(jù)如表1所示。表中灌溉種植面積和產(chǎn)量數(shù)據(jù)由查閱《民勤縣2015年國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計(jì)資料匯編》[37]和實(shí)地調(diào)研得到；作物單價(jià)以及成本數(shù)據(jù)由《2016年全國農(nóng)產(chǎn)品成本收益資料匯編》[38]獲得；單位面積灌溉水量數(shù)據(jù)由武威市水務(wù)局獲得。

表1 民勤6種作物基礎(chǔ)數(shù)據(jù)

注：傳統(tǒng)灌溉方式包括畦灌、溝灌、管灌和大水漫灌。

Note:Traditional irrigation methods include border irrigation, furrow irrigation, pipe irrigation and flood irrigation.

3 結(jié)果與分析

3.1 模型率定與驗(yàn)證

表2 模型校準(zhǔn)過程中得到的相關(guān)參數(shù)

注：差值為種植面積模擬值減去實(shí)際值。

Note: Difference is obtained by subtracting actual value from simulated value.

3.2 不同水價(jià)制度情景下的結(jié)果分析

3.2.1 單一水價(jià)情景

計(jì)量水價(jià)是指僅按水量計(jì)算的水價(jià)。在該情景中，不存在基本水價(jià)，也不考慮不同灌溉方式下的水價(jià)差別。則式（2）在該情景中轉(zhuǎn)化為

在單一水價(jià)情景下對農(nóng)業(yè)水價(jià)實(shí)證數(shù)學(xué)規(guī)劃模型進(jìn)行求解。6種作物總的和分項(xiàng)的種植面積與滴灌面積比例（即節(jié)水灌溉工程面積比例）隨水價(jià)上升的變化情況分別如圖2和圖3所示。

圖3 不同作物種植面積和滴灌比例隨水價(jià)上升變化情況

由圖2可知，總種植面積隨著水價(jià)的升高而減小，當(dāng)水價(jià)升高至6元/m3時(shí)，總種植面積減小至最小容許耕地面積28 640 hm2并保持不變。由圖3可知，隨著總種植面積的減少，各作物的種植面積均有減少。其中，小麥和玉米的種植面積隨水價(jià)上升會迅速下降，分別在水價(jià)為0.42元/m3和2.10元/m3時(shí)趨于穩(wěn)定；葵花和棉花的種植面積隨水價(jià)的上升逐漸下降，且均在水價(jià)達(dá)到6元/m3之后趨于穩(wěn)定；瓜類和蔬菜的種植面積幾乎保持不變。這說明為應(yīng)對水價(jià)的上升，農(nóng)戶會選擇先壓縮低收益的糧食作物的種植面積，再減少高收益的經(jīng)濟(jì)作物的種植面積。隨著水價(jià)的上升，各作物的滴灌面積比例均有上升。其中，棉花的滴灌比例增加最大，小麥、玉米次之，而葵花、蔬菜和瓜類最不明顯。這說明，為應(yīng)對水價(jià)的上升，農(nóng)戶會選擇增加節(jié)水灌溉工程面積比例來減少灌溉用水量。

農(nóng)戶收入、灌溉用水量以及單方水效益隨水價(jià)上升的變化情況如圖4所示。灌溉用水量隨水價(jià)上升而迅速減少，當(dāng)水價(jià)上升至2.10元/m3時(shí)，下降幅度達(dá)到17.5%；之后下降過程變得平緩，下降速度不超過之前的1/3；當(dāng)水價(jià)進(jìn)一步上升到6元/m3后，灌溉用水量基本維持不變，這與總種植面積下降到最小容許面積有關(guān)。農(nóng)戶收入隨水價(jià)的提升勻速下降，當(dāng)水價(jià)為2.10元/m3時(shí)，下降幅度達(dá)到13.0%。單方水效益隨著水價(jià)的上升先增加后減小，且在水價(jià)為2.10元/m3時(shí)達(dá)到最高值13.74元/m3，比水價(jià)為0.24元/m3時(shí)上升了5.5%，之后不斷下降。

圖4 農(nóng)戶收入、灌溉用水量以及單方水效益隨水價(jià)上升變化情況

結(jié)合圖2和圖4中指標(biāo)的變化趨勢和拐點(diǎn)看出，在水價(jià)達(dá)到2.10元/m3之前，隨著水價(jià)的上升，灌溉用水量的減少幅度比農(nóng)戶收入下降的幅度更大，且單方水效益不斷上升，占地面積不斷減小。而在水價(jià)上升至2.10元/m3之后，作物種植面積和灌溉用水量的下降速度顯著變低，同時(shí)農(nóng)戶收入和單方水效益均不斷減小。因此，在單一水價(jià)情景下，水價(jià)應(yīng)處于0.24～2.10元/m3之間。

3.2.2 兩部制水價(jià)情景

兩部制水價(jià)是指計(jì)量水價(jià)和基本水價(jià)相結(jié)合的方式。在該政策情景中，除單一水價(jià)情景中的計(jì)量水價(jià)外，還需按種植面積支付基本水價(jià)。其灌溉水費(fèi)計(jì)算公式（21）可轉(zhuǎn)化為

基本水價(jià)的加收，會使兩部制水價(jià)情景下的水價(jià)比單一水價(jià)情景更高，但由于目前民勤縣基本水價(jià)收費(fèi)標(biāo)準(zhǔn)偏低（僅為30元/hm2），其影響較小?？傮w來說，2種政策的實(shí)施效果無顯著差別。具體的，兩部制水價(jià)情景中，隨水價(jià)的上升，種植面積、滴灌比例、農(nóng)戶收入、灌溉用水量和單方水效益的變化趨勢和變化拐點(diǎn)與單一水價(jià)情景基本一致。因此，與單一水價(jià)情景類似，從占地面積小、灌溉用水省、單方水效益高、對農(nóng)戶收入影響小4個(gè)方面考慮，如果實(shí)施兩部制水價(jià)政策，上漲后的水價(jià)應(yīng)處于0.24～2.10元/m3之間。在下述的兩部制水價(jià)和差別水價(jià)情景中，滴灌比例也呈現(xiàn)出與單一水價(jià)情景相似的規(guī)律。由于基本水價(jià)的存在，當(dāng)計(jì)量部分水價(jià)相等時(shí)，兩部制水價(jià)情景滴灌比例會高于單一水價(jià)情景；差別水價(jià)情景因兼有基本水價(jià)和對滴灌作物的優(yōu)惠，其滴灌比例又高于兩部制水價(jià)。

雖然單一和兩部制水價(jià)情景下各指標(biāo)的變化趨勢基本一致，但具體的指標(biāo)值有一定差別（圖5）。除單方水效益以外，兩部制情景下的總種植面積、灌溉用水量、農(nóng)戶收入和滴灌比例均低于單一水價(jià)情景。在作物總種植面積和灌溉用水量方面，2種政策情景的差值呈臺階式變化（圖5a和5b）。其中，第1個(gè)差值突變臺階始于水價(jià)上升至0.42元/m3時(shí)，即2種政策下小麥種植面積均下降至穩(wěn)定水平的節(jié)點(diǎn)；第2個(gè)差值突變臺階出現(xiàn)在水價(jià)進(jìn)一步上升至2.10元/m3時(shí)，此時(shí)2種政策下的玉米種植面積均基本穩(wěn)定、不再下降；第3個(gè)突變臺階出現(xiàn)于水價(jià)提升至6元/m3時(shí)，即2種政策下的總種植面積均降至最小允許值后，差值為0。

圖5 兩部制水價(jià)情景中各指標(biāo)減去單一水價(jià)情景中各指標(biāo)的差值

這說明，在2種政策下，總種植面積和灌溉用水量隨水價(jià)上升的下降速度相同，其政策間差值在一般情況下不隨水價(jià)上升而改變。但在使作物種植面積趨于穩(wěn)定的水價(jià)拐點(diǎn)上，對應(yīng)作物的種植面積和灌溉用水量差值變?yōu)?，這時(shí)總種植面積和灌溉用水量差值會發(fā)生跳躍式變化，且水價(jià)越高、差值越小。

在農(nóng)戶收入方面（圖5c），兩部制水價(jià)情景略低于單一水價(jià)情景，但2情景的相差幅度（（兩部制水價(jià)的農(nóng)戶收入-單一水價(jià)的農(nóng)戶收入）/兩部制水價(jià)的農(nóng)戶收入）僅為-0.05%左右，其絕對值低于灌溉用水量和單方水效益相差幅度絕對值。在單方水效益方面（圖5c），當(dāng)水價(jià)小于2.10元/m3時(shí)，兩部制水價(jià)情景略高于單一水價(jià)情景，其差距在水價(jià)為0.36元/m3時(shí)出現(xiàn)峰值。然而，當(dāng)水價(jià)達(dá)到2.10元/m3后出現(xiàn)逆轉(zhuǎn)，單一水價(jià)情景的單方水效益將略高于兩部制水價(jià)情景，其優(yōu)勢在水價(jià)上漲到6元/m3后將進(jìn)一步加大。

就不同水價(jià)水平對應(yīng)的適宜政策而言，當(dāng)水價(jià)在2.10元/m3以下時(shí)，兩部制水價(jià)政策單方水效益高于單一水價(jià)政策，種植面積、灌溉用水量低于單一水價(jià)政策，農(nóng)戶收入方面雖低于單一水價(jià)政策，但其減少幅度低于灌溉用水量的減少幅度和單方水效益的增加幅度，因此推薦兩部制水價(jià)政策。若水價(jià)需要上漲到拐點(diǎn)2.10元/m3以上，單一水價(jià)政策的單方水效益和農(nóng)戶收入要比兩部制水價(jià)政策更高，且其比兩部制水價(jià)政策多占的地、多用的水急劇減少，因此推薦單一水價(jià)政策。

3.2.3 差別水價(jià)情景

該情景在兩部制水價(jià)政策的基礎(chǔ)上實(shí)行差別水價(jià)，水價(jià)隨作物灌溉方式和作物種類而不同。將滴灌灌溉的水價(jià)優(yōu)惠50%，將傳統(tǒng)種植方式下糧食作物的水價(jià)上浮50%。

差別水價(jià)情景中，隨水價(jià)的上升，總種植面積、滴灌比例、農(nóng)戶收入、灌溉用水量和單方水效益的變化趨勢與其他2種情景基本一致，但變化拐點(diǎn)不同，且模擬得到的各指標(biāo)具體值與其他2種情景有所差別。

作物總種植面積和分項(xiàng)種植面積隨水價(jià)變化情況如圖6a所示。在該情景中，小麥和玉米的種植面積分別在水價(jià)上升至0.3元/m3和1.50元/m3時(shí)下降至穩(wěn)定不變，與前2種政策情景相比有所提前；而其余作物的種植面積以及總種植面積隨水價(jià)下降的穩(wěn)態(tài)拐點(diǎn)則出現(xiàn)在水價(jià)高達(dá)6.80元/m3時(shí)，滯后于前2種政策情景。這主要是由于，差別水價(jià)政策對采用傳統(tǒng)灌溉方式的糧食作物加價(jià)，使得糧食作物種植面積易隨水價(jià)上漲而快速壓縮；對采取滴灌的經(jīng)濟(jì)作物種植實(shí)行優(yōu)惠政策，使其對水價(jià)上漲的響應(yīng)更為遲緩。

差別水價(jià)情景中單方水效益、農(nóng)戶收入和灌溉用水量隨水價(jià)上升的變化情況如圖6b所示。在該情景中，單方水效益在水價(jià)上升至1.50元/m3時(shí)達(dá)到峰值14.00元/m3，相比于起漲點(diǎn)的上升幅度為4.2%；作物總灌溉用水量隨水價(jià)上升而下降，在水價(jià)為1.50元/m3時(shí)下降幅度達(dá)到14.9%，之后下降速度明顯后變慢；農(nóng)戶收入隨水價(jià)的提升的下降速度波動(dòng)不大，其在水價(jià)上升至1.50元/m3時(shí)，下降幅度為9.2%。

圖6 差別水價(jià)情景中各指標(biāo)值隨水價(jià)上升變化情況

結(jié)合圖6a和6b中各指標(biāo)的變化趨勢和拐點(diǎn)可以看出，在水價(jià)達(dá)到1.50元/m3之前，隨著水價(jià)的上升，雖然農(nóng)戶收入有所下降，但不及灌溉用水量的減少幅度，且單方水效益不斷上升，占地面積也不斷減小。而在水價(jià)上升至1.50元/m3之后，隨著水價(jià)的上升，農(nóng)戶收入和單方水效益均不斷減小，同時(shí)作物種植面積和灌溉用水量的下降速度顯著放緩。因此，在差別水價(jià)情景下，從單方水效益高、占地面積小、灌溉用水省、農(nóng)戶收入影響小四方面考慮，上漲后的水價(jià)應(yīng)處于0.24～1.50元/m3之間。

圖7比較了差別水價(jià)情景與其他2種水價(jià)政策情景的模擬結(jié)果。較之其他2種政策情景，當(dāng)水價(jià)少于2.10元/m3時(shí)，差別水價(jià)情景下總種植面積更小、灌溉用水更省、單方水效益更高；當(dāng)計(jì)量部分水價(jià)為2.10～6元/m3時(shí)，差別水價(jià)情景灌溉用水更多、總種植面積更大、單方水效益和農(nóng)戶收入更低；水價(jià)高于6元/m3后，差別水價(jià)情景的單方水效益更高、灌溉用水更省、農(nóng)戶收入和種植面積總體來看更高。

注：差值①表示差別水價(jià)減去兩部制水價(jià)得到的各指標(biāo)差值，差值②表示差別水價(jià)減去單一水價(jià)得到的各指標(biāo)差值

綜合3種政策情景下的分析結(jié)果，當(dāng)計(jì)量部分水價(jià)低于2.10元/m3或高于6元/m3時(shí)，差別水價(jià)政策優(yōu)于其他2種水價(jià)政策；當(dāng)計(jì)量部分水價(jià)在2.10～6元/m3之間時(shí)，差別水價(jià)政策在3種水價(jià)政策中最不理想，而單一水價(jià)政策最為適宜。就單方水效益而言，差別水價(jià)政策在計(jì)量部分水價(jià)為1.50元/m3時(shí)所獲得的單方水效益是所有情景中最高的。

4 結(jié) 論

針對多水源多灌溉方式的農(nóng)業(yè)灌區(qū)，運(yùn)用實(shí)證數(shù)學(xué)規(guī)劃模型，構(gòu)建了農(nóng)業(yè)水價(jià)政策效應(yīng)模擬方法，以水價(jià)改革試點(diǎn)區(qū)甘肅省民勤縣為例，進(jìn)行實(shí)證研究。分別模擬了單一水價(jià)、兩部制水價(jià)以及差別水價(jià)3種政策情景之下，研究區(qū)農(nóng)戶收入、單方水效益、作物種植面積、灌溉用水量和節(jié)水灌溉工程面積比例隨水價(jià)上升的變化趨勢，并在其基礎(chǔ)上得到不同水價(jià)政策下的水價(jià)上漲閾值以及各種政策的適用范圍。研究得到的主要結(jié)論包括：

1）利用2015年的數(shù)據(jù)在現(xiàn)狀水價(jià)政策下對模型進(jìn)行校準(zhǔn)，校準(zhǔn)后的模型基準(zhǔn)期最優(yōu)解與實(shí)際種植面積的差值絕對值不超過1.03×10-3hm2，模型可用于水價(jià)政策效應(yīng)的模擬。根據(jù)模擬結(jié)果可得，隨著水價(jià)的上升，3種水價(jià)政策情景下各效應(yīng)指標(biāo)變化趨勢相似：農(nóng)戶收入不斷減少，作物灌溉用水量和作物種植面積先減少后趨于穩(wěn)定，節(jié)水灌溉工程面積比例增加，單方水效益先增加后減少。

2）3種水價(jià)政策情景下，為應(yīng)對水價(jià)的上升，農(nóng)戶普遍傾向先減少糧食作物的種植面積，然后再壓縮經(jīng)濟(jì)作物的面積，還會通過增加作物節(jié)水灌溉面積比例來減少灌溉用水量。

3）在單一水價(jià)政策和兩部制水價(jià)政策情景下，上漲后的水價(jià)宜處于0.24～2.10元/m3之間；在差別水價(jià)政策情景下，上漲后的水價(jià)宜處于0.24～1.50元/m3之間。當(dāng)計(jì)量部分水價(jià)低于2.10元/m3或高于6元/m3后，差別水價(jià)政策的綜合效益最高；當(dāng)計(jì)量部分水價(jià)位于2.10～6元/m3之間時(shí)，單一水價(jià)政策最為理想；實(shí)行差別水價(jià)政策且計(jì)量部分水價(jià)為1.50元/m3時(shí)取得的單方水效益在所有情景中最高。

本文證明了水價(jià)在閾值內(nèi)穩(wěn)步上漲，不僅不會降低用水的綜合效益，反而可以提升單方用水效益。同時(shí)，證明了民勤縣現(xiàn)行的差別水價(jià)政策是合理的。民勤縣目前的水價(jià)遠(yuǎn)低于閾值，因此，在未來中短期內(nèi)，應(yīng)推進(jìn)差別水價(jià)的實(shí)行，計(jì)量部分水價(jià)應(yīng)在0.24～1.50元/m3之間逐步上升。在實(shí)施差別水價(jià)政策的同時(shí)，應(yīng)注意加大對節(jié)水灌溉工程的投入，也可考慮制訂配套的補(bǔ)貼政策對農(nóng)戶進(jìn)行適當(dāng)補(bǔ)償。

在實(shí)際應(yīng)用中，可在本研究結(jié)論的基礎(chǔ)上，廣泛開展公眾參與，并綜合考慮當(dāng)?shù)噩F(xiàn)有水價(jià)、農(nóng)民支付意愿、政府財(cái)政補(bǔ)貼、供水工程收入等多種因素和評判指標(biāo)得到最終的水價(jià)具體實(shí)施方案。本研究的方法可廣泛應(yīng)用于其他多水源、多灌溉方式灌區(qū)，為其水價(jià)政策的制定提供科學(xué)依據(jù)。

[1] 彭致功，張寶忠，劉鈺，等. 基于灌溉制度優(yōu)化和種植結(jié)構(gòu)調(diào)整的用水總量控制[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(3)：103－109. Peng Zhigong, Zhang Baozhong, Liu Wei, et al. Constraint of total water consumption amount based on optimized irrigation schedule and planting structure ajustment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(3): 103－109. (in Chinese with English abstract)

[2] 姜秋香，周智美，王子龍，等.基于水土資源耦合的水資源短缺風(fēng)險(xiǎn)評價(jià)及優(yōu)化[J].農(nóng)業(yè)工程學(xué)報(bào)，2017，33(12)：136－143. Jiang Qiuxiang, Zhou Zhimei, Wang Zilong, et al. Risk assessment and optimization of water resources shortage based on water and land resources coupling [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 136－143. (in Chinese with English abstract)

[3] James L D，Lee R R. 水資源規(guī)劃經(jīng)濟(jì)學(xué)[M]. 北京：水利電力出版社，1984：358.

[4] 周春應(yīng)，章仁俊. 農(nóng)業(yè)需水價(jià)格彈性分析模型[J]. 節(jié)水灌溉，2005(6)：24－26. Zhou Chunying, Zhang Renjun. Analytical model for price elasticity of agricultural water demand[J]. Water Saving Irrigation, 2005(6): 24－26. (in Chinese with English abstract)

[5] 王密俠，汪志農(nóng)，尚虎軍，等.陜西關(guān)中灌區(qū)水費(fèi)計(jì)收方式與水資源合理利用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2005，21(增刊1)：11－14. Wang Mixia, Wang Zhinong, Shang Hujun, et al. Water fee collecting methods and reasonable utilization of water resources in guanzhong irrigation district of shaanxi province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005, 21(Supp.1): 11－14. (in Chinese with English abstract)

[6] 賈大林，姜文來. 農(nóng)業(yè)水價(jià)改革是促進(jìn)節(jié)水農(nóng)業(yè)發(fā)展的動(dòng)力[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，1999，18(5)：4－7. Jia Dalin, Jiang Wenlai. Agricultural water price reform is the driving force for promoting the development of water-saving agriculture[J]. Journal of Agrotechnical Economics, 1995, 18(5): 4－7. (in Chinese with English abstract)

[7] 杜榮江，張鈞. 水資源浪費(fèi)的經(jīng)濟(jì)學(xué)分析與控制對策[J]. 河海大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，35(6)：722－726. Du Rongjiang, Zhang Jun. Economics analysis of waste of water resources and countermeasures for its control[J]. Journal of Hohai University: Natural Sciences, 2007, 35(6): 722－726. (in Chinese with English abstract)

[8] 李晨洋，于偉銘，陳正銳，等.考慮生態(tài)的灌區(qū)水資源區(qū)間兩階段隨機(jī)規(guī)劃模型建立與應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(21)：105－114. Li Chenyang, Yu Weiming, Chen Zhengrui, et al. Two-stage interval parameters water resources model considering ecology and application in irrigation district[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(21): 105－114. (in Chinese with English abstract)

[9] 粟曉玲，宋悅，劉俊民，等. 耦合地下水模擬的渠井灌區(qū)水資源時(shí)空優(yōu)化配置[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(13)：43－51. Su Xiaoling, Song Yue, Liu Junmin, et al. Spatiotemporal optimization allocation of water resources coupling groundwater simulation model in canal-well irrigation district[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(13): 43－51. (in Chinese with English abstract)

[10] 李然，田代貴. 農(nóng)業(yè)水價(jià)的困境擺脫與當(dāng)下因應(yīng)[J]. 改革，2016，32(9)：107－114. Li Ran, Tian Daigui. The dilemma of agricultural water price and its countermeasures[J]. Reform, 2016, 32(9): 107－114. (in Chinese with English abstract)

[11] 馬超. 農(nóng)業(yè)水價(jià)需求彈性調(diào)查及政策建議[J]. 水利發(fā)展研究，2018，18(4)：7－9，13. Ma Chao. Agricultural water price demand flexibility survey and policy suggestions[J]. Water Resources Development Research, 2018, 18(4): 7－9, 13. (in Chinese with English abstract)

[12] Espey M, Espey J, Shaw W D. Price elasticity of residential demand for water: A meta-analysis[J]. Water Resources Research, 1997(33): 1369－1374.

[13] 裴源生，方玲，羅琳. 黃河流域農(nóng)業(yè)需水價(jià)格彈性研究[J]. 資源科學(xué)，2003，25(6)：25－30. Pei Yuansheng, Fang Ling, Luo ling. Price elasticity of agricultural water demand in China[J]. Resources Science, 2003, 25(6): 25－30. (in Chinese with English abstract)

[14] 嚴(yán)冬，周建中. 水價(jià)改革及其相關(guān)因素的一般均衡分析[J]. 水利學(xué)報(bào)，2010，41(10)：1220－1227. Yan Dong, Zhou Jianzhong. A general equilibrium analysis of water price reform and related factors[J]. Journal of Hydraulic Engineering, 2010, 41(10): 1220－1227. (in Chinese with English abstract)

[15] 趙永，竇身堂，賴瑞勛. 基于靜態(tài)多區(qū)域CGE模型的黃河流域灌溉水價(jià)研究[J]. 自然資源學(xué)報(bào)，2015，30(3)：433－445. Zhao Yong, Dou Shentang, Lai Ruixun. Research on agriculture water pricing in the yellow river basin based on static multi-region CGE model[J]. Journal of Natural Resources, 2015, 30(3): 433－445. (in Chinese with English abstract)

[16] Bartolini F, Bazzani G M, Gallerani V, et al. The impact of water and agriculture policy scenarios on irrigated farming systems in Italy: An analysis based on farm level multi-attribute linear programming models[J]. Agricultural Systems, 2007(93): 90－114.

[17] Howitt R E. Positive mathematical programming [J]. American Journal of Agricultural Economics, 1995, 77(2): 329－342.

[18] 王裕雄，肖海峰. 實(shí)證數(shù)學(xué)規(guī)劃模型在農(nóng)業(yè)政策分析中的應(yīng)用：兼與計(jì)量經(jīng)濟(jì)學(xué)模型的比較[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2012，31(7)：15－21. Wang Yuxiong, Xiao Haifeng. The application of positive mathematial programming model in agricultural policy analysis: comparison with econometric model[J]. Journal of Agrotechnical Economics, 2012, 31(7): 15－21.(in Chinese with English abstract)

[19] Heckelei T, Britz W. Models based on positive mathematical programming: State of the art and further extensions [C]. 89th EAAE Seminar Paper, Parma, 2005：48-73.

[20] Merel P, Howitt R. Theory and application of positive mathematical programming in agriculture and the environment [J]. Annual Review of Resource Economics, 2014(6): 451－470.

[21] Iglesias E, Blanco M. New directions in water resources management: The role of water pricing policies [J]. Water Resources Research, 2008, 44(6): 1－11.

[22] Gallego-Ayala J, Gomez-Limon J A. Irrigation water pricing instruments: A sustainability assessment[J]. Spanish Journal of Agricultural Research, 2011, 9(4): 981－999.

[23] Gallego-Ayala J. Selecting irrigation water pricing alternatives using a multi-methodological approach[J]. Mathematical and Computer Modelling, 2012, 55(3/4): 861－883.

[24] Watto M A, Mugera A W. Irrigation water demand and implications for groundwater pricing in Pakistan[J]. Water Policy, 2016, 18(3): 565－585.

[25] Chebil A, Frija A, Thabet C. Irrigation water pricing between governmental policies and farmers' perception: Implications for green-houses horticultural production in Teboulba (Tunisia)[J]. Agricultural Economics Review, 2010, 11(2): 44－54.

[26] Franco-Crespo C, Vi?as J M S. The impact of pricing policies on irrigation water for agro-food farms in Ecuador[J]. Sustainability, 2017 9(9): 1515.

[27] Heady E O, Egbert A C. Regional programming of efficient agricultural production patterns[J]. Econometrica, 1964, 32(3): 374－386.

[28] Martin L J, Meister A D, Chen C C, et al. Quadratic programming models applied to agricultural policies[J]. American Journal of Agricultural Economics, 1979, 61(3): 586－587.

[29] Arfini F, Donati M, Solazzo R. Positive Mathematical Programming[M]. Roma, Italy: Research and Economics, 2016.

[30] 田聰穎，肖海峰. 目標(biāo)價(jià)格補(bǔ)貼與生產(chǎn)者補(bǔ)貼的比較:對我國大豆直補(bǔ)方式選擇的思考[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題，2018，39(12)：107－117. Tian Chongying, Xiao Haifeng. Comparison between target price subsidy and producer subsidy: Reflection on the selection of China's soybean direct subsidy[J]. Issues in Agricultural Economy, 2018, 39(12): 107－117. (in Chinese with English abstract)

[31] 中國營養(yǎng)學(xué)會. 中國居民膳食指南（2016）[M]. 北京：人民衛(wèi)生出版社，2016：269－287.

[32] Howitt R E. A Calibration method for agricultural economic production model[J]. Journal of Agricultural Economics, 1995, 46(2): 147－159.

[33] 謝臻，張鳳榮，王瀚巍，等. 基于節(jié)水灌溉技術(shù)的民勤綠洲土地利用空間布局和利用方式調(diào)整[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2017，22(10)：65－74. Xie Zheng, Zhang Fengrong, Wang Hanwei, et al. Adjustment of land use spatial structure and use pattern based on modern agriculture in Minqin oasis[J]. Journal of China Agricultural University, 2017, 22(10): 65－74. (in Chinese with English abstract)

[34] 柯珊珊，趙偉，楊晴. 水利工程現(xiàn)行水價(jià)分析[J]. 水利規(guī)劃與設(shè)計(jì)，2017，30(7)：161－165. Ke Shanshan, Zhao Wei, Yang Qing. Current water price analysis of water conservancy projects[J]. Water Resources Planning and Design, 2017, 30(7): 161－165. (in Chinese with English abstract)

[35] 王春艷，張博，沈大軍. 淠河總干渠供水成本分析計(jì)算及兩部制水價(jià)設(shè)計(jì)[J]. 水利經(jīng)濟(jì)，2012，30(5)：30－33. Wang Chunyan, Zhang Bo, Shen Dajun. Calculation of water supply costs and design of two-part water price for main canal of Pihe river[J]. Journal of Economics of Water Resources, 2012, 30(5): 30－33. (in Chinese with English abstract)

[36] 張麗霞. 民勤縣農(nóng)業(yè)水價(jià)綜合改革現(xiàn)狀分析及思考[J]. 南方農(nóng)業(yè)，2017，11(30)：61－65. Zhang Lixia. Analysis and thinking on the current situation of comprehensive reform of agricultural water price in Minqin county[J]. South China Agriculture. 2017, 11(30): 61－65. (in Chinese with English abstract)

[37] 民勤縣統(tǒng)計(jì)局民勤調(diào)查隊(duì). 民勤縣2015年國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計(jì)資料匯編[M]. 甘肅：民勤縣統(tǒng)計(jì)局，2016：110－120.

[38] 國家發(fā)展和改革委員會價(jià)格司. 全國農(nóng)產(chǎn)品成本收益資料匯編2016[M]. 北京：中國統(tǒng)計(jì)出版社，2016：21－81.

Simulation of effects of agricultural water price policy based on positive mathematical programming

Tan Qian, Wang Shuping, Zhang Tianyuan

(100083,)

Rational water price policy is of great significance for the efficient use of water resources. However, development processes of water price policy are mostly subjective and few methods are suitable to quantify policy effects. In this study, we focused on irrigation areas with multiple water sources and irrigation techniques. A method based on positive mathematical programming (PMP) was proposed to simulate the impacts of agricultural water price policy on economic benefits of water use, irrigation water consumption, planting land occupation and water-saving irrigation area. This method was applied to address a case in Minqin County, Gansu Province. Under the current water price policy of the study area, the model was calibrated with the actual data in 2015, the absolute value of the difference between the optimal solution of the calibrated model and the actual planting area in the base year didn’t exceed 1.3×10-3hm2, indicating that the model was reliable to simulate the effect of water price policy. Subsequently, the model was used to simulate the change trend of index with the rise of water price under the volumetric water price policy, the 2-part water price policy and the differential water price policy. The simulation results showed that, with the rise of water price, the change trend of index under 3 water price policy was similar:the farmer income reduced, the irrigation water consumption and planting area reduced first and then stabilized, the water-saving irrigation area increased, and the benefit per cubic meter of water increased first and then decreased. Simultaneously, in response to rising water price, farmers tended to reduce the planting area of food crops first and then the area of cash crops to alleviate the loss of income. Based on the simulation results, the water price threshold under different policies and the scope of application of various policies were further analyzed, and the following conclusions were drawn. Under the volumetric water price and the 2-part water price policy, the water price should be between 0.24 yuan/m3-2.10 yuan/m3; under the differential water price policy, the water price should be between 0.24-1.50 yuan/m3; when the metering part of the water price was lower than 2.10 yuan/m3or higher than 6 yuan/m3, the comprehensive benefit of the differential water price policy was the highest; when the metering part of the water price was between 2.10 yuan/m3and 6 yuan/m3, the volumetric water price policy was the most ideal; when the differential water price policy was implemented and the metering part of the water price was 1.50 yuan/m3, the benefit per cubic meter of water was the highest. This paper confirmed that the steady increase of water price would not reduce the comprehensive benefits of water use, but increase the benefit per cubic meter of water. The differential water price policy currently implemented in Minqin County was appropriate. In the future, in Minqin County, the water price under the differential water price policy should be increased within the range of 0.24-1.50 yuan/m3. Moreover, during the course of water price rising in the study area, special attentions should be paid to increase investment in water-saving irrigation projects and farmer income subsidy policies should be established under local contexts. Based on the conclusions of this study, various local factors and evaluation indexes such as current water price, farmers' willingness to pay, government financial subsidies, and income from water supply projects, could be comprehensively considered to obtain detailed implementation plan of water price policy. The method established in this paper was also applicable to other areas with similar problems.

agriculture; water; models; positive mathematical programming; water price policy; economic benefits of water use

2019-03-10

2019-07-10

國家自然科學(xué)基金資助項(xiàng)目（51822905、51779255）；國家千人計(jì)劃青年項(xiàng)目

譚 倩，教授，博士生導(dǎo)師，主要從事水資源水環(huán)境系統(tǒng)分析方面的研究。Email：qian_tan@cau.edu.cn

10.11975/j.issn.1002-6819.2019.16.018

S607+.1

A

1002-6819(2019)-16-0161-11

譚 倩，王淑萍，張?zhí)镦?基于實(shí)證數(shù)學(xué)規(guī)劃模型的農(nóng)業(yè)水價(jià)政策效應(yīng)模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(16)：161－171. doi：10.11975/j.issn.1002-6819.2019.16.018 http://www.tcsae.org

Tan Qian, Wang Shuping, Zhang Tianyuan. Simulation of effects of agricultural water price policy based on positive mathematical programming[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(16): 161－171. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.16.018 http://www.tcsae.org