車國偉,王澤民,楊正麗
(1.武漢大學,湖北 武漢 430079;2.天津市測繪院,天津 300381;3.四川大學水利水電學院,四川 成都 610065; 4.四川大學水力學與山區河流開發保護國家重點實驗室,四川 成都 610065)
自1987年可持續發展的概念正式提出以來[1],可持續發展的理論與實踐研究受到各國政府和學術界的關注[2-4]。國內外有較多關于區域可持續發展評價的研究,比較常見的有DEA評價模型及其拓展模型(如CCR、BCC、兩階段DEA等),雖然較其他綜合評價方法具有很多無法替代的優勢,但是在評價決策單元的效率時仍然存在諸多弊端,如非完全排序問題、“自評自利”效應和權重方案缺陷。鑒于傳統DEA評價模型存在的弊端,國內外專家學者逐漸開始對傳統的DEA模型進行拓展及完善,其中最具代表性的研究成果之一就是DEA交叉效率評價方法[5-7]。
雖然交叉效率方法能夠綜合考慮自評與他評體系使得各個決策單元的測評結果可相互對比,還能得到所有決策單元的完全排序評價結果,但是,交叉效率模型也存在一些短板,如由于某些決策單元之間可能存在多組最優權重使最終交叉效率評價值不一致。正是這種評價結果的不唯一性嚴重制約了DEA交叉效率模型的穩定性和在現實生活中的實用性。因此,本文針對傳統的DEA交叉效率模型的缺陷,提出綜合協調交叉效率模型,并以我國31個省、自治區和直轄市(不含港、澳、臺)的數據為例,進行可持續發展效率評價。
眾所周知,數據包絡分析是目前測算被評價單元效率最基本、最有效的方法[8-10],但是DEA的方法通常只能對截面數據及時間序列進行評價,而區域發展不是一個靜止的狀態,而是一個連續的、長期的變化過程,即DEA測算方法缺乏對決策單元的面板數據評價處理能力。因此本文引用Malmquist指數評價模型,根據區域可持續發展指標的面板數據評價區域可持續發展效率變化情況。
Malmquist最早是在1953被提出,在1982年,Caves等將Malmquist指數方法用于評估生產效率變化情況。后由F?re等于1994年建立了基于DEA的Malmquist指數方法[11],將數據包絡分析方法與非參數線性規劃法結合起來,用以測算多投入產出角度的生產率變化情況。該方法可以用來比較兩個連續時期的全要素生產效率的增長情況,彌補了種種靜態數據包絡分析模型在對面板數據評估方面的缺陷,因此,Malmquist指數模型得到了學術界、產業界廣泛的認可和應用。
可持續發展問題可歸結為復雜的系統問題,故本文從區域可持續發展的多維度框架出發,將可持續發展的具體表現形式分為4個板塊:生態環境可持續、經濟發展可持續、社會福祉可持續及科技創新可持續[12-14]。
為了能夠描述與評價各個區域的環境可持續性、經濟可持續性、社會可持續性和科技可持續性,需要尋求具有以上意義、滿足上述要求的指標;同時,應從區域可持續發展系統結構出發,兼顧科學性、完備性、靈敏性、可靠性、動態性、協調性等。本文評價指標體系設計具體如下:首先針對本文研究的問題及遵循指標遴選的各項原則,理論遴選出40個區域可持續發展評價指標,在此基礎上,根據區域可持續發展的“四維交互結構”理論框架,從隸屬度、鑒別力和相關性3個方面對指標進行多重實證篩選,最后對指標體系的有效性進行檢驗。最終,保證本文所建立的指標體系是建立在科學創新的理論基礎之上,且經過多重的客觀定量梳理和檢驗得到的一套科學、有效的指標體系,主要包括19個指標(見表1)。

表1 區域可持續發展評價指標體系

據此,本文將上文篩選后的指標體系進行投入與產出的分類,詳見表1。需要注意的是其中節能環保支出、社會保障財政支出、R&D經費內部支出、科學技術財政支出、高技術產業R&D人員為逆向投入指標,即在區域可持續發展過程中,投入越多的財政支出進行節能環保事業、社會保障及科學技術研發越符合可持續發展的基本準則,人口密度為產出的逆向指標。
本文測評31個省域的可持續發展效率情況,因此即有31個決策單元(n=31),其中每個決策單元(DMUj,j=1,2,3,…,n)均有11個投入指標和8個產出指標(m=11,s=8)。即投入和產出變量分別為:xj=(x1j,x2j,x3j,…x11j),yj=(y1j,y2j,y3j,…y8j),則Edd代表被評價單元EMUd的效率值,通過DEA的CCR模型進行線性規劃可得
(1)
則任意一個“區”可表示為
(2)

(3)
同時,31個區域可持續發展評價的平衡型交叉效率評價模型為
(4)

(5)

(6)

