吳浩然+張強

摘要首先運用超效DEA模型測算全要素水資源利用效率，然后利用我國2005—2015年31個省市自治區(qū)的面板數(shù)據(jù)，構(gòu)建面板VAR模型，測度水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的動態(tài)關(guān)系.結(jié)果表明：水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間存在協(xié)整關(guān)系；水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的相互作用存在區(qū)域差異，主要體現(xiàn)在反應(yīng)時滯、作用力度等方面.

關(guān)鍵詞水資源保護(hù)；產(chǎn)業(yè)結(jié)構(gòu)；面板VAR模型

中圖分類號F224文獻(xiàn)標(biāo)識碼A

The Dynamic Relationship between Water Resource Utilization and Optimizing Industrial Structure under the Framework of Panel VAR Model

Haoran WU1，Qiang ZHANG2

（1.Department of Mathematics，Guangdong Polytechnic of Water Resources and Electric Engineering，Guangzhou 510925，

Guangdong，China； 2.School of Management and Economics，Beijing Institute of Technology，Beijing 100081，China）

AbstractBased on the panel data of 31 provinces from 2005 to 2015，this paper established a panel VAR model to measure the dynamic relationship between water resource utilization and optimizing industrial structure.The results show that there is longrun cointegration relationship between water resource utilization and optimizing industrial structure，and the dynamic relationship between water resource utilization and optimizing industrial structure varies with the regions from the perspective of response lag，action intensity and so on.

Key wordswater resource conservation； industrial structure； Panel VAR model

1引言

水是人類生存發(fā)展的基礎(chǔ)性自然資源，其開發(fā)利用與社會經(jīng)濟密切相關(guān).在當(dāng)前產(chǎn)業(yè)進(jìn)行升級性轉(zhuǎn)變的背景下，明晰水資源利用與產(chǎn)業(yè)結(jié)構(gòu)的關(guān)系，既可判斷產(chǎn)業(yè)結(jié)構(gòu)的資源使用效率，也是建立高效節(jié)水產(chǎn)業(yè)結(jié)構(gòu)的理論參考.

現(xiàn)階段，關(guān)于兩者關(guān)系的研究主要分為兩類.一類運用脫鉤分析法、環(huán)境庫茲涅茨曲線考查產(chǎn)業(yè)結(jié)構(gòu)對水資源利用的影響，較為典型的有：劉昌明等（2012）[1]總結(jié)了發(fā)達(dá)國家水資源利用“脫鉤”的顯著特點之一是第三產(chǎn)業(yè)占三次產(chǎn)業(yè)的比重達(dá)到60%左右，并指出通過最嚴(yán)格的需水管理，我國在2030年前后可望實現(xiàn)水資源利用“脫鉤”；賈紹鳳等（2004）[2]以發(fā)達(dá)國家為研究對象對工業(yè)用水庫茲涅茲曲線進(jìn)行研究，檢驗了工業(yè)用水庫茲涅茨曲線的存在性，并指出產(chǎn)業(yè)結(jié)構(gòu)調(diào)整是工業(yè)用水下降的來源之一；魯曉東等（2016）[3]利用我國2004—2013 年八大流域監(jiān)測斷面水文數(shù)據(jù)，依照環(huán)境庫茲涅茨曲線的經(jīng)典和擴展分析框架，發(fā)現(xiàn)產(chǎn)業(yè)特征會影響流域水體的質(zhì)量.另一類則反方向研究水資源利用對產(chǎn)業(yè)結(jié)構(gòu)的影響，如：尹偉強等（2007）[4]以山西省為例，研究了水資源供給對產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)的作用模式，指出在供水短缺壓力下，水資源作為產(chǎn)業(yè)演進(jìn)驅(qū)動力必然在推動區(qū)域產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)中產(chǎn)生重要影響；吳丹、王福林等（2009）[5-6]構(gòu)建了基于水資源優(yōu)化配置的區(qū)域產(chǎn)業(yè)結(jié)構(gòu)演化模型，采用遺傳算法對模型進(jìn)行求解，并結(jié)合實證分析驗證了模型的有效性.

