吳浩力，吳方力，黃 愉

(1.陜西省水利發(fā)展調(diào)查與引漢濟渭工程協(xié)調(diào)中心，陜西 西安 710004；2.陜西省地下水管理中心，陜西 西安 710003)

1 引言

截止到2017年末，全國耕地面積134881.2千hm2，耕地灌溉面積68272千hm2，耕地灌溉率僅為50.6%。陜西省的耕地面積3982.9千hm2，耕地灌溉面積1274.99千hm2，耕地灌溉率32.0%，明顯低于全國水平。

灌區(qū)是我國經(jīng)濟社會發(fā)展的重大公益性基礎(chǔ)設(shè)施，是維系國家糧食安全、水生態(tài)健康、水土流失和風沙治理的重要紐帶，同時是城市和工業(yè)以及生態(tài)環(huán)境供水的重要載體，也是山水林田湖草系統(tǒng)治理和鄉(xiāng)村振興的重要支撐。因此，如何高效、合理的評價灌區(qū)運管水平對灌區(qū)可持續(xù)發(fā)展具有重要意義。本文主要使用數(shù)據(jù)包絡(luò)分析法對陜西省值五大灌區(qū)進行評價，為現(xiàn)代化灌區(qū)評價提供借鑒。

2 數(shù)據(jù)包絡(luò)分析法(DEA，Data Envelop ment Analysis)

數(shù)據(jù)包絡(luò)分析法是運籌學、管理科學和數(shù)理經(jīng)濟學交叉研究的一個新的領(lǐng)域，也稱“平衡基準”。其第一個模型CCR由運籌學家Charnes，Cooper等人于1978年提出[1]。該方法以凸分析、線性規(guī)劃和對策論為基礎(chǔ)，使用輸入、輸出數(shù)據(jù)建立非參數(shù)的經(jīng)濟數(shù)學模型，評價具有多個輸入和輸出的，擁有相同目標的運營單位或“決策單元”的相對效率。其主要特點一是不依賴過多條件假設(shè)，簡化了指標間的關(guān)系分析；二是采用綜合指標，不必對數(shù)據(jù)進行規(guī)范化處理，實現(xiàn)了無量綱化；三是不需要人為干預，由模型自動確定最佳權(quán)重，實現(xiàn)了客觀分析；四是相關(guān)軟件成熟，借助電算操作便捷。

2.1 決策單元(DMU，Decision Making Units)

任何一個生產(chǎn)過程或企事業(yè)單位(或部門)都涉及投入和產(chǎn)出兩項基本活動，而從投入到產(chǎn)出的過程是通過一系列決策來實現(xiàn)的，因此，將這個生產(chǎn)過程或生產(chǎn)單位稱為一個決策單元DMU(即評價對象)。通常，這些DMU應(yīng)為同種類型或同屬于某個領(lǐng)域，具有類似的投入和產(chǎn)出，并且，可以認為這些活動的期望都是向著資源的高效利用，實現(xiàn)帕累托最優(yōu)狀態(tài)。以灌區(qū)為例，假設(shè)對A、B、C三個灌區(qū)進行相對有效性評價，則這每個灌區(qū)就可以看作是一個DMU。

2.2 生產(chǎn)可能集及其公理

設(shè)每個DMU在生產(chǎn)活動中有m項投入，n項產(chǎn)出，寫成向量形式分別為，；于是便可以用來表示該DMU的整個生產(chǎn)活動，并稱集合S為該生產(chǎn)活動的生產(chǎn)可能集。它反映的是多投入和產(chǎn)出條件下，兩者之間的技術(shù)關(guān)系(即經(jīng)濟學中的生產(chǎn)函數(shù))。由于DEA最早是應(yīng)用在基于收入和產(chǎn)出的經(jīng)濟學分析，為了合理描述生產(chǎn)活動的內(nèi)在邏輯關(guān)系，從中引入一些公理體系，包括平凡公理、凸性公理、無效性公理、錐性公理等[2-3]。在滿足上述部分公理的基礎(chǔ)上，又引入了三個取值0或1的參數(shù)。生產(chǎn)可能集S的一般形式如式(1)所示：

(1)

2.3 生產(chǎn)前沿面與DEA有效性的判斷

生產(chǎn)前沿面(又稱包絡(luò)面)是經(jīng)濟學中生產(chǎn)函數(shù)向多產(chǎn)出情況的推廣，有關(guān)其研究始于1957年，由經(jīng)濟學家Machaell Farrell基于線性規(guī)劃模型的生產(chǎn)效率測度思想所構(gòu)建[4]。國內(nèi)魏權(quán)齡、王金祥等人[5-6]對該理論有過系統(tǒng)深入的論述。簡而言之，生產(chǎn)前沿面指從DMU中觀測輸入和輸出數(shù)據(jù)包絡(luò)面的有效部分，可以理解為帕累托最優(yōu)點構(gòu)成的面。而DEA方法的本質(zhì)就是判斷DMU是否位于生產(chǎn)可能集的前沿面上。

