陸泉志，陸桂軍，范稚蓮，羅明智，莫良玉

(1.廣西大學(xué)農(nóng)學(xué)院，南寧 530004；2.廣西壯族自治區(qū)科學(xué)技術(shù)情報研究所，南寧 530022； 3.廣西壯族自治區(qū)蠶業(yè)科學(xué)研究院，南寧 530007)

農(nóng)業(yè)生產(chǎn)引發(fā)的水環(huán)境污染問題始終是綠色發(fā)展行動中的重要議題。根據(jù)《廣西壯族自治區(qū)環(huán)境統(tǒng)計年報》的數(shù)據(jù)，2015年廣西農(nóng)業(yè)廢水中的氨氮排放量占全部廢水中氨氮排放量的32%，種植業(yè)廢水中的氨氮排放量占農(nóng)業(yè)廢水氨氮排放量的38%，化肥淋失是種植業(yè)水污染的主要因素。盡管2011-2015年廣西化肥使用量年均增長率為1.72%，增幅開始呈現(xiàn)緩和態(tài)勢，但是短期來看，要實現(xiàn)化肥使用量大幅減少的現(xiàn)實可行性較低，化肥污染排放的“約束紅線”始終存在。因此在研究種植業(yè)用水效率的時候不能只考慮經(jīng)濟產(chǎn)出，也要考慮化肥污染排放所帶來的負(fù)面效應(yīng)，這對保障廣西種植業(yè)水資源可持續(xù)利用，建設(shè)生態(tài)環(huán)保型綠色種植業(yè)有重要意義。

國內(nèi)學(xué)者大多基于數(shù)據(jù)包絡(luò)分析方法對用水效率進行測算。許多學(xué)者基于經(jīng)濟產(chǎn)出角度對不同地區(qū)的農(nóng)業(yè)用水效率進行了測度[1-5]，其考慮的是如何以最小的水資源投入實現(xiàn)更大的經(jīng)濟產(chǎn)出，忽視了環(huán)境約束對用水效率的影響，這可能導(dǎo)致真實的用水績效發(fā)生扭曲。為此，一些學(xué)者開始關(guān)注農(nóng)業(yè)生產(chǎn)給水體生態(tài)帶來的負(fù)面效應(yīng)，將水體污染指標(biāo)作為非期望產(chǎn)出納入全要素用水效率分析框架，在同時考慮經(jīng)濟效益與生態(tài)效益的基礎(chǔ)上采用非徑向的SBM-Undesirable模型修正以往研究中用水效率的偏差估計[6-8]。但這些學(xué)者直接采用整個投入產(chǎn)出系統(tǒng)的效率值作為農(nóng)業(yè)用水效率，并沒有考慮水資源投入的松弛量，這是存在一定偏差的。而一些學(xué)者借鑒Hu等[9]提出的全要素用水效率思想，通過構(gòu)建包含水資源投入的全要素農(nóng)業(yè)生產(chǎn)系統(tǒng)，再基于Hu等“最優(yōu)用水量與實際用水量比值”的用水效率定義將水資源投入從整個農(nóng)業(yè)生產(chǎn)系統(tǒng)中剝離出來，得到更為精確的用水效率值[1-3,10]。在用水效率的空間格局研究方面，大多學(xué)者采用空間自相關(guān)分析方法對用水效率的時空格局及其演化規(guī)律進行了分析[11-13]。

綜合以上研究，本文基于全要素用水效率思想，將種植業(yè)生產(chǎn)中的化肥污染排放作為非期望產(chǎn)出納入評價體系中，運用超效率SBM-Undesirable模型對廣西區(qū)域內(nèi)14個地市的種植業(yè)用水效率進行測算，以求更準(zhǔn)確地把握各地市在經(jīng)濟與環(huán)境制衡下的真實用水績效。在此基礎(chǔ)上，運用ESDA方法對廣西種植業(yè)用水效率的動態(tài)演進趨勢及空間相關(guān)性進行分析，以期深化對廣西區(qū)域內(nèi)種植業(yè)用水效率的認(rèn)識和理解，為“十三五”期間廣西種植業(yè)水資源可持續(xù)利用前景和綠色區(qū)域發(fā)展的相關(guān)對策提供借鑒。