通過對文獻(xiàn)的梳理，可以發(fā)現(xiàn)水資源利用與產(chǎn)業(yè)結(jié)構(gòu)之間存在雙向關(guān)系：一方面，水作為人類發(fā)展不可缺少的自然資源，是產(chǎn)業(yè)進(jìn)行升級性轉(zhuǎn)變的支撐；另一方面，在產(chǎn)業(yè)部門水資源消耗強度不變的前提下，產(chǎn)業(yè)結(jié)構(gòu)的變動會影響水資源消耗結(jié)構(gòu)的變化.然而，現(xiàn)有絕大多數(shù)文獻(xiàn)僅分析了兩者之間的單向關(guān)系，可能導(dǎo)致變量的內(nèi)生性偏差，使結(jié)果出現(xiàn)偏誤.少數(shù)文獻(xiàn)即使探討了兩者之間的雙向關(guān)系，但采用的是靜態(tài)分析方法，只關(guān)注兩者之間的同期影響，不能捕捉面對沖擊時變量的動態(tài)響應(yīng).并且，大多數(shù)文獻(xiàn)以水資源消耗強度作為水資源利用的表征，令水資源投入與經(jīng)濟產(chǎn)出的聯(lián)系簡單化.

為克服上述不足，利用我國2005—2015年間的省際面板數(shù)據(jù)，在以超效DEA模型測算水資源利用效率的基礎(chǔ)上，采用面板向量自回歸模型（簡稱Panel VAR模型），研究水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的雙向、動態(tài)關(guān)系.文章的結(jié)構(gòu)如下：第二部分介紹超效DEA和Panel VAR模型，第三部分進(jìn)行實證分析，第四部分是研究結(jié)論.

2數(shù)據(jù)與模型

2.1模型設(shè)定

2.1.1超效DEA模型

DEA模型利用非參數(shù)方法測算決策單元（Decision Making Unit，簡稱DMU）的相對效率.根據(jù)規(guī)模報酬是否可變，DEA模型可分為兩類：一類是規(guī)模報酬不變的CCR模型，另外一類是規(guī)模報酬可變的BCC模型.本文采用投入導(dǎo)向下的CCR模型，測算我國31個省市自治區(qū)的水資源利用效率.設(shè)有n個DMU，每個DMU有m種投入和s種產(chǎn)出，xjk為第k個DMU的第j個投入變量，yjk為第k個DMU的第j個產(chǎn)出變量，則DMUk的效率值θ可以從以下線性規(guī)劃中獲得：endprint

式（1）中，xk=（x1k，x2k，…，xmk）T；yk=（y1k，y2k，…，ymk）T；λj為單位組合系數(shù)；θ為DMUk的效率值，滿足0≤θ≤1.當(dāng)θ=1時，DMUk是效率前沿面上的點，處于有效狀態(tài)；θ<1時，DMUk處于無效狀態(tài)，存在效率損失.

DEA模型的弱點是當(dāng)θ=1的有效DMU較多時，無法對有效DMU的效率水平進(jìn)一步區(qū)分[7-8].針對這一不足，Andersen和Petersen建立了超效DEA模型，能使有效DMU之間也能比較效率水平的高低，其思想是在評估DMUk的效率時，將DMUk本身排除在決策單元的集合之外：

利用式（2）評估DMU的效率水平時，對于無效DMU，其效率值與CCR模型相同；對于有效DMU，可以使其投入按比例增加，而效率值保持不變，其投入增加比例即超效率值，因此有效DMU的效率值θ≥1，并且θ越大，效率越高.

2.1.2Panel VAR模型