落在前沿面上的DMU被認為具有相對有效性，稱“DEA有效”。確切的說應(yīng)是“DEA弱有效”，因為DEA有效還需要滿足兩項權(quán)重系數(shù)u、v大于0的條件。相關(guān)數(shù)學表達如式(2)所示，其中hjo表示模型的相對效率值，簡稱效率指數(shù)，其經(jīng)濟學含義為相對生產(chǎn)率，根據(jù)相關(guān)性質(zhì)[7]，其與量綱選取無關(guān)，其值等于1即表示DMU位于生產(chǎn)前沿面上，此時即為DEA弱有效。相反的，沒有在前沿面上的DMU是相對無效率的，其值介于0～1之間；u、v表示權(quán)重系數(shù)，其并非由主觀的價值取向得出，而是由模型本身求得的最優(yōu)權(quán)重，這就較好的排除了主觀因素的影響，pj、qj分別表示投入和產(chǎn)出。

(2)

2.4 技術(shù)有效、規(guī)模收益和綜合技術(shù)效率

技術(shù)有效又稱純技術(shù)效率，反映的是DMU在最優(yōu)規(guī)模時投入要素的生產(chǎn)效率，可通過產(chǎn)出/投入的比值來測度，通常受技術(shù)和管理等因素制約。規(guī)模效率反映的是實際規(guī)模與最優(yōu)生產(chǎn)規(guī)模的差距，主要受規(guī)模因素影響。綜合技術(shù)效率是對DMU的資源配置能力、資源使用效率等多方面能力的綜合衡量與評價。通常認為，綜合技術(shù)效率=技術(shù)有效×規(guī)模效率。綜合技術(shù)效率=1，表示該決策單元的投入產(chǎn)出是綜合有效的，即同時技術(shù)有效和規(guī)模有效，當采用CCR模型時，就認為DMU同時為技術(shù)有效和規(guī)模有效。純技術(shù)效率=1，表示在目前的技術(shù)水平上，其投入資源的使用是有效率的，未能達到綜合有效的根本原因在于其規(guī)模無效，因此，其改革的重點在于如何更好地發(fā)揮其規(guī)模效益。

2.5 DEA方法原理及其經(jīng)典模型

DEA方法核心思想是使生產(chǎn)率越大越好。經(jīng)典的模型方法是基于DMU的投入、產(chǎn)出構(gòu)造生產(chǎn)函數(shù)和不失一般性的約束條件(為分式規(guī)劃)，利用分式規(guī)劃CC-變換[8]將模型中生產(chǎn)函數(shù)轉(zhuǎn)化為的線性規(guī)劃的目標函數(shù)，同相關(guān)約束條件一起合稱為原線性規(guī)劃(PCCR)，為了降低運算量，再利用線性規(guī)劃理論對其進行變換，將其轉(zhuǎn)化為對偶線性規(guī)劃(DCCR)，易知對偶規(guī)劃與原規(guī)劃是等價的，再求DCCR最優(yōu)解的過程。由于生產(chǎn)活動一般都伴隨著多項投入和產(chǎn)出，按照模型的要求，需要將這些投入和產(chǎn)出合成單一的投入和產(chǎn)出，而這種合成，就需要借助權(quán)重向量，這里權(quán)重向量正好揭示出DMU的重要性排序。利用約束條件求解線性規(guī)劃的最優(yōu)解即是計算權(quán)重向量。方法流程圖見圖1。

圖1 DEA方法流程圖

DEA方法自誕生以來，受到眾多學者的追捧，并針對不同類型的決策單元開發(fā)了大量模型，如CCR模型、BCC模型和交叉效率模型等。其中，最早同時也是普適度最高的當屬CCR模型，包括經(jīng)典CCR模型和含非Archimedes無窮小量的CCR模型。非Archimedes無窮小量是一個小于任何正數(shù)且大于零的數(shù)，在實際應(yīng)用中通常量級為1×10-5，過大有可能造成目標函數(shù)無邊界，過小則可能使對偶模型無解[9]。有關(guān)其引入模型的合理性仍存在爭議，相關(guān)文獻可參考[10]。

ε→0 Λε>0

CCR模型假定DMU的規(guī)模收益是不變的，也即是說投入量和產(chǎn)出量是成正比例關(guān)系，本文既采用該模型。其以投入為導向的分式規(guī)劃模型如(3)所示：

(3)

式中：vT、uT分別表示DMU輸入和輸出的權(quán)重系數(shù)向量；I表示為投入導向模型。

CCR的原線性規(guī)劃模型PCRR見式(4)：

(4)

CCR的對偶線性規(guī)劃模型DCRR見式(5)：

