疏勒河灌區(qū)水資源分配與種植結(jié)構(gòu)調(diào)整研究

虎海燕

(甘肅省疏勒河流域水資源管理局，甘肅 玉門 735211)

1 概況

1.1 水資源現(xiàn)狀

甘肅省疏勒河流域水資源管理局(以下簡稱省疏管局)轄昌馬、雙塔、花海三個灌區(qū)管理處和一個水庫電站管理處，疏勒河流域干流全長670km，流域面積17萬km2[1]。流域內(nèi)土地資源豐富，現(xiàn)有灌溉面積134.42萬畝，多年平均來水量10.31億m3，地表水不重復部分0.73億m3，流域內(nèi)水資源總量為11.55億m3。工農(nóng)業(yè)生產(chǎn)總供水量10.27億m3，其中：河道引水8.35億m3，占81%，提取地下水1.92億m3，占19%。2012年流域總耗水量14.37億m3，其中農(nóng)田灌溉耗水量12.07億m3，林牧漁畜耗水量1.24億m3，工業(yè)耗水量0.33億m3，城鎮(zhèn)生活耗水量0.09億m3，農(nóng)村生活耗水量0.06億m3，生態(tài)環(huán)境耗水量0.58億m3[2]。根據(jù)流域供需水實際情況，說明疏勒河流域農(nóng)田灌溉用水明顯偏高，流域工業(yè)及生活用水比例明顯偏低，水資源嚴重缺乏。

1.2 種植結(jié)構(gòu)情況

疏勒河流域2006—2016年調(diào)查數(shù)據(jù)顯示，主要種植作物為經(jīng)濟作物，其次為糧食作物，最后為林草地。其中：糧食作物主要種植在昌馬灌區(qū)，占總量的82.49%；經(jīng)濟作物主要種植在昌馬灌區(qū)和雙塔灌區(qū)，占總量的47.87%和43.44%；林草地在昌馬灌區(qū)、雙塔灌區(qū)、花海灌區(qū)種植占比均值分別為61.82%、30.22%和8.30%。在糧食作物中，三大灌區(qū)均以小麥和玉米為主要作物；在經(jīng)濟作物中，雙塔灌區(qū)以瓜菜和棉花為主導作物，兩者之和超過90%，其中棉花的占比均超過70%，昌馬和花海灌區(qū)經(jīng)濟作物以大麥、瓜菜、孜然、棉花和其它為主，其它類中主要包含食葵、洋蔥、茴香、枸杞、花卉、藥材等品種。

按照水資源總量分配情況和作物種植結(jié)構(gòu)，疏勒河流域農(nóng)業(yè)用水占總用水量的約90%左右，灌區(qū)內(nèi)種植作物多以高耗水作物為主，并且水分利用效率、綜合效益均比較低，嚴重制約了灌區(qū)農(nóng)業(yè)生產(chǎn)發(fā)展和水分效益提升。因此，為了改變畸形的用水結(jié)構(gòu)和流域生態(tài)惡化局面，亟需開展基于水資源分配的作物種植結(jié)構(gòu)模型優(yōu)化研究，調(diào)整流域內(nèi)現(xiàn)狀種植結(jié)構(gòu)，實現(xiàn)經(jīng)濟、社會、生態(tài)環(huán)境的協(xié)調(diào)可持續(xù)發(fā)展，將一定的水資源合理地分配到各種作物之間，使灌區(qū)獲得最大效益是灌溉管理者目前研究的主要課題之一。

2 構(gòu)建種植結(jié)構(gòu)模型

2.1 建立目標函數(shù)

本次研究所建立的種植結(jié)構(gòu)模型，基于水資源優(yōu)化配置，以保證灌區(qū)糧食產(chǎn)量、生態(tài)平衡及經(jīng)濟效益最大化為目標，綜合考慮了經(jīng)濟效益、社會效益和生態(tài)效益[3- 6]，目標函數(shù)為：

(1)

約束條件：

式中，N—選入目標函數(shù)中的種植作物種類；Aj—第j種糧食作物和經(jīng)濟作物的種植面積，畝；pj—第j種糧食作物和經(jīng)濟作物產(chǎn)品的單價，元/kg；yj—第j種糧食作物和經(jīng)濟作物灌溉條件下產(chǎn)量，kg/畝；Yjsp—根據(jù)社會人口、經(jīng)濟條件和人民生活水平而設(shè)定的灌區(qū)糧食產(chǎn)量閾值，kg。A3—林草地的種植面積，畝；A3sp—保障生態(tài)系統(tǒng)平衡而設(shè)定的林草地種植面積閾值，畝。