相比VAR模型，由Holtz等提出的Panel VAR模型放松了較大樣本觀測量的限制，在T≥2p+2（T為時間序列的長度，p為滯后項階數(shù)）的情況下能較穩(wěn)態(tài)地估計參數(shù).除此之外，Panel VAR模型允許樣本個體存在不可觀察的差異，可以捕捉到個體在橫截面上可能受到的共同沖擊.根據(jù)研究內(nèi)容，本文構(gòu)建的Panel VAR模型如下

式（3）中，i為地區(qū)，t為年份，yi，t為內(nèi)生變量向量，p為滯后項階數(shù)，αj為滯后內(nèi)生變量的回歸系數(shù)，ηi為個體固定效應(yīng)，i為時間效應(yīng)，εi，t為隨機擾動項.

2.2變量處理

2.2.1水資源利用效率（Wre）

水作為一種自然資源，必須和其他生產(chǎn)要素結(jié)合才能帶來經(jīng)濟產(chǎn)出，借助傳統(tǒng)生產(chǎn)函數(shù)Y=f（K，L，W），可以構(gòu)建包含水資源在內(nèi)的投入產(chǎn)出分析框架.其中，Y為經(jīng)濟產(chǎn)出，以各產(chǎn)業(yè)的總產(chǎn)值表示；K為資本投入，以各產(chǎn)業(yè)的固定資產(chǎn)投資表示；L為勞動投入，以各產(chǎn)業(yè)的從業(yè)人員數(shù)量表示；W為水資源投入，以三次產(chǎn)業(yè)的用水總量表示.

為消除物價波動的影響，各產(chǎn)業(yè)的經(jīng)濟產(chǎn)出Y和資本投入K，按照相應(yīng)年份的平減指數(shù)以2005年為基期計算得到；各產(chǎn)業(yè)的從業(yè)人員數(shù)量L取上一年年末從業(yè)人數(shù)與本年年末從業(yè)人數(shù)的平均數(shù)；第一、二產(chǎn)業(yè)用水量從歷年《中國統(tǒng)計年鑒》直接獲取，而第三產(chǎn)業(yè)用水量在年鑒中未直接給出，按照“第三產(chǎn)業(yè)用水量=生活用水量+生態(tài)用水量-居民家庭用水量”近似計算得到，其中，居民家庭用水量由居民生活用水定額（102.9升/人.天）和居民人口數(shù)加成折算.

2.2.2產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化

產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化指生產(chǎn)要素在經(jīng)濟各部門和不同產(chǎn)業(yè)之間的重新配置，主要表現(xiàn)為國民經(jīng)濟中產(chǎn)業(yè)結(jié)構(gòu)的合理化調(diào)整以及產(chǎn)業(yè)結(jié)構(gòu)的高級化[9].因此，為反映產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化的內(nèi)涵和本質(zhì)，引入兩個變量進(jìn)行描述：

產(chǎn)業(yè)結(jié)構(gòu)合理化（Ish）.產(chǎn)業(yè)結(jié)構(gòu)合理化著眼于要素資源在產(chǎn)業(yè)間的配置、協(xié)調(diào)和利用效率.采用國際基準(zhǔn)法，選擇區(qū)域產(chǎn)業(yè)結(jié)構(gòu)與錢納里三次產(chǎn)業(yè)結(jié)構(gòu)模式的Hamming貼近度，作為區(qū)域內(nèi)各單元產(chǎn)業(yè)結(jié)構(gòu)合理化的度量：

式（4）中，Ish為Hamming貼近度；Sdk為各產(chǎn)業(yè)產(chǎn)值比重；STk為錢納里三次產(chǎn)業(yè)結(jié)構(gòu)模式中各產(chǎn)業(yè)產(chǎn)值比重；k為產(chǎn)業(yè)類型.Ish的值越大，產(chǎn)業(yè)結(jié)構(gòu)越合理.