(5)

式中：θ、λ為DCRR的最優(yōu)解。θ表示DMU距離包絡(luò)面的徑向距離；λj可將各個有效點連起來形成包絡(luò)面。

(6)

式中：s+、s-分別為松弛變量與剩余變量，其作用是使包絡(luò)面沿著水平和垂直方向延伸。

2.6 陜西省直屬五大灌區(qū)基本情況

陜西省直屬灌區(qū)有寶雞峽引渭、涇惠渠、交口抽渭、石頭河水庫和桃曲坡水庫五大灌區(qū)，均為公益二類事業(yè)單位。五大灌區(qū)總設(shè)計灌溉面積429.06千hm2，總有效灌溉面積377.56千hm2，占全省灌區(qū)總有效灌溉面積的42.9%，總實際灌溉面積286.45千hm2，占自身有效灌溉面積的75.87%，累計節(jié)水灌溉面積205.09千hm2，占全省節(jié)水灌溉面積21.24%。涉及總灌溉人口492.01萬人，分別占全省總?cè)丝诤完P(guān)中地區(qū)總?cè)丝诘?2.69%和20.27%。詳細情況見圖2、圖3。

圖2 省直五大灌區(qū)灌溉面積情況

圖3 全省灌區(qū)有效灌溉面積占比情況

2.7 評價指標的選取

大型灌區(qū)的綜合運行狀況評價是一項復雜的系統(tǒng)工程，涉及項目目標評價、設(shè)施狀況和管理水平評價、項目效果和效益評價、資源及環(huán)境可持續(xù)性評價、社會評價等諸多方面。因此，需要大量的指標和數(shù)據(jù)的支撐，這些數(shù)據(jù)來源于不同的利益相關(guān)者，搜集難度較大，加之有些數(shù)據(jù)不全面、口徑不一致或更新不及時，不能客觀的反映評價對象的真實狀況。

本文經(jīng)過對省直五大灌區(qū)的現(xiàn)有數(shù)據(jù)資料的分析，指標體系整理結(jié)果見表1，共計28項。在此基礎(chǔ)上清除無效數(shù)據(jù)，見表2。

表1 灌區(qū)綜合評價指標體系

表2 指標數(shù)據(jù)

2.8 計算分析評價結(jié)果

Cooper[11]認為，DMU數(shù)量應(yīng)滿足式(7)的要求，而本例中只有5個DMU，故最多只應(yīng)有一項投入指標和一項產(chǎn)出指標，但若只選2項指標結(jié)果會以偏概全，使評價結(jié)果置信度低，為此，本文借助主成分分析法，使上述指標在降維的同時又保證了信息丟失最小化。降維后的投入和產(chǎn)出見表3、圖4。對于單投入和單產(chǎn)出的DEA模型則可以直接用產(chǎn)出比投入來表示：

圖4 五大灌區(qū)投入產(chǎn)出對數(shù)圖

表3 降維后指標

t≥max{m×n,(m+n)×3}

(7)

式中：t為DMU的個數(shù)，m、n分別為投入和產(chǎn)出指標的個數(shù)。

通過上述分析，可知目前交口抽渭的相對有效性最高，也即是說，在相同的生產(chǎn)條件下，其產(chǎn)出率最高，后面依次為涇惠渠、石頭河、寶雞峽和桃曲坡，作為五大灌區(qū)之首的寶雞峽，其規(guī)模最大，但其產(chǎn)出與其投入并不匹配，五大灌區(qū)產(chǎn)出投入比平均水平為9.99，寶雞峽為4.2，占平均水平的42%，和交口、涇惠渠還有一定差距，需要加強設(shè)施維護和提高管理水平；石頭河和桃曲坡排名靠后是由于其主要給城市供水，農(nóng)業(yè)供水相對薄弱，同時存在水資源短缺和灌區(qū)節(jié)水措施有待改進的問題。

3 結(jié)語

由于本文的分析主要是側(cè)重于農(nóng)水方向，同時也缺乏必要的數(shù)據(jù)，如綜合經(jīng)營情況、制度管理等，所以排名只反映灌區(qū)在支撐農(nóng)業(yè)方面的相對有效性。如石頭河雖然排名靠后，主要是石頭河主要是向城市供水，農(nóng)灌較少的緣故。要全面公允的評價五大灌區(qū)的綜合運行情況還需更多詳實數(shù)據(jù)資料，但這不影響DEA方法評價的有效性；再者，實例中的評價單元只有5個，對于DEA方法而言，往往評價單元越多，所構(gòu)建的前沿面就越精細，計算結(jié)果置信度就越高；此外，最后，DEA方法還有很多種模型，本文只拋磚引玉采用CCR模型予以評價，某些過程和假設(shè)或與現(xiàn)實情況有出入，值得進一步探索挖掘。