1 研究方法與數(shù)據(jù)

1.1 超效率SBM-Undesirable模型

s->0，sb>0，sg>0，λ>0

i=1，2，…,m；r=1，2，…，q；j=1，2，…，n(j≠0)

1.2 種植業(yè)全要素用水效率測度方法

在超效率SBM-Undesirable模型的基礎(chǔ)上，借鑒Hu等[9]提出的全要素用水效率的理論，該理論以產(chǎn)出不變作為基本假定，以“最優(yōu)用水量與實際用水量比值”對種植業(yè)用水效率進行衡量，最優(yōu)用水量即在目前種植業(yè)產(chǎn)出及生產(chǎn)技術(shù)水平不變條件下，所需要的最少或最優(yōu)的農(nóng)田水資源投入。此方法將農(nóng)田灌溉用水從整個種植業(yè)投入產(chǎn)出系統(tǒng)中分離出來，以期更精準(zhǔn)地測度廣西種植業(yè)用水效率。全要素用水效率(Total Factor Water Efficiency，TFWE)的具體測度方法如下：

式中：最優(yōu)用水量為種植業(yè)實際用水量減去超效率SBM-Undesirable模型結(jié)果中農(nóng)田水資源投入量的徑向調(diào)整值和松弛調(diào)整值之和。

1.3 空間自相關(guān)分析方法

空間自相關(guān)分析是探索性空間數(shù)據(jù)分析方法(Exploratory Spatial Data Analysis，ESDA)的主要內(nèi)容，旨在研究樣本對象的空間分布及其相關(guān)性，揭示樣本數(shù)據(jù)的空間依賴性與空間異質(zhì)性的可視現(xiàn)象?？臻g自相關(guān)系數(shù)主要包括全局自相關(guān)系數(shù)與局部自相關(guān)系數(shù)。全局自相關(guān)系數(shù)( Global Moran’sI)是從全局的角度來測算單元值的相關(guān)性強度，通過全局自相關(guān)指數(shù)檢驗廣西地區(qū)內(nèi)所有地市的種植業(yè)用水效率值是否出現(xiàn)了空間聚集以及是否存在穩(wěn)健的空間溢出效應(yīng)；局部自相關(guān)系數(shù)(Local Moran’sIi)是用來測度局部地區(qū)或相鄰地區(qū)間觀測值的相關(guān)性強度，檢測是否存在局部地區(qū)高值或者低值，一般采用 Moran’s散點圖表示。計算公式如下：

式中：xi、xj表示樣本的屬性觀測值；n表示樣本數(shù)；Wij表示鄰接空間權(quán)重，本文采用鄰接標(biāo)準(zhǔn)來計算各市的空間權(quán)重，若城市i與城市j相鄰，Wij為1，否則Wij為0；I表示全局自相關(guān)系數(shù)，其系數(shù)取值范圍為[-1，1]，當(dāng)I>0 時，表示研究區(qū)域存在空間正相關(guān)，當(dāng)I<0 時，表示研究區(qū)域存在空間負(fù)相關(guān)，當(dāng)I=0時，表示不存在空間自相關(guān)，服從隨機分布。

在 Moran’s 散點圖中分為4個象限，高高集聚區(qū)(第1象限)和低低集聚區(qū)(第3象限)表示廣西區(qū)域種植業(yè)用水效率空間正相關(guān)性較強，即具有均質(zhì)性；低高集聚區(qū)(第2象限)和高低集聚區(qū)(第4象限)表示空間負(fù)相關(guān)性較強，即具有異質(zhì)性。

1.4 指標(biāo)體系構(gòu)建與數(shù)據(jù)來源

