趙俊威，張惠琴，曾 晨，陳安航

（成都理工大學(xué)管理科學(xué)學(xué)院，四川 成都 610059）

0 引言

生態(tài)文明建設(shè)關(guān)乎人類未來，建設(shè)綠色家園是人類共同的夢想，為此，習(xí)近平總書記強(qiáng)調(diào)：“‘十四五’時期，我國生態(tài)文明建設(shè)進(jìn)入了實(shí)現(xiàn)生態(tài)環(huán)境質(zhì)量改善由量變到質(zhì)變的關(guān)鍵時期”。在此階段，水資源的高效利用便顯得至關(guān)重要。隨著經(jīng)濟(jì)社會的發(fā)展，我國水資源問題日益嚴(yán)重，一方面，人均水資源占有量不足，時空分布不均；另一方面，水污染嚴(yán)重，水資源利用效率低下使得水資源形勢進(jìn)一步惡化。因此如何高效可持續(xù)地利用水資源成為社會各界關(guān)注的重點(diǎn)，而探究水資源綠色效率的影響因素不僅是提高水資源利用效率直接有效的選擇，而且對協(xié)同碳排放有效控制和環(huán)境高水平保護(hù)具有重要意義。

近年來，關(guān)于水資源綠色效率影響因素的研究，國內(nèi)外學(xué)者均取得了一定的研究成果。如魏楚［1］等通過對文獻(xiàn)的集成式梳理和歸納，認(rèn)為技術(shù)水平、水資源價格等是影響水效率的主要因素；許朗［2］等利用Tobit模型考察了農(nóng)業(yè)收入占比和用水成本等因素對農(nóng)業(yè)灌溉用水效率的影響；孫思奧［3］等采用LMDI（Logarithmic Mean Divisia Index）模型發(fā)現(xiàn)人口與人均GDP是黃河流域用水量增加的主導(dǎo)因素；Shang［4］等研究表明技術(shù)進(jìn)步降低了天津市的工業(yè)用水量，而工業(yè)規(guī)模則會增加工業(yè)用水量；這些學(xué)者們的研究從不同角度探討了水資源綠色效率的影響因素，但均未考慮各地理單元之間的聯(lián)系，隨著空間計量分析技術(shù)的逐漸發(fā)展，學(xué)者們愈發(fā)看重空間因素對水資源綠色效率的影響，秦騰［5］、張峰［6］以及Zhao［7］等利用空間計量模型對水資源綠色效率進(jìn)行研究，發(fā)現(xiàn)不同地區(qū)水資源利用效率的驅(qū)動因素存在較大差異；由此可知水資源綠色效率是諸多因素共同作用的結(jié)果，且在不同的時間和空間下各因素的作用力大小和方向都不盡相同。鑒于此，地理加權(quán)回歸模型（Geographically weighted regression，GWR）由于可以識別空間非平穩(wěn)性的局部變異系數(shù)而得到了廣泛應(yīng)用［8］，孫才志［9］等基于GWR 模型證明了我國水資源綠色效率空間分布不均衡并且存在著顯著的空間正相關(guān)性。綜合來看，GWR模型能夠針對不同地區(qū)得出差異化的研究結(jié)論，有效克服了地理單元間的空間異質(zhì)性，但不足在于只能對截面數(shù)據(jù)進(jìn)行回歸分析，當(dāng)需要估計的參數(shù)過多時，將會極大損失參數(shù)估計的精確度［10］。于是Huang［11］等將時間效應(yīng)引入到GWR 模型中，提出了時空地理加權(quán)回歸模型（Geographically and Temporally Weighted Regression，GTWR），該模型能在時間和空間上對不同時空單元的參數(shù)變異情況進(jìn)行捕捉，有效彌補(bǔ)了GWR 模型的不足，在城市群發(fā)展規(guī)劃［12］、物流業(yè)發(fā)展［13］以及環(huán)境保護(hù)［14］等諸多背景下得到了廣泛應(yīng)用。