產(chǎn)業(yè)結(jié)構(gòu)高級化（Isc）.產(chǎn)業(yè)結(jié)構(gòu)高級化指產(chǎn)業(yè)結(jié)構(gòu)從勞動密集型到資本密集型、再到知識技術(shù)密集型產(chǎn)業(yè)的順次轉(zhuǎn)換.學(xué)界常以各產(chǎn)業(yè)比例關(guān)系定量產(chǎn)業(yè)結(jié)構(gòu)高度，本文以Clark指標(biāo)，即非農(nóng)產(chǎn)值比重作為產(chǎn)業(yè)高度的度量：

式（5）中，Sd2、Sd3分別為第二、三產(chǎn)業(yè)占比.Isc的值越大，產(chǎn)業(yè)結(jié)構(gòu)水平越高.

2.3數(shù)據(jù)來源及區(qū)域劃分

研究的時間跨度為2005—2015年，數(shù)據(jù)來源于《中國統(tǒng)計年鑒》.考慮到我國社會經(jīng)濟發(fā)展的區(qū)域差異，各區(qū)域水資源利用與產(chǎn)業(yè)結(jié)構(gòu)之間的關(guān)系未必會遵循同一經(jīng)驗規(guī)律，參照熊啟躍、張依茹（2012）[10]的研究方法，將我國分為經(jīng)濟發(fā)達(dá)地區(qū)、經(jīng)濟次發(fā)達(dá)地區(qū)和經(jīng)濟不發(fā)達(dá)地區(qū)分別考察.其中，經(jīng)濟發(fā)達(dá)地區(qū)包括北京、江蘇等7個省市，經(jīng)濟次發(fā)達(dá)地區(qū)包括河北、遼寧等13個省市自治區(qū)，經(jīng)濟不發(fā)達(dá)地區(qū)包括內(nèi)蒙古、吉林等11個省市自治區(qū).

3實證結(jié)果與分析

實證分析的思路是：首先利用（2）式對水資源利用效率進(jìn)行測算；然后進(jìn)行面板單位根檢驗和面板協(xié)整分析，以檢驗數(shù)據(jù)的平穩(wěn)性以及變量之間的長期均衡關(guān)系；最后運用Panel VAR模型，考察水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的雙向、動態(tài)關(guān)系.實證分析所采用的軟件為EMS1.3和stata12.0.

3.1面板單位根檢驗

為了避免“偽回歸”現(xiàn)象的出現(xiàn)，在進(jìn)行Panel VAR模型估計之前應(yīng)先進(jìn)行面板單位根檢驗.出于對結(jié)論穩(wěn)健性的考慮，采用同質(zhì)面板單位根LLC檢驗和異質(zhì)面板單位根IPS、ADF、PP檢驗等四種方法進(jìn)行檢驗，結(jié)果見表1.

由表1可知，對Wre、Ish和Isc的水平值進(jìn)行檢驗時，檢驗結(jié)果表明不能完全拒絕“存在單位根”的原假設(shè)，而對這三個變量的一階差分值進(jìn)行檢驗時，均顯著地拒絕“存在單位根”的原假設(shè).因此，三個地區(qū)的Wre、Ish和Isc均為一階單整序列，符合面板協(xié)整檢驗的要求.

3.2面板協(xié)整檢驗

面板單位根檢驗的結(jié)果表明Wre、Ish和Isc均為一階單整，因此，它們之間可能存在協(xié)整關(guān)系.根據(jù)Pedroni 提出的面板數(shù)據(jù)協(xié)整檢驗方法，以回歸殘差為基礎(chǔ)構(gòu)造的4個組內(nèi)統(tǒng)計量和3個組間統(tǒng)計量進(jìn)行面板協(xié)整檢驗，結(jié)果如表2所示.

