• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    新疆區(qū)域水資源對(duì)產(chǎn)業(yè)生態(tài)系統(tǒng)與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)聯(lián)
    ——基于VAR模型

    2017-10-13 11:29:05張振龍
    生態(tài)學(xué)報(bào) 2017年16期
    關(guān)鍵詞:貢獻(xiàn)度用水量總量

    張振龍,孫 慧,*

    1 新疆大學(xué)新疆創(chuàng)新管理研究中心, 烏魯木齊 830046 2 新疆大學(xué)經(jīng)濟(jì)與管理學(xué)院, 烏魯木齊 830046

    新疆區(qū)域水資源對(duì)產(chǎn)業(yè)生態(tài)系統(tǒng)與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)聯(lián)
    ——基于VAR模型

    張振龍1,2,孫 慧1,2,*

    1 新疆大學(xué)新疆創(chuàng)新管理研究中心, 烏魯木齊 830046 2 新疆大學(xué)經(jīng)濟(jì)與管理學(xué)院, 烏魯木齊 830046

    新疆正面臨“五化同步”對(duì)水資源的需求不斷增長與水資源開發(fā)利用過度的矛盾。水資源瓶頸制約已成為影響新疆經(jīng)濟(jì)可持續(xù)發(fā)展和長治久安的突出問題之一。運(yùn)用VAR模型,通過ADF檢驗(yàn)、脈沖響應(yīng)函數(shù)和方差貢獻(xiàn)度分解,對(duì)2000—2014年新疆耗水產(chǎn)業(yè)生態(tài)系統(tǒng)和經(jīng)濟(jì)增長的長期均衡關(guān)系進(jìn)行實(shí)證分析。結(jié)果表明:(1)經(jīng)濟(jì)增長與總用水量、工業(yè)用水量和農(nóng)業(yè)用水量之間均存在長期均衡關(guān)系;(2)經(jīng)濟(jì)發(fā)展對(duì)用水量產(chǎn)生負(fù)向沖擊,工業(yè)用水量和農(nóng)業(yè)用水量隨著經(jīng)濟(jì)發(fā)展出現(xiàn)正向沖擊效應(yīng)。(3)新疆經(jīng)濟(jì)快速增長伴隨著水資源的大力開發(fā)和過度利用。據(jù)此提出對(duì)策建議,通過實(shí)施嚴(yán)格的退地減水政策,明確用水總量控制和定額指標(biāo),加強(qiáng)跨流域調(diào)水工程建設(shè),防止浪費(fèi)等多途徑維持新疆水資源可持續(xù)利用和耗水產(chǎn)業(yè)健康發(fā)展。

    產(chǎn)業(yè)生態(tài)系統(tǒng);經(jīng)濟(jì)增長;動(dòng)態(tài)關(guān)系;VAR模型;新疆

    Abstract: Xinjiang is coping with the growing concern of increasing level of demand and overuse of water resources caused by 5-tion synchronization. Excessive utilization of water resources strongly influences local economic development and plays a vital role in accelerating ecological degradation. The bottleneck of water resources has become one of the outstanding issues that affects sustainable development and long-term stability of Xinjiang′s economy. To address these concerns, the government should strengthen the construction of water infrastructure for improving water use efficiency, and thereby increasing the quantity and frequency of irrigation during an ecological recovery. The VAR model can deal with conditions when some random parameters or variables appear in the modeling formulation of a program. Such random parameters are expressed as probability density functions (PDFs). This implies that an ADF test technique can be employed when the quality of uncertain information is comprehensive. The results can be interpreted under different levels of probabilities (or risks). Previously, a number of stochastic programming methods were developed to support decision making for water resources planning and management. Among them, the VAR model proved effective for problems where analysis of policy scenarios is desired, and the right-hand-side of the coefficients are random with known probability distributions. Therefore, as an extension of previous works, the objective of this study is to present practical research on the long-term equilibrium relationship between Xinjiang's utilization of water resources and economic growth using the VAR model, ADF test, and cointegration, accompanied by the decomposition of generalized impulse response and variance contribution degree. Firstly, the results show the existence of a long-run equilibrium relationship between economic growth and total water consumption, industrial water, and agricultural water. Secondly, rapid economic development has a negative impact on water use; yet, industrial and agricultural use of water has a positive impact. Finally, accompanied by the booming economy of Xinjiang, it carried the problem of over-exploitation and over-utilization of water resources. Based on the above factors, this article recommends several suggestions for the sustainable utilization of water resources and the robust development of Xinjiang’s economy, including implementing strict policies on conserving the depleting land and water resources, clarifying the amount of water use, strengthening the construction of inter-basin water transfer projects, and preventing excessive consumption, among others. The results obtained will help local decision makers to formulate policies that mandate effective water exploitation and allocation, and thus facilitate the local agricultural sustainability.

    KeyWords: water-intensive industrial ecosystem; economic growth; dynamic relationship; VAR model; Xinjiang

    生態(tài)系統(tǒng),水資源無疑是其中重要的組成部分,也是維系生態(tài)系統(tǒng)正常運(yùn)行的必要條件。對(duì)生態(tài)系統(tǒng)進(jìn)行系統(tǒng)化分析,通過調(diào)研或者模擬,合理配置水資源成為當(dāng)前研究的一個(gè)重要熱點(diǎn)。新疆地區(qū),以其特殊的自然環(huán)境、獨(dú)特的綠洲經(jīng)濟(jì)以及重要的國家戰(zhàn)略地位,決定了水資源合理配置是其其他戰(zhàn)略資源開發(fā)利用的基礎(chǔ)與前提[1]。新疆地域遼闊、資源豐富,但遠(yuǎn)離海洋,“三山夾兩盆”的地貌特征使其形成了典型的干旱氣候,生態(tài)環(huán)境十分脆弱,水資源相對(duì)緊缺。2014年全疆水資源總量726.93×108m3,其中,地表水資源量686.55×108m3,地下水資源量443.93×108m3,地表水與地下水之間重復(fù)量高達(dá)403.55×108m3。人均水資源量只有3130 m3/人。全疆多年平均入境水量為90×108m3,出境水量222×108m3。但是新疆水資源空間分布極不均衡,呈現(xiàn)出“北多南少,西多東少”的特征,主要體現(xiàn)在北疆面積占新疆的27.7%,但年徑流量則占到全疆的50%,南疆面積占新疆的72.3%,而年徑流量僅占50%,按單位面積徑流來看,北疆為南疆的3倍多。農(nóng)業(yè)用水比例偏大,工程性缺水問題突出,已成為制約新疆社會(huì)經(jīng)濟(jì)發(fā)展的主要因素之一。2014年全疆用水總量581.82×108m3,其中農(nóng)業(yè)用水量550.99×108m3,占全疆用水總量的94.70%,是新疆最主要的耗水部門[2]。在有限的水資源中,吐哈產(chǎn)業(yè)聚集區(qū)和天山南坡產(chǎn)業(yè)帶的水資源開發(fā)利用率已超過或接近100%。處理好經(jīng)濟(jì)發(fā)展與資源保護(hù)的關(guān)系,探索水資源可持續(xù)利用與主要耗水產(chǎn)業(yè)系統(tǒng)協(xié)調(diào)發(fā)展,是新疆社會(huì)各界高度重視的重大課題,也是探索生態(tài)與經(jīng)濟(jì)、人與自然協(xié)調(diào)發(fā)展新模式的重要方面。

    學(xué)者們針對(duì)水資源利用與區(qū)域復(fù)雜系統(tǒng)中經(jīng)濟(jì)發(fā)展與生態(tài)環(huán)境的相關(guān)性問題開展了相對(duì)豐富的研究,其研究領(lǐng)域涵蓋水資源供需均衡[3-7]、水資源承載力[8-9]、水資源配置[10- 13]、水資源持續(xù)利用[14-16]、生態(tài)需水[17-19]等方面。其中涉及到耗水產(chǎn)業(yè)生態(tài)系統(tǒng)與經(jīng)濟(jì)發(fā)展關(guān)系的研究主要集中在水資源承載力方面。但是,學(xué)者們關(guān)于水資源承載能力的研究多從承載力的內(nèi)涵及外延[20]、承載能力的衡量標(biāo)準(zhǔn)[21]等方面著手開展,并廣泛使用工程技術(shù)案例,而從生態(tài)經(jīng)濟(jì)、生態(tài)管理等角度的研究剛剛起步[22],尤其是就新疆這一特殊區(qū)域水資源利用進(jìn)行研究的文獻(xiàn)偏少。而現(xiàn)有的研究中,有研究對(duì)新疆水資源利用與經(jīng)濟(jì)增長進(jìn)行研究,主要就GDP總量與農(nóng)業(yè)用水、工業(yè)用水、生活用水的相關(guān)性進(jìn)行了分析[23];也有早期的研究以干旱地區(qū)塔里木河下游尉犁地區(qū)為例,選取生態(tài)-生產(chǎn)-生活系統(tǒng)承載力的綜合測度指標(biāo),綜合評(píng)估了生態(tài)-生產(chǎn)-生活系統(tǒng)承載力的現(xiàn)狀情勢,進(jìn)而預(yù)測未來30年塔里木河下游尉犁地區(qū)生態(tài)-生產(chǎn)-生活承載力變化的總體變化情景[24]。而事實(shí)上,不同產(chǎn)業(yè)水資源利用方式和效果差異很大,將新疆地區(qū)生產(chǎn)總值這一宏觀變量分別與各行業(yè)用水量進(jìn)行關(guān)聯(lián)性分析,容易忽略各產(chǎn)業(yè)的隱含信息,難以體現(xiàn)復(fù)雜生態(tài)巨系統(tǒng)各產(chǎn)業(yè)部門與水資源利用的微觀實(shí)際,而多產(chǎn)業(yè)變量間的相互作用也成為研究的難點(diǎn)。

    自從Sims[25]的開創(chuàng)性工作以來,矢量自回歸VAR (Vector AutoRegressive)模型已廣泛應(yīng)用于實(shí)證研究。VAR模型相對(duì)于單方程模型(Single Equation Models)的優(yōu)點(diǎn)是,可以通過脈沖響應(yīng)分析幾個(gè)變量之間的相互作用。VAR模型已廣泛用于時(shí)間序列的研究中,也有一些研究在面板數(shù)據(jù)上使用VAR模型,如基于工具變量(IVs)的估計(jì)[26]、偏差校正固定效應(yīng)(FE)估計(jì)[27]、準(zhǔn)最大似然(QML)和矩陣的廣義方法(GMM)估計(jì)[28]。此外,有幾個(gè)使用面板VAR模型的實(shí)證研究[29-31]?,F(xiàn)有研究分析針對(duì)全國水資源利用與經(jīng)濟(jì)增長關(guān)系的探討較多,但對(duì)區(qū)域水資源利用與經(jīng)濟(jì)增長和產(chǎn)業(yè)生態(tài)系統(tǒng)的動(dòng)態(tài)關(guān)系的研究相對(duì)較少。因此,本文將梳理新疆產(chǎn)業(yè)生態(tài)系統(tǒng)中不同組分與水資源利用之間的關(guān)聯(lián)性,采用VAR模型,分別就用水總量與新疆產(chǎn)業(yè)生態(tài)系統(tǒng)的經(jīng)濟(jì)增長、農(nóng)業(yè)用水量與農(nóng)業(yè)經(jīng)濟(jì)增長、工業(yè)用水量與工業(yè)經(jīng)濟(jì)增長三方面的關(guān)系進(jìn)行詳細(xì)探討,以期為新疆進(jìn)一步調(diào)整優(yōu)化產(chǎn)業(yè)結(jié)構(gòu),實(shí)現(xiàn)水資源持續(xù)利用以及新疆經(jīng)濟(jì)可持續(xù)發(fā)展提供理論參考。

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

    1.1 研究方法——VAR模型

    向量自回歸(VAR)模型,就是用所有當(dāng)期變量對(duì)其若干滯后期變量進(jìn)行回歸,通常用來估計(jì)相互聯(lián)系的時(shí)間序列系統(tǒng)以及分析隨機(jī)擾動(dòng)對(duì)變量系統(tǒng)的動(dòng)態(tài)關(guān)系[32],不需要提前設(shè)定任何約束條件。VAR方法通過把系統(tǒng)中每一個(gè)內(nèi)生變量作為系統(tǒng)中所有內(nèi)生變量的滯后值的函數(shù)來構(gòu)造模型,從而回避了結(jié)構(gòu)化模型的需要。一個(gè)VAR(p)模型的數(shù)學(xué)形式[33]:

    yt=A1yt-1+…+Apyt-p+Bxt+εt

    (1)

    式中,yt是內(nèi)生變量;yt-1, …,yt-p表示yt的滯后期;xt是外生變量;A1, …,Ap表示yt的待估系數(shù);B是xt待估系數(shù);εt是隨機(jī)擾動(dòng)項(xiàng)。VAR模型轉(zhuǎn)化為矩陣為:

    (2)

    式中,y1t, …,ykt作為內(nèi)生變量,可以同期相關(guān),而y1t-1, …,ykt-p作為滯后變量均在等號(hào)右邊,因此不會(huì)出現(xiàn)同期相關(guān)問題,OLS仍然是有效的。

    1.2 指標(biāo)選取與數(shù)據(jù)來源

    本文主要考察區(qū)域產(chǎn)業(yè)巨系統(tǒng)中產(chǎn)業(yè)用水與經(jīng)濟(jì)增長的互動(dòng)。在產(chǎn)業(yè)用水的指標(biāo)選擇中,我國水資源利用結(jié)構(gòu)中通常將用水分為農(nóng)業(yè)用水、工業(yè)用水、居民生活用水和生態(tài)環(huán)境用水[34]。由于數(shù)據(jù)的支持情況,本研究將主要研究產(chǎn)業(yè)用水與經(jīng)濟(jì)增長的關(guān)系,故居民生活用水和生態(tài)環(huán)境用水將在后續(xù)研究進(jìn)一步開展。因此,本研究將產(chǎn)業(yè)用水分解為3個(gè)指標(biāo),以用水總量(X1)表示產(chǎn)業(yè)用水總量,以農(nóng)業(yè)用水量(X2)表示第一產(chǎn)業(yè)的用水總量,以工業(yè)用水量(X3)表示第二產(chǎn)業(yè)的用水總量。在經(jīng)濟(jì)增長指標(biāo)選擇方面,由于GDP是用來衡量國家、地區(qū)經(jīng)濟(jì)活動(dòng)總量的國際通用指標(biāo),因此,本研究采用GDP作為度量經(jīng)濟(jì)增長的指標(biāo),并將其分解為三個(gè)指標(biāo),以新疆GDP總量(Y1)表示宏觀經(jīng)濟(jì)增長,以農(nóng)業(yè)產(chǎn)值(Y2)衡量農(nóng)業(yè)經(jīng)濟(jì),以工業(yè)產(chǎn)值(Y3)衡量工業(yè)經(jīng)濟(jì)。以用水總量和新疆GDP來考察產(chǎn)業(yè)用水與經(jīng)濟(jì)增長的總體關(guān)系;以農(nóng)業(yè)用水量和農(nóng)業(yè)產(chǎn)值來考察農(nóng)業(yè)用水與農(nóng)業(yè)經(jīng)濟(jì)的動(dòng)態(tài)關(guān)系;以工業(yè)用水量和工業(yè)產(chǎn)值來考察工業(yè)用水與工業(yè)經(jīng)濟(jì)的動(dòng)態(tài)關(guān)系。

    由于我國自1997年才開始正式編制《中國水資源公報(bào)》,而新疆作為邊遠(yuǎn)落后地區(qū),水資源統(tǒng)計(jì)工作開始更晚,再考慮到數(shù)據(jù)的可得性和可靠性,本文將研究的樣本區(qū)間確定為2000—2014年。模型所需的新疆歷年GDP、農(nóng)業(yè)產(chǎn)值、工業(yè)產(chǎn)值、用水總量、農(nóng)業(yè)用水量以及工業(yè)用水量均來源于《新疆統(tǒng)計(jì)年鑒》(2001—2015)。為避免數(shù)據(jù)的劇烈波動(dòng),消除可能存在的異方差,對(duì)各指標(biāo)進(jìn)行了對(duì)數(shù)化處理。

    2 模型結(jié)果與分析

    2.1 模型建立及基本檢驗(yàn)

    (1)單位根檢驗(yàn)

    要想建立VAR模型,必須滿足所有的序列都是平穩(wěn)的,為了避免偽回歸現(xiàn)象,保證結(jié)論的可靠性,首先對(duì)數(shù)據(jù)進(jìn)行ADF檢驗(yàn)(表1)。

    首先,對(duì)新疆GDP(LOGY1)和用水總量(LOGX1)進(jìn)行檢驗(yàn)時(shí),均不顯著,接受存在單位根的原假設(shè),說明這兩個(gè)變量的水平序列都是非平穩(wěn)序列;將LOGY1變量進(jìn)行一階差分(DLOGY1),LOGX1變量進(jìn)行二階差分(DDLOGX1),檢驗(yàn)結(jié)果顯示均顯著地拒絕原假設(shè),認(rèn)為兩個(gè)變量均為平穩(wěn)序列。但是,考慮到只有同階平穩(wěn)才能避免偽回歸,因此這里對(duì)LOGY1進(jìn)行二階差分,得到平穩(wěn)序列。由此可知,新疆GDP和用水總量均為二階單整序列。

    其次,對(duì)新疆農(nóng)業(yè)產(chǎn)值(LOGY2)、農(nóng)業(yè)用水量(LOGX2)、新疆工業(yè)產(chǎn)值(LOGY3)、工業(yè)用水量(LOGX3)進(jìn)行檢驗(yàn),4個(gè)變量均不顯著,都是非平穩(wěn)序列;分別對(duì)4個(gè)變量的一階差分進(jìn)行檢驗(yàn),均顯著地拒絕原假設(shè),認(rèn)為差分后的4個(gè)變量均為平穩(wěn)序列。由此可知,新疆農(nóng)業(yè)產(chǎn)值、農(nóng)業(yè)用水量、工業(yè)產(chǎn)值、工業(yè)用水量均為一階單整序列。

    表1 變量的ADF檢驗(yàn)結(jié)果

    ADF:單位根檢驗(yàn),Augmented Dickey-Fuller Test;LOGX1:對(duì)X1求對(duì)數(shù),log(x1);LOGX2:對(duì)X2求對(duì)數(shù),log(x2);LOGX3:對(duì)X3求對(duì)數(shù),log(x3);LOGY1:對(duì)Y1求對(duì)數(shù),log(Y1);LOGY2:對(duì)Y2求對(duì)數(shù),log(Y2);LOGY3:對(duì)Y3求對(duì)數(shù),log(Y3);DDLOGX1:對(duì)LOGX1取二階差分;DLOGX2:對(duì)LOGX2取一階差分;DLOGX3:對(duì)LOGX3取一階差分;DDLOGY1:對(duì)LOGY1取二階差分;DLOGY2:對(duì)LOGY2取一階差分;DLOGY3:對(duì)LOGY3取一階差分

    (2)滯后期的選擇

    滯后期的選擇是VAR模型建立的重要前提,確定了滯后期,才能確定VAR模型形式。這里分別從新疆經(jīng)濟(jì)增長與用水總量、農(nóng)業(yè)經(jīng)濟(jì)增長與農(nóng)業(yè)用水量、工業(yè)經(jīng)濟(jì)增長與工業(yè)用水量3個(gè)方面分析滯后期。檢驗(yàn)的依據(jù)是按照少數(shù)服從多數(shù)原則,根據(jù)LogL、LR、FPE、AIC、SC、HQ六種檢驗(yàn)法則的結(jié)果進(jìn)行綜合判斷。

    總體來看,新疆經(jīng)濟(jì)增長與用水總量的最佳滯后期為1;農(nóng)業(yè)經(jīng)濟(jì)增長與農(nóng)業(yè)用水量、工業(yè)經(jīng)濟(jì)增長與工業(yè)用水量3個(gè)VAR模型的最佳滯后期均為3。

    (3)VAR模型的建立

    根據(jù)前述單位根檢驗(yàn)及滯后期選擇結(jié)果,將檢驗(yàn)過的平穩(wěn)序列代入到模型中,利用最小二乘法估計(jì)模型中的參數(shù),擬合出VAR模型的系數(shù)矩陣:

    DLOGY1=-0.80 DDLOGY1(t-1)-0.30 DDLOGX1(t-1)-0.001 C

    (3)

    DLOGY2=-0.47 DLOGY2(t-1)-0.81 DLOGY2(t-2)-0.004 DLOGY2(t-3) -2.74DLOGX2(t-1)- 1.41DLOGX2(t-2)-3.18 DLOGX2(t-3) + 0.40 C

    (4)

    DLOGY3=-0.27 DLOGY3(t-1)-0.16 DLOGY3(t-2)-0.13 DLOGY3(t-3) + 0.18 DLOGX3(t-1) - 0.05DLOGX3(t-2)-0.33 DLOGX3(t-3) + 0.25 C

    (5)

    上述模型中,式(3)代表經(jīng)濟(jì)增長與用水總量的VAR模型;式(4)代表農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量的VAR模型;式(5)代表工業(yè)產(chǎn)值與工業(yè)用水量的VAR模型。

    表2 VAR模型最佳滯后階數(shù)檢驗(yàn)結(jié)果

    LR:似然比,Likelihood Ratio;FPE:最終預(yù)測誤差準(zhǔn)則,F(xiàn)inal Prediction Error;AIC:赤池信息準(zhǔn)則,Akaike Information Criterion;SC:施瓦茨準(zhǔn)則,Schwarz Criterion;HQ:漢南-奎因準(zhǔn)則Hannan-Quinn Criterion;NA:不適用,Not Applicable;*:表示對(duì)應(yīng)檢驗(yàn)的最佳滯后期選擇

    從表3可以看出,VAR模型的3個(gè)單方程檢驗(yàn)結(jié)果都比較顯著,AIC和SC值都很低,擬合效果較好。

    表3 VAR模型的整體檢驗(yàn)結(jié)果

    (4)VAR模型的穩(wěn)定性檢驗(yàn)

    由表4可知所有根的模全部小于1,說明本研究建立的新疆經(jīng)濟(jì)增長與用水總量、農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量、工業(yè)產(chǎn)值與工業(yè)用水量3個(gè)VAR模型是穩(wěn)定的,基于該VAR模型上的各種檢驗(yàn)是有效的。

    表4 VAR模型的穩(wěn)定性檢驗(yàn)結(jié)果

    i:表示虛根,Imaginary Root;-:表示不含有數(shù)據(jù)

    (5)Granger因果檢驗(yàn)

    因果關(guān)系檢驗(yàn)是用于檢驗(yàn)經(jīng)濟(jì)變量之間的時(shí)間先后順序,并不表示這真正存在因果關(guān)系,是否存在因果關(guān)系需要根據(jù)理論、經(jīng)驗(yàn)和模型進(jìn)行綜合判斷。為了分別考察經(jīng)濟(jì)增長與用水總量、農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量、工業(yè)產(chǎn)值與工業(yè)用水量之間的傳導(dǎo)方向,對(duì)三組序列分別做Granger因果關(guān)系檢驗(yàn)(表5)。

    通過Granger因果關(guān)系檢驗(yàn)結(jié)果可知,從DDLOGX1對(duì)DDLOGY1的因果關(guān)系來看,P值為0.0174,小于0.05,F統(tǒng)計(jì)量落在拒絕域,拒絕原假設(shè),接受備擇假設(shè),說明用DDLOGX1的平穩(wěn)序列能夠引起DLOGY1的平穩(wěn)變換,用水總量是GDP的Granger原因,說明水資源的開發(fā)利用極大地推動(dòng)了新疆GDP的快速增長。從DDLOGY1對(duì)DDLOGX1的因果關(guān)系來看,P值為0.0626,大于0.05而小于0.1,F統(tǒng)計(jì)量落在拒絕域,拒絕原假設(shè),接受備擇假設(shè),說明DDLOGY1的平穩(wěn)序列能夠引起DDLOGX1的平穩(wěn)變換,即DDLOGY1是DDLOGX1的Granger原因,說明GDP增長帶動(dòng)了新疆用水總量的不斷增加,是造成新疆水資源消耗的重要原因。綜上所述,GDP和用水總量之間存在雙向的因果關(guān)系,GDP的變動(dòng)會(huì)引起用水總量的變化。

    表5 Granger因果關(guān)系檢驗(yàn)結(jié)果

    通過Granger因果關(guān)系檢驗(yàn)結(jié)果可知,從DLOGX2對(duì)DLOGY2的因果關(guān)系來看,P值為0.3834,大于0.1,F統(tǒng)計(jì)量落在接受域,接受原假設(shè),即DLOGX2的平穩(wěn)序列不能夠引起DLOGY2的平穩(wěn)變換,說明農(nóng)業(yè)用水不是農(nóng)業(yè)產(chǎn)值的Granger原因。從DLOGY2對(duì)DLOGX2的因果關(guān)系來看,P值為0.1021,大于0.1,F統(tǒng)計(jì)量落在接受域,接受原假設(shè),拒絕備擇假設(shè),說明DLOGY2的平穩(wěn)序列不能夠引起DLOGX2的平穩(wěn)變換,即農(nóng)業(yè)產(chǎn)值不是農(nóng)業(yè)用水量的Granger原因。

    從DLOGX3對(duì)DLOGY3的因果關(guān)系來看,P值為0.8983,大于0.1,F統(tǒng)計(jì)量落在接受域,接受原假設(shè),說明DLOGX3的平穩(wěn)序列不能夠引起DLOGY3的平穩(wěn)變換,即工業(yè)用水總量不是工業(yè)產(chǎn)值的Granger原因。從DLOGY3對(duì)DLOGX3的因果關(guān)系來看,P值為0.3387,大于0.1,F統(tǒng)計(jì)量落在接受域,接受原假設(shè),拒絕備擇假設(shè),說明DLOGY3的平穩(wěn)序列不是DLOGX3的平穩(wěn)變換,即工業(yè)產(chǎn)值也不是工業(yè)用水總量的Granger原因。

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

    本部分采用脈沖響應(yīng)函數(shù)來分析經(jīng)濟(jì)增長與用水量之間的引致關(guān)系。脈沖響應(yīng)函數(shù)描述的是VAR模型中的一個(gè)內(nèi)生變量的沖擊給其他內(nèi)生變量所帶來的影響,是系統(tǒng)對(duì)某一變量擾動(dòng)的一個(gè)沖擊(或新息)所做出的動(dòng)態(tài)反應(yīng),并由此判斷各個(gè)變量間的時(shí)滯關(guān)系。

    2.2.1 經(jīng)濟(jì)增長與用水總量的脈沖響應(yīng)

    圖1表示了用水總量、GDP兩個(gè)指標(biāo)對(duì)自身及其相互的沖擊效應(yīng)以及系統(tǒng)的動(dòng)態(tài)影響。

    從DDLOGY1對(duì)DDLOGX1的一個(gè)標(biāo)準(zhǔn)差新息的沖擊,DDLOGX1在第1期為0,第2期為正值,第3期變?yōu)樨?fù)值,第4期和第5期為正值,正負(fù)影響也是交替出現(xiàn),但沖擊的幅度越來越小,并最終區(qū)域收斂。說明隨著時(shí)間的推移,新疆經(jīng)濟(jì)增長對(duì)水資源的依賴越來越弱。

    從DDLOGX1對(duì)DDLOGY1的一個(gè)標(biāo)準(zhǔn)差新息的沖擊,第1期就具有非常大的影響,接近0.06,但第2期就變?yōu)樨?fù)值,達(dá)-0.1,第3期為正值,第4期和5期又將為負(fù)值,呈現(xiàn)不規(guī)則的正負(fù)交替形態(tài),但沖擊響應(yīng)的絕對(duì)值越來越小,并最終收斂。說明伴隨著經(jīng)濟(jì)的不斷發(fā)展,水資源的進(jìn)一步開發(fā)對(duì)經(jīng)濟(jì)增長帶來的促進(jìn)作用越來越弱。

    圖1 DDLOGY1與DDLOGX1的脈沖響應(yīng)函數(shù)圖Fig.1 DDLOGY1 and DDLOGX1 impulse response function diagram

    2.2.2 農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量的脈沖響應(yīng)

    由圖2可知,DLOGY2對(duì)DLOGX2的一個(gè)標(biāo)準(zhǔn)差新息的沖擊,初期影響為0,第2期就變?yōu)樨?fù)值,且達(dá)到-0.04,第3期為0.001,第4期開始為正值,并持續(xù)到第6期,第7期開始又變?yōu)樨?fù)影響,并持續(xù)到第9期,第10期重又呈正值。由此可見,DLOGY2對(duì)DLOGX2的影響不平穩(wěn),正負(fù)影響交替出現(xiàn),呈現(xiàn)出“余弦型”波動(dòng)特征,但影響趨于越大,表明農(nóng)業(yè)增產(chǎn)增收目標(biāo)驅(qū)動(dòng)下,農(nóng)業(yè)用水將不斷得到強(qiáng)化。

    DLOGX2對(duì)DLOGY2的一個(gè)標(biāo)準(zhǔn)差新息的沖擊,第1期和第2期均為負(fù)值,也均為-0.01,第2期和第4期變?yōu)檎?第5—6期為負(fù)值,第7—8期又變?yōu)檎?正負(fù)影響交替出現(xiàn),但影響逐漸縮小,最終表現(xiàn)為正值,表明隨著農(nóng)業(yè)現(xiàn)代化水平的提高,水資源對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長的作用呈減弱趨勢。

    圖2 DLOGY2對(duì)DLOGX2的脈沖響應(yīng)函數(shù)圖Fig.2 DLOGY2 on the DLOGX2 impulse response function diagram

    圖3 DLOGY3對(duì)DLOGX3的脈沖響應(yīng)函數(shù)圖Fig.3 DLOGY3 on the DLOGX3 impulse response function diagram

    由圖3可知,DLOGY3對(duì)DLOGX3的一個(gè)標(biāo)準(zhǔn)差新息的沖擊,初期影響為0,第2期為正值,為0.005,第3期變?yōu)樨?fù)值,為-0.01,持續(xù)到第4期,第5期轉(zhuǎn)變?yōu)檎?第6期開始又變?yōu)樨?fù)值,交替出現(xiàn)至第10期,并最終表現(xiàn)為正影響。雖然DLOGY3對(duì)DLOGX3的影響呈現(xiàn)正負(fù)交替,但主要在零值附近波動(dòng),所產(chǎn)生的影響較小。表明工業(yè)產(chǎn)值的增加對(duì)工業(yè)用水量增加的帶動(dòng)作用較小,其原因可能同當(dāng)前新疆工業(yè)用水量長期不足有關(guān)。

    DLOGX3對(duì)DLOGY3一個(gè)標(biāo)準(zhǔn)差新息的沖擊,第1期為正值,且為最大值0.03,第2期變?yōu)樨?fù)值,為-0.02,之后交替正負(fù)影響至第4期,第5—7期均為正影響,第8期為負(fù)值,并交替至第10期,最終為負(fù)影響??傮w來看,DLOGX3對(duì)DLOGY3在初期就具有較大的正向沖擊,之后開始波動(dòng)下降,并最終收斂。表明水資源開發(fā)對(duì)經(jīng)濟(jì)增長的影響在中短期較為強(qiáng)烈,長期來看,其影響呈減弱趨勢,這可能同工業(yè)各行業(yè)的技術(shù)進(jìn)步有關(guān)。

    2.3 方差分解分析

    方差分解是一種描述系統(tǒng)動(dòng)態(tài)變化的方法,通過分析每一個(gè)結(jié)構(gòu)沖擊對(duì)內(nèi)生變量變化的貢獻(xiàn)度,通常以方差來度量,進(jìn)一步評(píng)價(jià)不同結(jié)構(gòu)沖擊的相對(duì)重要性。

    2.3.1 經(jīng)濟(jì)增長與用水總量預(yù)期誤差的方差分解

    由表6可知,GDP對(duì)水資源利用的方差分解的貢獻(xiàn)度較低,從首期來看,GDP對(duì)水資源利用方差分解的貢獻(xiàn)度為0,但第2期就增加到10.98%,并繼續(xù)增加,第3期開始保持在20%以上,并于第8期增至最大值21.93%。表明新疆經(jīng)濟(jì)快速增長的過程是對(duì)水資源的大力開發(fā)和過度利用的過程,這與新疆經(jīng)濟(jì)增長的現(xiàn)實(shí)是相符的。

    水資源利用對(duì)GDP增長的方差分解貢獻(xiàn)度很高,除前1期的貢獻(xiàn)度低于30%外,第2期開始,水資源利用對(duì)GDP增長方差分解的貢獻(xiàn)度高接近并超過60%,表明水資源利用對(duì)新疆經(jīng)濟(jì)增長起到較大促進(jìn)作用。相對(duì)于我國其他省區(qū)經(jīng)濟(jì)增長的多樣化因素而言,水資源對(duì)新疆經(jīng)濟(jì)可持續(xù)發(fā)展至關(guān)重要。

    2.3.2 農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量預(yù)期誤差的方差分解

    由表7可知,相較于農(nóng)業(yè)用水對(duì)農(nóng)業(yè)產(chǎn)值而言,農(nóng)業(yè)產(chǎn)值對(duì)農(nóng)業(yè)用水的方差分解的貢獻(xiàn)度較低,但總體水平仍處于較高水平,自第2期開始,農(nóng)業(yè)產(chǎn)值對(duì)農(nóng)業(yè)用水方差分解貢獻(xiàn)度均在20%左右,即農(nóng)業(yè)產(chǎn)值可以解釋農(nóng)業(yè)用水20%以上的方差,且貢獻(xiàn)度呈現(xiàn)波動(dòng)變化趨勢。一方面,伴隨著新疆農(nóng)業(yè)產(chǎn)值的快速提升,對(duì)前中期農(nóng)業(yè)用水造成了大量的引致需求,農(nóng)業(yè)用水不斷增加;另一方面,隨著新疆農(nóng)業(yè)生產(chǎn)節(jié)水灌溉的廣泛推廣及農(nóng)業(yè)現(xiàn)代化程度的提高,后期農(nóng)業(yè)產(chǎn)值對(duì)農(nóng)業(yè)用水的方差分解貢獻(xiàn)度呈下降趨勢。

    表6 經(jīng)濟(jì)增長與用水總量預(yù)期誤差的方差分解結(jié)果

    表7 農(nóng)業(yè)產(chǎn)值與農(nóng)業(yè)用水量預(yù)期誤差的方差分解結(jié)果

    水資源利用對(duì)農(nóng)業(yè)產(chǎn)值增長的方差分解貢獻(xiàn)度很高,除前2期的貢獻(xiàn)度處于30%以下外,第3期開始,水資源利用對(duì)GDP增長方差分解的貢獻(xiàn)度不斷提高,并于第6期達(dá)最大值,為71.21%??梢钥闯?水資源利用對(duì)新疆農(nóng)業(yè)的發(fā)展起到了至關(guān)重要的作用,水資源約束是主要的控制性因素,這與新疆干旱區(qū)綠洲農(nóng)業(yè)的特色非常吻合。

    2.3.3 工業(yè)產(chǎn)值與工業(yè)用水量預(yù)期誤差的方差分解

    由表8可知,相較于工業(yè)用水對(duì)工業(yè)產(chǎn)值而言,工業(yè)產(chǎn)值對(duì)工業(yè)用水的方差分解貢獻(xiàn)度非常低,自首期至末期的影響均不到1%,即工業(yè)產(chǎn)值所能解釋的農(nóng)業(yè)用水的方差還不到1%。表明伴隨著新疆工業(yè)產(chǎn)值的快速提升,并沒有帶動(dòng)新疆工業(yè)用水量的大幅增加,水資源可能在未來成為制約新疆工業(yè)發(fā)展的主要因素。

    工業(yè)用水對(duì)工業(yè)產(chǎn)值增長的方差分解貢獻(xiàn)度很高,首期貢獻(xiàn)度就高達(dá)74.97%,此后雖有波動(dòng),但貢獻(xiàn)度仍在65%以上,且自第4期開始,貢獻(xiàn)度一直在80%以上并不斷提高。由此可以看出,工業(yè)用水對(duì)新疆工業(yè)產(chǎn)值的增加具有非常大的促進(jìn)作用,其作用已經(jīng)達(dá)到了80%以上,工業(yè)用水與工業(yè)產(chǎn)值具有長期穩(wěn)定的均衡關(guān)系,這與當(dāng)前新疆工業(yè)發(fā)展中的水資源利用現(xiàn)狀是完全相符的。

    表8 工業(yè)產(chǎn)值與工業(yè)用水量預(yù)期誤差的方差分解結(jié)果

    3 結(jié)論及建議

    本文基于2000—2014年新疆水資源利用與經(jīng)濟(jì)增長的關(guān)聯(lián)視角,采用VAR模型,驗(yàn)證新疆大尺度產(chǎn)業(yè)生態(tài)系統(tǒng)中的用水與經(jīng)濟(jì)增長的協(xié)整關(guān)系,進(jìn)一步利用脈沖響應(yīng)和方差分解對(duì)新疆產(chǎn)業(yè)用水與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)系進(jìn)行擬合,得到如下結(jié)論:

    (1)用水總量與經(jīng)濟(jì)增長存在穩(wěn)定的協(xié)整關(guān)系,最佳滯后期為1期。格蘭杰因果關(guān)系檢驗(yàn)來看,GDP的平穩(wěn)序列能夠引起用水總量的平穩(wěn)變換,表明新疆經(jīng)濟(jì)增長過程中對(duì)用水量造成引致需求;但用水總量的平穩(wěn)序列不能夠引起GDP平穩(wěn)變換,用水量對(duì)經(jīng)濟(jì)增長的單一推動(dòng)效應(yīng)并不強(qiáng)烈。就脈沖響應(yīng)來看,經(jīng)濟(jì)增長對(duì)用水總量的沖擊響應(yīng)累計(jì)值為負(fù)值(-0.01),用水總量對(duì)經(jīng)濟(jì)增長的沖擊響應(yīng)累計(jì)值也為負(fù)值(-0.01)。GDP對(duì)水資源利用的方差分解的貢獻(xiàn)度較低,不到4%,水資源利用對(duì)GDP增長的方差分解貢獻(xiàn)度很高,超過90%。由此可知,隨著新疆經(jīng)濟(jì)的快速增長,對(duì)用水總量造成了巨大的引致需求,水資源消耗量不斷增加。表明水資源利用對(duì)新疆經(jīng)濟(jì)增長起到了非常大的促進(jìn)作用。相對(duì)于我國其他省區(qū)經(jīng)濟(jì)增長的多樣化因素而言,水資源對(duì)新疆經(jīng)濟(jì)可持續(xù)發(fā)展至關(guān)重要。未來一段時(shí)期,全面建成小康社會(huì)過程中,各方主體將會(huì)再次推動(dòng)新疆經(jīng)濟(jì)快速發(fā)展,這將對(duì)本就有限的水資源造成更大需求,水資源利用與經(jīng)濟(jì)增長之間的矛盾有可能進(jìn)一步加劇。

    (2)農(nóng)業(yè)生態(tài)系統(tǒng)用水量與農(nóng)業(yè)產(chǎn)值存在穩(wěn)定的協(xié)整關(guān)系,最佳滯后期為3期。兩者的平穩(wěn)序列均不能引起對(duì)方的平穩(wěn)變換,可能與新疆農(nóng)業(yè)用水在用水總量中的絕對(duì)地位導(dǎo)致其對(duì)農(nóng)業(yè)產(chǎn)值的影響不顯著。農(nóng)業(yè)產(chǎn)值對(duì)農(nóng)業(yè)用水量的沖擊響應(yīng)累計(jì)值為負(fù)值(-0.03),農(nóng)業(yè)用水量對(duì)農(nóng)業(yè)產(chǎn)值的沖擊響應(yīng)累計(jì)值為正值(0.01)。農(nóng)業(yè)產(chǎn)值對(duì)農(nóng)業(yè)用水的方差分解的貢獻(xiàn)度較低,在20%以上;水資源利用對(duì)農(nóng)業(yè)產(chǎn)值增長的方差分解貢獻(xiàn)度很高,均在60%以上。可以看出,水資源對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長具有較大的正向促進(jìn)作用,但農(nóng)業(yè)產(chǎn)值的增加對(duì)水資源利用構(gòu)成較大壓力。作為我國重要的“糧、棉、果、畜和特色農(nóng)業(yè)生產(chǎn)基地”,未來新疆農(nóng)業(yè)生產(chǎn)的強(qiáng)度仍將進(jìn)一步強(qiáng)化,這對(duì)水資源日益緊張的新疆影響重大,大力推行節(jié)水農(nóng)業(yè)、強(qiáng)化產(chǎn)業(yè)生態(tài)系統(tǒng)綜合用水平衡將是今后新疆農(nóng)業(yè)科技甚至農(nóng)業(yè)產(chǎn)業(yè)化的重點(diǎn)內(nèi)容,也將成為新疆綠洲農(nóng)業(yè)現(xiàn)代化的重要方向。

    (3)工業(yè)系統(tǒng)用水與工業(yè)產(chǎn)值存在穩(wěn)定的協(xié)整關(guān)系,最佳滯后期為3期。兩者的平穩(wěn)序列均不能引起對(duì)方的平穩(wěn)變換,這可能與新疆工業(yè)產(chǎn)業(yè)水資源供應(yīng)不足且空間分布不協(xié)調(diào)有關(guān)。工業(yè)產(chǎn)值對(duì)工業(yè)用水量的沖擊響應(yīng)累計(jì)值為負(fù)值(-0.001),工業(yè)用水量對(duì)工業(yè)產(chǎn)值的沖擊響應(yīng)累計(jì)值為正值(0.02)。工業(yè)產(chǎn)值對(duì)工業(yè)用水的方差分解貢獻(xiàn)度非常低,首期至末期均不到1%;工業(yè)用水對(duì)工業(yè)產(chǎn)值增長的方差分解貢獻(xiàn)度很高,達(dá)80%以上,可以看出,工業(yè)用水對(duì)工業(yè)產(chǎn)值的增加具有非常強(qiáng)的促進(jìn)作用,但工業(yè)產(chǎn)值對(duì)水資源也帶來了很大壓力,主要工業(yè)城市的水資源開發(fā)利用程度已接近或超過100%,地表水過度引用和地下水超采問題普遍存在[35]。面對(duì)發(fā)展日趨活躍,規(guī)模逐漸攀升,結(jié)構(gòu)更加完善的新疆工業(yè)體系,水資源條件與工業(yè)產(chǎn)業(yè)發(fā)展不匹配問題有可能成為未來新疆工業(yè)發(fā)展的重要瓶頸。

    由以上結(jié)論可知,新疆目前水資源利用與經(jīng)濟(jì)增長存在較強(qiáng)的相關(guān)性及動(dòng)態(tài)關(guān)系。伴隨著未來新疆經(jīng)濟(jì)的快速發(fā)展,用水需求將會(huì)繼續(xù)增加,對(duì)水資源的可持續(xù)利用將會(huì)造成巨大障礙,如何實(shí)現(xiàn)經(jīng)濟(jì)增長與水資源利用增加的脫鉤或弱化其強(qiáng)相關(guān)關(guān)系,是目前及未來新疆經(jīng)濟(jì)發(fā)展中要考慮的重要議題。要實(shí)現(xiàn)水依賴型區(qū)域生態(tài)系統(tǒng)的永續(xù)發(fā)展,從產(chǎn)業(yè)結(jié)構(gòu)入手考量其對(duì)生態(tài)環(huán)境的維護(hù)功能尤為重要,需要從以下方面著手:一是實(shí)施嚴(yán)格的退地減水政策。目前新疆農(nóng)業(yè)用水占比過高,大大超過用水總量的控制指標(biāo),如果不實(shí)施退地減水政策,生態(tài)用水將會(huì)繼續(xù)被大量擠占,并且使得工業(yè)發(fā)展無水可用。二是明確用水總量控制和定額指標(biāo),嚴(yán)格控制用水量,提高用水效率??傮w來看,目前新疆用水總量超標(biāo),主要由于用水效率偏低,應(yīng)加快制定農(nóng)業(yè)、工業(yè)用水效率標(biāo)準(zhǔn),將用水效率作為項(xiàng)目核準(zhǔn)和備案的依據(jù)。三是加強(qiáng)跨流域調(diào)水工程建設(shè)力度。新疆水資源區(qū)域分布不均衡是導(dǎo)致區(qū)域性缺水的重要原因之一,可以借鑒我國南水北調(diào)案例,加強(qiáng)跨流域調(diào)水,從而滿足農(nóng)業(yè)、工業(yè)生產(chǎn)用水。四是防止水資源浪費(fèi)。新疆地處內(nèi)陸,河流均為內(nèi)陸河,工農(nóng)業(yè)生產(chǎn)中所產(chǎn)生的污染物難以通過徑流向外轉(zhuǎn)移。因此需要加強(qiáng)水資源節(jié)約,加大治污和生態(tài)保護(hù)力度,維持新疆生態(tài)環(huán)境健康發(fā)展,并大力加強(qiáng)污水治理工作,防止因污水或其他水源浪費(fèi)而造成的人為缺水問題。

    [1] 羅蘭, 丁曉娟, 馬琴, 代燕. 新疆地州市水資源與經(jīng)濟(jì)社會(huì)可持續(xù)協(xié)調(diào)發(fā)展研究. 中國城市經(jīng)濟(jì), 2011, (30): 307-309.

    [2] 黃寶連, 黃海平, 王生貴. 新疆水資源及其產(chǎn)業(yè)間用水結(jié)構(gòu)特征研究. 石河子大學(xué)學(xué)報(bào): 哲學(xué)社會(huì)科學(xué)版, 2012, 26(5): 16-22.

    [3] 謝書玲, 王錚, 薛俊波. 中國經(jīng)濟(jì)發(fā)展中水土資源的“增長尾效”分析. 管理世界, 2005, (7): 22-25.

    [4] Liu K K, Li C H, Yang X L, Hu J, Xia X H. Water resources supply-consumption (demand) balance analyses in the Yellow River Basin in 2009. Procedia Environmental Sciences, 2012, 13: 1956-1965.

    [5] Guo Y, Shen Y J. Agricultural water supply/demand changes under projected future climate change in the arid region of northwestern China. Journal of Hydrology, 2016, 540: 257-273.

    [6] Shang Y Z, Lu S B, Li X F, Hei P F, Lei X H, Gong J G, Liu J H, Zhai J Q, Wang H. Balancing development of major coal bases with available water resources in China through 2020. Applied Energy, 2016, doi: 10.1016/j.apenergy.2016.07.002.

    [7] Ji X B, Kang E S, Chen R S, Zhao W Z, Xiao S C, Jin B W. Analysis of water resources supply and demand and security of water resources development in irrigation regions of the middle reaches of the Heihe River Basin, Northwest China. Agricultural Sciences in China, 2006, 5(2): 130-140.

    [8] 龍騰銳, 姜文超, 何強(qiáng). 水資源承載力內(nèi)涵的新認(rèn)識(shí). 水利學(xué)報(bào), 2004, 35(1): 38-45.

    [9] Ait-Aoudia M N, Berezowska-Azzag E. Water resources carrying capacity assessment: The case of Algeria's capital city. Habitat International, 2016, 58: 51-58

    [10] 姚治君, 王建華, 江東, 陳傳友. 區(qū)域水資源承載力的研究進(jìn)展及其理論探析. 水科學(xué)進(jìn)展, 2002, 13(1): 111-115.

    [11] Davijani M H, Banihabib M E, Anvar A N, Hashemi S R. Optimization model for the allocation of water resources based on the maximization of employment in the agriculture and industry sectors. Journal of Hydrology, 533: 430-438.

    [12] Hu Z N, Chen Y Z, Yao L M, Wei C T, Li C Z. Optimal allocation of regional water resources: From a perspective of equity-efficiency tradeoff. Resources, Conservation and Recycling, 109: 102-113.

    [13] Mianabadi H, Mostert E, Zarghami M, Van De Giesen N. A new bankruptcy method for conflict resolution in water resources allocation. Journal of Environmental Management, 144: 152-159.

    [14] Varis O, Vakkilainen P. China's challenges to water resources management in the first quarter of the 21st Century. Geomorphology, 2001, 41(2-3): 93-104.

    [15] Cai Y P, Yue W C, Xu L Y, Yang Z F, Rong Q Q. Sustainable urban water resources management considering life-cycle environmental impacts of water utilization under uncertainty. Resources, Conservation and Recycling, 2016, 108: 21-40.

    [16] Kotir J H, Smith C, Brown G, Marshall N, Johnstone R. A system dynamics simulation model for sustainable water resources management and agricultural development in the Volta River Basin, Ghana. Science of the Total Environment, 573: 444-457.

    [17] 張華, 張?zhí)m, 趙傳燕. 極端干旱區(qū)尾閭湖生態(tài)需水估算——以東居延海為例. 生態(tài)學(xué)報(bào), 2014, 34(8): 2102-2108.

    [18] 孫棟元, 楊俊, 胡想全, 金彥兆, 張?jiān)屏? 基于生態(tài)保護(hù)目標(biāo)的疏勒河中游綠洲生態(tài)環(huán)境需水研究. 生態(tài)學(xué)報(bào), 2017, 37(3).

    [19] Pan B Z, Wang H Z, Ban X, Yin X A. An exploratory analysis of ecological water requirements of macroinvertebrates in the Wuhan branch of the Yangtze River. Quaternary International, 2015, 380-381: 256-261.

    [20] 王薇, 雷學(xué)東, 余新曉, 陳麗華. 基于SD模型的水資源承載力計(jì)算理論研究——以青海共和盆地水資源承載力研究為例. 水資源與水工程學(xué)報(bào), 2005, 16(3): 11-15.

    [21] Howe C W. Water resources and regional economic growth in the United States, 1950-1960. Southern Economic Journal, 1968, 34(4): 477-489.

    [22] 路寧, 周海光. 中國城市經(jīng)濟(jì)與水資源利用壓力的關(guān)系研究. 中國人口·資源與環(huán)境, 2010, 20(S2): 48-50.

    [23] 李青, 陳紅梅, 王雅鵬. 基于面板VAR模型的新疆農(nóng)業(yè)用水與農(nóng)業(yè)經(jīng)濟(jì)增長的互動(dòng)效應(yīng)研究. 資源科學(xué), 2014, 36(8): 1679-1685.

    [24] 方創(chuàng)琳, 鮑超, 張傳國. 干旱地區(qū)生態(tài)-生產(chǎn)-生活承載力變化情勢與演變情景分析. 生態(tài)學(xué)報(bào), 2003, 23(9): 1915-1923.

    [25] Sim C A. Macroeconomics and reality. Econometrica, 1980, 48(1): 1-48.

    [26] Holtz-Eakin D, Newey W K, Rosen H S. Estimating vector autoregressions with panel data. Econometrica, 1988, 56(6): 1371-1395.

    [27] Hahn J, Kuersteiner G. Asymptotically unbiased inference for a dynamic panel model with fixed effects when bothnandTare large. Econometrica, 2002, 70(4): 1639-1657.

    [28] Binder M, Hsiao C, Pesaran M H. Estimation and inference in short panel vector autoregressions with unit roots and cointegration. Econometric Theory, 2005, 21(4): 795-837.

    [29] Chong A, Gradstein M. Inequality and institutions. The Review of Economics and Statistics, 2007, 89(3): 454-465.

    [30] Aspachs O, Goodhart C A E, Tsomocos D P, Zicchino L. Towards a measure of financial fragility. Annals Finance, 2007, 3(1): 37-74.

    [31] Huang B N, Hwang M J, Yang C W. Causal relationship between energy consumption and GDP growth revisited: A dynamic panel data approach. Ecological Economics, 2008, 67(1): 41-54.

    [32] 王舒健, 李釗. 中國地區(qū)經(jīng)濟(jì)增長互動(dòng)關(guān)系的脈沖響應(yīng)分析. 數(shù)理統(tǒng)計(jì)與管理, 2007, 26(3): 385-390.

    [33] 孫雪蓮, 鄧峰. 干旱區(qū)水資源對(duì)經(jīng)濟(jì)增長的約束作用實(shí)證分析——以新疆為例. 新疆社會(huì)科學(xué), 2013, (2): 43-47.

    [34] 汪奎, 邵東國, 顧文權(quán), 岑棟浩, 譚學(xué)智, 楊豐順. 中國用水量與經(jīng)濟(jì)增長的脫鉤分析. 灌溉排水學(xué)報(bào), 2011, 30(3): 34-38.

    [35] 鄧朝暉, 劉洋, 薛惠鋒. 基于VAR模型的水資源利用與經(jīng)濟(jì)增長動(dòng)態(tài)關(guān)系研究. 中國人口·資源環(huán)境, 2012, 22(6): 128-135.

    Researchonthedynamicrelationshipbetweenwater-intensiveindustrialecosystemandeconomicgrowthinXinjiangbasedonVARmodel

    ZHANG Zhenlong1,2, SUN Hui1,2,*

    1CenterforInnovationManagementResearchofXinjiangUniversity,Urumqi830046,China2SchoolofEconomicsandManagement,XinjiangUniversity,Urumqi830046,China

    國家自然科學(xué)基金項(xiàng)目(71463056,41461114);新疆維吾爾自治區(qū)自然科學(xué)基金項(xiàng)目(2016D01B054)

    2016- 05- 19; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 03- 27

    *通訊作者Corresponding author.E-mail: 18999921777@163.com

    10.5846/stxb201605190966

    張振龍,孫慧.新疆區(qū)域水資源對(duì)產(chǎn)業(yè)生態(tài)系統(tǒng)與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)聯(lián)——基于VAR模型.生態(tài)學(xué)報(bào),2017,37(16):5273- 5284.

    Zhang Z L, Sun H.Research on the dynamic relationship between water-intensive industrial ecosystem and economic growth in Xinjiang based on VAR model.Acta Ecologica Sinica,2017,37(16):5273- 5284.

    猜你喜歡
    貢獻(xiàn)度用水量總量
    “十三五”期間山西省與10省簽約糧食總量2230萬噸
    2020年全國農(nóng)民工總量比上年減少517萬人
    你的用水量是多少?
    澳大利亞研發(fā)出新型農(nóng)業(yè)傳感器可預(yù)測農(nóng)作物用水量
    充分把握教育對(duì)經(jīng)濟(jì)社會(huì)發(fā)展的貢獻(xiàn)度
    基于貢獻(xiàn)度排序的腎透明細(xì)胞癌串?dāng)_通路分析
    為何化肥淡儲(chǔ)總量再度增加
    需求側(cè)資源促進(jìn)可再生能源消納貢獻(xiàn)度綜合評(píng)價(jià)體系
    工業(yè)給排水用水量計(jì)算的重要性
    中國對(duì)外貿(mào)易對(duì)經(jīng)濟(jì)增長貢獻(xiàn)度探索——基于VAR模型的分析
    精品第一国产精品| 国内久久婷婷六月综合欲色啪| 亚洲av片天天在线观看| 黄色 视频免费看| 精品国产一区二区三区四区第35| 久久人妻福利社区极品人妻图片| e午夜精品久久久久久久| 91大片在线观看| 国产成+人综合+亚洲专区| 嫩草影视91久久| 亚洲熟妇中文字幕五十中出 | 欧美丝袜亚洲另类 | 精品久久蜜臀av无| 手机成人av网站| av中文乱码字幕在线| 捣出白浆h1v1| 视频区图区小说| 一级毛片女人18水好多| 少妇粗大呻吟视频| 国产精品.久久久| 日韩一卡2卡3卡4卡2021年| 狂野欧美激情性xxxx| 亚洲av片天天在线观看| 色在线成人网| 在线天堂中文资源库| 高清毛片免费观看视频网站 | 久久国产亚洲av麻豆专区| 欧美日韩黄片免| 亚洲第一青青草原| 一区二区三区精品91| 国产片内射在线| 美女国产高潮福利片在线看| 一区二区日韩欧美中文字幕| 91av网站免费观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 久久精品人人爽人人爽视色| 亚洲免费av在线视频| 在线国产一区二区在线| 亚洲美女黄片视频| 亚洲一码二码三码区别大吗| 高清在线国产一区| 美女视频免费永久观看网站| 午夜福利免费观看在线| 精品欧美一区二区三区在线| 成年版毛片免费区| 国产91精品成人一区二区三区| 18在线观看网站| 欧美精品av麻豆av| 亚洲成人免费电影在线观看| 十八禁网站免费在线| 中文欧美无线码| 久久国产精品人妻蜜桃| 999久久久国产精品视频| 99在线人妻在线中文字幕 | 亚洲色图av天堂| 丝袜美足系列| 精品一区二区三区视频在线观看免费 | 大型黄色视频在线免费观看| 丰满饥渴人妻一区二区三| 777久久人妻少妇嫩草av网站| 欧美日韩av久久| 99国产精品一区二区三区| 天天躁日日躁夜夜躁夜夜| 亚洲精品久久午夜乱码| 亚洲 欧美一区二区三区| 国精品久久久久久国模美| 真人做人爱边吃奶动态| 淫妇啪啪啪对白视频| 国产成人欧美在线观看 | 12—13女人毛片做爰片一| 欧美精品高潮呻吟av久久| 国产欧美亚洲国产| 1024视频免费在线观看| 精品久久久久久电影网| 成熟少妇高潮喷水视频| 在线观看免费午夜福利视频| 欧美一级毛片孕妇| 热re99久久国产66热| 国产亚洲av高清不卡| 免费在线观看日本一区| 国产精品98久久久久久宅男小说| 精品亚洲成a人片在线观看| 99精品欧美一区二区三区四区| 男女下面插进去视频免费观看| 国产麻豆69| 黄色片一级片一级黄色片| 久久这里只有精品19| 免费不卡黄色视频| 国产aⅴ精品一区二区三区波| 亚洲国产欧美网| 亚洲精品成人av观看孕妇| 国产精品影院久久| 黄片小视频在线播放| 亚洲欧美一区二区三区黑人| 在线观看日韩欧美| 午夜久久久在线观看| 国产亚洲精品一区二区www | 超色免费av| 精品无人区乱码1区二区| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 18禁裸乳无遮挡免费网站照片 | 国产蜜桃级精品一区二区三区 | 日本五十路高清| 18禁美女被吸乳视频| 亚洲综合色网址| 午夜激情av网站| 宅男免费午夜| x7x7x7水蜜桃| 99国产精品一区二区三区| 精品一品国产午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 成人精品一区二区免费| 国产一区有黄有色的免费视频| 男女免费视频国产| 老司机福利观看| 久久久国产一区二区| 婷婷丁香在线五月| 久久人人爽av亚洲精品天堂| 在线十欧美十亚洲十日本专区| 超碰97精品在线观看| av线在线观看网站| 日韩 欧美 亚洲 中文字幕| 欧美在线黄色| 高清欧美精品videossex| 黄网站色视频无遮挡免费观看| 999久久久国产精品视频| 亚洲,欧美精品.| 18禁美女被吸乳视频| 午夜老司机福利片| 久久人人爽av亚洲精品天堂| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 日韩熟女老妇一区二区性免费视频| a在线观看视频网站| 法律面前人人平等表现在哪些方面| 欧美日韩国产mv在线观看视频| 女人被躁到高潮嗷嗷叫费观| 国产精品国产av在线观看| 欧美激情高清一区二区三区| 亚洲熟女精品中文字幕| 久久久久视频综合| a级毛片黄视频| 宅男免费午夜| 亚洲成国产人片在线观看| 国产乱人伦免费视频| 高清视频免费观看一区二区| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲熟妇少妇任你| 一二三四社区在线视频社区8| 女警被强在线播放| 又紧又爽又黄一区二区| 亚洲人成77777在线视频| 99久久综合精品五月天人人| 飞空精品影院首页| 国产有黄有色有爽视频| 国内久久婷婷六月综合欲色啪| 丝袜在线中文字幕| 免费人成视频x8x8入口观看| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 久99久视频精品免费| 不卡av一区二区三区| 欧美午夜高清在线| 最新在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 91成年电影在线观看| 精品午夜福利视频在线观看一区| 亚洲专区中文字幕在线| 最新的欧美精品一区二区| 久久 成人 亚洲| 日本欧美视频一区| 午夜精品久久久久久毛片777| 国产在视频线精品| 国产精品国产av在线观看| 欧美日韩一级在线毛片| 侵犯人妻中文字幕一二三四区| 亚洲熟妇熟女久久| 日韩欧美在线二视频 | 国产伦人伦偷精品视频| 久久99一区二区三区| 亚洲国产精品一区二区三区在线| 日韩人妻精品一区2区三区| 制服诱惑二区| e午夜精品久久久久久久| 在线av久久热| 国产1区2区3区精品| 黄色a级毛片大全视频| 757午夜福利合集在线观看| 97人妻天天添夜夜摸| 亚洲 国产 在线| 国产精品影院久久| 日本精品一区二区三区蜜桃| 在线永久观看黄色视频| 国产亚洲精品第一综合不卡| 91成年电影在线观看| 日本一区二区免费在线视频| 国产成人精品久久二区二区免费| 成人国产一区最新在线观看| 亚洲人成77777在线视频| 亚洲人成伊人成综合网2020| 日日摸夜夜添夜夜添小说| 欧美精品亚洲一区二区| 巨乳人妻的诱惑在线观看| 免费高清在线观看日韩| 久久亚洲真实| 亚洲人成伊人成综合网2020| 免费一级毛片在线播放高清视频 | 国产av精品麻豆| 成人手机av| 中文字幕色久视频| 日韩欧美在线二视频 | 亚洲全国av大片| 不卡一级毛片| 亚洲,欧美精品.| 99国产精品一区二区蜜桃av | 亚洲全国av大片| 久久久久久久午夜电影 | 国产男女内射视频| 亚洲va日本ⅴa欧美va伊人久久| a级片在线免费高清观看视频| 色婷婷av一区二区三区视频| 美女视频免费永久观看网站| 久久久久国内视频| 久久精品亚洲精品国产色婷小说| 国产亚洲欧美在线一区二区| 黑人欧美特级aaaaaa片| 一区福利在线观看| 国产精品亚洲av一区麻豆| 动漫黄色视频在线观看| 老鸭窝网址在线观看| xxx96com| 午夜福利一区二区在线看| 精品一区二区三区四区五区乱码| 老汉色∧v一级毛片| 日本a在线网址| 国产精品电影一区二区三区 | 国精品久久久久久国模美| 欧美黑人精品巨大| 成年女人毛片免费观看观看9 | 欧美亚洲 丝袜 人妻 在线| 91九色精品人成在线观看| 男女床上黄色一级片免费看| 一级片'在线观看视频| 久久久久久久国产电影| 国产精品偷伦视频观看了| 99国产精品一区二区三区| 国产三级黄色录像| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 欧美黄色淫秽网站| 纯流量卡能插随身wifi吗| a级毛片在线看网站| 搡老熟女国产l中国老女人| 国产高清激情床上av| 黄色 视频免费看| 日韩制服丝袜自拍偷拍| 午夜福利欧美成人| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| a级毛片黄视频| av在线播放免费不卡| 国产伦人伦偷精品视频| 亚洲色图 男人天堂 中文字幕| 国产高清国产精品国产三级| 欧美成人午夜精品| 日本精品一区二区三区蜜桃| 一区在线观看完整版| 黄色毛片三级朝国网站| 精品亚洲成a人片在线观看| 母亲3免费完整高清在线观看| 法律面前人人平等表现在哪些方面| 午夜免费成人在线视频| 精品熟女少妇八av免费久了| 国产亚洲一区二区精品| 亚洲一区二区三区欧美精品| 色婷婷久久久亚洲欧美| 国产精品自产拍在线观看55亚洲 | 久久香蕉国产精品| 亚洲在线自拍视频| 欧美黄色片欧美黄色片| 久久影院123| 欧美在线黄色| 99久久综合精品五月天人人| 国产野战对白在线观看| aaaaa片日本免费| 中文亚洲av片在线观看爽 | 久久精品亚洲熟妇少妇任你| 国产精品.久久久| 久久久久久久久免费视频了| 日韩欧美三级三区| 欧美大码av| av欧美777| 欧美黄色淫秽网站| 国产在线一区二区三区精| 国产熟女午夜一区二区三区| 热re99久久精品国产66热6| av网站免费在线观看视频| 一本综合久久免费| 精品国产乱码久久久久久男人| 黄色a级毛片大全视频| 少妇被粗大的猛进出69影院| 亚洲国产中文字幕在线视频| 啦啦啦免费观看视频1| 国产成人免费观看mmmm| 国产精品影院久久| 日韩一卡2卡3卡4卡2021年| 老熟妇仑乱视频hdxx| 国产精品偷伦视频观看了| 黄色成人免费大全| 亚洲一码二码三码区别大吗| 成年人午夜在线观看视频| 老司机影院毛片| 黄色毛片三级朝国网站| 捣出白浆h1v1| 亚洲色图综合在线观看| 欧美激情极品国产一区二区三区| 中文字幕av电影在线播放| 99re6热这里在线精品视频| 正在播放国产对白刺激| 热99久久久久精品小说推荐| 久久ye,这里只有精品| 国产高清激情床上av| 99久久综合精品五月天人人| 精品一区二区三卡| 国产精品久久电影中文字幕 | 日本撒尿小便嘘嘘汇集6| 亚洲全国av大片| 一级作爱视频免费观看| 女人被躁到高潮嗷嗷叫费观| 99热国产这里只有精品6| 午夜激情av网站| 国产精品一区二区在线观看99| 久久久久视频综合| 一本综合久久免费| 国产男靠女视频免费网站| 大陆偷拍与自拍| 黄色片一级片一级黄色片| 久久精品亚洲精品国产色婷小说| 精品一区二区三卡| 亚洲黑人精品在线| 亚洲avbb在线观看| 亚洲中文字幕日韩| 国产精品九九99| 波多野结衣一区麻豆| 亚洲国产欧美网| 天天躁狠狠躁夜夜躁狠狠躁| 老鸭窝网址在线观看| 精品第一国产精品| 国产欧美日韩一区二区三区在线| 久久中文字幕人妻熟女| av网站在线播放免费| 欧美 亚洲 国产 日韩一| 波多野结衣av一区二区av| 免费观看a级毛片全部| 欧美日韩瑟瑟在线播放| 亚洲精品美女久久久久99蜜臀| 99久久综合精品五月天人人| 国产精品久久久久久精品古装| 亚洲七黄色美女视频| 黄色毛片三级朝国网站| 国产欧美日韩一区二区三| 精品国产一区二区久久| 成人精品一区二区免费| 久久国产精品男人的天堂亚洲| 一区二区三区精品91| 国产成人免费无遮挡视频| 少妇被粗大的猛进出69影院| 国产97色在线日韩免费| ponron亚洲| 国产黄色免费在线视频| 搡老乐熟女国产| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区黑人| 国产激情欧美一区二区| 久久精品aⅴ一区二区三区四区| x7x7x7水蜜桃| 日本黄色视频三级网站网址 | 国产亚洲精品一区二区www | 热re99久久国产66热| 大香蕉久久网| 色播在线永久视频| 黄频高清免费视频| 国产精品一区二区精品视频观看| 午夜福利在线观看吧| 一级a爱视频在线免费观看| 在线观看午夜福利视频| 身体一侧抽搐| 久久久精品国产亚洲av高清涩受| tube8黄色片| 老司机福利观看| av欧美777| 国内久久婷婷六月综合欲色啪| 国产欧美日韩一区二区三| 精品第一国产精品| 久久久久精品人妻al黑| 亚洲国产欧美网| 日韩欧美免费精品| 国产精品国产高清国产av | 999久久久国产精品视频| 国产成人一区二区三区免费视频网站| 色在线成人网| 国产成人啪精品午夜网站| 午夜福利乱码中文字幕| 两个人看的免费小视频| 女性生殖器流出的白浆| 精品久久久精品久久久| 国产精品久久久久久人妻精品电影| 成年版毛片免费区| 成人手机av| 一进一出抽搐gif免费好疼 | 又黄又粗又硬又大视频| 老鸭窝网址在线观看| 国产亚洲精品久久久久久毛片 | 十八禁高潮呻吟视频| 首页视频小说图片口味搜索| 欧美黑人精品巨大| 精品久久久久久久久久免费视频 | 一夜夜www| 久久久久国产精品人妻aⅴ院 | 欧洲精品卡2卡3卡4卡5卡区| 欧美午夜高清在线| 国精品久久久久久国模美| 午夜福利在线免费观看网站| 色婷婷久久久亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 久久精品成人免费网站| 欧美性长视频在线观看| 成在线人永久免费视频| 免费在线观看亚洲国产| 亚洲五月色婷婷综合| 性色av乱码一区二区三区2| 色尼玛亚洲综合影院| 成人国产一区最新在线观看| 视频在线观看一区二区三区| 怎么达到女性高潮| 亚洲精品在线观看二区| 中亚洲国语对白在线视频| 看黄色毛片网站| 午夜免费观看网址| 美女扒开内裤让男人捅视频| 亚洲精品在线观看二区| 亚洲五月天丁香| 国产精品乱码一区二三区的特点 | 老熟妇乱子伦视频在线观看| 欧美成狂野欧美在线观看| 午夜精品国产一区二区电影| 久久草成人影院| 丰满的人妻完整版| 亚洲在线自拍视频| 亚洲精品美女久久av网站| 亚洲国产欧美网| 满18在线观看网站| 午夜两性在线视频| 女警被强在线播放| 99re在线观看精品视频| 91老司机精品| 亚洲av电影在线进入| 高清欧美精品videossex| a级毛片在线看网站| 身体一侧抽搐| 欧美丝袜亚洲另类 | 国产亚洲精品第一综合不卡| 国产精品 欧美亚洲| 在线永久观看黄色视频| av国产精品久久久久影院| 美女高潮喷水抽搐中文字幕| 久久国产亚洲av麻豆专区| 亚洲九九香蕉| 又紧又爽又黄一区二区| 黄色片一级片一级黄色片| 精品亚洲成a人片在线观看| 法律面前人人平等表现在哪些方面| 99精品欧美一区二区三区四区| 99精品在免费线老司机午夜| 日韩欧美免费精品| 久久香蕉精品热| 无人区码免费观看不卡| avwww免费| 怎么达到女性高潮| 国产伦人伦偷精品视频| 亚洲中文字幕日韩| 精品国产一区二区久久| 精品国产超薄肉色丝袜足j| 麻豆国产av国片精品| 婷婷精品国产亚洲av在线 | 久久天堂一区二区三区四区| 黄色毛片三级朝国网站| 天天操日日干夜夜撸| 成人特级黄色片久久久久久久| 日本wwww免费看| 99re在线观看精品视频| 免费观看精品视频网站| 欧美日韩亚洲综合一区二区三区_| 一二三四在线观看免费中文在| 久9热在线精品视频| 免费看a级黄色片| 人妻久久中文字幕网| 亚洲av美国av| 亚洲欧美一区二区三区久久| 免费看十八禁软件| 99国产极品粉嫩在线观看| svipshipincom国产片| 欧美日本中文国产一区发布| 在线十欧美十亚洲十日本专区| 青草久久国产| 天天躁狠狠躁夜夜躁狠狠躁| 丰满迷人的少妇在线观看| 亚洲人成77777在线视频| 国产精品1区2区在线观看. | 免费在线观看日本一区| 久久久久久久午夜电影 | 很黄的视频免费| 国产亚洲精品久久久久久毛片 | 9191精品国产免费久久| 欧美黄色片欧美黄色片| www.精华液| 午夜影院日韩av| 亚洲人成电影免费在线| 精品国产亚洲在线| 欧美精品一区二区免费开放| 女同久久另类99精品国产91| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕av电影在线播放| 宅男免费午夜| 久久久国产一区二区| 久久国产乱子伦精品免费另类| 国产成人免费观看mmmm| 日韩制服丝袜自拍偷拍| 免费观看a级毛片全部| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯 | 男女床上黄色一级片免费看| 国产精品一区二区免费欧美| 亚洲精品美女久久av网站| av国产精品久久久久影院| 成年人免费黄色播放视频| 日日摸夜夜添夜夜添小说| 女性生殖器流出的白浆| 丝袜美足系列| 黄频高清免费视频| 亚洲欧美精品综合一区二区三区| 国产成人影院久久av| 成人特级黄色片久久久久久久| 欧美日韩成人在线一区二区| 777久久人妻少妇嫩草av网站| 亚洲七黄色美女视频| 欧美日韩精品网址| 免费人成视频x8x8入口观看| 老司机深夜福利视频在线观看| 大型av网站在线播放| 成人永久免费在线观看视频| 国产亚洲欧美在线一区二区| 王馨瑶露胸无遮挡在线观看| а√天堂www在线а√下载 | 黄频高清免费视频| 久热爱精品视频在线9| 在线观看日韩欧美| 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| 午夜福利视频在线观看免费| 久久这里只有精品19| 涩涩av久久男人的天堂| 亚洲欧美激情综合另类| 欧美在线黄色| 久久久久久久久免费视频了| 岛国毛片在线播放| 91在线观看av| 精品福利永久在线观看| 久久婷婷成人综合色麻豆| 国产精品亚洲av一区麻豆| 欧美老熟妇乱子伦牲交| 国产成人精品久久二区二区免费| 午夜亚洲福利在线播放| 麻豆乱淫一区二区| 日本wwww免费看| 久久亚洲精品不卡| 人妻久久中文字幕网| 欧美午夜高清在线| 国产高清videossex| 老汉色av国产亚洲站长工具| 久久 成人 亚洲| 亚洲九九香蕉| 午夜激情av网站| 欧美成狂野欧美在线观看| 亚洲av成人不卡在线观看播放网| 中文字幕色久视频| 丰满的人妻完整版| 国精品久久久久久国模美| 亚洲九九香蕉| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频| 波多野结衣av一区二区av| 免费观看精品视频网站| 麻豆av在线久日| 在线视频色国产色| 久久国产亚洲av麻豆专区| 18在线观看网站| 啦啦啦免费观看视频1| 国产成人精品久久二区二区免费| 两性夫妻黄色片| 日韩 欧美 亚洲 中文字幕| 亚洲av片天天在线观看| tocl精华| 丰满的人妻完整版| 后天国语完整版免费观看| 久久精品国产亚洲av香蕉五月 | 亚洲国产欧美一区二区综合| 精品亚洲成a人片在线观看| 日日夜夜操网爽| 欧美日韩成人在线一区二区| 99re在线观看精品视频|