綜上所述，同時考慮時間因素和空間因素的GTWR 模型具備優(yōu)良的統(tǒng)計性質(zhì)［15］，將其用于分析水資源綠色效率的影響因素，能夠更好地厘清各因素的作用機(jī)制及強(qiáng)度，準(zhǔn)確揭示區(qū)域水資源綠色效率的發(fā)展差異。但目前國內(nèi)外利用GTWR 模型來分析水資源綠色效率的相關(guān)文獻(xiàn)較少，而且在已有的流域探究中，大多以長江、黃河為例，鮮有研究關(guān)注沱江流域水資源利用情況。沱江作為長江的一級支流，不僅是四川省經(jīng)濟(jì)產(chǎn)出的重要支撐，同時也是長江上游的綠色屏障，沱江流域的高質(zhì)量發(fā)展關(guān)乎到四川省乃至整個西南地區(qū)人民綠色生活的健康。有鑒于此，本文以沱江流域?yàn)槔?，運(yùn)用超效率SBM（Slacks-Based Measure）模型測度水資源綠色效率水平，并結(jié)合GTWR模型從時空異質(zhì)性視角探究水資源綠色效率的影響因素，進(jìn)而推進(jìn)不同空間范圍內(nèi)水資源綠色效率的治理機(jī)制研究，為我國區(qū)域水資源綠色效率治理提供有效依據(jù)和科學(xué)參考。

1 研究區(qū)域概況與方法

1.1 研究區(qū)域概況

沱江是四川省腹部地區(qū)的重要河流之一，河流全長627.4 km，是四川省工業(yè)和人口密度最高的區(qū)域，承擔(dān)了四川省內(nèi)絕大多數(shù)居民的日常生活用水和工業(yè)、農(nóng)業(yè)發(fā)展用水。受地質(zhì)、地貌、氣候等因素影響，沱江流域降水量和徑流量季節(jié)差異明顯，降水量主要集中在每年的6-9 月，全年降水量呈北部多、東部少的分布趨勢。由此可見，沱江流域的水資源分布存在顯著的時空差異。隨著城市化集聚擴(kuò)張和工業(yè)增速的加快，沱江流域不同區(qū)域的水資源綠色效率也呈現(xiàn)出明顯的差異。因此，以沱江流域?yàn)槔芯繀^(qū)域水資源綠色效率的影響因素極具代表性。

由沱江流域高質(zhì)量發(fā)展研究中心（四川省社會科學(xué)重點(diǎn)研究基地）可知，全流域流經(jīng)四川省10 個市（州），36 個縣（區(qū)），如表1 所示。需要說明的是，本文中沱江流域與從水系分區(qū)角度的沱江流域并不一致。為保證數(shù)據(jù)來源的可靠性以及水資源綠色效率測度的有效性，本文將水資源綠色效率測度的對象選定為除阿壩藏族羌族自治州外的9 個地級市，選取2007-2017年的面板數(shù)據(jù)，利用時空地理加權(quán)回歸模型探究不同時間點(diǎn)、不同區(qū)域的水資源綠色效率的分異性，因地制宜地制定水資源管理措施，促進(jìn)流域內(nèi)社會、經(jīng)濟(jì)以及環(huán)境的可持續(xù)發(fā)展。

表1 四川境內(nèi)各行政流域面積統(tǒng)計表Tab.1 Statistical table of the area of each administrative watershed in Sichuan

1.2 研究方法

1.2.1 超效率SBM模型

綠色發(fā)展的本質(zhì)在于減少環(huán)境污染，降低資源消耗。依照魏楚［1］、孫才志［9］等對水資源效率的界定，結(jié)合實(shí)際發(fā)展情況及綠色發(fā)展的理念，水資源綠色效率是指水資源等生產(chǎn)要素投入和帶來的經(jīng)濟(jì)、社會和生態(tài)環(huán)境的產(chǎn)出的比率。對水資源綠色效率的測度是一個復(fù)雜的投入產(chǎn)出過程，勞動、資本、能源等生產(chǎn)要素的投入不僅產(chǎn)生期望的經(jīng)濟(jì)效益，也會帶來非期望的環(huán)境損害。超效率SBM 模型不僅能恰當(dāng)?shù)靥幚砹朔瞧谕a(chǎn)出，而且消除了松弛變量對測度值的影響［16］，使得水資源利用效率的測度值更加準(zhǔn)確。具體公式如下所示：

