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

    基于VAR模型的水資源利用與經(jīng)濟(jì)增長動(dòng)態(tài)關(guān)系研究

    2012-11-29 05:35:12鄧朝暉薛惠鋒
    中國人口·資源與環(huán)境 2012年6期
    關(guān)鍵詞:脈沖響應(yīng)協(xié)整用水量

    鄧朝暉 劉 洋 薛惠鋒

    (1.西安理工大學(xué)經(jīng)濟(jì)與管理學(xué)院,陜西西安710054;2.西北工業(yè)大學(xué)資源與環(huán)境信息化工程研究所,陜西西安710072)

    基于VAR模型的水資源利用與經(jīng)濟(jì)增長動(dòng)態(tài)關(guān)系研究

    鄧朝暉1,2劉 洋1薛惠鋒2

    (1.西安理工大學(xué)經(jīng)濟(jì)與管理學(xué)院,陜西西安710054;2.西北工業(yè)大學(xué)資源與環(huán)境信息化工程研究所,陜西西安710072)

    基于VAR模型,通過變量平穩(wěn)性檢驗(yàn)和協(xié)整分析,廣義脈沖響應(yīng)和預(yù)測方差分解分析,利用中國1980-2007年水資源利用和經(jīng)濟(jì)增長相關(guān)數(shù)據(jù),對中國經(jīng)濟(jì)增長與水資源利用的長期均衡關(guān)系及其動(dòng)態(tài)性進(jìn)行了實(shí)證分析。研究結(jié)果表明:①經(jīng)濟(jì)增長與總用水量、工業(yè)用水量和生活用水量之間存在長期均衡關(guān)系,而農(nóng)業(yè)用水與經(jīng)濟(jì)增長之間不存在長期均衡關(guān)系,這與研究期間中國處于工業(yè)化中期階段的事實(shí)相吻合;②經(jīng)濟(jì)增長對水資源利用的沖擊響應(yīng)的滯后期短且是非漸進(jìn)的,而水資源對經(jīng)濟(jì)增長產(chǎn)生顯著影響的滯后期較長且是非漸進(jìn)的,我國經(jīng)濟(jì)發(fā)展中工業(yè)用水、生活用水量增加趨勢明顯,農(nóng)業(yè)用水量隨著經(jīng)濟(jì)發(fā)展出現(xiàn)零增長和負(fù)增長;③經(jīng)濟(jì)增長對水資源利用的預(yù)測方差起著重要作用,而水資源利用對經(jīng)濟(jì)增長的預(yù)測方差的貢獻(xiàn)度較小。建議加強(qiáng)工業(yè)用水和生活用水的節(jié)水措施,減少工業(yè)和生活用水量,實(shí)現(xiàn)水資源可持續(xù)利用。

    水資源利用;經(jīng)濟(jì)增長;VAR模型;脈沖響應(yīng)

    水資源是基礎(chǔ)性自然資源和戰(zhàn)略性經(jīng)濟(jì)資源,是生態(tài)環(huán)境的控制性要素,是一個(gè)國家綜合國力的有機(jī)組成部分。進(jìn)入21世紀(jì)以來,隨著水資源利用量的急劇上升,越來越多的國家出現(xiàn)了水危機(jī),不僅嚴(yán)重制約著經(jīng)濟(jì)發(fā)展,而且威脅著人類生存安全。

    國外關(guān)于可持續(xù)水資源管理的研究始于上世紀(jì)90年代,以聯(lián)合國可持續(xù)發(fā)展委員會(huì)(UNCSD)開展的研究工作和國外學(xué)者在眾多國際會(huì)議上發(fā)表的論文為標(biāo)志開始興起,已由定性研究為主逐步過渡到以定量研究為主;國內(nèi)有關(guān)可持續(xù)水資源管理的研究也可追溯到上世紀(jì)90年代,許多學(xué)者及水利工作者都對該問題進(jìn)行了早期探索,如馮尚友等學(xué)者介紹了水資源系統(tǒng)的概念,并開展了一系列研究工作;胡治郡等學(xué)者則進(jìn)一步把可持續(xù)水資源管理的概念、研究方法、進(jìn)展等引入我國水文水資源學(xué)術(shù)界,并提出了“水資源可持續(xù)管理”的熱點(diǎn)問題和研究方向[1]。

    國內(nèi)外對水資源管理的研究主要涉及水資源供需分析、水資源承載力、水資源配置、水資源可持續(xù)利用等層面和角度展開。其中,在水資源供需分析中,主要就水資源供需管理的意義、內(nèi)涵、實(shí)施策略、轉(zhuǎn)變方式等問題進(jìn)行了研究[2];水資源承載能力研究則主要從概念、內(nèi)涵、意義[3-5]及評(píng)價(jià)模型[6]等方面入手進(jìn)行了研究;此外有關(guān)水資源配置的研究也較多,但主要是從工程技術(shù)角度進(jìn)行研究,在社會(huì)科學(xué)、經(jīng)濟(jì)學(xué)、政治學(xué)、管理學(xué)等領(lǐng)域的研究還較為滯后。關(guān)于水資源利用與經(jīng)濟(jì)增長之間關(guān)系的研究還很少,美國學(xué)者Charles定性分析了水資源利用對經(jīng)濟(jì)增長的促進(jìn)作用和靜態(tài)關(guān)系[7]。路寧等利用中國52個(gè)城市的截面數(shù)據(jù),通過建立計(jì)量模型對經(jīng)濟(jì)增長與水資源利用壓力之間的關(guān)系進(jìn)行了實(shí)證檢驗(yàn)[8],李周利用GDP和用水量變化的時(shí)間序列證實(shí)了中國經(jīng)濟(jì)增長和水資源利用之間庫茲涅茨曲線的存在性[9],但二者都沒有給出水資源利用與經(jīng)濟(jì)增長之間的動(dòng)態(tài)作用關(guān)系。

    基于此,本文擬采用已廣泛應(yīng)用于地區(qū)經(jīng)濟(jì)[10]、環(huán)境污染與經(jīng)濟(jì)增長[11-14]、產(chǎn)業(yè)結(jié)構(gòu)[15]等關(guān)系研究中并證明取得良好效果的VAR模型計(jì)量經(jīng)濟(jì)學(xué)分析方法,對水資源利用與經(jīng)濟(jì)增長之間的關(guān)系進(jìn)行研究,利用1980-2007年中國主要用水指標(biāo)和GDP數(shù)據(jù),對中國經(jīng)濟(jì)增長與水資源利用之間的關(guān)系進(jìn)行研究,通過對二者之間的協(xié)整檢驗(yàn)、廣義脈沖響應(yīng)和預(yù)測方差分解分析,揭示水資源利用與經(jīng)濟(jì)增長之間的長期動(dòng)態(tài)變化關(guān)系,為認(rèn)識(shí)和解決經(jīng)濟(jì)社會(huì)發(fā)展中水資源短缺問題提供科學(xué)依據(jù)。

    1 中國水資源利用發(fā)展變化

    中國是一個(gè)干旱缺水嚴(yán)重的國家。淡水資源總量為28 000 m3,占全球水資源的6%,僅次于巴西、俄羅斯和加拿大,居世界第四位,但人均只有2 300 m3,僅為世界平均水平的1/4、美國的1/5,在世界上名列121位,是全球13個(gè)人均水資源最貧乏的國家之一。隨著經(jīng)濟(jì)的快速增長,用水量不斷增加,水資源短缺已成為經(jīng)濟(jì)、社會(huì)、環(huán)境協(xié)調(diào)發(fā)展的主要制約因素。

    圖1 1949-2007年中國用水量變化示意圖Fig.1 Change of water consumption in China(1949-2007)數(shù)據(jù)來源:中國可持續(xù)發(fā)展水資源戰(zhàn)略研究報(bào)告;《中國水資源公報(bào)》(1997-2007)

    多年來,全國用水量增長迅速,如圖1所示。1949年估計(jì)為1 031億m3,到1959年翻了一番;1980年全國用水總量達(dá)4 437億m3,與1965年相比,年均增長約為3.3%;1993年與1980年相比,全國用水量年均增長59億m3,增長率為1.2%;1997年全國總用水量為5 566億 m3,與1993年相比,年均增長率為1.7%,年均增幅為92億m3。但從1997年開始,總用水量下降并保持平穩(wěn)趨勢,2004年開始出現(xiàn)上升趨勢,但是上升趨勢較為緩慢,表明隨著經(jīng)濟(jì)增長總用水量增長率逐漸減小,這與世界發(fā)達(dá)國家用水量已達(dá)到零增長甚至負(fù)增長的發(fā)展趨勢是一致的,我國通過提高水資源利用、控制人口增長等措施,實(shí)現(xiàn)我國用水量的零增長是完全可以實(shí)現(xiàn)的。

    從用水結(jié)構(gòu)來看,我國農(nóng)業(yè)用水占全部用水的比重較高,但始終呈遞減趨勢:由1949年的97.1%下降到1980年的88.2%、1993年的78.1%、1997年的75.4%和2007年的61.9%。工業(yè)和生活用水快速上升,工業(yè)用水占總用水量的比重在1949年僅為2.3%,而1997年上升為20%,到2007年為24.1%。生活用水增長更加迅速,1980-1997年間的年均增長率為7.9%。隨著工業(yè)化進(jìn)程的加快和城鎮(zhèn)化水平的提高,如果不加以管理和控制,這種趨勢仍將持續(xù)下去。

    2 數(shù)據(jù)制備和研究方法

    2.1 變量選取與數(shù)據(jù)制備

    在本研究中,對經(jīng)濟(jì)增長和水資源利用兩方面選取具有代表性指標(biāo)變量進(jìn)行計(jì)量分析。GDP是以一個(gè)國家或地區(qū)所有常住經(jīng)濟(jì)單位的生產(chǎn)成果為對象進(jìn)行核算,覆蓋國民經(jīng)濟(jì)所有行業(yè),并具有國際上通用的核算原則與方法,是衡量國家之間、地區(qū)之間經(jīng)濟(jì)活動(dòng)總量的國際通用指標(biāo)。因此,本研究選用GDP作為度量經(jīng)濟(jì)增長的指標(biāo),單位為億元。按照水資源使用結(jié)構(gòu),將用水分為農(nóng)業(yè)用水、工業(yè)用水、生活用水和生態(tài)用水,由于生態(tài)用水量較小,加之生態(tài)用水是近幾年開始作為單獨(dú)的指標(biāo)進(jìn)行統(tǒng)計(jì),為了研究方便,本文將其歸入生活用水中計(jì)算。因此,本文選取總用水量、農(nóng)業(yè)用水量、工業(yè)用水量和生活用水量作為水資源利用指標(biāo)(單位:億m3)。

    水資源專門統(tǒng)計(jì)工作開始相對較晚,中國自1997年開始正式編制《中國水資源公報(bào)》,這給水資源統(tǒng)計(jì)研究本身帶來一定困難。鑒于數(shù)據(jù)的可靠性和可得性,本研究樣本區(qū)間確定為1980-2007年。GDP數(shù)據(jù)源自《中國統(tǒng)計(jì)年鑒》(1980-2007年);1997-2007年水資源數(shù)據(jù)源自《中國水資源公報(bào)》,1980-1997年水資源數(shù)據(jù)在查閱《中國統(tǒng)計(jì)年鑒》、《城市建設(shè)統(tǒng)計(jì)年報(bào)》、《中國城市年鑒》、《中國水利年鑒》、《中國農(nóng)業(yè)年鑒》、《中國工業(yè)年鑒》、《中國環(huán)境年鑒》和相關(guān)研究成果基礎(chǔ)上,進(jìn)行計(jì)算和統(tǒng)計(jì),對部分缺失的數(shù)據(jù)采用灰度預(yù)測和專家評(píng)估等方法進(jìn)行數(shù)據(jù)填充。

    為避免數(shù)據(jù)的劇烈波動(dòng),消除可能存在的異方差,考慮到對時(shí)間序列數(shù)據(jù)進(jìn)行對數(shù)化后容易得到平穩(wěn)序列,并且還不會(huì)改變時(shí)序數(shù)據(jù)的特征,本文對GDP、總用水量、農(nóng)業(yè)用水量、工業(yè)用水量、生活用水量數(shù)據(jù)序列進(jìn)行對數(shù)化處理,分別命名為 LNGDP、LNTAL、LNAGR、LNIND、LNLIV。

    2.2 研究方法

    傳統(tǒng)的經(jīng)濟(jì)計(jì)量方法是以經(jīng)濟(jì)理論為基礎(chǔ)來描述變量關(guān)系的模型。但是,經(jīng)濟(jì)理論通常不足以對變量之間的動(dòng)態(tài)聯(lián)系提供一個(gè)嚴(yán)密的說明,而且內(nèi)生變量既可以出現(xiàn)在方程左端又可出現(xiàn)在方程右端,使得估計(jì)和推斷變得更加復(fù)雜。VAR就是一種用非結(jié)構(gòu)性方法來建立各個(gè)變量之間關(guān)系的模型,其表示如式(1)所示。

    其中:Xt為時(shí)間序列構(gòu)成的向量;c為常數(shù)項(xiàng);p為自回歸滯后階數(shù);Aj為時(shí)間序列系數(shù)矩陣;εt為白噪聲序列向量,滿足:①E(εt)=0,誤差項(xiàng)的均值為 0;②E(εtε't)= Ω,誤差項(xiàng)的協(xié)方差矩陣為 Ω;③E(εtε't-k)=0,誤差項(xiàng)不存在自相關(guān)。

    本文建立由4個(gè)水資源指標(biāo)(總用水量、農(nóng)業(yè)用水量、工業(yè)用水量、生活用水量)與經(jīng)濟(jì)增長指標(biāo)(GDP)組成的雙變量VAR模型,并通過協(xié)整分析和脈沖響應(yīng)分析,對中國經(jīng)濟(jì)增長與水資源利用之間的關(guān)系進(jìn)行實(shí)證研究。在協(xié)整分析及脈沖響應(yīng)函數(shù)分析之前,需要對變量的時(shí)序數(shù)據(jù)進(jìn)行ADF(Augmented Dickey Fuller)平穩(wěn)性檢驗(yàn)。

    協(xié)整理論是2003年諾貝爾經(jīng)濟(jì)學(xué)獎(jiǎng)得主恩格爾(R.F.Engle)和格蘭杰(C.W.J.Granger)在 1987 年首先提出的[17]。所謂協(xié)整是指兩個(gè)或多個(gè)非平穩(wěn)的變量序列某個(gè)線性組合后的序列呈平穩(wěn)性。經(jīng)濟(jì)意義在于兩個(gè)變量雖然具有各自的長期波動(dòng)規(guī)律,但如果是協(xié)整的,那么它們之間存在著一個(gè)長期穩(wěn)定的比例關(guān)系,反之,如果兩個(gè)變量具有各自的長期波動(dòng)規(guī)律,但如果不是協(xié)整的,它們之間就不存在一個(gè)長期穩(wěn)定的關(guān)系。

    為研究經(jīng)濟(jì)增長與水資源利用的長期動(dòng)態(tài)關(guān)系,本文擬采用脈沖響應(yīng)函數(shù)法來刻畫水資源利用與經(jīng)濟(jì)增長兩類變量之間的長期的相互動(dòng)態(tài)作用。脈沖響應(yīng)函數(shù)(IRF)是描述一個(gè)內(nèi)生變量對誤差的反應(yīng),也即在擾動(dòng)項(xiàng)上加一個(gè)標(biāo)準(zhǔn)差大小的新息(Innovation)沖擊對內(nèi)生變量的當(dāng)前值和未來值的影響,其定義為:

    式中,δk代表來自第k個(gè)變量的沖擊,n是沖擊響應(yīng)時(shí)期數(shù),t-1代表沖擊發(fā)生時(shí)所有可獲得的信息。要求n期沖擊的IRF值,即考慮δk沖擊對 xt+n期望值所導(dǎo)致的差異。

    與脈沖響應(yīng)函數(shù)方法不同,VAR預(yù)測方差分解法能給出隨機(jī)信息的相對重要性。其主要思想是,把系統(tǒng)中每個(gè)內(nèi)生變量的預(yù)測均方誤差(Mean Square Error,MSE)按其成因分解為與各方程相關(guān)聯(lián)的m個(gè)部分,從而了解各新息對模型的內(nèi)生變量的相對重要性。VAR(p)模型的s步預(yù)測誤差為:

    它的均方誤差(MSE)為:

    式中,pp'=Ω,根據(jù)式(6)可以將任意1個(gè)內(nèi)生變量的預(yù)測均方誤差分解成系統(tǒng)中的各變量的沖擊貢獻(xiàn)值,然后計(jì)算每個(gè)變量沖擊的相對重要性,即變量的貢獻(xiàn)占總貢獻(xiàn)的比例。本文擬運(yùn)用VAR預(yù)測方差分解法來考察經(jīng)濟(jì)增長與水資源利用之間的相互影響程度[4]。

    3 中國水資源利用與經(jīng)濟(jì)增長關(guān)系實(shí)證

    3.1 VAR模型的建立

    本文分析的VAR模型為中國經(jīng)濟(jì)增長與水資源利用指標(biāo)之間的雙變量系統(tǒng),需要建立GDP與總用水量、GDP與農(nóng)業(yè)用水量、GDP與工業(yè)用水量、GDP與生活用水量4個(gè)雙變量 VAR模型。根據(jù)上述的數(shù)據(jù)分析,利用EVIEWS5.1對動(dòng)態(tài)方程的參數(shù)進(jìn)行估計(jì)見表1。模型中各變量均經(jīng)過取對數(shù)處理,從方程的擬合度和系數(shù)的顯著性,以及滯后階數(shù)判斷的AIC準(zhǔn)則綜合考慮,取各變量的最大滯后階數(shù)為2。

    表1 GDP、總用水量、農(nóng)業(yè)用水量、工業(yè)用水量、生活用水量向量自回歸方程參數(shù)估計(jì)Tab.1 Auto regression parameter estimates of GDP,total water consumption,agricultural water consumption,industrial waterz consumption and living water consumption vector

    對于VAR模型而言,如果VAR模型所有根模的倒數(shù)小于1,即位于單位圓內(nèi),則VAR模型是穩(wěn)定的。如果模型不穩(wěn)定,則某些結(jié)果不是有效的,例如脈沖響應(yīng)函數(shù)的標(biāo)準(zhǔn)差。通過檢驗(yàn),表2和圖2可以判定VAR模型是穩(wěn)定的,可以進(jìn)行脈沖響應(yīng)分析。

    3.2 ADF檢驗(yàn)和協(xié)整檢驗(yàn)

    在進(jìn)行協(xié)整分析之前,首先需要檢驗(yàn)被分析序列變量是否平穩(wěn),即單位根檢驗(yàn)。常用的單位根檢驗(yàn)方法DF檢驗(yàn)由于不能保證方程中的殘差項(xiàng)是白噪聲(white noise),所以Dickey和Fuller對DF檢驗(yàn)法進(jìn)行了擴(kuò)充,形成ADF(Augented Dickey-Fuller Test)檢驗(yàn),這是目前普遍應(yīng)用的單位根檢驗(yàn)方法[16]。本文選用ADF法對水資源利用與經(jīng)濟(jì)增長各變量時(shí)間序列的平穩(wěn)性檢驗(yàn)如表3所示,由于樣本容量的限制,最大滯后階數(shù)取3。

    表2 VAR模型滯后結(jié)構(gòu)檢驗(yàn)Tab.2 Lag structure test of VAR model

    圖2 VAR特征多項(xiàng)式根模倒數(shù)Fig.2 VAR model characteristic polynomial roots reciprocal

    由檢驗(yàn)結(jié)果表明,樣本區(qū)間在5%的顯著水平下,接受所有變量序列水平值有單位根的假設(shè),拒絕所有變量一階差分存在單位根的假設(shè)。檢驗(yàn)結(jié)果說明,1980-2007年的LNGDP、LNTAL、LNAGR、LNIND 和 LNLIV 序列一階差分都是平穩(wěn)的,表明GDP與用水總量、工業(yè)用水量、農(nóng)業(yè)用水量、生活用水量之間可能存在協(xié)整關(guān)系,可以進(jìn)一步檢驗(yàn)其協(xié)整性。

    檢驗(yàn)變量間協(xié)整關(guān)系的方法有EG兩步法和Johansen極大似然法兩種[17]。EG兩步法是 Engle和 Granger于1987年提出的,用來檢驗(yàn)兩個(gè)變量之間協(xié)整關(guān)系的一種方法[18]。本文采用這種簡便的方法來檢驗(yàn)GDP與工業(yè)用水量、GDP與生活用水量之間的協(xié)整關(guān)系,具體檢驗(yàn)步驟為:

    表3 變量序列的單位根檢驗(yàn)(ADF)結(jié)果Tab.3 Unit Root Test(ADF)results of variable sequence

    (1)用 OLS分別對 LNGDP和 LNGYS、LNGDP和LNSHS進(jìn)行靜態(tài)回歸,回歸方程如下:

    (2)分別檢驗(yàn)四個(gè)殘差序列的單整階數(shù),方法與檢驗(yàn)GDP序列平穩(wěn)性的方法相同,結(jié)果如表4所示。

    由檢驗(yàn)結(jié)果可知,回歸方程的殘差序列的ADF檢驗(yàn)值都小于顯著水平的5%時(shí)的臨界值,即殘差序列μ^1t、μ^2t、μ^3t和 μ^4t是平穩(wěn)序列。說明 LNGDP 與 LNTAL、LNGDP 與LNAGR、LNGDP與LNIND、LNGDP與LNLIV之間存在協(xié)整關(guān)系,因此,水資源與經(jīng)濟(jì)增長之間具有長期的均衡關(guān)系。

    3.3 廣義脈沖響應(yīng)分析

    由于VAR模型各個(gè)估計(jì)方程擾動(dòng)項(xiàng)的方差協(xié)方差矩陣不是對焦矩陣,因此必須首先進(jìn)行正交處理得到對角化矩陣,正交化處理常用的是喬利斯基分解(Cholesky)。喬利斯基分解為VAR模型的變量增加一個(gè)次序,并將所有影響變量的公共因素歸結(jié)到VAR模型中第一次出現(xiàn)的變量上,并且如果改變變量的次序,將會(huì)明顯改變變量的響應(yīng)結(jié)果[19]。由于喬利斯基分解依賴次序的缺陷,1998年P(guān)esaran和Shin提出了廣義脈沖響應(yīng)分析[20],這種分析方法不依賴VAR模型中變量次序的正交的殘差矩陣,可提高估計(jì)結(jié)果的穩(wěn)定性與可靠性。為了分析經(jīng)濟(jì)增長與水資源利用之間的動(dòng)態(tài)影響關(guān)系,本文運(yùn)用廣義脈沖響應(yīng)函數(shù)分析二者之間的沖擊響應(yīng)[21],這里將沖擊響應(yīng)期設(shè)定為10期,分析結(jié)果見表5。

    3.3.1 總用水量與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)系

    總用水量與經(jīng)濟(jì)增長的脈沖響應(yīng)分析結(jié)果如表5和圖3。就總用水量對 GDP一個(gè)單位沖擊的響應(yīng)來看,LNTAL當(dāng)期反應(yīng)為負(fù)值(-0.000 49),下一期反應(yīng)上升為正值(0.004 31),然后開始上升,至第三期為最高值(0.006 8),然后開始平穩(wěn)下降,在整個(gè)分析期內(nèi)的LNTAL對LNGDP的累積響應(yīng)值為 0.041,即當(dāng)期 LNGDP對LNTAL的總體影響為正,表明隨著GDP的增長總用水量在增加,但從第三個(gè)周期開始具有下降趨勢。就GDP對總用水量一個(gè)單位沖擊的響應(yīng)來看,LNGDP的當(dāng)前反應(yīng)為負(fù)值(-0.001 66),然后一直下降,整個(gè)分析期內(nèi)的沖擊反應(yīng)均為負(fù)值,累積響應(yīng)值為-0.338,表明總用水量變動(dòng)對經(jīng)濟(jì)增長產(chǎn)生負(fù)面效應(yīng),水資源對經(jīng)濟(jì)發(fā)展具有約束作用。

    3.3.2 農(nóng)業(yè)用水量與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)系

    由表5和圖4可知,LNAGR對LNGDP的單位新息沖擊的響應(yīng)曲線大致為N型,農(nóng)業(yè)用水對GDP一個(gè)單位沖擊響應(yīng),LNAGR當(dāng)期反應(yīng)為負(fù)值(-0.000 19),到第三期上升為最高值(0.004 2),隨后開始下降。單位LNGDP新息沖擊對LNAGR的累積響應(yīng)值為-0.005 6。表明隨著經(jīng)濟(jì)的發(fā)展,農(nóng)業(yè)用水量出現(xiàn)減小趨勢。而LNGDP對LNAGR的沖擊反應(yīng)曲線大致為U曲線,當(dāng)GDP對農(nóng)業(yè)用水量一個(gè)單位沖擊響應(yīng),LNGDP的當(dāng)期反應(yīng)為負(fù)值(-0.000 43),然后開始下降,到第 6期降至最低(-0.032),從第7期開始上升。整個(gè)分析期內(nèi),LNGDP對LNAGR的累積響應(yīng)值為-0.209,表明農(nóng)業(yè)用水量的變動(dòng)對經(jīng)濟(jì)增長產(chǎn)生負(fù)面效應(yīng)。

    圖3 總用水量與經(jīng)濟(jì)增長脈沖響應(yīng)曲線Fig.3 Impulse response curves of the total water consumption and economic growth

    表5 廣義脈沖響應(yīng)分析結(jié)果Tab.5 Results of the generalized impulse response analysis

    圖4 農(nóng)業(yè)用水量與經(jīng)濟(jì)增長脈沖響應(yīng)曲線Fig.4 Impulse response curves of agricultural water consumption and economic growth

    3.3.3 工業(yè)用水量與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)系

    由表5和圖5可知,就工業(yè)用水量對GDP一個(gè)單位沖擊的響應(yīng)來看,LNIND的當(dāng)期反應(yīng)為負(fù)值(-0.001 4),第2期上升為正值(0.011 8),上升到第6期(0.011 8)后開始平穩(wěn)上升,在整個(gè)分析期內(nèi)LNIND對LNGDP的累積響應(yīng)值為0.096,表明GDP的增加導(dǎo)致工業(yè)用水量的增加。而LNGDP對LNIND的一個(gè)單位沖擊的響應(yīng),整個(gè)分析周期內(nèi)均為負(fù)值,并呈現(xiàn)下降趨勢,至第7期降至最低(-0.048)。LNGDP對LNIND累積響應(yīng)值為-0.368,表明工業(yè)用水量變化對經(jīng)濟(jì)增長產(chǎn)生負(fù)面效應(yīng)。

    圖5 工業(yè)用水量與經(jīng)濟(jì)增長脈沖響應(yīng)曲線Fig.5 Impulse response curves of industrial water consumption and economic growth

    圖6 生活用水量與經(jīng)濟(jì)增長脈沖響應(yīng)曲線Fig.6 Impulse response curves of living water consumption and economic growth

    3.3.4 生活用水量與經(jīng)濟(jì)增長的動(dòng)態(tài)關(guān)系

    由表5和圖6可知,就生活用水量對GDP一個(gè)單位沖擊的響應(yīng)來看,LNLIV的當(dāng)期反應(yīng)為負(fù)值(-0.005 6),第2期上升為正值(0.011 8),在第3周期和第4周期發(fā)生突變,至第8期達(dá)到最大值(0.038 2),在整個(gè)分析期內(nèi)LNLIV對LNGDP的累積響應(yīng)值為0.263,表明GDP的增加導(dǎo)致生活用水量的增加。就LNGDP對LNLIV的一個(gè)單位沖擊的響應(yīng)來看,整個(gè)分析周期內(nèi)均為負(fù)值,并呈現(xiàn)下降趨勢,累積響應(yīng)值為-0.121 1,表明生活用水量變化對經(jīng)濟(jì)增長產(chǎn)生負(fù)面效應(yīng)。

    3.4 水資源利用與經(jīng)濟(jì)增長的預(yù)測方差分解

    由水資源利用指標(biāo)與GDP的預(yù)測方差分解結(jié)果(表6)可知,就總體而言,GDP解釋各用水量指標(biāo)的預(yù)測方差分解的貢獻(xiàn)度較高,GDP解釋了總用水量、工業(yè)用水量、生活用水量三變量15%以上的方差,其中對總用水量的預(yù)測方差高達(dá)24.89%,GDP對農(nóng)業(yè)用水量的預(yù)測方差較小(4.6577%)。此分析刻畫了自20世紀(jì)80年代以來,中國水資源利用與經(jīng)濟(jì)增長之間的變化關(guān)系:經(jīng)濟(jì)增長、工業(yè)化、城市化進(jìn)程的加快伴隨著對水資源的過度開發(fā)利用與水資源浪費(fèi),工業(yè)用水和生活用水增加是總用水量增加的主要原因。相比而言,水資源利用對GDP的預(yù)測方差的解釋貢獻(xiàn)度較小,三類主要用水指標(biāo)對GDP的預(yù)測方差的解釋貢獻(xiàn)度均低于2%,尤其是工業(yè)用水量對GDP的方差分解平均貢獻(xiàn)度僅為0.428%,幾乎可以忽略??傆盟繉?jīng)濟(jì)增長的方差分解平均貢獻(xiàn)度也僅僅為7.07%,遠(yuǎn)遠(yuǎn)低于GDP對水資源利用預(yù)測方差的貢獻(xiàn)度。說明引起經(jīng)濟(jì)發(fā)展變化的原因是多方面的,而水資源對經(jīng)濟(jì)增長的影響作用僅僅是一個(gè)方面,這與中國當(dāng)前經(jīng)濟(jì)發(fā)展中水資源利用的現(xiàn)狀是完全相符的。

    表6 水資源利用與GDP的預(yù)測方差分解平均值Tab.6 Forecast variance decomposition average of water resources use and GDP

    4 結(jié)論與建議

    本文基于1980年-2007年中國水資源利用與GDP的時(shí)間序列數(shù)據(jù),建立了中國水資源利用與經(jīng)濟(jì)增長之間的VAR模型,檢驗(yàn)了水資源利用與經(jīng)濟(jì)增長之間的協(xié)整關(guān)系,運(yùn)用廣義脈沖響應(yīng)函數(shù)和預(yù)測方差分析了水資源利用與經(jīng)濟(jì)增長之間的動(dòng)態(tài)響應(yīng)關(guān)系,得到以下結(jié)論和建議:

    (1)研究期間,中國經(jīng)濟(jì)增長與總用水量、工業(yè)用水量、生活用水量之間存在協(xié)整關(guān)系,而農(nóng)業(yè)用水量與經(jīng)濟(jì)增長之間不具有協(xié)整關(guān)系。也就是說,除農(nóng)業(yè)用水外,中國水資源利用與經(jīng)濟(jì)增長之間存在長期穩(wěn)定的均衡關(guān)系。說明隨著經(jīng)濟(jì)增長,中國農(nóng)業(yè)用水量變化基本保持平穩(wěn)狀態(tài),出現(xiàn)零增長甚至負(fù)增長,農(nóng)業(yè)節(jié)水初見成效。但是總用水量、工業(yè)用水量和生活用水量仍保持著較快的增長趨勢,經(jīng)濟(jì)增長對降低水資源使用量的作用不夠明顯,尤其是工業(yè)用水量和生活用水量的增加沒有得到有效控制,這與我國目前處于工業(yè)化中期階段的事實(shí)相符。建議通過實(shí)施中水利用、提高水資源利用率、倡導(dǎo)生活節(jié)水等戰(zhàn)略措施和政策,降低工業(yè)用水量和生活用水量,以實(shí)現(xiàn)中國經(jīng)濟(jì)發(fā)展用水量的零增長乃至負(fù)增長的目標(biāo)。

    (2)經(jīng)濟(jì)增長對水資源利用的沖擊響應(yīng)滯后期短(3年左右)且是非漸進(jìn)的,而水資源對經(jīng)濟(jì)增長產(chǎn)生顯著影響的滯后期較長(5年左右)且是非漸進(jìn)的。經(jīng)濟(jì)增長對總用水量、農(nóng)業(yè)用水量、工業(yè)用水量、生活用水量的單位沖擊響應(yīng)累積值均為負(fù)值 (-0.338,-0.209,-0.368,-0.121 1),而總用水量、工業(yè)用水量、生活用水量對經(jīng)濟(jì)增長的單位沖擊響應(yīng)累積值均為正值(0.041,0.096,0.263),農(nóng)業(yè)用水量對經(jīng)濟(jì)增長的單位沖擊響應(yīng)累積值為負(fù)(-0.005 6)。以上結(jié)果表明,經(jīng)濟(jì)增長帶來總用水量、工業(yè)用水量、生活用水量的增加,農(nóng)業(yè)用水量隨著經(jīng)濟(jì)增長出現(xiàn)零增長和負(fù)增長的趨勢;水資源的減少對經(jīng)濟(jì)發(fā)展具有約束作用。說明我國經(jīng)濟(jì)發(fā)展中工業(yè)用水、生活用水量增加趨勢明顯,水資源對經(jīng)濟(jì)增長的約束作用較為明顯。建議在經(jīng)濟(jì)發(fā)展開發(fā)利用水資源的同時(shí),必須加大對水利工程、節(jié)水新技術(shù)、水資源管理等的投入,充分發(fā)揮經(jīng)濟(jì)對水資源利用的積極促進(jìn)作用,減緩工業(yè)用水和生活用水的增長速度。

    (3)經(jīng)濟(jì)增長對水資源利用的預(yù)測方差起著重要作用,而水資源利用對經(jīng)濟(jì)增長的預(yù)測方差的貢獻(xiàn)度較小。當(dāng)前,一方面要緩解經(jīng)濟(jì)發(fā)展帶來的用水量增加的壓力,另一方面要重視水資源短缺對經(jīng)濟(jì)發(fā)展帶來的潛在反作用。要通過政策和措施減少工業(yè)用水量和生活用水量,建立起有效的水資源保護(hù)體系和虛擬水交易機(jī)制,以形成水資源對經(jīng)濟(jì)發(fā)展長效良好反饋機(jī)制。

    Reference)

    [1]胡治郡.水資源可持續(xù)管理研究進(jìn)展綜述[J].廣東科技,2005,(5):57.[ Hu Zhijun.A Summary of Sustainable Water Resources Management Research[J].Guangdong Science & Technology,2005,(5):57.]

    [2]姜東暉.農(nóng)用水資源需求管理研究:一個(gè)綜述[J].山東農(nóng)業(yè)大學(xué)學(xué)報(bào):社會(huì)科學(xué)版,2011,(1):14 -17.[Jiang Donghui.An Overview of the Research on Management of Agricultural Water Demand[J].Journal of Shandong Agricultural University:Social Science Edition,2011,(1):14 -17.]

    [3]龍騰銳,姜文超,何強(qiáng).水資源承載力內(nèi)涵的新認(rèn)識(shí)[J].水利學(xué)報(bào),2004,35(1):38 - 45.[Long Tengrui,Jiang Wenchao,He Qiang.Water Resources Carrying Capacity:New Perspectives based on Eco-economic Analysis and Sustainable Development[J].Journal of Hydraulic Engineering,2004,35(1):38 -45.]

    [4]姚治君,王建華,江東,等.區(qū)域水資源承載力的研究進(jìn)展及其理論探析[J].水科學(xué)進(jìn)展,2002,13(1):111 -115.[Yao Zhijun,Wang Jianhua,Jiang Dong,et al.Advances in Study on Regional Water Resources Carrying Capacity and Research on Its Theory[J].Advances in Water Science,2002,13(1):111 -115.]

    [5]Varis O,Vakkilainen P.China's Challenges to Water Resources Management in the First Quarterofthe 21stCentury [J].Geomorphology,2001,41(2 -3):93 -104.

    [6]王薇,雷學(xué)東,余新曉,等.基于SD模型的水資源承載力計(jì)算理論研究[J].水資源與水工程學(xué)報(bào),2005,16(3):11 -15.[Wang Wei,Lei Xuedong,Yu Xinxiao,et al.Study on the Region Carrying Capacity of Water Resources based on System Dynamics(SD)Model[J].Journal of Water Resources and Water Engineering,2005,16(3):11 -15.]

    [7]W Howe C.Water Resources and Regional Economic Growth in the United States,1950 -1960[J].Southern Economic Journal,1968,34(4):477-489.

    [8]路寧,周海光.中國城市經(jīng)濟(jì)與水資源利用壓力的關(guān)系研究[J].中國人口·資源與環(huán)境,2010,20(5):48-50.[Lu Ning,Zhou Haiguang.Study on the Relationship Between Water-use-press and Economic Growth of Cities in China[J]. China Population,Resources and Environment,2010,20(5):48 -50.]

    [9]李周.中國農(nóng)村發(fā)展研究報(bào)告No.6[M].北京:社會(huì)科學(xué)文獻(xiàn)出版社,2008:2.[Li Zhou.Report of China Rural Development Research No.6[M].Beijing:Social Sciences Academic Press,2008:2.]

    [10]王舒健,李釗.中國地區(qū)經(jīng)濟(jì)增長互動(dòng)關(guān)系的脈沖響應(yīng)分析[J].數(shù)理統(tǒng)計(jì)與管理,2007,26(3):385 -390.[Wang Shujian,Li Zhao.A Pulse-responding Analysis on Mutual Effect of Regional Economic Growth in China[J].Application of Statistics and Management,2007,26(3):385 - 390.]

    [11]劉坤,劉賢趙,常文靜.煙臺(tái)市經(jīng)濟(jì)增長與環(huán)境污染關(guān)系實(shí)證研究:基于VAR計(jì)量技術(shù)的檢驗(yàn)分析[J].環(huán)境科學(xué)學(xué)報(bào),2007,27(11):1929 -1936.[Liu Kun,Liu Xianzhao,Chang Wenjing.An Empirical Study on the Interactions of Economic Growth and Environmental Pollution in Yantai:Based on the Unrestricted VAR Method[J].Acta Scientiae Circumstantiae,2007,27(11):1929 -1936.]

    [12]李慶華,鄧萍萍,宋琴.中國環(huán)境污染與經(jīng)濟(jì)增長關(guān)系的實(shí)證分析[J].資源開發(fā)與市場,2011,(2):131 -134.[Li Qinhua,Deng Pingping,Song Qin.Empirical Study on Interactions of Chinese EnvironmentalPollution and EconomicGrowth[J].Resource Development& Market,2011,(2):131 -134.]

    [13]李琳,朱金山,高潤霞.基于VAR模型的重慶市經(jīng)濟(jì)增長與環(huán)境污染的關(guān)系研究[J].西南大學(xué)學(xué)報(bào):自然科學(xué)版.2009,31(11):92 - 96.[Li Lin,Zhu Jinshan,Gao Runxia.The Study on Relationship Between Economic Growth and Environmental Pollution in Chongqing Based on VAR Model[J].Journal of Southwest University:Natural Sciences Edition,2009,31(11):92 -96.]

    [14]吳丹,吳仁海.不同地區(qū)經(jīng)濟(jì)增長與環(huán)境污染關(guān)系的VAR模型分析:基于廣州、佛山、肇慶經(jīng)濟(jì)圈的實(shí)證研究[J].環(huán)境科學(xué)學(xué)報(bào),2011,31(4):880 - 888.[Wu Dan,Wu Renhai.VAR Model for the Analysis of Dynamic Interactions between Economic Growth and Environmental Pollution:A Case Study on the Economic Circle of the Cities of Guangzhou,F(xiàn)oshan and Zhaoqing[J].Acta Scientiae Circumstantiae,2011,31(4):880 -888.]

    [15]趙春艷.我國經(jīng)濟(jì)增長與產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)關(guān)系的研究[J].數(shù)理統(tǒng)計(jì)與管理,2008,27(5):487 -492.[Zhao Chunyan.The Research on the Relationship between Economic Growth and the Industrial Structure Evolvement of Our Country:Demonstration Investigation based on the Panel Data Model[J].Application of Statistics and Management,2008,27(5):487 - 492.]

    [16]李子奈,葉阿忠.高等計(jì)量經(jīng)濟(jì)學(xué)[M].北京:清華大學(xué)出版社,2000.[Li Zinai,Ye Azhong.Advanced Econometrics[M].Bei Jing:Tsinghua University Press,2000.]

    [17]EngerR, GrangerC. Co-integration and ErrorCorrection:Representation,Estimation and Testing[J].Econometrics,1987,35:391-407.

    [18]Engle R F,Yoo B S.Forecasting and Testing in Cointegrated Systems[J].Journal of Econometrics,1987,35:143 -159.

    [19]Granger C W J.Some Properties of Time Series Data and Their Use in Econometric Model Specification[J].Journal of Econometrics,1981,16:120 -130.

    [20]Koop G,Pesaran M ,Potter S.Impulse Response Analysis in Nonlinear Multivariate Models[J].Journal of Econometrics,1996,74(1):119-147.

    [21]李國柱.經(jīng)濟(jì)增長與環(huán)境協(xié)調(diào)發(fā)展的計(jì)量分析[M].中國經(jīng)濟(jì)出版社,2007:9.[LiGuozhu. TheEconometricAnalysisof Coordinated Development of Economic Growth and Environment[M].China Economic Publishing House,2007:9.]

    Study on the Dynamic Relationship Between Economic Growth and Water Resources-Use Based on the VAR Model

    DENG Zhao-hui1,2LIU Yang2XUE Hui-feng1
    (1.College of Business,Xi'an University of Technology,Xi'an Shaanxi 710054,China;2.Information Engineering Institute of Resources and Environment,Northwestern Polytechnic University,Xi'an Shaanxi 710072,China)

    Based on the VAR model,the long term dynamic relationship between Chinese economic growth and water resources use is analyzed by using the variable stationary test and co-integration analysis,impulse response and variance decomposition analysis,and the related data of water resources and GDP indices in China from 1980 to 2007.The results show that economic growth is extendedly and dynamically related to total water resources consumption,industrial water consumption and domestic water consumption,but it is not the case with agricultural water consumption.This accords with the fact that China is in the mid-term of industrialization.Moreover,the lag phase of the shock response of economic growth to water resources-use is short and non-asymptotic,whereas the lag phase of the impact of water resources-use on economic growth is long and non-asymptotic.With the development of Chinese economy,industrial and living water consumption is obviously increasing while agricultural water consumption shows the trends of zero growth and negative growth.Additionally,variance decomposition analysis indicates that economic growth is the main variable to the forecasting mean square error of water resources while it is not so vice versa.Therefore,water-saving measures of industrial and living consumption should be strengthened so as to reduce water consumption and achieve sustainable use of water resources.

    VAR model;water resources;economic growth;impulse response

    C812

    A

    1002-2104(2012)06-0128-08

    10.3969/j.issn.1002-2104.2012.06.021

    (編輯:王愛萍)

    2011-11-19

    徐鶴,博士生,主要研究方向?yàn)樗乃Y源。

    國家自然科學(xué)基金重點(diǎn)項(xiàng)目“社會(huì)水循環(huán)系統(tǒng)演化機(jī)制與過程模擬研究——以海河流域?yàn)槔?編號(hào):40830637);國家自然科學(xué)基金創(chuàng)新研究群體基金項(xiàng)目“流域水循環(huán)模擬與調(diào)控”(編號(hào):50721006)。

    猜你喜歡
    脈沖響應(yīng)協(xié)整用水量
    你的用水量是多少?
    基于重復(fù)脈沖響應(yīng)的發(fā)電機(jī)轉(zhuǎn)子繞組匝間短路檢測技術(shù)的研究與應(yīng)用
    澳大利亞研發(fā)出新型農(nóng)業(yè)傳感器可預(yù)測農(nóng)作物用水量
    外商直接投資對我國進(jìn)出口貿(mào)易影響的協(xié)整分析
    河南金融發(fā)展和城鄉(xiāng)居民收入差距的協(xié)整分析
    脈沖響應(yīng)函數(shù)下的我國貨幣需求變動(dòng)與決定
    基于有限元素法的室內(nèi)脈沖響應(yīng)的仿真
    電大理工(2015年3期)2015-12-03 11:34:12
    工業(yè)給排水用水量計(jì)算的重要性
    玻璃氣體放電管與陶瓷氣體放電管的納秒脈沖響應(yīng)特性比較
    中國居民消費(fèi)與經(jīng)濟(jì)增長的協(xié)整關(guān)系檢驗(yàn)
    又黄又爽又刺激的免费视频.| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av男天堂| 一二三四中文在线观看免费高清| 免费看av在线观看网站| 国产伦理片在线播放av一区| 国产精品无大码| 国产精品日韩av在线免费观看| 禁无遮挡网站| 欧美一级a爱片免费观看看| 免费av不卡在线播放| 五月伊人婷婷丁香| 一个人观看的视频www高清免费观看| 久久精品久久久久久久性| 波多野结衣巨乳人妻| 91久久精品电影网| 国产欧美另类精品又又久久亚洲欧美| 一本一本综合久久| 精品久久久噜噜| h日本视频在线播放| eeuss影院久久| 午夜亚洲福利在线播放| 中文字幕亚洲精品专区| 夫妻性生交免费视频一级片| 亚洲国产精品专区欧美| 久久久久九九精品影院| 国产亚洲av片在线观看秒播厂 | 2022亚洲国产成人精品| 国产高清国产精品国产三级 | 一区二区三区乱码不卡18| 欧美成人午夜免费资源| 国产精品一区二区在线观看99 | 麻豆精品久久久久久蜜桃| 欧美一级a爱片免费观看看| 免费观看在线日韩| 日韩av免费高清视频| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 亚洲欧美精品自产自拍| 国产国拍精品亚洲av在线观看| 午夜福利视频1000在线观看| 国产色婷婷99| 日韩av不卡免费在线播放| 少妇的逼好多水| 一本一本综合久久| 中文字幕久久专区| 身体一侧抽搐| 夫妻午夜视频| 国产麻豆成人av免费视频| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 成人漫画全彩无遮挡| 三级国产精品片| 日韩在线高清观看一区二区三区| 久久99蜜桃精品久久| 我要看日韩黄色一级片| 免费看日本二区| 久久精品国产亚洲av涩爱| freevideosex欧美| 日本免费a在线| 久久精品久久久久久久性| 亚洲,欧美,日韩| av播播在线观看一区| 国产高清国产精品国产三级 | 久99久视频精品免费| 秋霞伦理黄片| 国产极品天堂在线| 成人漫画全彩无遮挡| 日本黄色片子视频| 日韩成人伦理影院| 亚洲成人精品中文字幕电影| 久久97久久精品| 中文字幕免费在线视频6| 欧美激情久久久久久爽电影| 亚洲av男天堂| 在现免费观看毛片| 天堂av国产一区二区熟女人妻| 22中文网久久字幕| 偷拍熟女少妇极品色| 亚洲精品中文字幕在线视频 | 成人漫画全彩无遮挡| 欧美丝袜亚洲另类| 久久久久免费精品人妻一区二区| 亚洲精品成人久久久久久| 国产高清有码在线观看视频| 国产在视频线在精品| 九九在线视频观看精品| 日韩大片免费观看网站| av女优亚洲男人天堂| 日韩一区二区视频免费看| 91午夜精品亚洲一区二区三区| 丰满乱子伦码专区| 亚洲精品456在线播放app| 欧美激情久久久久久爽电影| 免费电影在线观看免费观看| 最近2019中文字幕mv第一页| 天天一区二区日本电影三级| 七月丁香在线播放| 亚洲成人久久爱视频| 午夜免费男女啪啪视频观看| 亚洲熟女精品中文字幕| 大片免费播放器 马上看| 乱码一卡2卡4卡精品| 日本-黄色视频高清免费观看| 国产亚洲91精品色在线| 99热6这里只有精品| 尾随美女入室| 热99在线观看视频| 在线免费十八禁| 国产午夜精品一二区理论片| 亚洲国产精品sss在线观看| 国产黄色小视频在线观看| 晚上一个人看的免费电影| 成人亚洲精品一区在线观看 | 亚洲av日韩在线播放| 有码 亚洲区| 97精品久久久久久久久久精品| 国产精品不卡视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 最后的刺客免费高清国语| 中文字幕久久专区| 亚洲国产精品sss在线观看| 国产成人a∨麻豆精品| 老师上课跳d突然被开到最大视频| 永久网站在线| 青春草亚洲视频在线观看| 免费看av在线观看网站| 老司机影院毛片| 亚洲18禁久久av| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲av涩爱| 99久久精品热视频| 亚洲成人久久爱视频| 人妻一区二区av| 男女下面进入的视频免费午夜| 亚洲婷婷狠狠爱综合网| 国产激情偷乱视频一区二区| 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 色5月婷婷丁香| 久久精品夜夜夜夜夜久久蜜豆| 日本一二三区视频观看| kizo精华| 精品久久久久久久久亚洲| 成人无遮挡网站| 免费人成在线观看视频色| 日本与韩国留学比较| 国产成人a∨麻豆精品| 国产成人a区在线观看| 男女边摸边吃奶| 观看美女的网站| 亚洲精品乱码久久久v下载方式| 最近中文字幕2019免费版| 九色成人免费人妻av| 成年免费大片在线观看| 欧美日韩综合久久久久久| 国产av不卡久久| 天堂俺去俺来也www色官网 | 日本av手机在线免费观看| 亚洲不卡免费看| 国产精品一区二区性色av| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| av网站免费在线观看视频 | 欧美变态另类bdsm刘玥| 亚洲精品日韩av片在线观看| 午夜精品在线福利| 日韩国内少妇激情av| 99久久精品热视频| 免费在线观看成人毛片| 大陆偷拍与自拍| 有码 亚洲区| 欧美成人精品欧美一级黄| 插逼视频在线观看| 成人漫画全彩无遮挡| 国产在线一区二区三区精| 永久免费av网站大全| 深爱激情五月婷婷| 丝袜喷水一区| 精品人妻视频免费看| 丰满少妇做爰视频| 波多野结衣巨乳人妻| 男女下面进入的视频免费午夜| 久久久色成人| 欧美激情国产日韩精品一区| 中文精品一卡2卡3卡4更新| 建设人人有责人人尽责人人享有的 | 国产 一区精品| 丰满少妇做爰视频| 国产伦一二天堂av在线观看| 国产成人免费观看mmmm| 极品教师在线视频| 麻豆久久精品国产亚洲av| 国产精品女同一区二区软件| 日本一本二区三区精品| 日本色播在线视频| 简卡轻食公司| 精品国产一区二区三区久久久樱花 | 国产成人freesex在线| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲网站| 精品一区二区三区人妻视频| 国产精品久久久久久av不卡| 亚洲国产日韩欧美精品在线观看| 日韩欧美三级三区| 欧美日韩亚洲高清精品| 久久热精品热| 欧美日韩国产mv在线观看视频 | 赤兔流量卡办理| av国产免费在线观看| 国内精品宾馆在线| 国产白丝娇喘喷水9色精品| 日本猛色少妇xxxxx猛交久久| 亚洲国产成人一精品久久久| 纵有疾风起免费观看全集完整版 | 亚洲欧美成人精品一区二区| 99久久精品热视频| 99热全是精品| 国产一级毛片七仙女欲春2| 欧美日韩视频高清一区二区三区二| 日日啪夜夜爽| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区国产| av天堂中文字幕网| 最近中文字幕高清免费大全6| 欧美一级a爱片免费观看看| 美女主播在线视频| 久久精品国产亚洲网站| 成人毛片60女人毛片免费| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 久久久久久伊人网av| 美女高潮的动态| 欧美xxⅹ黑人| 18禁在线无遮挡免费观看视频| 国产成人午夜福利电影在线观看| 国产一区亚洲一区在线观看| 亚洲美女搞黄在线观看| 国产精品无大码| 久久国内精品自在自线图片| 伦理电影大哥的女人| 日韩伦理黄色片| 一区二区三区免费毛片| 看免费成人av毛片| 国产成人91sexporn| 男人和女人高潮做爰伦理| 精品酒店卫生间| 激情五月婷婷亚洲| 午夜精品一区二区三区免费看| 看非洲黑人一级黄片| 水蜜桃什么品种好| 人体艺术视频欧美日本| 卡戴珊不雅视频在线播放| 日韩欧美一区视频在线观看 | 波多野结衣巨乳人妻| 七月丁香在线播放| 亚洲图色成人| 国产三级在线视频| 少妇人妻精品综合一区二区| 亚洲av日韩在线播放| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| 国产精品一区二区三区四区久久| 尤物成人国产欧美一区二区三区| 黄片wwwwww| 国产中年淑女户外野战色| 插逼视频在线观看| 精品久久久久久久末码| 看免费成人av毛片| 日韩成人av中文字幕在线观看| 日韩一区二区视频免费看| 国产精品一及| 国产乱人偷精品视频| 亚洲欧美日韩无卡精品| 青春草国产在线视频| 国产一级毛片七仙女欲春2| 一二三四中文在线观看免费高清| 一个人看视频在线观看www免费| 国产av不卡久久| 国产乱人偷精品视频| 日本色播在线视频| 干丝袜人妻中文字幕| 国产乱人视频| 成人性生交大片免费视频hd| 三级国产精品欧美在线观看| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 偷拍熟女少妇极品色| 久久久精品免费免费高清| 久久久久久久久大av| 国产精品一区www在线观看| 久久久久久久久久成人| 精华霜和精华液先用哪个| 在线天堂最新版资源| 日韩av免费高清视频| 狂野欧美白嫩少妇大欣赏| 中国美白少妇内射xxxbb| 天天一区二区日本电影三级| 久久久精品免费免费高清| 汤姆久久久久久久影院中文字幕 | 色播亚洲综合网| 国产精品人妻久久久影院| 国产在视频线精品| 国产色婷婷99| av在线观看视频网站免费| 免费大片18禁| 亚洲精品乱码久久久v下载方式| 肉色欧美久久久久久久蜜桃 | 99热网站在线观看| 国产精品美女特级片免费视频播放器| 国产伦理片在线播放av一区| 成人国产麻豆网| 麻豆成人午夜福利视频| 搞女人的毛片| 国产成人91sexporn| 久99久视频精品免费| 日韩国内少妇激情av| 国产精品国产三级国产av玫瑰| 亚洲精品成人av观看孕妇| 色综合色国产| 国产亚洲91精品色在线| 蜜桃亚洲精品一区二区三区| 精品熟女少妇av免费看| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 午夜福利在线在线| videossex国产| 亚洲精华国产精华液的使用体验| 亚洲av.av天堂| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 别揉我奶头 嗯啊视频| 美女主播在线视频| 精品99又大又爽又粗少妇毛片| 美女大奶头视频| 麻豆精品久久久久久蜜桃| 日韩精品青青久久久久久| 国产成人aa在线观看| 1000部很黄的大片| 国产精品美女特级片免费视频播放器| 男女边摸边吃奶| 国产成人aa在线观看| 成人毛片a级毛片在线播放| 久久久色成人| av网站免费在线观看视频 | 欧美日韩综合久久久久久| 毛片一级片免费看久久久久| 久久久久久久久久人人人人人人| videossex国产| av福利片在线观看| av国产久精品久网站免费入址| 看十八女毛片水多多多| 欧美日韩国产mv在线观看视频 | 黄色配什么色好看| 国产免费视频播放在线视频 | 男人爽女人下面视频在线观看| 亚洲国产高清在线一区二区三| 亚洲va在线va天堂va国产| 大陆偷拍与自拍| 97热精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 好男人视频免费观看在线| 国产精品1区2区在线观看.| 国产精品久久久久久av不卡| 国产老妇女一区| 性色avwww在线观看| 国产淫片久久久久久久久| 欧美最新免费一区二区三区| 插阴视频在线观看视频| 亚洲欧美一区二区三区黑人 | 欧美日韩精品成人综合77777| 婷婷色av中文字幕| 久久精品国产自在天天线| 18禁动态无遮挡网站| 久久精品国产自在天天线| 欧美日韩亚洲高清精品| 男女边摸边吃奶| 欧美成人a在线观看| 国产69精品久久久久777片| 久久精品熟女亚洲av麻豆精品 | 国国产精品蜜臀av免费| 国产色爽女视频免费观看| 2021天堂中文幕一二区在线观| 插阴视频在线观看视频| 亚洲欧美日韩东京热| 我要看日韩黄色一级片| 可以在线观看毛片的网站| 美女黄网站色视频| 高清毛片免费看| 男女边摸边吃奶| 久久久久九九精品影院| 十八禁网站网址无遮挡 | 在线观看免费高清a一片| 天堂av国产一区二区熟女人妻| 亚洲国产av新网站| 天堂俺去俺来也www色官网 | 淫秽高清视频在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲精品影视一区二区三区av| 欧美高清成人免费视频www| 一级毛片 在线播放| 蜜桃亚洲精品一区二区三区| 在线观看av片永久免费下载| 一区二区三区免费毛片| 高清视频免费观看一区二区 | 国产又色又爽无遮挡免| 一级毛片久久久久久久久女| 国产精品.久久久| 精品久久久久久久末码| 大话2 男鬼变身卡| 在线免费观看不下载黄p国产| 午夜福利在线在线| 亚洲av不卡在线观看| 九九久久精品国产亚洲av麻豆| videos熟女内射| 亚洲自偷自拍三级| 国产成人午夜福利电影在线观看| 最近中文字幕高清免费大全6| 国产精品.久久久| 丰满少妇做爰视频| 中文天堂在线官网| 欧美3d第一页| 久久99热这里只有精品18| videos熟女内射| 久久人人爽人人爽人人片va| 日韩在线高清观看一区二区三区| 国产av国产精品国产| av天堂中文字幕网| 亚洲国产色片| 黄片无遮挡物在线观看| 精品99又大又爽又粗少妇毛片| 亚洲综合色惰| 日韩一区二区视频免费看| 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 好男人视频免费观看在线| 成年女人看的毛片在线观看| 国产中年淑女户外野战色| 97人妻精品一区二区三区麻豆| 99热这里只有精品一区| 国产精品人妻久久久影院| 成人午夜精彩视频在线观看| 久久99热这里只频精品6学生| 欧美zozozo另类| 精品酒店卫生间| 18禁裸乳无遮挡免费网站照片| 1000部很黄的大片| 亚洲一区高清亚洲精品| 亚洲怡红院男人天堂| 亚洲在线观看片| 97在线视频观看| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看吧| 亚洲欧美一区二区三区黑人 | 色吧在线观看| 久久久久久久久中文| 综合色丁香网| 精品国产露脸久久av麻豆 | 免费观看av网站的网址| 免费黄色在线免费观看| 26uuu在线亚洲综合色| 国产大屁股一区二区在线视频| 日本猛色少妇xxxxx猛交久久| 日日摸夜夜添夜夜添av毛片| 性插视频无遮挡在线免费观看| 午夜福利在线观看免费完整高清在| 最近的中文字幕免费完整| 亚洲成人中文字幕在线播放| 亚洲天堂国产精品一区在线| av一本久久久久| 成人高潮视频无遮挡免费网站| 插阴视频在线观看视频| 亚洲精品日韩在线中文字幕| 亚洲人成网站在线播| 爱豆传媒免费全集在线观看| 亚洲欧美成人精品一区二区| 亚洲精品影视一区二区三区av| 亚洲av日韩在线播放| 国产精品一二三区在线看| 日韩伦理黄色片| 国产成人精品一,二区| a级毛片免费高清观看在线播放| 成人高潮视频无遮挡免费网站| 亚洲欧美清纯卡通| 国产亚洲5aaaaa淫片| 蜜桃久久精品国产亚洲av| 黑人高潮一二区| 免费av观看视频| av卡一久久| 精品一区二区三区人妻视频| 婷婷色av中文字幕| 国产欧美另类精品又又久久亚洲欧美| 亚洲av国产av综合av卡| 如何舔出高潮| www.色视频.com| 久久久精品免费免费高清| 国产av国产精品国产| 亚洲欧美精品自产自拍| 亚洲成人一二三区av| 色吧在线观看| 国产在线男女| 毛片一级片免费看久久久久| 久久热精品热| 高清毛片免费看| 国产精品三级大全| 少妇高潮的动态图| 久99久视频精品免费| 3wmmmm亚洲av在线观看| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 亚洲第一区二区三区不卡| 天堂俺去俺来也www色官网 | 高清欧美精品videossex| 蜜臀久久99精品久久宅男| 高清毛片免费看| 国国产精品蜜臀av免费| 成人亚洲精品一区在线观看 | 青春草亚洲视频在线观看| 老司机影院毛片| 在线免费观看不下载黄p国产| 亚洲av成人精品一二三区| 亚洲人成网站高清观看| 日本免费a在线| 欧美三级亚洲精品| 一级av片app| 日韩三级伦理在线观看| 亚洲,欧美,日韩| 毛片女人毛片| 免费观看a级毛片全部| 十八禁国产超污无遮挡网站| 午夜亚洲福利在线播放| 精品一区二区三区人妻视频| 18禁动态无遮挡网站| 一边亲一边摸免费视频| 国内精品一区二区在线观看| 国产淫片久久久久久久久| 看非洲黑人一级黄片| 国产欧美日韩精品一区二区| 久久韩国三级中文字幕| 男人舔奶头视频| 国产在线一区二区三区精| 国产亚洲最大av| 在线a可以看的网站| 免费黄色在线免费观看| 99热网站在线观看| 人妻一区二区av| 99久国产av精品国产电影| 国产 一区精品| 中文字幕久久专区| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 色尼玛亚洲综合影院| 激情 狠狠 欧美| 一级毛片黄色毛片免费观看视频| 婷婷色av中文字幕| 亚洲内射少妇av| 亚洲成色77777| 秋霞在线观看毛片| 97超视频在线观看视频| 久久精品国产亚洲av涩爱| 午夜精品在线福利| 少妇丰满av| 小蜜桃在线观看免费完整版高清| 能在线免费看毛片的网站| 人人妻人人澡欧美一区二区| 久久久久久久亚洲中文字幕| 69人妻影院| 中文字幕av成人在线电影| 高清在线视频一区二区三区| 久久久精品欧美日韩精品| 熟女电影av网| 国产v大片淫在线免费观看| 麻豆国产97在线/欧美| 免费av不卡在线播放| 三级国产精品片| xxx大片免费视频| 欧美不卡视频在线免费观看| 免费黄色在线免费观看| 久久久久久久久久黄片| 精品酒店卫生间| 久久人人爽人人片av| 97超视频在线观看视频| www.av在线官网国产| 午夜激情欧美在线| 爱豆传媒免费全集在线观看| 欧美激情久久久久久爽电影| 国产中年淑女户外野战色| 日韩欧美三级三区| 欧美激情久久久久久爽电影| 午夜激情福利司机影院| 亚洲精品色激情综合| 欧美不卡视频在线免费观看| 三级国产精品片| 欧美精品国产亚洲| 国产精品一区二区在线观看99 | 久久热精品热| 国产综合懂色| 最近中文字幕高清免费大全6| 国产永久视频网站| 精品久久久久久久久久久久久| 日本免费a在线| 久久久亚洲精品成人影院| 中文字幕人妻熟人妻熟丝袜美| 国产黄色小视频在线观看| 国产乱人视频| 男女那种视频在线观看| 亚洲国产精品成人综合色| 人人妻人人澡欧美一区二区| 啦啦啦韩国在线观看视频| 一级毛片我不卡| 黄片wwwwww| 久久久久免费精品人妻一区二区| 国产淫片久久久久久久久| 亚洲av.av天堂| 九草在线视频观看|