本文通過構(gòu)建包含水資源投入及化肥污染產(chǎn)出的種植業(yè)全要素投入產(chǎn)出系統(tǒng)來測度種植業(yè)用水效率。根據(jù)農(nóng)業(yè)經(jīng)濟學(xué)理論，農(nóng)業(yè)生產(chǎn)要素的投入主要包含土地、勞動力和資本等指標(biāo)。本文依照科學(xué)性、合理性和數(shù)據(jù)可獲得性的原則，在投入要素方面選取農(nóng)作物總播種面積、種植業(yè)勞動力投入、化肥折純量、農(nóng)田灌溉用水量、農(nóng)業(yè)機械總動力5個指標(biāo)。由于種植業(yè)勞動力數(shù)據(jù)缺失，根據(jù)戴俊[16]的研究，種植業(yè)勞動人數(shù)平均占農(nóng)村從業(yè)人員的59.38%。因此本文選取廣西農(nóng)村人口數(shù)的59.38%近似替代種植業(yè)勞動力數(shù)量。在產(chǎn)出指標(biāo)方面，由于化肥的施用會造成氨氮淋失，本文選取種植業(yè)廢水中的氨氮排放量作為非期望產(chǎn)出。選取種植業(yè)生產(chǎn)總值作為期望產(chǎn)出，并以2011年為基期按城市所在省份的CPI指數(shù)進行消脹處理。由于運用DEA模型時各項投入與產(chǎn)出之間必須滿足“同向性”假設(shè)，即投入量增加時產(chǎn)出不得減少。本文運用皮爾遜(Pearson)相關(guān)性檢驗方法對各項指標(biāo)進行檢驗，結(jié)果表明各投入項與產(chǎn)出項之間均能在1%的置信水平下通過雙尾檢驗且相關(guān)系數(shù)為正，各投入產(chǎn)出指標(biāo)符合“同向性”原則，具有合理性，可作進一步分析。

由于廣西環(huán)境保護廳僅公布了“十二五”期間的水污染排放數(shù)據(jù)，因此研究樣本數(shù)據(jù)跨度相應(yīng)受到限制。另一方面，研究“十二五”期間廣西各地區(qū)種植業(yè)用水效率時空分異演化特征，對“十三五”期間廣西種植業(yè)水資源管理工作的優(yōu)化改進具有重要的現(xiàn)實指導(dǎo)意義。因此，本文研究對象為廣西所轄的14個地市，時間段為2011-2015年，觀測值共70個，本文所選取的種植業(yè)廢水中的氨氮排放量數(shù)據(jù)來源于《廣西壯族自治區(qū)環(huán)境統(tǒng)計年報》(2011-2015年)；其他指標(biāo)數(shù)據(jù)均來源于《廣西壯族自治區(qū)統(tǒng)計年鑒》(2012-2016年)。

2 結(jié)果與討論

2.1 廣西種植業(yè)用水效率測算結(jié)果分析

運用MaxDEA.Ultra7.5軟件，基于投入導(dǎo)向的超效率SBM-Undesirable模型和全要素用水效率測度方法計算出2011-2015年廣西各地市種植業(yè)用水效率值，并以各地市種植業(yè)用水效率均值進行了排名，具體結(jié)果見表1。

表1 2011-2015年種植業(yè)用水效率及排名Tab.1 Water efficiency and ranking of planting industry from 2011 to 2015

由表1可知，從廣西總體用水效率均值來看，2011-2015年廣西歷年種植業(yè)用水效率均值均超過1，整體并未表現(xiàn)出明顯的增長或降低趨勢，而是呈現(xiàn)小幅波動的均衡走勢。但需要注意的是，廣西歷年用水效率平均值較高主要是由于個別地市的超效率值過高引起的，效率均值大于1的城市只有3個，剩余的11個地市均未達到有效前沿面。本文將超效率均值過高的防城港市和崇左市去除后，發(fā)現(xiàn)研究期內(nèi)廣西種植業(yè)用水效率均值為0.83，距離有效前沿面還有一定距離，廣西種植業(yè)還有較大的節(jié)水潛力和用水效率提升空間。從各地市的效率均值及排名來看，防城港、崇左和南寧的用水效率均值最高，呈“三足鼎立”的態(tài)勢。但在研究期內(nèi)這3個地市的用水效率呈現(xiàn)不斷回落走勢，表明其近年來未能較好地平衡經(jīng)濟產(chǎn)出與污染排放之間的關(guān)系，污染約束壓力有所增強。貴港、河池、來賓、北海的用水效率均值都較低，一方面主要是因為其水資源投入的冗余量(可節(jié)約用水量)較大，粗放式灌溉用水是主要原因，也從側(cè)面反映了這些地市在種植業(yè)水資源管理方面表現(xiàn)不佳；另一方面這些地市的經(jīng)濟產(chǎn)出松弛量都為0，但污染排放的松弛調(diào)整值較大，這也反映了化肥過量投入導(dǎo)致的污染排放在一定程度上降低了用水效率。