式中：pSE為效率值；x和y表示投入和產(chǎn)出的要素；m和s表示投入指標(biāo)和產(chǎn)出指標(biāo)的個數(shù)；i和r代表投入和產(chǎn)出的決策單元；s+和s-分別代表投入和產(chǎn)出的松弛量；λ為權(quán)重向量。pSE≥1時，代表決策單元相對有效；pSE≤1時，代表決策單元相對無效。

1.2.2 時空地理加權(quán)回歸模型

GTWR 模型嵌入了研究數(shù)據(jù)的時空特性，不僅能夠較好地體現(xiàn)各驅(qū)動要素的時空差異性，而且使得模型的參數(shù)估計和統(tǒng)計檢驗(yàn)結(jié)果都更加顯著。此外，相較于一般的回歸模型，GTWR模型還可利用ArcGIS 等軟件將各樣本點(diǎn)的回歸參數(shù)在空間上可視化地呈現(xiàn)，使得模型結(jié)果的表現(xiàn)更為直觀。因此，本文引入了GTWR模型，模型結(jié)構(gòu)如下：

式中：yi為第i個樣本點(diǎn)的因變量；xik為第k個自變量在第i個樣本點(diǎn)處的觀測值；n為樣本點(diǎn)的個數(shù)；(ui，vi，ti)是第i個樣本點(diǎn)的時空坐標(biāo)位置；是第k個自變量在第i個樣本點(diǎn)處的回歸系數(shù)；β0(ui，vi，ti)是第i個樣本點(diǎn)的時空截距；εi為殘差。

2 水資源綠色效率測度及時空差異性

2.1 水資源綠色效率指標(biāo)體系的構(gòu)建及數(shù)據(jù)處理

水資源綠色效率測度指標(biāo)體系的構(gòu)建要綜合考慮各方面因素。本文在柯布道格拉斯生產(chǎn)函數(shù)的基礎(chǔ)之上，結(jié)合目前已有的相關(guān)研究［17-23］以及沱江流域區(qū)域特性，構(gòu)建出如下的水資源綠色效率評價體系。

表2 水資源綠色效率評價體系Tab.2 Evaluation System of Green Efficiency of Water Resources

本文選取2007-2017 年沱江流域9 個地級市（受數(shù)據(jù)缺失影響，不含阿壩藏族羌族自治州）的水資源相關(guān)數(shù)據(jù)，各項(xiàng)數(shù)據(jù)主要來源于2007-2017年《中國城市統(tǒng)計年鑒》、《四川省統(tǒng)計年鑒》、《中國農(nóng)村統(tǒng)計年鑒》、《全國污染源普查：農(nóng)業(yè)污染源系數(shù)手冊》以及各地級市水資源公報。具體指標(biāo)說明如下：

（1）勞動力投入：用年均從業(yè)人數(shù)來衡量實(shí)際勞動力投入。

（2）資本投入：本文選擇資本存量作為資本投入。采用張軍［21］等計算資本存量的方法，使用永續(xù)盤存法估算沱江流域9個地級市的資本存量。

（3）水資源投入：本文選擇總用水量作為水資源投入。其中總用水量為生活用水總量、農(nóng)業(yè)用水總量和工業(yè)用水總量之和。