從表2可知，發(fā)達(dá)地區(qū)7個統(tǒng)計量均通過顯著性檢驗，所以，該地區(qū)的Wre、Ish和Isc存在面板協(xié)整關(guān)系.次發(fā)達(dá)、不發(fā)達(dá)地區(qū)分別有Panel rho和Panel v統(tǒng)計量沒有通過顯著性檢驗，但Pedroni的Monte Carlo模擬實驗結(jié)果顯示，對于小樣本（T<20，T為時間序列的長度）來說， Panel ADF和Group ADF統(tǒng)計量最有效力，其次是Panel PP和Group PP統(tǒng)計量，其他則最差[11].所以Panel rho和Panel v統(tǒng)計量沒有通過顯著性檢驗對次發(fā)達(dá)、不發(fā)達(dá)地區(qū)的三個變量之間存在協(xié)整關(guān)系的結(jié)論沒有影響.endprint

3.3Panel VAR估計

Panel VAR估計主要由三個部分構(gòu)成：一是GMM估計，用于說明變量之間的回歸關(guān)系；二是脈沖響應(yīng)分析，用于觀察變量面對沖擊的反應(yīng)情況；三是誤差項的方差分解，用于分析誤差項的影響因素大小.在估計前，綜合考慮BIC信息準(zhǔn)則和脈沖響應(yīng)函數(shù)的收斂情況，發(fā)達(dá)地區(qū)、次發(fā)達(dá)地區(qū)和不發(fā)達(dá)地區(qū)的最優(yōu)滯后階數(shù)p分別取2、2和1（由于篇幅限制，略去各地區(qū)的最優(yōu)滯后階數(shù)檢驗結(jié)果）.

3.3.1GMM估計

由于模型（3）中包含了個體固定效應(yīng)ηi和時間效應(yīng)i，可能造成估計結(jié)果的偏誤.因此，首先采用前向均值差分法剔除個體固定效應(yīng)，采用橫截面上的均值差分法剔除時間效應(yīng)，之后利用GMM估計得到Wre、Ish和Isc之間的Panel VAR模型參數(shù)估計結(jié)果，如表3所示.

就發(fā)達(dá)地區(qū)而言，對于Wre方程，滯后一期Ish的估計系數(shù)顯著為負(fù)，表明上一年Ish的提高會顯著引起下一年Wre的降低，說明產(chǎn)業(yè)結(jié)構(gòu)的合理化對水資源利用起到了阻礙作用；滯后二期Isc的估計系數(shù)顯著為正，說明產(chǎn)業(yè)結(jié)構(gòu)的高級化對水資源利用起到了促進(jìn)作用.結(jié)合其它兩個方程，可以發(fā)現(xiàn)：滯后二期Wre對Ish有顯著的負(fù)向影響；雖然無論滯后一期還是滯后二期，Wre對Isc均無顯著影響，但Wre通過顯著影響Ish而間接引起Isc的變化.

考察其它地區(qū)待估參數(shù)的符號和顯著性，可知：在次發(fā)達(dá)地區(qū)，滯后二期Isc對Wre的影響顯著為正，Ish通過顯著影響Isc而間接引起Wre的變化，滯后二期Wre對Isc的影響顯著為正，Wre通過顯著影響Isc而間接引起Ish的變化；在不發(fā)達(dá)地區(qū)，滯后一期Isc對Wre的影響顯著為正.

由于Panel VAR模型引入了較多變量的滯后項，會導(dǎo)致某些估計參數(shù)不顯著，因此，在實證中往往還需要運用脈沖響應(yīng)和方差分解作進(jìn)行進(jìn)一步分析.

3.3.2脈沖響應(yīng)分析

脈沖響應(yīng)函數(shù)描述在隨機擾動項εi，t上施加一個標(biāo)準(zhǔn)差大小的沖擊對內(nèi)生變量當(dāng)前值和未來值的影響，能直觀地刻畫出變量之間的動態(tài)關(guān)系[12].取沖擊響應(yīng)期為10 期，通過Monte Carlo實驗?zāi)M500次后得到我國三個地區(qū)的Wre、Ish和Isc之間的脈沖響應(yīng)函數(shù)，如圖1所示.

