• <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)地正外部性價值量估算
    專利
    91久久精品国产一区二区成人| 久久午夜福利片| 欧美在线一区亚洲| 久久精品国产亚洲网站| 国产色婷婷99| 一级黄色大片毛片| 久久中文看片网| 亚洲最大成人av| 99久国产av精品国产电影| 欧美成人免费av一区二区三区| 99热精品在线国产| 精品熟女少妇av免费看| 亚洲av成人av| 亚洲欧美日韩高清专用| 国产91av在线免费观看| 亚洲经典国产精华液单| 黄片wwwwww| 观看美女的网站| 久久久久性生活片| 国产精品亚洲美女久久久| 中文字幕av成人在线电影| 长腿黑丝高跟| 国产精品美女特级片免费视频播放器| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久久毛片| 国产精品福利在线免费观看| 日韩欧美免费精品| 少妇人妻精品综合一区二区 | 国产在视频线在精品| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| 国产精品久久视频播放| 国产在视频线在精品| 亚洲专区国产一区二区| 国产精品免费一区二区三区在线| 18禁在线无遮挡免费观看视频 | 少妇丰满av| 欧美日韩乱码在线| 特大巨黑吊av在线直播| 黄色日韩在线| 国产精品免费一区二区三区在线| 99久国产av精品国产电影| 国产高清三级在线| 国产高潮美女av| 波多野结衣高清作品| 男人的好看免费观看在线视频| 成熟少妇高潮喷水视频| 亚洲av二区三区四区| 男人舔女人下体高潮全视频| 黄片wwwwww| 日韩高清综合在线| 综合色丁香网| 搡女人真爽免费视频火全软件 | 国产精品不卡视频一区二区| 日本三级黄在线观看| 在线免费观看不下载黄p国产| 在线观看午夜福利视频| 一本一本综合久久| 久久草成人影院| 精品一区二区三区av网在线观看| 小蜜桃在线观看免费完整版高清| 99热精品在线国产| 精品一区二区三区视频在线| 欧美成人免费av一区二区三区| 久久99热6这里只有精品| 男女视频在线观看网站免费| 国产淫片久久久久久久久| 九九爱精品视频在线观看| 欧美潮喷喷水| 成人美女网站在线观看视频| 赤兔流量卡办理| 国内精品一区二区在线观看| 一个人看的www免费观看视频| 极品教师在线视频| 亚洲人成网站在线观看播放| 美女大奶头视频| 久久久午夜欧美精品| 亚洲自偷自拍三级| 男人和女人高潮做爰伦理| 淫妇啪啪啪对白视频| 永久网站在线| 久久久久久久午夜电影| 久久草成人影院| 日韩国内少妇激情av| 亚洲自拍偷在线| 综合色丁香网| 在线播放国产精品三级| 午夜福利在线观看免费完整高清在 | 日韩精品中文字幕看吧| 国产一区二区在线观看日韩| 久久久久九九精品影院| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久久电影| 国产黄a三级三级三级人| 国产精品美女特级片免费视频播放器| 尾随美女入室| 草草在线视频免费看| 一级毛片我不卡| 国产 一区 欧美 日韩| 国产午夜福利久久久久久| 亚洲精品色激情综合| 亚洲av免费高清在线观看| 午夜精品一区二区三区免费看| 中文亚洲av片在线观看爽| 日韩欧美三级三区| 精品无人区乱码1区二区| 亚洲欧美成人精品一区二区| 精品国产三级普通话版| 亚洲欧美日韩高清专用| 久久久久免费精品人妻一区二区| 在线观看一区二区三区| 国内精品宾馆在线| 国产精品电影一区二区三区| 国产真实伦视频高清在线观看| 精品免费久久久久久久清纯| 观看免费一级毛片| 国产精品福利在线免费观看| 丝袜美腿在线中文| 亚洲欧美日韩高清在线视频| 亚洲乱码一区二区免费版| 最近中文字幕高清免费大全6| 久久精品国产99精品国产亚洲性色| 亚洲无线观看免费| 免费av毛片视频| 春色校园在线视频观看| 亚洲欧美成人精品一区二区| 久久久久久伊人网av| 一本久久中文字幕| 国产精品99久久久久久久久| 午夜a级毛片| 在线观看午夜福利视频| 国产亚洲精品av在线| 美女大奶头视频| 中国美女看黄片| 亚洲人成网站高清观看| 午夜福利在线在线| 18禁在线无遮挡免费观看视频 | 99在线人妻在线中文字幕| 一进一出抽搐动态| 欧美最新免费一区二区三区| 中文字幕熟女人妻在线| 91精品国产九色| 一个人看的www免费观看视频| 一进一出好大好爽视频| 最近2019中文字幕mv第一页| 熟妇人妻久久中文字幕3abv| 国产精品国产三级国产av玫瑰| 日本黄大片高清| 成人国产麻豆网| 久久久久国内视频| 国产色婷婷99| 波多野结衣高清无吗| av在线播放精品| 久久精品国产鲁丝片午夜精品| av天堂中文字幕网| 亚洲最大成人av| 欧美极品一区二区三区四区| 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 老司机福利观看| 欧美三级亚洲精品| 日韩欧美免费精品| 国产色婷婷99| 日韩精品有码人妻一区| 国产黄色视频一区二区在线观看 | 日本三级黄在线观看| 美女被艹到高潮喷水动态| 一个人看的www免费观看视频| 日本黄色视频三级网站网址| 国产 一区精品| 九九热线精品视视频播放| 毛片女人毛片| 中文在线观看免费www的网站| 国产伦一二天堂av在线观看| 一区二区三区高清视频在线| 亚洲国产欧洲综合997久久,| 欧美成人一区二区免费高清观看| 国产精品一二三区在线看| av专区在线播放| 国产视频一区二区在线看| 精品国产三级普通话版| 日日干狠狠操夜夜爽| 秋霞在线观看毛片| 免费人成在线观看视频色| 国产老妇女一区| 国产免费一级a男人的天堂| 久久精品人妻少妇| 超碰av人人做人人爽久久| 日本在线视频免费播放| av中文乱码字幕在线| 毛片一级片免费看久久久久| 中文字幕av在线有码专区| 日韩中字成人| 人妻少妇偷人精品九色| 欧美最新免费一区二区三区| 久久综合国产亚洲精品| 99热这里只有是精品50| 国内揄拍国产精品人妻在线| 99热这里只有是精品在线观看| 少妇熟女aⅴ在线视频| 国产精品国产三级国产av玫瑰| av女优亚洲男人天堂| 国产免费男女视频| 变态另类丝袜制服| 久久久久久久久久久丰满| 天天躁夜夜躁狠狠久久av| 岛国在线免费视频观看| 午夜视频国产福利| 久久中文看片网| 啦啦啦啦在线视频资源| 国内少妇人妻偷人精品xxx网站| 欧美精品国产亚洲| 成年女人毛片免费观看观看9| 看十八女毛片水多多多| 日韩欧美免费精品| 精品久久久久久久人妻蜜臀av| 亚洲成人中文字幕在线播放| 亚洲国产精品成人久久小说 | 亚洲精品一区av在线观看| 欧美zozozo另类| 人妻久久中文字幕网| 午夜视频国产福利| 精品熟女少妇av免费看| 亚洲美女视频黄频| 99在线人妻在线中文字幕| 麻豆成人午夜福利视频| 一级黄色大片毛片| 岛国在线免费视频观看| 久久久久免费精品人妻一区二区| 免费av毛片视频| 麻豆精品久久久久久蜜桃| 在线a可以看的网站| 亚洲精品国产成人久久av| 如何舔出高潮| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 中文资源天堂在线| 99在线人妻在线中文字幕| 成年av动漫网址| 日韩一本色道免费dvd| 毛片一级片免费看久久久久| av在线老鸭窝| 国产精品人妻久久久久久| 欧美在线一区亚洲| 精品人妻熟女av久视频| 免费观看精品视频网站| 亚洲精品一区av在线观看| 亚洲五月天丁香| 国产乱人偷精品视频| 老熟妇乱子伦视频在线观看| 久久国产乱子免费精品| 国产精品久久久久久av不卡| 国产色婷婷99| 人人妻人人澡欧美一区二区| 久久久久久国产a免费观看| 嫩草影院入口| 午夜精品国产一区二区电影 | 丝袜喷水一区| 日韩在线高清观看一区二区三区| 热99在线观看视频| 干丝袜人妻中文字幕| 国产极品精品免费视频能看的| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 日日啪夜夜撸| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| av在线天堂中文字幕| 大又大粗又爽又黄少妇毛片口| 一进一出抽搐gif免费好疼| 色哟哟·www| 综合色av麻豆| 国产精品嫩草影院av在线观看| 亚洲天堂国产精品一区在线| 国产免费一级a男人的天堂| 国产精品女同一区二区软件| 亚洲不卡免费看| 不卡一级毛片| 综合色av麻豆| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 天堂√8在线中文| 成人鲁丝片一二三区免费| 精品福利观看| 久久久久久久久大av| 一级黄片播放器| 国产免费男女视频| 国产私拍福利视频在线观看| 国产人妻一区二区三区在| 亚洲欧美日韩东京热| 十八禁网站免费在线| 久久久国产成人精品二区| 免费观看在线日韩| 久久久久久国产a免费观看| 高清毛片免费看| 久久热精品热| av卡一久久| 精品日产1卡2卡| 国产视频一区二区在线看| 如何舔出高潮| 真人做人爱边吃奶动态| 亚洲成人精品中文字幕电影| 嫩草影院入口| 青春草视频在线免费观看| 3wmmmm亚洲av在线观看| 午夜福利视频1000在线观看| 亚洲精华国产精华液的使用体验 | 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 综合色丁香网| 国产aⅴ精品一区二区三区波| 亚洲第一区二区三区不卡| 国产精品久久久久久精品电影| 日本精品一区二区三区蜜桃| 国产真实乱freesex| 国产成人一区二区在线| 国产精品久久久久久亚洲av鲁大| 亚洲最大成人av| 极品教师在线视频| 国产成人精品久久久久久| 亚洲国产欧洲综合997久久,| 1024手机看黄色片| 99久久九九国产精品国产免费| 无遮挡黄片免费观看| 精品久久久久久久末码| 成人综合一区亚洲| 村上凉子中文字幕在线| 午夜精品在线福利| 国产精品乱码一区二三区的特点| 男女下面进入的视频免费午夜| 禁无遮挡网站| 日韩成人伦理影院| 草草在线视频免费看| 国产成人aa在线观看| 白带黄色成豆腐渣| 国产男人的电影天堂91| 久99久视频精品免费| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| 欧美激情国产日韩精品一区| 亚洲av美国av| 久久精品国产亚洲网站| 婷婷亚洲欧美| 国产高清视频在线观看网站| 乱系列少妇在线播放| 成年版毛片免费区| 久久久精品欧美日韩精品| 午夜精品一区二区三区免费看| 免费观看人在逋| 午夜福利成人在线免费观看| av天堂中文字幕网| 日韩强制内射视频| 日本在线视频免费播放| 蜜臀久久99精品久久宅男| 91久久精品国产一区二区成人| 午夜福利成人在线免费观看| 悠悠久久av| 综合色av麻豆| 成熟少妇高潮喷水视频| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 国产aⅴ精品一区二区三区波| 国产女主播在线喷水免费视频网站 | 12—13女人毛片做爰片一| 综合色av麻豆| 久久久久精品国产欧美久久久| 国产毛片a区久久久久| 国产一区亚洲一区在线观看| 国产精品一二三区在线看| 国产精品1区2区在线观看.| 十八禁网站免费在线| 晚上一个人看的免费电影| 日本免费一区二区三区高清不卡| 99九九线精品视频在线观看视频| 日韩大尺度精品在线看网址| 国产亚洲91精品色在线| 成人精品一区二区免费| 中文字幕熟女人妻在线| 精品少妇黑人巨大在线播放 | 日韩中字成人| 99热这里只有精品一区| 成人国产麻豆网| 免费在线观看成人毛片| 午夜福利18| 最新在线观看一区二区三区| 午夜影院日韩av| 桃色一区二区三区在线观看| www日本黄色视频网| 国产视频内射| 五月伊人婷婷丁香| 亚洲国产精品国产精品| 免费电影在线观看免费观看| 成人美女网站在线观看视频| 色综合站精品国产| 亚洲av中文字字幕乱码综合| 欧美一区二区亚洲| 国产精品乱码一区二三区的特点| 搡老妇女老女人老熟妇| 少妇高潮的动态图| 中国美女看黄片| 你懂的网址亚洲精品在线观看 | 乱人视频在线观看| 欧美+日韩+精品| 欧美性感艳星| 久久鲁丝午夜福利片| 精品免费久久久久久久清纯| 亚洲欧美日韩东京热| av视频在线观看入口| 少妇高潮的动态图| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 97在线视频观看| 国产精品不卡视频一区二区| 男女之事视频高清在线观看| 国产精品一区二区三区四区免费观看 | aaaaa片日本免费| 国内精品美女久久久久久| 久99久视频精品免费| 欧美一区二区精品小视频在线| 最后的刺客免费高清国语| 亚洲中文字幕日韩| 久久久色成人| 婷婷精品国产亚洲av在线| 久久久久性生活片| 欧美日韩精品成人综合77777| 成人二区视频| 精品国产三级普通话版| 丝袜喷水一区| 国产一区二区在线av高清观看| 国产精品一区二区免费欧美| АⅤ资源中文在线天堂| 国产激情偷乱视频一区二区| 国产精品女同一区二区软件| 啦啦啦观看免费观看视频高清| 91av网一区二区| 亚洲国产精品成人久久小说 | 永久网站在线| 在线免费观看不下载黄p国产| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| av国产免费在线观看| 亚洲国产精品国产精品| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 观看美女的网站| 人妻制服诱惑在线中文字幕| 一本久久中文字幕| 男人舔奶头视频| 欧美丝袜亚洲另类| 少妇人妻一区二区三区视频| 中国美白少妇内射xxxbb| 我要搜黄色片| 久久久久国产精品人妻aⅴ院| 少妇猛男粗大的猛烈进出视频 | 精品一区二区免费观看| 久久久精品大字幕| 女生性感内裤真人,穿戴方法视频| 亚洲三级黄色毛片| 天堂√8在线中文| 亚洲成av人片在线播放无| 成人午夜高清在线视频| 在线观看66精品国产| 美女cb高潮喷水在线观看| 亚洲中文字幕日韩| 99在线视频只有这里精品首页| 天天躁夜夜躁狠狠久久av| 成人二区视频| 国产蜜桃级精品一区二区三区| 国产大屁股一区二区在线视频| 亚洲av免费在线观看| 中文字幕av在线有码专区| av在线播放精品| 夜夜看夜夜爽夜夜摸| 国产精品一区www在线观看| 亚洲精品一区av在线观看| 精品久久久久久成人av| 亚州av有码| 麻豆一二三区av精品| 亚洲av第一区精品v没综合| 国产老妇女一区| 国产精品美女特级片免费视频播放器| 蜜臀久久99精品久久宅男| 亚洲久久久久久中文字幕| 国产精品人妻久久久影院| 91在线观看av| 卡戴珊不雅视频在线播放| 身体一侧抽搐| 欧美丝袜亚洲另类| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看| 亚洲熟妇中文字幕五十中出| 搡老熟女国产l中国老女人| 亚洲婷婷狠狠爱综合网| 寂寞人妻少妇视频99o| 最近2019中文字幕mv第一页| 婷婷色综合大香蕉| 亚洲人成网站高清观看| 色综合站精品国产| 黄色一级大片看看| 日韩高清综合在线| 日韩精品青青久久久久久| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久精品电影| 亚洲欧美日韩高清在线视频| 18禁在线播放成人免费| 日本成人三级电影网站| 午夜爱爱视频在线播放| 少妇人妻一区二区三区视频| 日本三级黄在线观看| 国产精品三级大全| 精品少妇黑人巨大在线播放 | 91在线精品国自产拍蜜月| 亚洲电影在线观看av| 少妇人妻精品综合一区二区 | 天堂动漫精品| 国产精品亚洲美女久久久| 国产成人影院久久av| 性色avwww在线观看| 日韩人妻高清精品专区| 老司机午夜福利在线观看视频| 黄片wwwwww| 成人综合一区亚洲| 一夜夜www| 日本-黄色视频高清免费观看| 午夜激情福利司机影院| 婷婷六月久久综合丁香| 国产伦精品一区二区三区四那| 日韩成人av中文字幕在线观看 | 最近的中文字幕免费完整| av国产免费在线观看| 成人综合一区亚洲| 狠狠狠狠99中文字幕| 在线观看66精品国产| 久久人人爽人人片av| 九九久久精品国产亚洲av麻豆| 精品国内亚洲2022精品成人| 俄罗斯特黄特色一大片| 欧美激情久久久久久爽电影| 亚洲av免费高清在线观看| 亚洲国产精品sss在线观看| 精品久久久久久久久久免费视频| 欧美性猛交黑人性爽| av女优亚洲男人天堂| 国产精品福利在线免费观看| 真实男女啪啪啪动态图| 亚洲av一区综合| 精品熟女少妇av免费看| 成人二区视频| 欧美精品国产亚洲| 久久久久国产网址| 国产精品嫩草影院av在线观看| av中文乱码字幕在线| 国产av一区在线观看免费| 嫩草影院精品99| 高清午夜精品一区二区三区 | aaaaa片日本免费| 精品99又大又爽又粗少妇毛片| 熟女人妻精品中文字幕| 永久网站在线| 久久久国产成人免费| 12—13女人毛片做爰片一| 简卡轻食公司| 日日啪夜夜撸| 亚洲欧美精品综合久久99| 日韩欧美一区二区三区在线观看| 九九爱精品视频在线观看| 亚洲成人精品中文字幕电影| 校园春色视频在线观看| 一进一出好大好爽视频| 又黄又爽又刺激的免费视频.| 丰满人妻一区二区三区视频av| 观看免费一级毛片| 麻豆国产97在线/欧美| 中文资源天堂在线| 欧美zozozo另类| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 又粗又爽又猛毛片免费看| 欧美中文日本在线观看视频| 少妇丰满av| 国产色爽女视频免费观看| 久久精品综合一区二区三区| 欧美人与善性xxx| 精品免费久久久久久久清纯| 波野结衣二区三区在线| 久久久久九九精品影院| 日产精品乱码卡一卡2卡三| 99热这里只有是精品50| 日日摸夜夜添夜夜爱| 大型黄色视频在线免费观看| 亚洲国产精品国产精品| 深爱激情五月婷婷| 成人av一区二区三区在线看| 别揉我奶头 嗯啊视频| 久久综合国产亚洲精品| 国产高清视频在线播放一区| 高清毛片免费看| 国产精品女同一区二区软件| 国产高潮美女av| 尤物成人国产欧美一区二区三区| 听说在线观看完整版免费高清| 中文亚洲av片在线观看爽| 99久国产av精品| 男女边吃奶边做爰视频| 一级黄色大片毛片| 亚洲精品456在线播放app| 久久久久久久午夜电影| av.在线天堂| 女的被弄到高潮叫床怎么办| 色综合站精品国产| 国产午夜福利久久久久久| 亚洲精品乱码久久久v下载方式| 九色成人免费人妻av| 在线观看免费视频日本深夜|