（4）期望產(chǎn)出。本文以綠色發(fā)展和可持續(xù)發(fā)展理念為出發(fā)點(diǎn)，結(jié)合沱江流域的區(qū)域特性，將經(jīng)濟(jì)效益和社會效益納入指標(biāo)體系來衡量沱江流域的期望產(chǎn)出。以沱江流域各地級市生產(chǎn)總值來衡量經(jīng)濟(jì)效益，引入社會環(huán)境指數(shù)（Socio-Environmen?tal Index，SEI）來衡量社會和環(huán)境效益。其中，社會環(huán)境指數(shù)在結(jié)合現(xiàn)有研究的基礎(chǔ)上參考中國統(tǒng)計學(xué)會提出的綜合發(fā)展指數(shù)［22，23］，從產(chǎn)業(yè)發(fā)展、政府環(huán)境管制力度以及科技創(chuàng)新3個層面選取7 個指標(biāo)來衡量。通過AHP 專家打分法，對產(chǎn)業(yè)發(fā)展、政府環(huán)境管制、科技創(chuàng)新3個層進(jìn)行賦權(quán)，并運(yùn)用主成分分析法對7個指標(biāo)進(jìn)行降維，從而計算得到社會環(huán)境指數(shù)，具體測算過程如下：

①產(chǎn)業(yè)發(fā)展（B1）。

② 政府環(huán)境管制（B2）。

③科技創(chuàng)新（B3）。

（5）非期望產(chǎn)出。參考已有的研究，多數(shù)學(xué)者采用生活污水排放總量和工業(yè)污水排放總量之和來衡量非期望產(chǎn)出［24，25］，但沱江流域內(nèi)農(nóng)業(yè)用水量占比較高，且農(nóng)業(yè)面源污染嚴(yán)重，因此我們引入農(nóng)業(yè)污水排放總量，用三者共同來表示非期望產(chǎn)出。但由于農(nóng)業(yè)污水排放總量不好量化測算，因此我們用農(nóng)業(yè)面源污染排放量來替代，本文參考陳敏鵬等［26］的研究，用如下公式進(jìn)行測算：

式中：E為農(nóng)業(yè)面源污染排放量；i表示污染來源，鑒于數(shù)據(jù)的可獲得性，本文將化肥、農(nóng)藥和農(nóng)膜作為農(nóng)業(yè)面源污染的主要來源；EUi為i的規(guī)模；ρi、σi和Ci分別為i的產(chǎn)污強(qiáng)度系數(shù)、利用效率系數(shù)以及污染物排放系數(shù)。

2.2 水資源綠色效率的時空差異性分析

根據(jù)沱江流域高質(zhì)量發(fā)展中心對沱江流域的描述，本文對沱江流域進(jìn)行劃分：德陽和成都以上為沱江上游，眉山市、資陽市、樂山市為沱江中游，內(nèi)江市、自貢市、宜賓市、瀘州市為沱江下游，利用超效率SBM 模型對水資源綠色效率進(jìn)行測度。為了更加直觀的展示出沱江流域水資源綠色效率的差異性，采用自然斷點(diǎn)法對水資源綠色效率進(jìn)行階段劃分［27］，如圖1所示。

圖1 沱江流域水資源綠色效率時空分異圖Fig.1 Spatial and temporal differentiation of green efficiency of water resources in the Tuojiang River Basin

由圖1 可知，沱江流域各地級市的水資源綠色效率存在顯著差異。從時間維度上看水資源綠色效率水平波動較大但總體呈上升趨勢，空間維度上表現(xiàn)為從南向北逐漸遞減。其中資陽、成都、自貢以及內(nèi)江的水資源綠色效率水平較高，但近年來增長速度緩慢；眉山、樂山、德陽雖然水資源綠色效率水平較低，但增速較快?？傮w而言，沱江流域水資源綠色效率呈現(xiàn)上游強(qiáng)，中游次之，下游弱的發(fā)展趨勢，各地級市水資源綠色效率的初始水平、發(fā)展路徑、發(fā)展速度均存在明顯差異。

3 水資源綠色效率影響因素的時空異質(zhì)性分析

3.1 變量選取與指標(biāo)說明

水資源綠色效率是一個與自然、社會、環(huán)境3個系統(tǒng)密切相關(guān)的概念。因此，在變量選取時應(yīng)緊靠這3個系統(tǒng)，實(shí)現(xiàn)水資源的可持續(xù)性利用。本文在現(xiàn)有相關(guān)研究的基礎(chǔ)之上，選取了以下變量來分析沱江流域水資源綠色效率的影響因素：

