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

    污染減排、外部性與環(huán)境技術(shù)創(chuàng)新:來自省級環(huán)境專利的證據(jù)

    2014-05-12 10:23:22王文普
    統(tǒng)計與信息論壇 2014年7期
    關(guān)鍵詞:因變量外部性專利

    王文普

    (南通大學(xué) 商學(xué)院,江蘇 南通 226019)

    污染減排、外部性與環(huán)境技術(shù)創(chuàng)新:來自省級環(huán)境專利的證據(jù)

    王文普

    (南通大學(xué) 商學(xué)院,江蘇 南通 226019)

    基于污染外部性機(jī)制,推導(dǎo)出含有污染減排外部性的空間Tobit設(shè)定。利用中國31個省級環(huán)境專利數(shù)據(jù),在控制技術(shù)外部性的情況下,分離出污染減排外部性對環(huán)境技術(shù)創(chuàng)新的影響。貝葉斯MCMC方法估計結(jié)果顯示,污染減排對環(huán)境技術(shù)創(chuàng)新有顯著的正的直接效應(yīng),而且還存在顯著的負(fù)外部效應(yīng),從而污染減排的總效應(yīng)為不顯著的正影響,表明忽略污染減排外部性對環(huán)境技術(shù)創(chuàng)新有重要影響。進(jìn)而探討避免污染減排外部性的政策含義。

    環(huán)境技術(shù)創(chuàng)新;污染減排;外部性;貝葉斯MCMC方法

    一、引 言

    作為世界第二大經(jīng)濟(jì)體,中國在創(chuàng)造經(jīng)濟(jì)奇跡的同時,也帶來了嚴(yán)峻的環(huán)境污染問題。特別是2013年12月持續(xù)的霧霾侵襲了中國東中部地區(qū),全國20個省100多個城市空氣質(zhì)量達(dá)到嚴(yán)重污染。不斷出現(xiàn)的嚴(yán)重霧霾事件表明,中國的環(huán)境污染不是局部,而是全國性的。環(huán)境污染成為中國經(jīng)濟(jì)可持續(xù)發(fā)展的一大障礙。

    緩解經(jīng)濟(jì)發(fā)展對資源的過度依賴和對環(huán)境的損害,科技進(jìn)步與創(chuàng)新是一個有效的途徑。特別是環(huán)境技術(shù),包括降低末端污染排放的技術(shù),也包括生產(chǎn)工藝的技術(shù),都能夠減輕對環(huán)境的不利影響。本文試圖通過設(shè)定帶有外部性的空間Tobit模型,在控制技術(shù)外部性的條件下,考察污染減排對環(huán)境技術(shù)創(chuàng)新的影響,特別是污染減排外部性的影響,對于完善中國的環(huán)境政策,提高技術(shù)自主創(chuàng)新能力,建設(shè)生態(tài)文明社會,推動可持續(xù)發(fā)展,具有重要的現(xiàn)實意義。

    二、文獻(xiàn)評述

    有關(guān)污染減排政策如何引起技術(shù)創(chuàng)新,文獻(xiàn)中至少有兩種基本觀點(diǎn)。一是需求拉動說。該觀點(diǎn)認(rèn)為污染減排政策能夠使企業(yè)的環(huán)境成本內(nèi)部化,引起投入要素的相對價格發(fā)生變化,刺激企業(yè)對低成本減污技術(shù)的需求,進(jìn)而激勵企業(yè)從事環(huán)境技術(shù)的開發(fā)與應(yīng)用,然而由于技術(shù)外部性的存在,使創(chuàng)新者僅能取得一部分的研發(fā)回報,導(dǎo)致企業(yè)技術(shù)創(chuàng)新的動力不足。另一個是政策推力說。它從供給的角度來解釋污染減排政策的技術(shù)創(chuàng)新效應(yīng)。作為一種外部推動力,污染減排政策將促使企業(yè)采用新的更清潔技術(shù),或激發(fā)企業(yè)從事科技研發(fā)的內(nèi)在動力。Taylor等人的觀點(diǎn)是,需求拉動可能會促進(jìn)更多的環(huán)境技術(shù)進(jìn)步與創(chuàng)新[1-2]。

    污染減排的技術(shù)創(chuàng)新效應(yīng)也獲得一定的經(jīng)驗證據(jù)。Popp發(fā)現(xiàn),美國的專利申請數(shù)對污染減排所引起的相對價格變化有顯著的反應(yīng),特別在1972-1994年期間,美國的專利申請數(shù)對投入價格變化有靈敏的反應(yīng);Popp還發(fā)現(xiàn)不斷上升的化石燃料價格,尤其是石油和天然氣價格的上升,提高了這類能源的使用成本,從而導(dǎo)致能源節(jié)約技術(shù)專利數(shù)的增加[3]。Brunnermeir和Cohen也發(fā)現(xiàn)了類似的證據(jù),他們發(fā)現(xiàn)減排支出對美國環(huán)境專利數(shù)有正的影響;Rehfeld等利用德國制造部門的公司水平數(shù)據(jù),通過離散選擇模型,發(fā)現(xiàn)污染減排政策對環(huán)境產(chǎn)品的創(chuàng)新存在顯著的正向影響;Noailly等以1979-2006年荷蘭建筑業(yè)中9個領(lǐng)域的專利數(shù),分析污染減排與建筑業(yè)能源效率改進(jìn)的技術(shù)創(chuàng)新之間關(guān)系,發(fā)現(xiàn)20世紀(jì)90年代中期的嚴(yán)格污染減排政策對建筑企業(yè)的創(chuàng)新有正向激勵作用[4-6]。

    國內(nèi)一些研究也表明,中國的污染減排政策對技術(shù)創(chuàng)新有一定的激勵作用。黃德春等將技術(shù)系數(shù)引入Robert模型,發(fā)現(xiàn)污染減排政策給企業(yè)帶來直接成本的同時,也激發(fā)了技術(shù)創(chuàng)新,從而部分或全部抵消污染減排政策所引起的成本;李斌等使用污染排放水平作為衡量治污技術(shù)創(chuàng)新的指標(biāo),動態(tài)分析顯示,污染減排政策嚴(yán)厲性與治污技術(shù)創(chuàng)新的關(guān)系呈U形,然而,污染減排政策對當(dāng)期治污技術(shù)創(chuàng)新卻有不利影響[7-8]。

    許多研究分析忽略了污染減排所引起的外部效應(yīng)。由于技術(shù)創(chuàng)新也存在外部效應(yīng),這兩種外部性相互交織在一起使得污染減排政策與技術(shù)創(chuàng)新的關(guān)系更為復(fù)雜。本文關(guān)注兩個問題:第一,基于污染外部性機(jī)制,將污染減排的外部性視作影響技術(shù)創(chuàng)新的一個重要因素,從理論上分析污染減排外部性對技術(shù)創(chuàng)新的影響機(jī)制。第二,運(yùn)用中國省級環(huán)境專利數(shù)據(jù)和空間Tobit經(jīng)驗設(shè)定,在控制技術(shù)外部性的條件下,測算中國污染減排及其外部性的技術(shù)創(chuàng)新效應(yīng)。此外,借鑒Levinson的方法,構(gòu)造一個改進(jìn)的污染減排指標(biāo)[9-10]。

    三、模型推導(dǎo)與經(jīng)驗設(shè)定

    根據(jù)政策推力假說,將污染減排視作一種外部驅(qū)動因素,從理論上解釋技術(shù)創(chuàng)新函數(shù)包括污染減排外部效應(yīng)的動機(jī),設(shè)定環(huán)境技術(shù)創(chuàng)新的計量模型。

    (一)含有外部性的模型推導(dǎo)

    污染減排的外部性是由污染溢出引起的,因為許多污染物在地區(qū)之間是可轉(zhuǎn)移的,從而一個地區(qū)的環(huán)境質(zhì)量不僅取決于本地區(qū)污染政策,而且也受其他地區(qū)污染排放政策的影響。根據(jù)政策推力假說,并參照Embora等的研究方法,假定污染減排S及其外部性S*是影響技術(shù)創(chuàng)新的決定因素[9],因而,技術(shù)創(chuàng)新函數(shù)被定義為:

    這里,yit表示i地區(qū)t期的技術(shù)創(chuàng)新指標(biāo);Sit為i地區(qū)t年的污染減排變量表示其他地區(qū)j污染減排所產(chǎn)生的外部效應(yīng)。

    由于污染減排的外部效應(yīng),使得式(1)中S和S*在地區(qū)間相互影響。為了衡量這種相互作用,再假設(shè)S和S*具有式(2)、(3)和(4)所表示的地區(qū)間相互聯(lián)系,并引入具有零均值和常方差的擾動項u、v和e,則有:

    其中,b1=b+γ,θ1=-ρb-λγ。式(5)描述了污染減排外部性引起的地區(qū)之間的相互聯(lián)系,文獻(xiàn)稱之空間Durbin模型(SDM)。SDM是空間模型設(shè)定的一般形式,其包含著空間滯后因變量模型(SAR)和空間誤差自回歸模型(SEM)。特別地,當(dāng)γ=0時,這意味著對影響S的隨機(jī)沖擊與影響S*的隨機(jī)沖擊之間沒有相關(guān)性,SDM可簡化為SEM;當(dāng)θ1=-ρb1成立時,SDM也可退化為SEM;當(dāng)θ1=0時,SDM就變?yōu)镾AR。

    上述推導(dǎo)過程說明了技術(shù)創(chuàng)新函數(shù)中包含空間滯后因變量和滯后解釋變量,有助于解釋污染減排外部性在地區(qū)間的相互作用。

    (二)經(jīng)驗設(shè)定與相關(guān)估計問題

    本文選取環(huán)境專利授權(quán)數(shù)作為環(huán)境技術(shù)創(chuàng)新的代理變量,環(huán)境專利數(shù)是離散的,特別在所使用的樣本(NT=279)中,有180個零值,其占樣本總數(shù)的比例近66%。LeSage和Pace指出,因變量在零處截斷有一定的風(fēng)險。如果忽略零值,將有可能造成一些信息損失,因而建議使用Tobit模型。該模型不僅考慮零值因變量的非正態(tài)分布特征,而且將因變量視作潛變量,以表示與專利決策有關(guān)的不利因素[11]。

    考慮到技術(shù)外部效應(yīng),為了分離出污染減排外部效應(yīng),需要控制技術(shù)外部性。因而,將技術(shù)及其外部性直接添加到式(5)之中,從而得到經(jīng)驗設(shè)定:

    上述設(shè)定稱之空間Tobit模型,它是空間Durbin模型的一個變體,因變量是離散變量而非連續(xù)變量。如果y*it$0,則yit=0;否則=y(tǒng)it。因變量為各省的環(huán)境專利授權(quán)數(shù)。X = (LEV,KLEV,ES),這里L(fēng)EV表示各省的污染減排政策指標(biāo);KLEV代表各省的技術(shù)水平,以捕捉技術(shù)水平對環(huán)境專利的貢獻(xiàn);ES為規(guī)??刂谱兞?。b和θ為待估參數(shù)。假設(shè)隨機(jī)擾動項e~ N(0,σ2V),其中V =diag(v1,v2,…,vn)。參數(shù)ρ表示地區(qū)間的聯(lián)系強(qiáng)度。

    經(jīng)驗設(shè)定中包含滯后因變量和滯后解釋變量,以解釋污染減排和技術(shù)的外部性。如果空間權(quán)重矩陣W 被行標(biāo)準(zhǔn)化,則W·y*表示“相鄰”省的平均技術(shù)創(chuàng)新水平,W·LEV和W·KLEV分別為“相鄰”省的污染減排和技術(shù)外部性的平均水平。此外,回歸中包含滯后因變量以反映可能與外部性有關(guān)的遺漏變量問題。在估計之前,對有關(guān)的估計問題做簡要說明。

    第一,估計方法。由于因變量是受限制的,故不能直接使用解析法。這里采用貝葉斯法來估計受限因變量。與EM法相比,貝葉斯法將零觀察處截斷的不可觀測因變量視為待估參數(shù),并通過MCMC程序來估計模型的參數(shù)。為了獲得觀察處截斷的不可觀測因變量,貝葉斯法使用m步Gibbs抽樣,從多變量的截斷正態(tài)分布中生成隨機(jī)樣本。在獲得不可觀測因變量估計值后,將估計問題轉(zhuǎn)換為非截斷數(shù)據(jù)的估計問題,再使用最大似然法。此外,貝葉斯法也很容易處理數(shù)據(jù)中異方差和極端值問題,且計算簡便,因而選擇LeSage、LeSage和Pace提出的貝葉斯法來估計上述空間Tobit模型。具體來看,在超參數(shù)q滿足q/vi~i.d.χ2(q)/q(其中q~Gamma(m,k)),1/σ2~ Gamma(μ,d0)和ρ~U(min,max)的假設(shè)下①根據(jù)LeSage(2010)提供的Matlab估計程序,對于分布中的參數(shù),先驗設(shè)定q=4,m=0;μ=0,d0=0;min=0,max=1。,貝葉斯法的步驟如下:(i)利用參數(shù)的任意初始值②在LeSage(2010)提供的 Matlab估計程序中,取ρ=0.5,b=0,σ2 =1,vi=0。,記ρ0,b0,v0i,從一個隨機(jī)χ2分布中抽樣來確定σ估計值,記σ1。進(jìn)而計算σ的條件概率p(σ|ρ,b,V)。(ii)利用ρ0,b0,σ1來確定b的均值和方差。根據(jù)這個均值與方差進(jìn)行多變量正態(tài)隨機(jī)抽樣來確定b估計量,記b1。(iii)根據(jù)ρ0,b1,σ1計算vi的條件后驗分布,進(jìn)而通過隨機(jī)χ2(q+1)分布抽取n個矢量來確定vi,記v1i(i=1,2,…,n)。(iv)給定b1,σ1,v1i值,根據(jù)ρ 的 條 件 概 率 公 式 p(ρ|b,σ,V),使 用Metropolis抽樣確定ρ,記ρ1。(v)利用ρ1,b1,σ1,v1i所確定的計算均值和方差的截斷正態(tài)分布中抽取截尾yi觀察值。上述步驟構(gòu)成一次Gibbs抽樣過程,經(jīng)過多次重復(fù)抽樣,在舍去抽樣初始估計值后,收集每次抽樣中的參數(shù),進(jìn)而得到后驗參數(shù)估計及其后驗分布[11-12]。

    四、變量與數(shù)據(jù)

    本文所使用的數(shù)據(jù)是2001-2009年中國31個省區(qū)的面板數(shù)據(jù),數(shù)據(jù)來源于各年的《中國環(huán)境年鑒》、《中國統(tǒng)計年鑒》和《中國工業(yè)統(tǒng)計年鑒》。

    省區(qū)環(huán)境技術(shù)創(chuàng)新是因變量y。文獻(xiàn)中測度技術(shù)創(chuàng)新常用指標(biāo)有研發(fā)支出和專利數(shù)兩個指標(biāo),由于《中國統(tǒng)計年鑒》和《環(huán)境年鑒》都沒有報告各省的環(huán)境研發(fā)支出數(shù)據(jù),環(huán)境年鑒從2002年起報告了各省的環(huán)境專利授權(quán)數(shù),有關(guān)研究表明,企業(yè)的專利活動類似于研發(fā)支出,而且它還可能與外部事件(如政府政策)有關(guān),因而,選取環(huán)境專利數(shù)來衡量環(huán)境技術(shù)創(chuàng)新,以捕捉各省的環(huán)境技術(shù)創(chuàng)新活動。

    污染減排指標(biāo)LEV,用來捕捉各省的污染減排政策的嚴(yán)厲性。污染減排政策的嚴(yán)厲性可以從多個角度來衡量,例如,中國科學(xué)研究院可持續(xù)發(fā)展戰(zhàn)略研究組從污染排放量、環(huán)境治理強(qiáng)度、氣候變化、土壤侵蝕等不同方面來衡量環(huán)境政策。早期研究使用污染減排支出除以總產(chǎn)出來構(gòu)建污染減排政策嚴(yán)厲性指標(biāo)。因缺少省區(qū)各工業(yè)行業(yè)部門的污染減排支出數(shù)據(jù),借鑒Levinson的方法來構(gòu)建各省的污染減排強(qiáng)度[10][13]133-134[14]。先計算單位產(chǎn)出的污染減排支出Sit=Pit/Qit,其中Pit為i省t年的污染治理支出,用各省的工業(yè)污染治理年投資完成額、工業(yè)廢水和工業(yè)廢氣設(shè)施運(yùn)行費(fèi)用三項之和來衡量;Qit為i省t年的工業(yè)增加值。由于省際間的工業(yè)構(gòu)成是不同的,因而,預(yù)測的單位產(chǎn)出污染減排支出為=,這里Qt為t年的全國工業(yè)增加值,n為省區(qū)總數(shù),從而得到經(jīng)工業(yè)構(gòu)成調(diào)整的污染減排指標(biāo)LEVit=Sit/。如果LEVit>1,則表明i省t年比其他省有更多的環(huán)境支出。LEV值越大,則該省的污染減排政策越嚴(yán)厲,從而產(chǎn)生較低的污染排放。計算該指標(biāo)時,分別用工業(yè)增加值指數(shù)和地區(qū)GDP平減指數(shù)換算為1995年價格。

    技術(shù)水平KLEV。使用各省的國內(nèi)三種專利授權(quán)數(shù)的對數(shù)來表示技術(shù)水平,各省專利授權(quán)數(shù)由國內(nèi)的發(fā)明專利、實用新型專利和外觀設(shè)計專利三項構(gòu)成。專利一定程度上反映了發(fā)明活動,與研發(fā)總支出相比,專利總數(shù)被視為與技術(shù)知識水平的關(guān)系更為密切,因而,用各省的國內(nèi)三種專利數(shù)(包括環(huán)境技術(shù)專利數(shù))來表示各省技術(shù)水平,它從供給層面上捕捉“知識池”對環(huán)境技術(shù)創(chuàng)新的影響[15]。

    控制變量ES。工業(yè)是一個地區(qū)經(jīng)濟(jì)發(fā)展的主要貢獻(xiàn)部門,也是污染減排的重點(diǎn)對象,因而選取工業(yè)增加值作為各省的規(guī)模變量[16]。ES為各省人均工業(yè)增加值。ES變量不僅捕捉了各省工業(yè)規(guī)模的影響,而且,根據(jù)需求拉動說的觀點(diǎn),一個省區(qū)的工業(yè)規(guī)模越大,其污染減排的壓力就越大,該省對環(huán)境技術(shù)創(chuàng)新需求越強(qiáng)烈。所有解釋變量均取自然對數(shù)。

    五、結(jié)果與分析

    (一)估計結(jié)果

    在估計之前,檢驗?zāi)P椭泻诵慕忉屪兞縇EV和KLEV是否具有空間自相關(guān),進(jìn)行空間自回歸檢驗,設(shè)定LEV=a+λW·LEV+u。檢驗結(jié)果表明W·LEV系數(shù)的t統(tǒng)計量為2.33,在5%水平上顯著;W·KLEV系數(shù)的t統(tǒng)計量為5.45,在1%水平上高度顯著。與上文假設(shè)的一樣,污染減排變量和技術(shù)變量存在空間相關(guān)性,初步表明經(jīng)驗設(shè)定中包含滯后解釋變量是合適的。

    表1 空間Tobit估計表

    表1報告了空間Tobit模型的估計結(jié)果。其來自于空間Tobit模型合并面板數(shù)據(jù)的貝葉斯估計,它是在Gibbs抽樣中使用一步且在MCMC中使用1 500次抽樣①為了避免估計參數(shù)的相關(guān)性,舍去(burn-in)前500個參數(shù)估計值。匿名審稿人建議做收斂性檢驗,但`Raftery等(1992)指出,為了在累積后驗分布中達(dá)到期望的精確水平,提供了一種確定抽樣次數(shù)的方法來滿足這一要求;根據(jù)這一建議,LeSage(1997)發(fā)現(xiàn)在通常的空間回歸模型應(yīng)用中需要大約進(jìn)行1 100次抽樣就達(dá)到所設(shè)定收斂性的要求(見LeSage 2000)。因而,我們通過設(shè)定抽樣次數(shù)來保證收斂性要求。。由于所有解釋變量都取自然對數(shù),系數(shù)可理解為半彈性。

    表中最大似然估計為未考慮專利計數(shù)因零值所引起的數(shù)據(jù)非正態(tài)分布。與空間Tobit估計結(jié)果相比(見表1中第3和4列),ML估計中大多數(shù)解釋變量的系數(shù)絕對值小,表明忽略零觀察將導(dǎo)致參數(shù)估計向下偏倚,即在ML估計中,專利對解釋變量的變化有較小的反應(yīng),意味著忽略零值觀察對參數(shù)的顯著性有較大影響。

    Tobit1中僅考慮污染減排外部性對環(huán)境技術(shù)創(chuàng)新的影響,但不包括技術(shù)外部性。滯后因變量系數(shù)ρ為0.139,且顯著性較低(顯著水平僅為15%)。Tobit2估計同時考慮了污染減排和技術(shù)的外部效應(yīng),與Tobit1的估計結(jié)果相比較,LEV和ES的系數(shù)值有很大差異,表明遺漏技術(shù)外部性有可能放大污染減排的影響,因而Tobit2的結(jié)果可能比較真實地反映污染減排對環(huán)境技術(shù)創(chuàng)新的影響。

    在兩個Tobit估計中,參數(shù)ρ均為不顯著的正號,表明地區(qū)之間存在不顯著的正的相關(guān)關(guān)系。此外,回歸方程中包含滯后因變量和滯后解釋變量,解釋變量的系數(shù)并不直接對應(yīng)于所關(guān)注的直接和外部效應(yīng)。因而,使用LeSage和Pace提出的方法,計算解釋變量的直接效應(yīng)和外部效應(yīng),有助于進(jìn)行更合理的判斷[11]。下面集中于討論 Tobit2的估計結(jié)果。

    從Tobit2中可以發(fā)現(xiàn),參數(shù)ρ在常規(guī)顯著水平上不顯著,LEV的系數(shù)及其滯后項系數(shù)卻都是高度顯著的,從而系數(shù)限制θ1=-ρb1不成立。然而,KLEV的系數(shù)b2是顯著的,W·KLEV的參數(shù)θ2在統(tǒng)計上則不顯著,故θ2=-ρb2有可能成立。類似地,對于規(guī)模變量ES來說,ES的系數(shù)b3和W·ES的系數(shù)θ3都顯著大于0,從而系數(shù)約束θ3=-ρb3也將不成立。表明數(shù)據(jù)中存在污染減排的外部性,而技術(shù)外部性似乎并不明顯,進(jìn)一步表明空間Tobit設(shè)定是合理的。

    (二)穩(wěn)健性分析

    Tobit2的估計結(jié)果是否可靠,進(jìn)行了兩種不同的敏感性檢驗。一個檢驗是,使用不同的空間權(quán)重矩陣來衡量地區(qū)間不同的聯(lián)系程度。選取經(jīng)濟(jì)加權(quán)矩陣和4階權(quán)重矩陣來檢驗結(jié)果的穩(wěn)健性。表2是使用不同空間權(quán)重的估計結(jié)果。結(jié)果表明:系數(shù)ρ仍為不顯著的正值。LEV和KLEV的系數(shù)符號和顯著性都沒有發(fā)生變化。在經(jīng)濟(jì)加權(quán)矩陣中,ES的系數(shù)顯著性有所下降,但W·LEV和W·KLEV的系數(shù)符號和顯著性均沒有出現(xiàn)逆轉(zhuǎn)。然而,在4階權(quán)重矩陣中,ES的系數(shù)顯著變得不顯著;W·LEV和W·KLEV的系數(shù)符號仍保持不變,但W·LEV的系數(shù)變得不顯著。表明Tobit2中的結(jié)果是比較可靠的。

    表2 不同權(quán)重的空間設(shè)定檢驗表

    另一個檢驗是選擇不同的抽樣方法。具體地講,逐步增加Gibbs抽樣中不同的步數(shù)m和MCMC程序中的抽樣次數(shù)n。其暗含的基本假設(shè)是,如果推斷是相同的,那么估計量可視為是一致的。改變Gibbs抽樣步數(shù)m,目的是檢驗替代不可觀測因變量(y*$0)參數(shù)矢量的計算精確性;其二是改變抽樣次數(shù),為了檢驗待估參數(shù)值的穩(wěn)定性,以便比較不同抽樣次數(shù)(這里取n=1 500和n=15 000)的估計結(jié)果[11]。

    表3給出了不同抽樣的估計結(jié)果。分別進(jìn)行1步、5步和10步的Gibbs抽樣,而且抽樣數(shù)為1 500次和15 000次的估計結(jié)果。系數(shù)ρ仍為正,LEV和KLEV的系數(shù)仍是高度顯著的,這進(jìn)一步表明基準(zhǔn)估計結(jié)果(見表1中的第4列)是可靠的。

    (三)外部效應(yīng)的量化

    參數(shù)ρ不僅不顯著,而且解釋變量的系數(shù)值并不反映直接效應(yīng)和外部效應(yīng),所以有必要使用LeSage和Pace提供的新方法來計算它們的大小和顯著性,以量化地區(qū)間的聯(lián)系強(qiáng)度[11]121-125[15]。

    與線性回歸估計不同,在帶有滯后因變量的空間設(shè)定中,系數(shù)值并不直接表示為半彈性。系數(shù)的正確解釋應(yīng)考慮模型中非線性關(guān)系中自偏導(dǎo)數(shù)(直接效應(yīng))和交叉偏導(dǎo)數(shù)(外部效應(yīng))。事實上,經(jīng)驗設(shè)定中隱含著一種非線性關(guān)系,空間Tobit設(shè)定可改寫為:

    其中,X = (LEV,KLEV,ES)′,b= (b0,b1,b2),θ=(θ1,θ2),In為NT×NT 的單位矩陣??梢姡瑈*/x′≠b,而是y*/x′= (In-ρW)-1(Inbk+Wθk),這是因為“相鄰”地區(qū)的解釋變量變化所引起的地區(qū)間交互作用而產(chǎn)生的。為了解釋設(shè)定中變量的影響,需要計算變量的直接效應(yīng)、外部效應(yīng)和總效應(yīng)。事實上,地區(qū)i解釋變量的變化不僅影響該地區(qū)的因變量y*,而且通過外部性影響到“相鄰”地區(qū)的因變量y*。具體計算詳見LeSage和Pace[11]121-125。

    表3 不同抽樣方法的空間設(shè)定檢驗表

    表4 直接效應(yīng)、外部效應(yīng)和總效應(yīng)的估計表

    表4計算了基于Tobit2(見表2中第4列)中解釋變量變化引起的直接效應(yīng)、外部效應(yīng)和總效應(yīng)。除各種效應(yīng)的平均值外,還報告了下限0.01和上限0.99的置信區(qū)間。這些效應(yīng)的t值和置信區(qū)間是從MCMC抽樣(取1 000次)得到的后驗分布中構(gòu)建的。從表4中可以看出:

    第一,直接效應(yīng),LEV的直接效應(yīng)為10.82,且在1%水平上高度顯著。當(dāng)一個省的污染減排強(qiáng)度每增加1%,則該省的環(huán)境專利數(shù)將增加10.82個單位。此外,污染減排的直接效應(yīng)并不等于表1中LEV的估計系數(shù)為10.735(見表1中第4列),二者之差為0.085(占直接效應(yīng)的0.79%)。這捕捉了來自于省際間相互作用所引起的反饋效應(yīng),這種反饋效應(yīng)部分來自于滯后因變量,部分源自于解釋變量的滯后項WX。從置信區(qū)間來看,LEV的直接效應(yīng)總為正的,表明提高污染減排支出對本省區(qū)的環(huán)境技術(shù)創(chuàng)新具有直接的正向激勵,這與政策推力說和需求拉動說的預(yù)測相一致,也與很多經(jīng)驗分析的結(jié)果相同。

    KLEV的直接效應(yīng)為5.792,是LEV直接效應(yīng)的近50%,且統(tǒng)計上高度顯著,表明技術(shù)知識的總水平對環(huán)境技術(shù)創(chuàng)新產(chǎn)生了顯著的積極作用。技術(shù)水平每提高1%,則該省區(qū)環(huán)境專利數(shù)將增加5.8個單位。技術(shù)的反饋效應(yīng)等于-0.224,是其直接效應(yīng)的-3.87%。從99%的置信區(qū)間來看,技術(shù)水平對本省區(qū)的環(huán)境技術(shù)創(chuàng)新的影響總是正的。因而,一個地區(qū)如果實施有利于提高污染減排和促進(jìn)技術(shù)進(jìn)步的公共政策,那么將對該地區(qū)的環(huán)境技術(shù)創(chuàng)新產(chǎn)生正影響。

    ES的直接效應(yīng)顯著為負(fù)。ES每增加1%,環(huán)境專利數(shù)將減少9.3個單位,這與需求拉動說的預(yù)期不同。也就是說,工業(yè)規(guī)模越大,越不利于該地區(qū)的環(huán)境技術(shù)創(chuàng)新。一個可能的解釋是,工業(yè)作為一個地區(qū)經(jīng)濟(jì)增長的重要貢獻(xiàn)部門,在追求經(jīng)濟(jì)增長的驅(qū)動下,放松污染減排政策的執(zhí)行,導(dǎo)致企業(yè)主動減排的動力不足,從而削弱企業(yè)技術(shù)創(chuàng)新的激勵。

    第二,外部效應(yīng),它則衡量了對“相鄰”地區(qū)影響的累積程度。LEV的外部效應(yīng)為-10.716,是其直接效應(yīng)的99%,且統(tǒng)計上高度顯著。如果一個省區(qū)提高污染減排強(qiáng)度,不僅影響本地區(qū)的環(huán)境專利數(shù),而且也影響“鄰居”省區(qū)的環(huán)境專利數(shù);后者與前者變動之比接近1∶-1,表明提高污染減排強(qiáng)度,促進(jìn)了本地區(qū)的環(huán)境技術(shù)創(chuàng)新,卻對相鄰地區(qū)的環(huán)境技術(shù)創(chuàng)新產(chǎn)生不利影響。一個比較合理的解釋是,當(dāng)一個省區(qū)比鄰近省區(qū)實施更嚴(yán)厲的污染減排政策時,推動了本地區(qū)內(nèi)企業(yè)的技術(shù)創(chuàng)新活動,然而,相鄰地區(qū)中有相對較低的污染減排成本,特別對那些減污成本較高的企業(yè)來說,將使企業(yè)考慮重新選址,從而削弱企業(yè)從事環(huán)境技術(shù)創(chuàng)新的內(nèi)在動力。

    技術(shù)外部效應(yīng)為不顯著的正數(shù)0.75,正的技術(shù)外部效應(yīng)也被很多技術(shù)創(chuàng)新文獻(xiàn)研究所支持。統(tǒng)計上不顯著的一個原因是,目前中國的技術(shù)以模仿、吸收消化者為主,自主創(chuàng)新者不多,從而使技術(shù)外部效應(yīng)不明顯[17]。

    ES變量的外部效應(yīng)約為24.5,且統(tǒng)計上高度顯著。即工業(yè)規(guī)模對“鄰居”地區(qū)產(chǎn)生了正的外部性,因為一個地區(qū)規(guī)模越大,對相鄰地區(qū)的工業(yè)構(gòu)成一定的威脅,加劇了省際間企業(yè)的競爭,從而刺激它們進(jìn)行環(huán)境技術(shù)創(chuàng)新的積極性。

    此外,盡管參數(shù)ρ是不顯著的,但污染減排和規(guī)模變量的外部效應(yīng)都是高度顯著的。因而,僅由參數(shù)ρ的顯著性來判斷是否存在空間作用,有可能做出誤導(dǎo)性推斷。

    第三,總效應(yīng),反映解釋變量的變化對環(huán)境技術(shù)創(chuàng)新的總體影響。污染減排的總效應(yīng)為0.104,且統(tǒng)計上不顯著。主要原因在于,污染減排對環(huán)境技術(shù)創(chuàng)新存在相互抵消的兩種效應(yīng),二者作用的結(jié)果使得總效應(yīng)數(shù)值較小,且統(tǒng)計上也不顯著。這也表明,污染減排在地區(qū)間相互影響的重要性。如果不控制污染減排外部性,可能會對污染減排所引起的技術(shù)創(chuàng)新效應(yīng)做出不合理的推斷。

    技術(shù)水平的總效應(yīng)為6.543,且在5%水平上高度顯著。這意味著,當(dāng)一個省區(qū)的技術(shù)水平每增加1%,該省區(qū)的環(huán)境專利數(shù)將增加6.543個單位。主要原因在于,技術(shù)的直接效應(yīng)占據(jù)支配地位,使KLEV的總效應(yīng)具有較大的顯著的正作用。

    規(guī)模變量對環(huán)境技術(shù)創(chuàng)新的總效應(yīng)顯著。盡管工業(yè)規(guī)模的直接效應(yīng)為負(fù),但其外部效應(yīng)顯著為正,并且遠(yuǎn)遠(yuǎn)大于其直接效應(yīng),使其總效應(yīng)仍是正的,表明工業(yè)規(guī)模總體上對環(huán)境技術(shù)創(chuàng)新具有顯著的正向作用。

    六、結(jié)論與政策啟示

    基于污染外部性機(jī)制,構(gòu)建包含污染減排和技術(shù)兩種外部性的空間Tobit經(jīng)驗設(shè)定,在控制技術(shù)外部性的條件下,測算出污染減排及其外部性對環(huán)境技術(shù)創(chuàng)新的影響。結(jié)果顯示:第一,一個地區(qū)的環(huán)境技術(shù)創(chuàng)新不僅取決于本地區(qū)的污染減排強(qiáng)度,而且也受鄰近地區(qū)污染減排政策的影響,這意味著污染減排外部性成為地區(qū)間相互影響的一個重要機(jī)制。第二,污染減排在地區(qū)間引起相當(dāng)大的負(fù)外部性。污染減排的外部效應(yīng)大小在數(shù)值上幾乎相當(dāng)于其直接效應(yīng),且二者的作用方向截然相反,使得污染減排的總效應(yīng)在統(tǒng)計上并不顯著,進(jìn)一步表明如果忽略污染減排的外部性,有可能對污染減排的作用做出不恰當(dāng)?shù)耐茢?。第三,在環(huán)境技術(shù)創(chuàng)新過程中,污染減排和技術(shù)等公共政策對環(huán)境技術(shù)創(chuàng)新具有顯著的直接的積極影響。第四,滯后因變量系數(shù)統(tǒng)計上并不顯著,但不一定意味著空間外部性也不顯著,因而,僅由參數(shù)ρ的顯著性判斷地區(qū)間相互作用是否存在,有可能會做出不合理的推斷。

    估計結(jié)果具有明顯的政策含義。為了充分發(fā)揮污染減排在技術(shù)創(chuàng)新中的積極作用,應(yīng)強(qiáng)化和完善污染減排機(jī)制:第一,加大與環(huán)境有關(guān)的區(qū)域間財政轉(zhuǎn)移支付力度,中央政府應(yīng)盡快完善省際間環(huán)境補(bǔ)償?shù)牧⒎ê蜆?biāo)準(zhǔn),以協(xié)調(diào)并指導(dǎo)省際間和區(qū)域間的環(huán)境合作,以減輕甚至消除污染減排的負(fù)外部效應(yīng)。第二,強(qiáng)化污染減排政策的執(zhí)行力度,嚴(yán)格污染減排政策的監(jiān)督和實施,促進(jìn)環(huán)境技術(shù)創(chuàng)新需求機(jī)制的形成。第三,提高對技術(shù)領(lǐng)域,特別是環(huán)境領(lǐng)域的公共研究支持力度,增強(qiáng)企業(yè)的自主創(chuàng)新能力,提高技術(shù)“知識池”對環(huán)境技術(shù)創(chuàng)新的“輸血”和“造血”功能。

    [1] Jaffe A B,Newell R G,Stavins T N.A Tale of Two Market Failures:Technology and Environmental Policy[J].Ecological Economics,2005(2).

    [2] Taylor M R,Rubin E S,Nemet G F.The Role of Technological Innovation in Meeting California’s Greenhouse Gas Emission Targets[C].Chapter Three in Hanemann,M.& A.Farrell(eds.).Managing Greenhouse Gases in California.Report Prepared for the Energy Foundation and the Hewlett Foundation.2006.

    [3] Popp D.Induced Innovation and Energy Prices[J].American Economic Review,2002(1).

    [4] Brunnermeier S B,Cohen M A.Determinants of Environmental Innovation in theU.S.Manufacturing Industries[J].Journal of Environmental Economics and Management,2003(2).

    [5] Rehfeld K M,Rennings M,Ziegler A.Integrated Product Policy and Environmental Product Innovations:An Empirical Analysis[J].Ecological Economics,2007,61(1).

    [6] Noailly J,S Batrakova R Lukach.The Impact of Environmental Policy on Energy-efficient Innovations in Buildings[C].C PB.Working Paper,2008.

    [7] 黃德春,劉志彪.環(huán)境規(guī)制與企業(yè)自主創(chuàng)新:基于波特假設(shè)的企業(yè)競爭優(yōu)勢構(gòu)建[J].中國工業(yè)經(jīng)濟(jì),2006(3).

    [8] 李斌,彭星,陳柱華.環(huán)境規(guī)制、FDI與中國治污技術(shù)創(chuàng)新:基于省際動態(tài)面板數(shù)據(jù)的分析[J].財經(jīng)研究,2011(10).

    [9] Embora N,T P.Mamuneas T Stengos Air Pollution Spillovers and U.S.state Productivity Growth.Working paper.University of Cyprus.2008(8).

    [10]Levinson A.An Industry-Adjusted Index of State Environmental Compliance Costs[R].NBER working paper,2001(7297).

    [11]LeSage J,Pace R K.Introduction to Spatial Econometrics[M].Chapman and Hall/CRC .New York.2009.

    [12]LeSage J P.Bayesian Estimation of Limited Dependent Variable Spatial Autoregressive Models[J].Geographical Analysis,2000(1).

    [13]中國科學(xué)研究院可持續(xù)發(fā)展戰(zhàn)略研究組.2007中國可持續(xù)發(fā)展戰(zhàn)略報告——水:治理與創(chuàng)新[M].北京:科學(xué)出版社,2007.

    [14]Konisky D M.Woods N D.Measuring State Environmental Policy[C].Working paper.Georgetown University.2012(1).

    [15]LeSage J P,F(xiàn)ischer M M.Estimates of the Impact of Static and Dynamic Knowledge Spillovers on Regional Factor Productivity[J].International Regional Science Review,2012(1).

    [16]Sauquet A,Marchand A,F(xiàn)eres G.Ecological Fiscal Incentives and Spatial Strategic Interactions :the Case of the ICMS-E inParana[J].CERDI Working Paper,2012(5).

    [17]王家庭、張俊韜.中國城市環(huán)保行業(yè)的技術(shù)效率研究.統(tǒng)計與信息論壇,2012(12).

    Pollution Abatement,Externalities and Environmental Technology Innovation:Evidences from Provincial Environmental Patents Data

    WANG Wen-pu
    (Business of School,Nantong University,Nantong 226019,China)

    Based on pollution externality,the paper derives the spatial Tobit model with pollution abatement externalities.Moreover,controlling the externality of technology,we isolate the impacts of pollution abatement externality on environmental technology innovation using the data of Chinese thirtyone provincial environmental patents.The estimation of the new approach of Bayesian Markov chain Monte Carlo shows that not only is there a significant positive direct effect of pollution abatement on environmental technology innovation,but does a significant negative externality as well,and the total impact of pollution abatement is a non-significantly positive effect,which indicates that pollution abatement externality exists an important influence on the environmental technology innovation.Furthermore we explore policy implications to avoid the externalities of pollution abatement.

    environmental technology innovation;pollution abatement;externalities;Bayesian MCMC approach

    F062;X196

    A

    1007-3116(2014)07-0095-08

    2014-03-05;修復(fù)日期:2014-05-07

    南通大學(xué)博士基金項目《中國工業(yè)綠色轉(zhuǎn)型路徑研究》(03080587)。

    王文普,男,安徽六安人,經(jīng)濟(jì)學(xué)博士,副教授,研究方向:公共經(jīng)濟(jì)與公共政策、經(jīng)濟(jì)增長。

    (責(zé)任編輯:張愛婷)

    猜你喜歡
    因變量外部性專利
    采暖區(qū)既有建筑節(jié)能改造外部性分析與應(yīng)用研究
    專利
    調(diào)整有限因變量混合模型在藥物經(jīng)濟(jì)學(xué)健康效用量表映射中的運(yùn)用
    中國藥房(2022年7期)2022-04-14 00:34:30
    環(huán)境負(fù)外部性的環(huán)境法新解析
    發(fā)明與專利
    傳感器世界(2019年4期)2019-06-26 09:58:44
    適應(yīng)性回歸分析(Ⅳ)
    ——與非適應(yīng)性回歸分析的比較
    偏最小二乘回歸方法
    既有商品住宅節(jié)能改造外部性及其測算研究
    湘潭市農(nóng)地正外部性價值量估算
    專利
    99在线视频只有这里精品首页| 久久午夜综合久久蜜桃| 人人澡人人妻人| 一级毛片女人18水好多| 国产伦一二天堂av在线观看| 亚洲国产精品久久男人天堂| 亚洲国产精品合色在线| 18禁观看日本| 久久婷婷人人爽人人干人人爱 | 日本五十路高清| 亚洲第一av免费看| 黄色视频,在线免费观看| 国产伦人伦偷精品视频| 一边摸一边抽搐一进一小说| 欧美成人免费av一区二区三区| www国产在线视频色| 国产成人精品久久二区二区免费| 97人妻天天添夜夜摸| 每晚都被弄得嗷嗷叫到高潮| 男女下面进入的视频免费午夜 | 巨乳人妻的诱惑在线观看| 欧美在线黄色| 国产精品一区二区三区四区久久 | 成人免费观看视频高清| 免费搜索国产男女视频| 中出人妻视频一区二区| 亚洲天堂国产精品一区在线| 亚洲va日本ⅴa欧美va伊人久久| 国产精品影院久久| 日本精品一区二区三区蜜桃| 人成视频在线观看免费观看| 美女扒开内裤让男人捅视频| 最新美女视频免费是黄的| 精品国产国语对白av| 99国产精品99久久久久| 国产精品99久久99久久久不卡| 久久久久亚洲av毛片大全| 在线观看免费视频网站a站| 午夜福利,免费看| 黄片小视频在线播放| 一级片免费观看大全| 亚洲无线在线观看| 免费在线观看影片大全网站| 级片在线观看| 久久这里只有精品19| 亚洲精品久久国产高清桃花| 丝袜美腿诱惑在线| 国产一区二区在线av高清观看| 午夜福利一区二区在线看| 在线视频色国产色| 国产精品影院久久| 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女| 亚洲熟妇中文字幕五十中出| 女生性感内裤真人,穿戴方法视频| 久久久久久久精品吃奶| www国产在线视频色| 一级毛片精品| 久热爱精品视频在线9| 国产主播在线观看一区二区| 日韩av在线大香蕉| 97超级碰碰碰精品色视频在线观看| 一区二区三区国产精品乱码| 亚洲成人久久性| 涩涩av久久男人的天堂| 国产蜜桃级精品一区二区三区| √禁漫天堂资源中文www| 制服诱惑二区| 天堂√8在线中文| 老熟妇乱子伦视频在线观看| 美女扒开内裤让男人捅视频| 如日韩欧美国产精品一区二区三区| 视频在线观看一区二区三区| 免费女性裸体啪啪无遮挡网站| 精品国产乱码久久久久久男人| 国产一区二区三区在线臀色熟女| 国产精品九九99| 琪琪午夜伦伦电影理论片6080| 99国产精品免费福利视频| 99国产精品免费福利视频| 黄片小视频在线播放| 天堂动漫精品| av在线天堂中文字幕| 非洲黑人性xxxx精品又粗又长| svipshipincom国产片| 国产亚洲精品一区二区www| 亚洲一区二区三区色噜噜| 999久久久国产精品视频| 亚洲熟妇中文字幕五十中出| 欧美成人一区二区免费高清观看 | 禁无遮挡网站| 中文字幕av电影在线播放| 久久精品亚洲熟妇少妇任你| 搡老岳熟女国产| av视频在线观看入口| 午夜老司机福利片| 国产99久久九九免费精品| 欧美精品啪啪一区二区三区| 丰满的人妻完整版| 一级,二级,三级黄色视频| 女人被狂操c到高潮| 午夜免费鲁丝| 日韩精品中文字幕看吧| 久久精品国产清高在天天线| 免费看十八禁软件| 国产激情久久老熟女| av免费在线观看网站| 亚洲欧美激情在线| 亚洲精品久久成人aⅴ小说| 深夜精品福利| 欧美成狂野欧美在线观看| 亚洲成国产人片在线观看| 国产精品乱码一区二三区的特点 | a在线观看视频网站| 最新美女视频免费是黄的| 欧美绝顶高潮抽搐喷水| 日韩精品中文字幕看吧| 在线观看一区二区三区| 中国美女看黄片| 日本vs欧美在线观看视频| 啦啦啦观看免费观看视频高清 | av超薄肉色丝袜交足视频| 成年女人毛片免费观看观看9| 搡老熟女国产l中国老女人| 伊人久久大香线蕉亚洲五| 午夜福利,免费看| 中亚洲国语对白在线视频| 午夜老司机福利片| 一二三四在线观看免费中文在| 国产又爽黄色视频| 国产成人精品久久二区二区91| 非洲黑人性xxxx精品又粗又长| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品999在线| 久久久久久大精品| 欧美另类亚洲清纯唯美| 韩国av一区二区三区四区| av天堂在线播放| 亚洲中文av在线| 午夜免费激情av| 国产精品电影一区二区三区| 亚洲一区二区三区色噜噜| 亚洲成人国产一区在线观看| 日韩大码丰满熟妇| 麻豆一二三区av精品| 亚洲欧美激情在线| 亚洲精品在线观看二区| 一二三四在线观看免费中文在| 国产激情久久老熟女| 少妇的丰满在线观看| 在线国产一区二区在线| 久久欧美精品欧美久久欧美| 久久香蕉国产精品| 国产成人啪精品午夜网站| 岛国在线观看网站| 日韩 欧美 亚洲 中文字幕| 黄色视频,在线免费观看| 国产精品秋霞免费鲁丝片| 国产99久久九九免费精品| 亚洲精品中文字幕在线视频| 国产成人精品在线电影| 国产麻豆成人av免费视频| 99精品在免费线老司机午夜| 亚洲最大成人中文| 在线免费观看的www视频| 在线观看午夜福利视频| 黑人欧美特级aaaaaa片| 精品电影一区二区在线| 中文字幕另类日韩欧美亚洲嫩草| 国产区一区二久久| 夜夜躁狠狠躁天天躁| 久久欧美精品欧美久久欧美| 日本a在线网址| 日本vs欧美在线观看视频| 咕卡用的链子| 亚洲伊人色综图| 国产一区二区激情短视频| 真人一进一出gif抽搐免费| 久久久久久久久免费视频了| 精品国产国语对白av| 欧美性长视频在线观看| 美女免费视频网站| 中文字幕精品免费在线观看视频| 欧美国产精品va在线观看不卡| 手机成人av网站| 精品国产乱子伦一区二区三区| 欧美成狂野欧美在线观看| 亚洲狠狠婷婷综合久久图片| 高清黄色对白视频在线免费看| 国产精品一区二区三区四区久久 | 一二三四在线观看免费中文在| 成人av一区二区三区在线看| 女同久久另类99精品国产91| 曰老女人黄片| 一本大道久久a久久精品| 色综合婷婷激情| 国产精品亚洲一级av第二区| 国产欧美日韩精品亚洲av| 午夜精品在线福利| 中文字幕久久专区| 91成人精品电影| 99国产综合亚洲精品| 国产97色在线日韩免费| 亚洲专区字幕在线| 日韩欧美三级三区| 一边摸一边抽搐一进一小说| 日日夜夜操网爽| 国产xxxxx性猛交| 男女下面插进去视频免费观看| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区激情视频| 男女之事视频高清在线观看| 久久国产精品人妻蜜桃| 天天躁夜夜躁狠狠躁躁| 亚洲国产欧美一区二区综合| 91大片在线观看| 亚洲欧美日韩另类电影网站| 国产私拍福利视频在线观看| 国产精品亚洲av一区麻豆| 天天一区二区日本电影三级 | 久久午夜综合久久蜜桃| 在线观看www视频免费| 亚洲成av片中文字幕在线观看| 久9热在线精品视频| 9191精品国产免费久久| av欧美777| 麻豆av在线久日| 欧美一级a爱片免费观看看 | 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲av高清不卡| 中文字幕人妻丝袜一区二区| 69精品国产乱码久久久| 亚洲 欧美一区二区三区| 欧美在线黄色| 制服诱惑二区| 中文字幕精品免费在线观看视频| 黄色视频不卡| 波多野结衣一区麻豆| 国产免费男女视频| 免费高清在线观看日韩| 亚洲色图综合在线观看| 亚洲第一青青草原| 亚洲欧美精品综合久久99| 女人被躁到高潮嗷嗷叫费观| 亚洲av熟女| 国内精品久久久久久久电影| 日本黄色视频三级网站网址| 亚洲人成77777在线视频| 午夜老司机福利片| 动漫黄色视频在线观看| 亚洲狠狠婷婷综合久久图片| 99久久精品国产亚洲精品| 午夜两性在线视频| 国产精品自产拍在线观看55亚洲| www.精华液| 国产av在哪里看| 欧美成人性av电影在线观看| 黄片播放在线免费| 嫩草影视91久久| 日韩欧美免费精品| 色尼玛亚洲综合影院| 老汉色∧v一级毛片| 亚洲国产精品成人综合色| 国产片内射在线| 久久久久久免费高清国产稀缺| 91麻豆精品激情在线观看国产| 性欧美人与动物交配| 老司机在亚洲福利影院| 啦啦啦 在线观看视频| 国产色视频综合| 热re99久久国产66热| 亚洲第一av免费看| 亚洲五月天丁香| 老司机午夜福利在线观看视频| 欧美在线一区亚洲| 91精品三级在线观看| 日韩av在线大香蕉| 久久亚洲真实| 国产成人av教育| 亚洲在线自拍视频| 老司机靠b影院| 丰满的人妻完整版| 一夜夜www| 精品一区二区三区四区五区乱码| 成人18禁高潮啪啪吃奶动态图| 色婷婷久久久亚洲欧美| 天天躁夜夜躁狠狠躁躁| 午夜福利,免费看| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美日韩高清在线视频| av欧美777| 久9热在线精品视频| 中文字幕精品免费在线观看视频| 啦啦啦免费观看视频1| 亚洲人成伊人成综合网2020| 老司机在亚洲福利影院| 国产一区二区激情短视频| 久久人妻av系列| 99国产极品粉嫩在线观看| 老鸭窝网址在线观看| 免费久久久久久久精品成人欧美视频| 乱人伦中国视频| 色播亚洲综合网| 亚洲在线自拍视频| 精品国产一区二区三区四区第35| aaaaa片日本免费| 激情在线观看视频在线高清| 啦啦啦 在线观看视频| 无人区码免费观看不卡| 久久 成人 亚洲| 99久久99久久久精品蜜桃| 别揉我奶头~嗯~啊~动态视频| 亚洲免费av在线视频| 亚洲第一电影网av| 99国产精品一区二区三区| 国产欧美日韩一区二区三区在线| 男女床上黄色一级片免费看| 黄频高清免费视频| 三级毛片av免费| 高潮久久久久久久久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 波多野结衣高清无吗| 国产私拍福利视频在线观看| 欧美在线一区亚洲| 女人被狂操c到高潮| 国产成人免费无遮挡视频| 757午夜福利合集在线观看| 国产精品一区二区三区四区久久 | 国产一区二区三区在线臀色熟女| 亚洲欧洲精品一区二区精品久久久| 黑人巨大精品欧美一区二区mp4| 成人亚洲精品一区在线观看| 国产99久久九九免费精品| 99re在线观看精品视频| 国产成人一区二区三区免费视频网站| 99久久久亚洲精品蜜臀av| 国产高清有码在线观看视频 | 涩涩av久久男人的天堂| 18禁黄网站禁片午夜丰满| 国产精品99久久99久久久不卡| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| 亚洲国产精品999在线| 久久午夜亚洲精品久久| 超碰成人久久| 叶爱在线成人免费视频播放| 久久青草综合色| 免费少妇av软件| 亚洲一区高清亚洲精品| 国产精品1区2区在线观看.| 91成年电影在线观看| 欧美色欧美亚洲另类二区 | 高清黄色对白视频在线免费看| 最新美女视频免费是黄的| 12—13女人毛片做爰片一| 国产成人精品久久二区二区91| 欧美中文综合在线视频| 日韩 欧美 亚洲 中文字幕| 欧美日韩福利视频一区二区| 国产亚洲精品久久久久久毛片| 久久亚洲真实| 国产伦一二天堂av在线观看| 有码 亚洲区| 联通29元200g的流量卡| 精品一区二区免费观看| 女生性感内裤真人,穿戴方法视频| 欧美最黄视频在线播放免费| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 亚洲精品一区av在线观看| 欧美成人一区二区免费高清观看| 亚洲最大成人中文| 国产精品女同一区二区软件 | 久久精品夜夜夜夜夜久久蜜豆| 欧美另类亚洲清纯唯美| 国产av在哪里看| 少妇猛男粗大的猛烈进出视频 | 天堂动漫精品| 亚洲美女搞黄在线观看 | 麻豆久久精品国产亚洲av| 午夜福利视频1000在线观看| 国产精品1区2区在线观看.| 露出奶头的视频| 午夜日韩欧美国产| 欧美黑人欧美精品刺激| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 欧美中文日本在线观看视频| 国产高潮美女av| 国产成人aa在线观看| 欧美最新免费一区二区三区| 亚洲av一区综合| 国产av麻豆久久久久久久| 亚洲av.av天堂| 欧美色视频一区免费| 一个人看的www免费观看视频| a级毛片免费高清观看在线播放| 成人亚洲精品av一区二区| 在线a可以看的网站| 十八禁国产超污无遮挡网站| 精品午夜福利视频在线观看一区| 亚洲成av人片在线播放无| 黄片wwwwww| 免费人成视频x8x8入口观看| 最新在线观看一区二区三区| 搡老熟女国产l中国老女人| 国产 一区精品| 99久久无色码亚洲精品果冻| 尤物成人国产欧美一区二区三区| 成人特级黄色片久久久久久久| 色在线成人网| 欧美一区二区亚洲| 真实男女啪啪啪动态图| 一进一出抽搐gif免费好疼| 国产精品日韩av在线免费观看| 1024手机看黄色片| 亚洲成av人片在线播放无| 国产精品无大码| 久久天躁狠狠躁夜夜2o2o| 成人一区二区视频在线观看| 淫秽高清视频在线观看| 成人一区二区视频在线观看| 亚洲一级一片aⅴ在线观看| 啦啦啦观看免费观看视频高清| 精品一区二区免费观看| 亚洲在线自拍视频| 亚洲自偷自拍三级| 亚洲图色成人| 免费不卡的大黄色大毛片视频在线观看 | 性欧美人与动物交配| 又粗又爽又猛毛片免费看| 又粗又爽又猛毛片免费看| 欧美精品啪啪一区二区三区| 最新中文字幕久久久久| 最后的刺客免费高清国语| 国产色婷婷99| 97碰自拍视频| a级毛片a级免费在线| 欧美丝袜亚洲另类 | 一级黄色大片毛片| 国产精品精品国产色婷婷| 伊人久久精品亚洲午夜| 偷拍熟女少妇极品色| 99热这里只有是精品50| 九色成人免费人妻av| 国产 一区 欧美 日韩| 搡女人真爽免费视频火全软件 | 麻豆国产av国片精品| 欧美人与善性xxx| 丝袜美腿在线中文| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品亚洲美女久久久| 露出奶头的视频| 啦啦啦观看免费观看视频高清| 亚洲最大成人av| 真人做人爱边吃奶动态| 在线免费观看的www视频| 国产精品野战在线观看| 国产高清视频在线观看网站| 久久久久久久久久久丰满 | 18禁黄网站禁片午夜丰满| 精品一区二区三区人妻视频| 一区二区三区免费毛片| 12—13女人毛片做爰片一| 亚洲av免费高清在线观看| 一级黄色大片毛片| 国产精品美女特级片免费视频播放器| 一进一出抽搐gif免费好疼| 久久精品国产亚洲网站| 2021天堂中文幕一二区在线观| 一级黄片播放器| 亚洲成人中文字幕在线播放| 午夜福利成人在线免费观看| 搡老岳熟女国产| 国产大屁股一区二区在线视频| 99久久精品热视频| 日韩欧美在线二视频| 亚洲av成人av| 啦啦啦韩国在线观看视频| 欧美性猛交╳xxx乱大交人| 国产亚洲精品av在线| or卡值多少钱| 色综合色国产| 欧美在线一区亚洲| 久久久久国产精品人妻aⅴ院| 性色avwww在线观看| 日本成人三级电影网站| 亚洲五月天丁香| 一级黄色大片毛片| 天天躁日日操中文字幕| 国产女主播在线喷水免费视频网站 | 日韩大尺度精品在线看网址| 中文字幕熟女人妻在线| 亚洲av二区三区四区| 精品欧美国产一区二区三| 日韩欧美国产一区二区入口| 蜜桃久久精品国产亚洲av| 国产男人的电影天堂91| 欧美激情国产日韩精品一区| 一本久久中文字幕| 日本五十路高清| 亚洲国产精品久久男人天堂| 麻豆国产97在线/欧美| 欧美潮喷喷水| 国产伦人伦偷精品视频| 国产精品久久久久久久电影| 欧美xxxx性猛交bbbb| av国产免费在线观看| 又黄又爽又刺激的免费视频.| 久久久精品欧美日韩精品| 色综合色国产| 午夜爱爱视频在线播放| 美女大奶头视频| 日韩欧美国产在线观看| 久久久久久久亚洲中文字幕| 午夜福利18| 午夜精品在线福利| 亚洲国产精品sss在线观看| 精品无人区乱码1区二区| 搡女人真爽免费视频火全软件 | 九色国产91popny在线| 亚洲电影在线观看av| 99国产极品粉嫩在线观看| 国产欧美日韩一区二区精品| 亚洲成人免费电影在线观看| 性色avwww在线观看| 精品福利观看| 岛国在线免费视频观看| 99久久九九国产精品国产免费| 亚洲av第一区精品v没综合| 伊人久久精品亚洲午夜| 狠狠狠狠99中文字幕| 国产毛片a区久久久久| 夜夜爽天天搞| 欧美最新免费一区二区三区| 99riav亚洲国产免费| av福利片在线观看| www.www免费av| 免费观看在线日韩| 国产淫片久久久久久久久| av在线蜜桃| 国产精品美女特级片免费视频播放器| 亚洲精品亚洲一区二区| 欧美色欧美亚洲另类二区| 亚洲自拍偷在线| 极品教师在线视频| 久久亚洲精品不卡| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级国产av玫瑰| 亚洲在线自拍视频| 97超视频在线观看视频| 欧美高清成人免费视频www| 欧美成人一区二区免费高清观看| 精品人妻1区二区| 中文字幕人妻熟人妻熟丝袜美| 俄罗斯特黄特色一大片| 永久网站在线| 91精品国产九色| 国产高清视频在线观看网站| 亚洲国产精品合色在线| 一进一出抽搐动态| 日日撸夜夜添| 搞女人的毛片| 成人av在线播放网站| 热99在线观看视频| 日韩在线高清观看一区二区三区 | 国产成人a区在线观看| 69人妻影院| 国产爱豆传媒在线观看| 国产精品永久免费网站| 欧美性感艳星| 欧美三级亚洲精品| 亚洲av不卡在线观看| 成年女人永久免费观看视频| 变态另类成人亚洲欧美熟女| 日本黄色片子视频| 天堂√8在线中文| 网址你懂的国产日韩在线| x7x7x7水蜜桃| 亚洲精品色激情综合| 黄色配什么色好看| 三级毛片av免费| 99久久中文字幕三级久久日本| 欧美激情国产日韩精品一区| 国产真实乱freesex| 99久久精品热视频| 成人永久免费在线观看视频| av女优亚洲男人天堂| 久久精品国产亚洲av香蕉五月| 中国美白少妇内射xxxbb| 真人一进一出gif抽搐免费| 久久精品夜夜夜夜夜久久蜜豆| 亚洲午夜理论影院| 日本欧美国产在线视频| 日韩精品有码人妻一区| 亚洲国产日韩欧美精品在线观看| 久久久色成人| 男人舔女人下体高潮全视频| 国产亚洲精品综合一区在线观看| 亚洲最大成人手机在线| 黄片wwwwww| 日韩在线高清观看一区二区三区 | 看十八女毛片水多多多| 99热这里只有是精品在线观看| 国产亚洲91精品色在线| 少妇高潮的动态图| 97超级碰碰碰精品色视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 欧美日韩国产亚洲二区| 国产淫片久久久久久久久| 丰满的人妻完整版| 日韩国内少妇激情av| 热99re8久久精品国产| 亚洲av免费在线观看|