根據(jù)表1，利用Geoda1.6.7軟件的分位功能繪制了各地市平均用水效率的地理分布圖(見圖1)，可以直觀地看出，污染排放約束下的廣西種植業(yè)用水效率在地理分布上呈現(xiàn)明顯的空間非均衡性特征。由圖1可知，廣西種植業(yè)用水效率較高的地區(qū)是南寧、防城港、崇左、百色、梧州，其中有4個地市位于廣西西南部且在地理空間上相鄰，呈現(xiàn)一定的空間集聚格局。南寧、防城港、崇左同屬于北部灣經(jīng)濟區(qū)，受經(jīng)濟區(qū)政策支持因素及首府經(jīng)濟圈“空間溢出效應(yīng)”影響，這些地區(qū)在經(jīng)濟、科技、資源等方面交互融合程度較高，社會協(xié)同發(fā)展能力較強，其經(jīng)濟發(fā)展水平與用水效率較高密切相關(guān)。王瑩[17]和王潔萍[18]等的研究結(jié)果也印證了社會經(jīng)濟發(fā)展水平對農(nóng)業(yè)水資源利用效率具有顯著的正向影響這一觀點。用水效率值也較高的百色屬于干季易旱地區(qū)，可能是由于常年在水資源稀缺的高壓約束下反而強化了農(nóng)民的節(jié)水意識，更重視水資源的高效利用，這與孫才志等[11]的研究結(jié)論相類似。廣西種植業(yè)用水效率較低的地區(qū)是來賓、賀州、河池、貴港等地，同樣這些地區(qū)在空間分布上趨于相鄰，表現(xiàn)出一定的地帶性特征。因此從廣西種植業(yè)用水效率的區(qū)域分布情況來看，高效率值地區(qū)與低效率值地區(qū)可能存在一定程度的空間關(guān)聯(lián)性，可運用空間自相關(guān)分析方法對廣西種植業(yè)用水效率的空間自相關(guān)性和空間異質(zhì)性作進一步探討。

圖1 廣西種植業(yè)用水效率的地理分布Fig.1 Geographical distribution map of water use efficiency in Guangxi

2.2 廣西種植業(yè)用水效率空間相關(guān)性分析

2.2.1 廣西種植業(yè)用水效率全局空間自相關(guān)分析

以不同年份的廣西種植業(yè)用水效率為統(tǒng)計變量，以地市為基礎(chǔ)單元，運用Geoda1.6.7軟件，全局 Moran’sI指數(shù)的計算結(jié)果如表2所示。由表2可知，在正態(tài)分布假設(shè)下，研究期內(nèi)所有年份的全局 Moran’sI指數(shù)均為正值，均通過5%水平顯著性檢驗，說明廣西種植業(yè)用水效率具有顯著的空間正相關(guān)關(guān)系，存在效率高值地市與效率高值地市相鄰、效率低值地市與效率低值地市相鄰的空間集聚格局。從動態(tài)分析可知，研究期內(nèi)全局 Moran’sI指數(shù)呈逐年下降趨勢，表明廣西種植業(yè)用水效率呈現(xiàn)分散態(tài)勢，空間集聚效應(yīng)隨著時間的推移在不斷減弱。按照此趨勢，“十三五”期間廣西種植業(yè)用水效率空間集聚效應(yīng)有可能會消失。表2中的全局 Moran’sI指數(shù)無法揭示具體是哪些區(qū)域產(chǎn)生了空間集聚效應(yīng)，無法對局部區(qū)域空間自相關(guān)關(guān)系進行刻畫，因此需進一步檢驗廣西種植業(yè)用水效率的局部空間相關(guān)性。