（1）技術(shù)進(jìn)步。技術(shù)進(jìn)步能夠反映沱江流域各城市的研發(fā)投入強(qiáng)度。選取R&D 經(jīng)費(fèi)內(nèi)部支出與地區(qū)生產(chǎn)總值的比例來衡量技術(shù)進(jìn)步。

（2）對外開放水平。本文借鑒秦騰［28］等的研究中對外開放水平的測量，選取沱江流域各城市外商直接投資與地區(qū)成產(chǎn)總值的占比作為測量依據(jù)。

（3）政府管制力度。用地方財政支出占地區(qū)生產(chǎn)總值的比重來衡量。

（4）水資源稟賦。水資源綠色效率是經(jīng)濟(jì)、社會、環(huán)境協(xié)調(diào)統(tǒng)一的總體表現(xiàn)，這里選擇人均水資源占有量來衡量水資源稟賦。

（5）產(chǎn)業(yè)結(jié)構(gòu)。產(chǎn)業(yè)結(jié)構(gòu)在一定程度上反映了水資源在各個產(chǎn)業(yè)的資源消耗程度，沱江流域內(nèi)大部分城市都是以第三產(chǎn)業(yè)為主。為此，本文用第三產(chǎn)業(yè)比重來衡量用水結(jié)構(gòu)。

（6）用水結(jié)構(gòu)。本文選擇生活用水占總用水量百分比來衡量用水結(jié)構(gòu)。

（7）經(jīng)濟(jì)發(fā)展水平。本文采用人均GDP 來衡量沱江流域經(jīng)濟(jì)發(fā)展水平。

3.2 共線性檢驗(yàn)及相關(guān)性分析

為了避免偽回歸情況的出現(xiàn)，基于變量獨(dú)立性原則，利用stata15.0 對各變量進(jìn)行了多重共線性檢驗(yàn)，結(jié)果如表3 所示。各變量的VIF 平均值為2.51，最大值為4.8，均小于5，表明各變量之間不存在多重共線性，變量選取合理。

表3 共線性檢驗(yàn)Tab.3 Collinearity test

進(jìn)一步對所選變量進(jìn)行相關(guān)性分析，各變量的相關(guān)系數(shù)如表4 所示。在5%及以下顯著性水平條件下，技術(shù)進(jìn)步、外開放水平、產(chǎn)業(yè)結(jié)構(gòu)以及經(jīng)濟(jì)發(fā)展水平與水資源綠色效率具有顯著相關(guān)性。在1%及以下顯著水平條件下，政府管制力度、水資源稟賦、用水結(jié)構(gòu)與水資源綠色效率具有顯著相關(guān)性，7個變量均可納入回歸模型進(jìn)行回歸分析。

表4 相關(guān)性檢驗(yàn)Tab.4 Correlation test

3.3 OLS模型和GTWR模型的回歸結(jié)果

GTWR 模型的參數(shù)估計會隨著時空的演變而不同，有效揭示了各影響要素的時空異質(zhì)性。但為了保證回歸結(jié)果的穩(wěn)健性，在進(jìn)行GTWR 回歸之前，需要先對各影響因素做OLS 回歸分析，如表5所示。

由表5可知，對外開放水平、技術(shù)進(jìn)步、用水結(jié)構(gòu)、政府管制力度以及水資源稟賦對沱江流域水資源綠色效率有著重要影響，其中，對外開放水平影響程度最大，技術(shù)進(jìn)步次之，用水結(jié)構(gòu)，政府管制力度和水資源稟賦的影響程度則較小，產(chǎn)業(yè)結(jié)構(gòu)和經(jīng)濟(jì)發(fā)展水平?jīng)]有影響。因此，本文根據(jù)OLS 的回歸結(jié)果，在剔除了產(chǎn)業(yè)結(jié)構(gòu)和經(jīng)濟(jì)發(fā)展水平兩個不顯著變量的基礎(chǔ)上，從時間、空間的角度構(gòu)建局部加權(quán)回歸模型進(jìn)行參數(shù)估計，結(jié)果如表6、7 所示?？梢钥闯鯣TWR 調(diào)整后的R2為93.8%，遠(yuǎn)高于OLS模型的43.9%，說明GTWR模型的擬合優(yōu)度更好。