從脈沖響應(yīng)圖可以看出，Wre、Ish和Isc的脈沖響應(yīng)經(jīng)過10期后均呈現(xiàn)收斂趨勢，表明本文所構(gòu)建的Panel VAR模型是穩(wěn)健的.

就發(fā)達(dá)地區(qū)而言，面對Ish的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，在第一期有一個較大的負(fù)響應(yīng)，在第二期呈現(xiàn)出微小的正響應(yīng)，之后逐步收斂于負(fù)的均衡值.表明在短期內(nèi)Ish對Wre的影響具有不確定性，但長期來看，兩者的關(guān)系趨于穩(wěn)定，即Ish的提高會引起Wre的降低.面對Wre的一個正交化沖擊，Ish在當(dāng)期出現(xiàn)一個負(fù)響應(yīng)，在第一期呈現(xiàn)出微小的正響應(yīng)，之后在第二期變?yōu)檩^大的負(fù)響應(yīng)并且逐步收斂于負(fù)的均衡值.表明盡管在短期Wre的提高會引起Ish的提高，但這種正效應(yīng)并不具有持續(xù)性.面對Isc的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，在第一期呈現(xiàn)較小的負(fù)響應(yīng)，第二期變?yōu)檎憫?yīng)，之后逐漸下降直到收斂于零，表明Isc的提高會引起Wre的提高，但存在兩個周期的滯后.面對Wre的一個正交化沖擊，Isc始終表現(xiàn)為正響應(yīng)，表明無論短期還是長期，Wre的提高均會引起Isc的提高.

就次發(fā)達(dá)地區(qū)而言，面對Ish的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，在第一期達(dá)到峰值后開始下降，并且逐步收斂于正的均衡值.面對Wre的一個正交化沖擊，Ish的響應(yīng)始終為負(fù)；面對Isc的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，在第一期開始呈現(xiàn)正的響應(yīng)，這種正響應(yīng)有所上升，在第二期達(dá)到峰值，之后逐漸下降直到收斂于零.面對Wre的一個正交化沖擊，Isc始終表現(xiàn)為正響應(yīng).

就不發(fā)達(dá)地區(qū)而言，面對Ish的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，但之后始終保持在一種負(fù)的響應(yīng)狀態(tài).面對Wre的一個正交化沖擊，Ish在當(dāng)期出現(xiàn)一個較大的負(fù)響應(yīng)，但第一期之后即收斂于零，呈現(xiàn)基本沒有響應(yīng)的狀態(tài).面對Isc的一個正交化沖擊，Wre在當(dāng)期沒有響應(yīng)，在第一期達(dá)到峰值后開始下降，并且逐步收斂于正的均衡值.面對Wre的一個正交化沖擊，Isc始終表現(xiàn)為正響應(yīng).

3.3.3方差分解

為了更深入地考察Wre、Ish和Isc之間的相互影響程度，通過方差分解的方法，獲得不同Panel VAR方程的沖擊反應(yīng)對內(nèi)生變量波動的貢獻(xiàn)度，如表5所示.

從表5可知：①選取10個預(yù)測期與選取20個預(yù)測期對方程分析的影響不大，說明經(jīng)過10個預(yù)測期以后，系統(tǒng)趨于穩(wěn)定狀態(tài).②Wre、Ish和Isc的波動主要來自于自身，三者對自身波動的貢獻(xiàn)率均在75%以上.③Wre對Ish的影響在1%-13%之間，其中發(fā)達(dá)地區(qū)所受影響最大，次發(fā)達(dá)和不發(fā)達(dá)地區(qū)次之；Wre對Isc的影響在2%-13%之間，其中不發(fā)達(dá)地區(qū)所受影響最大，發(fā)達(dá)地區(qū)次之，次發(fā)達(dá)地區(qū)相對較小.④Ish對Wre的影響在9%-12%之間，其中次發(fā)達(dá)地區(qū)所受影響最大，發(fā)達(dá)和不發(fā)達(dá)地區(qū)次之；Isc對Wre的影響在2%-15%之間，其中不發(fā)達(dá)地區(qū)所受影響最大，次發(fā)達(dá)地區(qū)次之，發(fā)達(dá)地區(qū)相對較小.