表2 廣西種植業(yè)用水效率全局 Moran’s I指數(shù)Tab.2 Global Moran’s I of planting industry water use efficiency in Guangxi

注：**和***分別表示在5%和1%的水平上顯著。

2.2.2 廣西種植業(yè)用水效率局部空間自相關(guān)分析

在全局自相關(guān)分析的基礎(chǔ)上，本文采用 Moran’s散點圖(見圖2)對2011年、2013年、2015年這3個代表性年份進行局部空間自相關(guān)分析，檢驗在統(tǒng)計顯著下，是否存在局部地區(qū)高值或者低值。

圖2 廣西種植業(yè)用水效率的 Moran’s散點圖Fig.2 Moran’s index scatter plot of planting industry water use efficiency in Guangxi

由圖2可知，Moran’sIi均大于0，3個代表性年份中大部分地市都處于第1象限與第3象限中，屬于高值簇和低值簇的現(xiàn)象，其空間均質(zhì)性顯著，表現(xiàn)出較強的正向空間自相關(guān)特征。從動態(tài)演變來看，廣西14個地市種植業(yè)用水效率的空間關(guān)聯(lián)類型變化較小，僅有1個地市在2015年移動到了第4象限，一定程度上反映了局部空間集聚現(xiàn)象存在減弱趨勢。具體來看，種植業(yè)用水效率穩(wěn)定存在高高(HH)集聚的地區(qū)是南寧、崇左、防城港，表明這些地區(qū)自身和周邊地區(qū)用水效率都較高，其空間溢出效應(yīng)顯著可歸因于這些地區(qū)水資源豐富，政治、經(jīng)濟地理位置優(yōu)越，北部灣經(jīng)濟區(qū)社會綜合發(fā)展能力較強，防城港、崇左與首府南寧在農(nóng)業(yè)科技創(chuàng)新、資本投資、農(nóng)產(chǎn)品貿(mào)易、社會教育等方面形成聯(lián)動效應(yīng)，提升了整個地區(qū)的用水效率。種植業(yè)用水效率穩(wěn)定存在低高(LH)集聚的地區(qū)是百色、欽州，表明這些地區(qū)自身用水效率較低而周邊地區(qū)較高，空間異質(zhì)性突出，高值簇對百色與欽州種植業(yè)用水效率的增長拉動能力較弱，未能對其產(chǎn)生空間溢出效應(yīng)。種植業(yè)用水效率穩(wěn)定存在低低(LL)集聚的地區(qū)是桂林、柳州、賀州、來賓、河池、貴港、梧州、玉林、北海，這些地區(qū)占有的比例較大，可以說是廣西種植業(yè)用水效率整體水平的縮影。與周邊西部省份一樣，廣西大部分地區(qū)社會發(fā)展綜合水平較低，大量水利基礎(chǔ)設(shè)施存在年久失修老化嚴(yán)重的問題，季節(jié)性缺水和工程性缺水并存，加之喀斯特地貌廣布，地表水流失較為嚴(yán)重，用水環(huán)境并不是十分理想[19]。另一方面，由于農(nóng)民教育功能缺失，其環(huán)保意識普遍不強，較難平衡好經(jīng)濟產(chǎn)出與化肥污染排放之間的關(guān)系，研究期內(nèi)廣西化肥施用量年均增幅1.72%，高于同期全國平均水平，其污染排放也抑制了用水效率的提升。一般來說，由于低值地區(qū)自身缺乏活力，規(guī)模性低值集聚將加劇用水環(huán)境惡化并易于形成惡性循環(huán)累積，最終對該區(qū)用水效率的快速發(fā)展起到抑制作用。在種植業(yè)用水效率高低(HL)集聚區(qū)方面，僅在2015年出現(xiàn)了梧州市，高低(HL)集聚表明梧州自身用水效率較高而周邊地區(qū)較低，說明梧州市開始打破“冷點屏障”，種植業(yè)用水效率出現(xiàn)回升。這可能是由于梧州與較為發(fā)達的廣東地區(qū)地緣鄰接，通過與廣東地區(qū)在資源、科技、貿(mào)易等方面交互融合從而不斷集聚能量，產(chǎn)生跨區(qū)域協(xié)同提升效應(yīng)。但梧州自身用水效率增長并未對周圍地區(qū)產(chǎn)生正向的空間交互影響，呈現(xiàn)一定程度的極化現(xiàn)象。