表5 最小二乘法模型回歸結(jié)果Tab.5 Least squares model regression results

表6 OLS模型與GTWR模型參數(shù)對比Tab.6 Comparison of parameters between OLS model and GTWR model

由表7 可知GTWR 模型測算的結(jié)果變化幅度較大，并且各因素的回歸系數(shù)有正有負(fù)。其結(jié)果表明OLS 回歸分析的參數(shù)值僅代表沱江流域整體的平均特征，從而掩蓋了沱江流域9 個地級市的局部系數(shù)特征，這表明沱江流域水資源綠色效率和各變量之間不僅是簡單的線性關(guān)系，各因素呈現(xiàn)出強(qiáng)烈的空間不穩(wěn)定性，并隨著時間而發(fā)生變化。

表7 GTWR模型回歸結(jié)果Tab.7 GTWR model regression results

3.4 水資源綠色效率影響因素分析

為進(jìn)一步探究各因素的空間變異特征，本文利用ArcGis10.5 對回歸結(jié)果進(jìn)行可視化表達(dá)，從而更加直觀地對各影響因素進(jìn)行空間變異分析，結(jié)果如圖2、3所示。

圖2 2007年沱江流域水資源綠色效率各影響因素回歸系數(shù)空間分布Fig.2 Spatial distribution of regression coefficients of factors influencing green efficiency of water resources in the Tuojiang River Basin in 2007

（1）技術(shù)進(jìn)步對沱江流域水資源綠色效率的影響。研究期初（2007 年），技術(shù)進(jìn)步回歸系數(shù)的高值主要分布在資陽市、內(nèi)江市等地，總體表現(xiàn)為上游、中游地區(qū)要大于下游地區(qū)。研究期末（2017 年），這種趨勢有所變化，總體表現(xiàn)為由西向東遞減的趨勢。從正值區(qū)數(shù)量來看，研究期內(nèi)回歸系數(shù)正值區(qū)數(shù)量有所增加，在空間上主要分布在流域的西部地區(qū)，但仍有部分地區(qū)促進(jìn)作用并不明顯，究其原因可能在于流域內(nèi)多數(shù)城市以工業(yè)發(fā)展為主，科技研發(fā)能力相對較弱，R&D 經(jīng)費(fèi)投入還未轉(zhuǎn)化為實(shí)質(zhì)的技術(shù)進(jìn)步，進(jìn)而未能有效改善水環(huán)境。

圖3 2017年沱江流域水資源綠色效率各影響因素回歸系數(shù)的空間分布Fig.3 Spatial distribution of regression coefficients of factors influencing green efficiency of water resources in the Tuojiang River Basin in 2017

（2）對外開放水平對沱江流域水資源綠色效率的影響。研究期初（2007年），沱江流域?qū)ν忾_放水平回歸系數(shù)的空間分布為“西高東低”，正值分布區(qū)域主要集中在流域的西部；研究期末（2017年），沱江流域?qū)ν忾_放水平回歸系數(shù)的空間分布呈現(xiàn)出由北向南逐步遞減的趨勢?？傮w而言，對外開放促進(jìn)水資源綠色效率提升的區(qū)域主要集中在流域內(nèi)經(jīng)濟(jì)比較發(fā)達(dá)的地區(qū)，如成都市、德陽市。這些地區(qū)擁有較為完整的產(chǎn)業(yè)體系，外商投資能夠更為有效地促進(jìn)各產(chǎn)業(yè)的發(fā)展，在帶來經(jīng)濟(jì)效益的同時也產(chǎn)生了“潛在的”技術(shù)效益，為水污染的綠色治理提供了充足的資金支持和技術(shù)保障。