2.2 確定模型參數(shù)

2.2.1 決策變量

根據(jù)種植結(jié)構(gòu)優(yōu)化模型目標函數(shù)，模型的決策變量為糧食作物、經(jīng)濟作物和林草地的種植面積，模型函數(shù)中還包括糧食作物產(chǎn)量閾值、林草地種植面積閾值、疏勒河流域種植總面積閾值等參數(shù)[7- 8]，見表1。

表1 疏勒河流域種植結(jié)構(gòu)優(yōu)化決策變量(一)

當優(yōu)化模型以糧食作物和經(jīng)濟作物的子類種植面積作為優(yōu)化目的時，決策變量由選定的作物種類決定，疏勒河流域現(xiàn)狀糧食作物主要以小麥和玉米為主，被調(diào)查的經(jīng)濟作物包括油料、啤酒大麥、籽瓜、啤酒花、甜菜、瓜菜、果園、棉花、孜然和其他，根據(jù)前面對疏勒河主要作物的綜合篩選，經(jīng)濟作物主要以棉花、蜜瓜、洋蔥等為主。見表2。

2.2.2 閾值等參數(shù)確定

根據(jù)疏勒河流域2006—2016年糧食作物產(chǎn)量、草地種植面積和總種植面積調(diào)查數(shù)據(jù)資料，分別對并進行甄別篩選，修正數(shù)據(jù)之后繪制繪制了疏勒河流域2006—2016各參數(shù)數(shù)據(jù)圖，如圖1—3所示。

表2 疏勒河流域種植結(jié)構(gòu)優(yōu)化決策變量(二)

圖1 疏勒河糧食作物產(chǎn)量/萬kg

圖2 疏勒河林草地種植面積/萬畝

圖3 疏勒河種植面積/萬畝

根據(jù)圖1—3中數(shù)據(jù)趨勢可知近年來糧食作物產(chǎn)量、林草地面積和種植總面積呈現(xiàn)出一定的趨勢，并有一定相關(guān)性。為保障灌區(qū)糧食產(chǎn)量和生態(tài)平衡，實現(xiàn)灌區(qū)糧食安全和社會穩(wěn)定，在種植結(jié)構(gòu)調(diào)整中應設(shè)定相應閾值，從圖中數(shù)據(jù)資料信息可知，可將糧食作物閾值設(shè)定為最小值Yjsp=8000萬kg和多年均值9737萬kg，林草地種植面積閾值A(chǔ)3sp閾值設(shè)定為多年均值，即A3sp=13.00萬畝，疏勒河流域種植面積圖中數(shù)據(jù)趨勢可知近年來總體平穩(wěn)，即A=130萬畝。

3 模型優(yōu)化設(shè)計與分析

3.1 優(yōu)化模型

在確定了優(yōu)化模型的目標函數(shù)、約束條件和相關(guān)參數(shù)之后，可知決定優(yōu)化模型最終結(jié)果的是畝均產(chǎn)量和市場價格，本次優(yōu)化通過對作物產(chǎn)量和市場價格的歷史數(shù)據(jù)進行分析，并結(jié)合上節(jié)約束條件和相關(guān)參數(shù)進行模型優(yōu)化設(shè)計，如圖4—7所示。

圖4 疏勒河糧食作物畝均產(chǎn)量/kg

圖5 疏勒河糧食作物平均價格/(元/kg)

圖6 疏勒河經(jīng)濟作物畝均產(chǎn)量/kg

圖7 疏勒河經(jīng)濟作物平均價格/(元/kg)

根據(jù)上述圖標分析，以畝均產(chǎn)量和平均價格作為變化因子，設(shè)定模型優(yōu)化方案。根據(jù)糧食、經(jīng)濟作物產(chǎn)量曲線可知，產(chǎn)量變幅整體呈現(xiàn)較小且增長趨勢，因此選定2012年為現(xiàn)狀年，在現(xiàn)狀年的基礎(chǔ)上認為產(chǎn)量仍然呈增長趨勢，近期2015年產(chǎn)量增加以5%計算，遠期2020年產(chǎn)量增加以10%計算。根據(jù)糧食、經(jīng)濟作物平均價格的曲線趨勢和現(xiàn)狀分析可知，平均價格的波動較大，且整體呈現(xiàn)出周期趨勢，因此以2009—2013年五年內(nèi)的最大和最小值代表現(xiàn)狀可能會出現(xiàn)的市場價格極值，未來的平均價格范圍為現(xiàn)狀極值的110%，并對糧食作物和經(jīng)濟作物的價格極值進行組合，共生成四種未來可能的市場平均價格，見表3—4。