3 結(jié) 論

綜上所述，可以得出以下結(jié)論。

(1)2011-2015年，污染排放約束下的廣西各地市種植業(yè)用水效率呈現(xiàn)明顯的非均衡特征，種植業(yè)用水效率整體不高，只有南寧、防城港、崇左達到有效前沿面，形成“三足鼎立”態(tài)勢。

(2)研究期內(nèi)，污染排放約束下的廣西種植業(yè)用水效率呈現(xiàn)顯著正向的全局空間自相關(guān)性，用水效率相似值(高高或低低)呈現(xiàn)空間集聚特征，但其空間溢出效應(yīng)逐年減弱。不同地市種植業(yè)用水效率的空間關(guān)聯(lián)類型變化較小，廣西有一半以上地市穩(wěn)定存在低低(LL)集聚區(qū)；高高(HH)集聚區(qū)呈現(xiàn)由南寧、防城港、崇左組成的三核心熱點集聚模式。

4 建 議

(1)用水效率與農(nóng)業(yè)經(jīng)濟發(fā)展水平密切相關(guān)，在全要素用水效率理論下，應(yīng)構(gòu)建與水資源相匹配的種植業(yè)生產(chǎn)布局，優(yōu)化各投入要素的資源配置，不斷提高種植業(yè)經(jīng)濟產(chǎn)出，縮小各地區(qū)農(nóng)業(yè)經(jīng)濟發(fā)展水平的差距，為進一步提升用水效率打下堅實經(jīng)濟基礎(chǔ)。此外需加大對農(nóng)田水利基礎(chǔ)設(shè)施與節(jié)水技術(shù)創(chuàng)新的資金投入，對用水效率低低集聚區(qū)采取適度的政策傾斜，尤其是農(nóng)民教育與科技政策，以促進低效率地區(qū)對高效率地區(qū)產(chǎn)生追趕效應(yīng)。

(2)提倡種植業(yè)綠色生態(tài)用水理念，穩(wěn)步推進化肥使用量零增長行動。加強綠色發(fā)展理念的宣傳力度，提高農(nóng)民生態(tài)意識與節(jié)水意識，在經(jīng)濟產(chǎn)出最大化和環(huán)境污染最小化的杠桿上保持平衡。同時要切實完成好農(nóng)業(yè)部要求的“到2020年化肥使用量零增長”任務(wù)。要因地制宜、循序漸進地優(yōu)化施肥結(jié)構(gòu)，通過政策與科技的配合，大力推廣有機肥替代化肥、高效新型肥料、水肥一體化技術(shù)，提高肥料利用效率，減少肥料揮發(fā)淋失，使其對水資源污染最小化，切實保障經(jīng)濟與環(huán)境良性協(xié)同發(fā)展。

(3)準(zhǔn)確把握廣西種植業(yè)用水效率的空間差異性，加強政府宏觀引導(dǎo)，實施區(qū)域協(xié)同發(fā)展戰(zhàn)略，促進各地區(qū)在節(jié)水科技、水資源管理、流域監(jiān)管等多方面交流融合，相互借鑒學(xué)習(xí)以實現(xiàn)資源優(yōu)勢互補，從而產(chǎn)生空間聯(lián)動效應(yīng)。也可通過高效率地區(qū)對口幫扶、資本融合、經(jīng)驗指導(dǎo)等方式對低效率集聚區(qū)實施“富鄰”政策，創(chuàng)造更多的空間溢出通道以拉動其用水效率增長。此外還應(yīng)針對廣西種植業(yè)用水效率空間自相關(guān)性逐年減弱的現(xiàn)象，深入研究其影響因素，以強化其空間溢出效應(yīng)。

□