（3）政府管制力度對水資源綠色效率的影響。研究期初（2007 年），政府管制力度的回歸系數(shù)呈自北向南遞減趨勢，高值區(qū)域主要集中在沱江流域的上游和中游地區(qū)；研究期末（2017 年），負(fù)值區(qū)域明顯增加，多數(shù)城市的回歸系數(shù)由正變?yōu)樨?fù)。原因可能在于這些城市的早期發(fā)展存在諸多不合理之處，政府適當(dāng)參與調(diào)控更有利于水資源綠色效率的提升，隨著經(jīng)濟(jì)發(fā)展的逐漸合理，過度的政府干預(yù)可能會提高企業(yè)的生產(chǎn)成本，減少研發(fā)投入。此外，流域內(nèi)絕大多數(shù)城市在污染治理上的投資較低，污染治理效果微乎其微，城市污水治理能力跟不上污水排放水平。為此，提升污水治理能力，完善污水管控體系是解決目前沱江流域大部分城市面臨的困境和難題的重要舉措。

（4）水資源稟賦對沱江流域水資源綠色效率的影響。研究期初（2007 年），各城市水資源稟賦回歸系數(shù)均為負(fù)值，總體呈自北向南遞減的趨勢。研究期末（2017年），水資源稟賦回歸系數(shù)的正值分布區(qū)域增加明顯，主要集中在沱江流域的中游。上下游地區(qū)水資源豐富，其回歸系數(shù)仍表現(xiàn)為負(fù)，表明隨著人均水資源量的增多，反而不利于水資源綠色效率的提升。這與孫才志等［9］的研究結(jié)論相符，即水資源越豐富的地區(qū)隨著人均水資源占有量的增加反而不利于水資源綠色效率的提升。究其原因可能在水資源豐富，取水容易的地區(qū)，節(jié)水意識普遍不強(qiáng)，錯誤的用水習(xí)慣可能導(dǎo)致水資源的大量浪費(fèi)。

（5）用水結(jié)構(gòu)對沱江流域水資源綠色效率的影響。研究期初（2007年），用水結(jié)構(gòu)回歸系數(shù)的正值主要分布在沱江流域的中下游，負(fù)值區(qū)域則集中在上游，總體呈現(xiàn)出由北向南逐漸增加的趨勢。研究期末（2017年），用水結(jié)構(gòu)的回歸系數(shù)分布有了較大變化，僅有內(nèi)江市、資陽市用水結(jié)構(gòu)回歸系數(shù)為負(fù)。由此可見，流域內(nèi)大部分地區(qū)的用水結(jié)構(gòu)對水資源綠色效率表現(xiàn)為正向促進(jìn)作用。這是因?yàn)橛盟盟Y(jié)構(gòu)反映了沱江流域城市水資源使用在工業(yè)、農(nóng)業(yè)以及生活中的比例，本文選擇生活用水占比來衡量用水結(jié)構(gòu)，當(dāng)生活用水占比越高，水資源投入在工業(yè)和農(nóng)業(yè)灌溉的總量就會降低，而沱江流域的污染主要以工業(yè)污染和農(nóng)藥中的磷、氨、氮等化學(xué)成分為主，因此生活用水占比越高，越有利于沱江流域水資源綠色效率的提升。

3.5 影響因素的趨勢分析

為進(jìn)一步研究沱江流域水資源綠色效率影響因素的空間特征，本文繪制了以各影響因素參數(shù)估計的平均值為高度的三維透視圖，并將點(diǎn)（默認(rèn)北和西）投影于地圖平面兩個方向的平面上，如圖4所示，以檢查沱江流域水資源綠色效率的影響因素的全局趨勢。

圖4 各影響因素回歸系數(shù)估計的趨勢變化圖Fig.4 Trend chart of estimated regression coefficients for each influencing factor

從X方向上看，各因素對沱江流域水資源綠色效率的影響呈現(xiàn)出不同的特征。技術(shù)進(jìn)步對沱江流域水資源綠色效率的影響自西向東呈下降趨勢；水資源稟賦對沱江流域水資源綠色效率自西向東呈增加趨勢；對外開放水平對沱江流域水資源綠色效率的影響自西向東呈倒“U”形；政府管制力度和用水結(jié)構(gòu)對沱江流域水資源綠色效率的影響自西向東呈正“U”形。