(7)
將Malmquist指數引入區域可持續發展評價,不僅可以考慮各區域可持續發展的技術進步變化(TEch)情況,而且能夠在此基礎上,將區域可持續發展的相對效率指數(EFFch)進一步分解為區域可持續發展純效率變化(PEch)和區域可持續發展規模效率變化(SEch),從而得出區域可持續全要素發展效率變化(TFPch)情況。Malmquist指數是通過對兩個連續時期的TFP取幾何平均值得到的,具體如下
TFPk=Effchk×Techk=(Pechk×Sechk)×Techk
(8)

(9)
對式(9)求幾何平均值即可得到某區域從第t年到t+1年的綜合Malmquist指數,如式(10)所示,并且經過RD分解可得到其詳細求解過程
(10)
本文數據來源于2009—2014年間的《中國統計年鑒》《中國高技術產業統計年鑒》《中國環境統計年鑒》《中國能源統計年鑒》中相關統計數據。通過上述統計數據,得到2008—2013年我國31個省、自治區和直轄市(不含港、澳、臺)區域可持續發展相關指標的統計值,即2009—2014年區域可持續發展的面板數據。面板數據也稱為平行數據,面板數據包含時間與截面2個維度,本文分析的是31個區域的可持續發展情況,則本文的數據具有樣本、指標及時間3個維度,利用該三維數據,既可以將每年不同的省域的可持續發展效率值進行比較(橫向比較),同時又可以將每個區域2008—2013年不同時期的可持續發展和演化動態(縱向比較)。即通過橫向縱向雙向剖析及深刻研究區域可持續發展的效率和演化情況。
利用Matlab R2015b軟件對傳統的DEA交叉效率評價方法和綜合協調評價方法進行公式編程,將指標統計數據代入,評價出2008—2013年的可持續發展效率。具體的區域可持續發展評價的詳細結果和分區域結果對比分析見下文。
圖1反映的是用綜合交叉效率集結評價模型得出的我國31個省、自治區和直轄市的區域可持續發展情況。從圖1可以看出,綜合協調交叉效率評價模型評價的結果波動相對明顯,說明在考慮自評與他評相結合及指標權重的分配等問題后,本文構建的區域可持續發展的綜合協調交叉效率評價模型的評價功能更為細致和精確。
本文采用DEAP version2.1進行計算,采用上文提到的2009—2014年區域可持續發展的面板數據,可以算出第k個區域在第t年到第t+1年的Malmquist指數值,可以衡量每個區域連續的、變化的發展狀況。根據計算,假設某一區域在2008—2009年的可持續發展Malmquist指數大于1,則表明該區域的可持續綜合發展水平提高;若小于1,則表明該區域的可持續綜合發展水平降低;等于1則表明區域可持續綜合發展水平不變。結果表明,除個別區域外,大部分區域可持續發展Malmquist指數一直大于1,表明可持續發展效率一直在平緩增長,如圖2所示。

圖2 2008—2013年區域可持續發展綜合交叉效率Malmquist指數
從圖3(a)可以看出我國東部沿海各區域最近的連續6年的區域可持續發展情況和走勢。不難發現,從北至南的東部10省黑龍江、吉林、遼寧、山東、江蘇、上海、浙江、福建、廣東、海南區域可持續發展呈現一定的規律。在以往價值理念中,區域可持續發展效率評估值可能會很高的上海、浙江、福建等區域的區域可持續發展效率并未達到人們“預期”的理想高度。如圖3(a)所示,廣東及海南的可持續發展效率值一直處于0.9~1,接近1或等于1,明顯高出其他幾個區域,浙江、上海、江蘇、福建、山東地區則比較一致地徘徊在0.8附近;黑龍江、吉林、遼寧較其他東部沿海省域相對得分較低,處于及格線的邊緣。根據圖3(b)顯示,東部沿海六省區域的可持續發展Malmquist指數一直大于1,說明其可持續發展效率一直在平緩增長。

圖3 東部沿海10省2008—2013年區域可持續發展綜合交叉效率評價和Malmquist指數
針對基于數據包絡評價模型在自評體系下的揚長避短、完全優劣排序、極端權重等問題,本文提出了綜合協調交叉效率評價模型。該模型包含19個評價指標,其中11個投入、8個產出,通過放松一般的交叉效率模型的平均假設性,克服平均交叉效率模型的多解性和不穩定性等問題,從而能夠減少區域可持續發展評價的加權投入與加權產出之間的離差。以我國31個省、自治區和直轄市的區域可持續發展為例,結果表明,考慮自評與他評相結合及指標權重的分配等問題后,本文構建的區域可持續發展的綜合協調交叉效率評價模型的評價功能更為細致和精確。同時,構建基于DEA的區域可持續發展Malmquist指數模型,結果表明,全國可持續發展Malmquist指數大部分大于1,說明其可持續發展效率一直在平緩增長。以我國東部沿海10省區域為例,有效地評價了其可持續發展效率和動態變化情況。