4結(jié)論

利用我國11年間的省際面板數(shù)據(jù)，構(gòu)建Panel VAR模型，對我國三個經(jīng)濟區(qū)域的水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的雙向、動態(tài)關(guān)系進(jìn)行研究，發(fā)現(xiàn)以下三點.

（1）水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間存在長期均衡關(guān)系.

（2）水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間具有雙向關(guān)系，但水資源利用與產(chǎn)業(yè)結(jié)構(gòu)的合理化、高級化之間的關(guān)系并不一致：水資源利用與產(chǎn)業(yè)結(jié)構(gòu)的合理化之間以負(fù)向關(guān)系為主，而水資源利用與產(chǎn)業(yè)結(jié)構(gòu)的高級化之間則為正向關(guān)系.

（3）水資源利用與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化之間的關(guān)系具有明顯的區(qū)域差異，主要體現(xiàn)在反應(yīng)時滯、作用力度等方面.

參考文獻(xiàn)

[1]劉昌明，趙彥琦.中國實現(xiàn)水需求零增長的可能性探討[J].中國科學(xué)院院刊，2012，27（4）：439-446.

[2]賈紹鳳，張士鋒，楊紅，夏軍.工業(yè)用水與經(jīng)濟發(fā)展的關(guān)系-用水庫茲涅茨曲線[J].自然資源學(xué)報，2004，19（3）：279-284.

[3]魯曉東，許羅丹，熊瑩.水資源環(huán)境與經(jīng)濟增長：EKC 假說在中國八大流域的表現(xiàn)[J].經(jīng)濟管理，2016，38（1）：20-29.

[4]尹偉強，司春霞.水資源與山西省產(chǎn)業(yè)結(jié)構(gòu)、空間結(jié)構(gòu)演進(jìn)分析[J].水資源與水工程學(xué)報，2007，18（5）：52-55.

[5]吳丹，吳鳳平.基于水資源環(huán)境綜合承載力的區(qū)域產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化研究[J].統(tǒng)計與決策，2009（22）：100-102.

[6]王福林，吳丹.基于水資源優(yōu)化配置的區(qū)域產(chǎn)業(yè)結(jié)構(gòu)動態(tài)演化模型[J].軟科學(xué)，2009，23（5）：92-96.

[7]袁海，吳振榮.中國省域文化產(chǎn)業(yè)效率測算及影響因素實證分析[J].軟科學(xué)，2012，26（3）：72-77.

[8]劉殿國，劉淼，張婕.基于超效模型的泛珠三角區(qū)第三產(chǎn)業(yè)生產(chǎn)效率評價[J].數(shù)學(xué)的實踐與認(rèn)識，2014，44（22）：8-12.

[9]韓永輝，黃亮雄，王賢彬.產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級改進(jìn)生態(tài)效率了嗎？[J].數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究，2016，33（4）：40-59.

[10]熊啟躍，張依茹.貨幣政策信貸渠道的經(jīng)濟區(qū)域效應(yīng)研究-基于我國31個省際面板數(shù)據(jù)的經(jīng)驗證據(jù)[J].投資研究，2012，31（7）：78-89.

[11]潘丹，應(yīng)瑞瑤.中國水資源與農(nóng)業(yè)經(jīng)濟增長關(guān)系研究[J].中國人口·資源與環(huán)境，2012，22（1）：161-166.

[12]黃可人，韋廷柒.經(jīng)濟增長、產(chǎn)業(yè)結(jié)構(gòu)變遷與城鄉(xiāng)居民收入差距-基于PVAR模型的動態(tài)分析[J].工業(yè)技術(shù)經(jīng)濟，2016，35（4）：145-152.endprint