從Y方向上看，水資源稟賦對沱江流域水資源綠色效率自北向南的影響呈增加趨勢；技術(shù)進(jìn)步、對外開放水平、用水結(jié)構(gòu)對沱江流域水資源綠色效率的影響自北向南呈倒“U”形，其中，技術(shù)進(jìn)步對沱江流域水資源綠色效率的影響的倒“U”形特征不明顯，敏感度不高；政府管制力度對沱江流域水資源綠色效率的影響自北向南呈正“U”形。綜合來看，各因素對沱江流域水資源綠色效率的影響程度在東西部地區(qū)的變化差異要強(qiáng)于南北地區(qū)。

4 結(jié)論與建議

基于2007-2017 年沱江流域9 個地級市的面板數(shù)據(jù)，在全面考慮水資源使用過程中的社會、經(jīng)濟(jì)以及環(huán)境效益的基礎(chǔ)上，運(yùn)用考慮非期望產(chǎn)出的超效率SBM 模型測度水資源綠色效率水平，并結(jié)合GTWR 模型從時空差異性的視角探究水資源綠色效率的影響因素。結(jié)果表明：沱江流域水資源綠色效率具有明顯的時空差異性，時間維度上水資源綠色效率水平波動較大但仍有上升，空間維度上呈現(xiàn)從南向北逐漸遞減趨勢；OLS模型和GTWR 模型均表明對外開放水平對水資源綠色效率的影響程度最大，技術(shù)進(jìn)步次之，用水結(jié)構(gòu)，政府規(guī)制力度和水資源稟賦的影響程度較??；此外，相較于OLS 模型，GTWR 模型能夠反映各影響因素在時空演變中的局部效應(yīng)，有效彌補(bǔ)了前者的不足。GTWR 模型的回歸結(jié)果顯示沱江流域水資源綠色效率的影響因素呈現(xiàn)出明顯的時空異質(zhì)性，不同年份，不同地區(qū)各影響因素的作用方向及強(qiáng)度都大不相同。

基于上述分析，本文將從治理能力、創(chuàng)新驅(qū)動、政策規(guī)制及用水結(jié)構(gòu)四個方面提出對策建議，以期促進(jìn)區(qū)域水資源的協(xié)調(diào)健康發(fā)展，破解沱江流域水資源發(fā)展失衡困境，為水資源綠色高效發(fā)展保駕護(hù)航。①提升沱江流域協(xié)同治理的能力。各地級市應(yīng)該以“成德眉資”同城化為指導(dǎo)，破除城市之間封閉發(fā)展，營造分工合作的良好氛圍，因地制宜地制定社會經(jīng)濟(jì)發(fā)展和綠色水資源保護(hù)的協(xié)同共生策略。②加強(qiáng)沱江流域創(chuàng)新內(nèi)驅(qū)力。從“創(chuàng)新”角度出發(fā)，提升各城市的自主創(chuàng)新內(nèi)驅(qū)力。改善管理創(chuàng)新“軟實(shí)力”，嚴(yán)格落實(shí)國務(wù)院提出水資源管理規(guī)章制度。③完善沱江流域環(huán)境規(guī)制政策。要根據(jù)各城市的具體污染情況，科學(xué)全面的調(diào)節(jié)水污染治理費(fèi)用，避免嚴(yán)苛的規(guī)制力度對經(jīng)濟(jì)增速的放緩，完善現(xiàn)有污染處理的排放與監(jiān)管機(jī)制，為工業(yè)企業(yè)發(fā)展提供良好的政策環(huán)境。④優(yōu)化沱江流域用水結(jié)構(gòu)。合理配置生產(chǎn)生活用水，建立有效的供給體制，實(shí)行階段性收費(fèi)標(biāo)準(zhǔn)，強(qiáng)化生產(chǎn)生活中節(jié)約用水意識，以實(shí)現(xiàn)沱江流域經(jīng)濟(jì)、社會和生態(tài)效益最大化。