表3 產(chǎn)量預測方案

表4 平均價格預測方案

3.2 模型優(yōu)化分析

根據(jù)產(chǎn)量和平均價格的未來預測變化，并以Yjsp=8000萬kg和9737萬kg，A3sp=13.00萬畝，A=130萬畝，共可設(shè)計16組優(yōu)化方案，對優(yōu)化模型和擬定預測方案利用規(guī)劃求解算法進行求解，得到種植面積和總產(chǎn)值，并對方案優(yōu)化結(jié)果進行分析總結(jié)，給出產(chǎn)量和價格變化情況下的疏勒河流域種植結(jié)構(gòu)優(yōu)化方案。

根據(jù)糧食產(chǎn)量閾值選擇的不同可以分為兩個情景：

情景一：Yjsp=8000萬kg，A3sp=13.00萬畝，A=130萬畝；

情景二：Yjsp=9737萬kg，A3sp=13.00萬畝，A=130萬畝。

在上述兩情景下，得出的各方案優(yōu)化后產(chǎn)值結(jié)果見表5—6。

表5 情景一下各優(yōu)化模型方案產(chǎn)值單位：億元

表6 情景二下各優(yōu)化模型方案產(chǎn)值單位：億元

根據(jù)上節(jié)產(chǎn)量及價格預測方案可以得到畝均產(chǎn)值預測方案，見表7，由表中數(shù)據(jù)可知，在不同方案設(shè)計下糧食作物的畝均產(chǎn)值均小于經(jīng)濟作物的畝均產(chǎn)值，因此在以總產(chǎn)值最大為目標的優(yōu)化方案設(shè)計和計算過程中，計算結(jié)果由畝均產(chǎn)值作為約束條件來決定。

表7 畝均產(chǎn)值預測方案 單位：元/畝

根據(jù)對優(yōu)化計算結(jié)果的分析總結(jié)可知，優(yōu)化結(jié)果根據(jù)糧食畝均產(chǎn)量和糧食產(chǎn)量閾值組合可以分為4類(此時A3sp=13.00萬畝，A=130萬畝)，見表8。

表8 畝均產(chǎn)量與產(chǎn)量閾值組合方案

根據(jù)對上述4種方案的總結(jié)，最終得到的種植結(jié)構(gòu)優(yōu)化結(jié)果，見表9。

表9 2015年種植結(jié)構(gòu)優(yōu)化結(jié)果表單位：萬畝

從表9的計算結(jié)果可知，糧食作物的種植面積約占13%～16%，經(jīng)濟作物的種植面積占總面積的74%～77%，林草地的種植面積在約束條件的邊界達到，占總種植面積的10%。糧食作物、經(jīng)濟作物和林草地的種植面積由糧食作物的畝均產(chǎn)量和糧食產(chǎn)量的閾值確定。由于經(jīng)濟作物的畝均產(chǎn)值高于糧食作物的畝均產(chǎn)值，因此經(jīng)濟作物仍然是疏勒河流域的主導作物，通過培育新品種，改善種植方式等提高糧食作物和經(jīng)濟作物的畝均產(chǎn)量，對于保障疏勒河糧食安全、提高流域農(nóng)業(yè)總產(chǎn)值及保障生態(tài)安全具有重要意義。

4 結(jié)論與討論

本文基于水資源合理分配，開展了農(nóng)業(yè)種植結(jié)構(gòu)優(yōu)化方案研究，以糧食畝均產(chǎn)量和糧食、經(jīng)濟作物單價作為變量進行模型優(yōu)化計算，結(jié)果顯示產(chǎn)值總量由市場價格和畝均產(chǎn)量共同決定，但市場價格作用更為明顯，并且優(yōu)化結(jié)果與糧食作物畝均產(chǎn)量、糧食作物產(chǎn)量和林草地種植面積閾值有關(guān)，且最優(yōu)解均為滿足約束條件邊界處，且計算分析，糧食作物種植面積應占13%～16%，經(jīng)濟作物的種植面積占總面積的74%～77%，林草地的種植面積占總種植面積的10%，這個種植比例既能保證水資源合理分配，又能實現(xiàn)產(chǎn)量總值最大化。目前，灌區(qū)作物種植面積和種類在很大程度上受主客觀因素的影響，在一定程度上還不能按照優(yōu)化模型來生產(chǎn)，所有還需要進一步調(diào)查，形成更加易于推廣的模型研究。