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

    京津冀PM2.5的主要影響因素及內(nèi)在關系研究

    2017-04-26 07:36:55周曙東??歐陽緯清葛繼紅
    中國人口·資源與環(huán)境 2017年4期
    關鍵詞:灰霾排放量大氣

    周曙東??歐陽緯清 葛繼紅

    摘要大氣污染物的源排放是形成灰霾天氣的內(nèi)因,氣象條件是形成灰霾天氣的外因。本研究通過構建PM2.5濃度的兩段式分布滯后模型,結(jié)合自然環(huán)境因素及經(jīng)濟因素對PM2.5的影響因素進行了綜合分析。在第一段模型中構建了PM2.5和大氣污染物排放量的分布滯后模型,第二段模型中構建了不同的大氣污染源對大氣污染物排放量的影響因素模型。大氣污染物排放源主要包括工業(yè)源、生活源、機動車源、集中式污染治理設施源。在工業(yè)源中,工業(yè)廢氣重度污染行業(yè)是大氣污染物排放主要的貢獻者;在生活源中,燃煤消費量對大氣污染物排放影響很大,這也是冬季供暖期間PM2.5劇增的原因;在機動車源中,盡管黃標車的保有量僅占汽車保有量的10%左右,但卻占據(jù)了顆粒物排放量的絕大部分。利用京津冀代表性城市PM2.5日度數(shù)據(jù)研究得出平均氣溫、平均風速、日照時數(shù)、平均氣壓、降雨量、平均相對濕度、沙塵暴等因素對PM2.5濃度的負向與正向作用。研究發(fā)現(xiàn),大氣污染物排放量對PM2.5濃度具有聚集的滯后效應,當期大氣污染物排放量、滯后一期、滯后兩期、滯后三期大氣污染物對PM2.5濃度具有顯著的正向作用,且影響依次遞減。構建的大氣污染物排放量的污染源影響因素模型揭示一個地區(qū)煤炭消費量、工業(yè)廢氣重度污染行業(yè)工業(yè)增加值、黃標車保有量對該地區(qū)大氣污染物排放量具有顯著影響。本研究對優(yōu)化能源消費結(jié)構和產(chǎn)業(yè)結(jié)構,減少空氣污染物排放提出了對策建議。

    關鍵詞PM2.5;大氣污染物排放;兩段式分布滯后模型;京津冀

    中圖分類號X513;F421文獻標識碼A文章編號1002-2104(2017)04-0102-08doi:10.12062/cpre.20170343

    PM2.5的空氣動力學含義是指直接小于等于2.5μm的顆粒物,也稱為細顆粒物。PM2.5能較長時間懸浮于空氣中,對空氣質(zhì)量和能見度等有重要的影響?;姻蔡鞖獾谋举|(zhì)是大氣顆粒物PM2.5的污染[1]。與PM10相比,PM2.5粒徑小,擴散面積大,活性強,易附帶有毒、有害物質(zhì),且在大氣中停留時間長,因而對人體健康和大氣環(huán)境質(zhì)量影響更大,灰霾污染中的PM2.5容易誘發(fā)呼吸系統(tǒng)疾病。2013年,灰霾成為年度關鍵詞,2013年1月,4次灰霾籠罩全國30個省。京津冀過去冬春季的灰霾比較嚴重,而最近幾年四季都會發(fā)生灰霾天氣。因此,研究PM2.5的主要影響因素,理順其內(nèi)在關系具有很迫切的現(xiàn)實意義,可以為相關研究和政策制定提供理論依據(jù)。

    京津冀都市圈是我國三大城市群之一,在經(jīng)濟發(fā)展中占據(jù)重要地位,而在空氣質(zhì)量最差的10個城市中,7個城市在京津冀地區(qū),京津冀已經(jīng)成為我國空氣污染最嚴重的地區(qū)之一。2015年4月30日,中共中央政治局審議通過的《京津冀協(xié)調(diào)發(fā)展規(guī)劃綱要》指出,要在京津冀生態(tài)環(huán)境保護、產(chǎn)業(yè)結(jié)構升級等重點領域率先取得突破。根據(jù)《2015中國環(huán)境狀況公報》顯示,京津冀PM2.5平均濃度為77 ug/m3(超過國家二級標準1.20倍),遠高于珠三角地區(qū)的34 ug/m3和長三角地區(qū)的53 ug/m3;京津冀區(qū)域13個城市達標天數(shù)比例在21.9%—86.4%,平均為42.8%,重度及以上污染天數(shù)比例為17%,而超標天數(shù)中以PM2.5為首要污染物天數(shù)最多。本研究將實證分析經(jīng)濟增長與大氣污染的關系,同時基于兩段式的分布滯后模型,結(jié)合自然環(huán)境因素和經(jīng)濟因素,對灰霾污染嚴重地區(qū)京津冀的PM2.5影響因素進行綜合分析。

    1文獻綜述

    對于可吸入顆粒物(PM)的最早研究,可追溯到20世紀80年代。1985年,美國將總懸浮顆粒物修改為PM10,隨后不久將PM2.5納入到大氣監(jiān)測范圍中[2]。2012年2月29日,中國環(huán)境保護部將環(huán)境空氣質(zhì)量標準(GB-2012)作為國家質(zhì)量標準,正式取代了環(huán)境空氣質(zhì)量標準(GB-1996),PM2.5開始納入到空氣質(zhì)量監(jiān)測中來。盡管PM2.5在我國監(jiān)測時間比較短,但由于PM2.5所引發(fā)的灰霾問題頻發(fā),引發(fā)人們的普遍關注。國內(nèi)外學者對PM2.5進行的廣泛深入的研究,主要集中在以下兩個方面:

    中國人口·資源與環(huán)境2017年第4期

    一是PM2.5的污染特征、排放清單、源解析等的物理化學過程、PM2.5與其他顆粒物直接的關系。梁延剛等利用2005年香港能見度與PM2.5濃度的每小時數(shù)據(jù),研究顯示PM2.5濃度與能見度在5%的顯著水平上呈負相關,相關系數(shù)為80%[3]。Chang W等對美國南加州的顆粒物進行源解析,將PM2.5分為機動車、道路揚塵、柴油車、海鹽等主要源類[4]。鄭玫等將PM2.5源種類分為塵類源、燃燒類源、二次元、其他源[5]。周靜、劉松華等通過對蘇州10個區(qū)15個工業(yè)行業(yè)現(xiàn)場咨詢、發(fā)放調(diào)查表等形式,研究得出工藝過程源和固定燃燒源分別占蘇州工業(yè)源PM2.5排放量的94%和6%[6]。朱增銀、郭清彬、于建華等研究PM2.5濃度與PM10濃度占比,分析得出PM2.5濃度與PM10濃度占比通常在29.10%—94.76%之間,呈現(xiàn)出季節(jié)差異,并且這一比例夏季要明顯低于秋季[7-9]。

    二是PM2.5與氣象因子及工業(yè)化、城市化、環(huán)境規(guī)制等因素之間的關系。張淑平等研究發(fā)現(xiàn)氣象條件與PM2.5日濃度顯著相關;劉瑞婷等研究表明高濃度的大氣顆粒物和特殊的氣象條件是灰霾形成的有利條件,小風或靜風容易造成PM2.5濃度增加,能見度降低;孟曉艷等研究得出較小的地面風速、較高的相對濕度以及穩(wěn)定的大氣層增加了PM2.5濃度的累積,導致空氣質(zhì)量惡化[10-12]。何楓等應用Tobit模型發(fā)現(xiàn)城市工業(yè)增加值占GDP比提高1個百分比,灰霾天數(shù)平均增加約4天 [13]。馬麗梅、張曉運用Donkelaar et al制作的首張PM2.5全球地圖,分析得出灰霾污染與GDP并不存在倒“U”型關系,灰霾污染與產(chǎn)業(yè)結(jié)構以及能源結(jié)構存在很強的正相關關系[14]。

    基于上述文獻回顧,我們可以發(fā)現(xiàn),現(xiàn)有研究主要有以下幾個方面不足:

    第一,灰霾的主要污染物PM2.5的形成是一個復合的大氣環(huán)境的機理形成過程,不僅與溫度、相對濕度、風速、壓強、降水等氣象因素影響密切,同時也受機動車、燃煤、工業(yè)生產(chǎn)等經(jīng)濟因素的影響,但絕大部分的研究要么只考慮自然環(huán)境因素,要么只考慮經(jīng)濟因素,即將燃煤、工業(yè)生產(chǎn)等經(jīng)濟因素直接與PM2.5進行回歸,很少有將自然環(huán)境因素和機動車、燃煤、工業(yè)生產(chǎn)等經(jīng)濟因素綜合起來進行研究。

    第二,由于我國將PM2.5納入空氣質(zhì)量檢測標準范圍時間較短,我國正式統(tǒng)計PM2.5開始于2013年,故經(jīng)濟學研究領域受限于數(shù)據(jù)的連續(xù)性和系統(tǒng)性,大部分研究采用PM10的年度數(shù)據(jù)或者模擬數(shù)據(jù)。PM2.5日度數(shù)據(jù)變化幅度很大,但年度平均濃度數(shù)據(jù)的變化幅度相對于日度數(shù)據(jù)并不是很大,因為平均數(shù)將劇烈變化拉平了,所以年度數(shù)據(jù)不能很好反映大氣污染物日度變化規(guī)律。另一些研究采用美國研究機構研制的根據(jù)氣象條件模擬出的PM2.5年度平均數(shù)據(jù),如馬麗梅等使用哥倫比亞大學國際科學信息網(wǎng)絡中心與巴特爾研究中心的全球2001—2010年PM2.5平均值的模擬數(shù)據(jù),但模擬數(shù)據(jù)與實際檢測數(shù)據(jù)存在誤差。

    本研究使用中國空氣質(zhì)量在線監(jiān)測分析平臺公布的京津冀2013年1月18日至2014年12月31日的日度數(shù)據(jù),結(jié)合影響PM2.5濃度的氣象因素和燃煤、機動車、工業(yè)產(chǎn)值等經(jīng)濟因素,對京津冀地區(qū)PM2.5的影響因素進行綜合分析。

    2模型構建

    2.1分析框架

    影響空氣質(zhì)量的主要因素有兩個:一個是污染源,還有一個是氣象條件。吳兌分析得出,大氣污染物的源排放是形成灰霾天氣的內(nèi)因,氣象條件是形成灰霾天氣的外因[1]。因此,一個城市的空氣質(zhì)量受自然條件和人類活動兩方面影響,自然條件包括風速、降雨、溫度等氣象條件以及地形地貌等地理條件和區(qū)域生態(tài)環(huán)境變化,人類活動包括城市工業(yè)發(fā)展、區(qū)域能源消耗。2014年《北京市環(huán)境狀況公報》對PM2.5進行源解析顯示,機動車源、燃煤、工業(yè)生產(chǎn)、揚塵對北京的PM2.5貢獻率分別達到了31.1%、22、4%、18.1%、14.3%。2014年《天津市環(huán)境狀況公報》對PM2.5進行源解析顯示,揚塵、燃煤、機動車、工業(yè)對天津的PM2.5貢獻率分別達到了30%、27%、20%、17%。

    我們對PM2.5濃度主要影響因素及其內(nèi)在關系的研究分為兩個階段,第一階段分析形成灰霾的外部因素,綜合考慮溫度、濕度、日照、風力、氣壓、降水等氣象因素。第二階段分析不同污染源對大氣污染物排放量的影響,主要包括四個方面:工業(yè)源、生活源、機動車源、集中式污染源。大氣污染物是造成PM2.5濃度上升的主要內(nèi)因。

    2.2模型設定

    通過以上對自然環(huán)境科學和環(huán)境經(jīng)濟學文獻綜述可以得出,PM2.5的影響因素主要有氣象因素,包括溫度、相對濕度、日照時數(shù)、氣壓、風力、降雨等氣象因子,以及煤炭燃燒產(chǎn)生的煙塵及一氧化碳、二氧化硫、氮氧化物及揮發(fā)性有機物等、工業(yè)源排放、機動車尾氣排放等因素。本文模型設定立足于自然環(huán)境科學基礎,并且將煤炭燃燒煙粉塵、機動車尾氣、工業(yè)源排放轉(zhuǎn)化為經(jīng)濟學上容易衡量的煤炭消費量、汽車保有量、工業(yè)增加值。一個地區(qū)污染物排放量的多少跟該地區(qū)的工業(yè)廢氣污染密集型產(chǎn)值有著密切的關系。本文借鑒張家峰對工業(yè)廢氣污染密集型行業(yè)的劃分方法[15],選取全國38個工業(yè)行業(yè)萬元產(chǎn)出的不同污染物排放量即排放強度(Emission Intensity,EI)大于等于0.75的五個重污染行業(yè),即選取電力和熱力的生產(chǎn)工業(yè)、非金屬礦物制品業(yè)、燃氣生產(chǎn)和供應業(yè)、黑色金屬冶煉及壓延加工、木材加工以及木竹藤棕草制品業(yè)五個工業(yè)廢氣重度污染行業(yè)。選取這五個重度污染行業(yè)工業(yè)增加值,研究工業(yè)廢氣重度污染行業(yè)工業(yè)增加值對PM2.5濃度的影響。本研究通過構建兩段式的分布滯后模型來研究京津冀PM2.5的影響因素,第一段模型是研究氣象因子、大氣污染物排放量與PM2.5的關系,第二段模型是研究不同的污染物排放源與大氣污染物排放量的關系。由于大氣污染物具有聚集的滯后效應,如果直接將氣象因子、經(jīng)濟因素放在一個模型里而忽略大氣污染物排放量的作用,勢必會造成模型估計結(jié)果的不精準。因此,本文的模型設定如下:

    PM2.5=β1+β2RHUit+β3TEMit+β4SSDit+β5PREit+β6Windit+β7Rainit+β8Dit+δ0Pollutionit+δ1Pollutionit-1+δ2Pollutionit-2+…+δkPollutionit-k+μit (1)

    lnPollution=α1+α2lnPIIit+α3lnCoalit+α4lnY_carit+

    瘙 窞 it (2)

    模型(1)中,PM2.5為PM2.5日濃度數(shù)值,RHUit為i城市t日的平均相對濕度,TEMit為i城市t日的平均溫度,SSDit為i城市t日的日照時數(shù),PREit為i城市t日的平均氣壓,Windit為i城市t日的平均風速,Rainit為i城市t日的降雨量,Pollutionit為i城市t日根據(jù)等距平滑法折算出來的大氣污染物排放量,Pollutionit-1、Pollutionit-2、…Pollutionit-k分別是Pollution的滯后1期、滯后2期、…、滯后k項,Dit為i城市t日沙塵暴天氣虛擬變量,β1代表總體均值截距項;μit代表隨機誤差項。模型(2)中PIIit為工業(yè)廢氣重度污染行業(yè)工業(yè)增加值,Coalit為煤炭消費量,Y_carit為機動車保有量,α1代表總體均值截距項, 瘙 窞 it代表隨機誤差項。

    2.3變量與數(shù)據(jù)

    PM2.5數(shù)據(jù)來源于中國空氣質(zhì)量在線監(jiān)測分析平臺,部分數(shù)據(jù)來源于青悅開放環(huán)境數(shù)據(jù)中心,選取京津冀代表性城市北京、天津、石家莊的PM2.5日度數(shù)據(jù)從2013年1月18日至2014年12月31日。相對濕度、氣溫、日照時數(shù)、壓強、風力、降雨等數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng),沙塵暴天氣數(shù)據(jù)來源于中國天氣網(wǎng)。大氣污染物排放量來源于《中國環(huán)境統(tǒng)計年鑒》《石家莊統(tǒng)計年鑒》。模型(1)中被解釋變量PM2.5濃度數(shù)據(jù)和相對濕度、氣溫、日照時數(shù)、壓強、風力、降雨等數(shù)據(jù)均是按天數(shù)據(jù),因此,為了統(tǒng)一數(shù)據(jù)口徑,本文首先根據(jù)《中國環(huán)境統(tǒng)計年鑒》中計算大氣污染物排放量方法,即根據(jù)公式大氣污染物排放量=∑工業(yè)源+生活源+機動車源+集中式污染治理設施源。其中生活源、機動車源、集中式污染治理設施源比較穩(wěn)定,可直接將生活源、機動車源、集中式污染治理設施源除于每年的天數(shù),得到按天的生活源、機動車源、集中式污染治理設施源的大氣污染物排放量。對于工業(yè)源而言,根據(jù)工業(yè)源折算系數(shù)=工業(yè)源大氣污染物排放量/工業(yè)增加值,得到工業(yè)源折算系數(shù),并通過月度工業(yè)源大氣污染物排放量=工業(yè)源折算系數(shù)×月度工業(yè)增加值,計算得到月度工業(yè)源大氣污染物排放量。第二,根據(jù)t月度工業(yè)源大氣污染物排放量增長率=t+1月工業(yè)源大氣污染物排放量/t月工業(yè)源大氣污染物排放量,得到該月工業(yè)源大氣污染物排放量增長率,并將該增長率根據(jù)等距平滑法平攤到每一天的大氣污染物排放量中去,得到按天的工業(yè)源大氣污染物排放量。最后,根據(jù)t日大氣污染物排放量=t日工業(yè)源大氣污染物排放量+t日生活源、機動車源、集中式污染治理設施源大氣污染物排放量計算出t日的大氣污染物排放總量。

    模型(2)是為了研究大氣污染物排放量同工業(yè)廢氣重度污染工業(yè)增加值、燃煤消費量、機動車保有量三者之間的關系,大氣污染物排放量數(shù)據(jù)計算方法跟模型1計算方法一樣。為了更好地體現(xiàn)京津冀大氣污染物排放量的演變規(guī)律,模型2中時間從2011年1月1日至2014年12月31日,按天的大氣污染物排放量計算方法跟模型1中一致。工業(yè)廢氣重度污染行業(yè)工業(yè)增加值根據(jù)前面五個重度污染產(chǎn)值之和計算得到,根據(jù)北京統(tǒng)計局、天津統(tǒng)計局、石家莊統(tǒng)計局統(tǒng)計的月度工業(yè)增加值,同樣根據(jù)上文中的等距平滑法,將其增長率平攤到每一天。燃煤消費量數(shù)據(jù)來源《中國能源統(tǒng)計年鑒》得到北京、天津、石家莊的年度煤炭消費量,同時根據(jù)張洪潮等的BP濾波分析得到我國2008—2013年月度煤炭修正系數(shù)[6],可得到月度煤炭消費量,同樣根據(jù)上文中的等距平滑法,將其增長率平攤到每一天得到按天的煤炭消費量數(shù)據(jù)。根據(jù)《2013年機動車污染防治年報》顯示,全國“黃標車”保有量僅占汽車保有量的10.7%,但其排放的一氧化碳(CO)、碳氫化合物(HC)、氮氧化物(NOX)、顆粒物(PM)分別占汽車排放量的49.0%、52.9%、52.4%、78.8%,因此黃標車占據(jù)了顆粒物(PM)排放量的絕大部分。由于無法獲得機動車保有量的數(shù)據(jù),本文將選擇黃標車保有量作為機動車的保有量的代理變量來研究。黃標車數(shù)據(jù)來源于《2011—2014年機動車污染防治年報》,通過Get Data Graph Digitizer軟件抓取各省市2011—2014年的數(shù)據(jù),同樣根據(jù)上文中的等距平滑法將其平攤計算得到按天的黃標車保有量。變量設定見表1。

    3實證檢驗結(jié)果及分析

    3.1PM2.5濃度描述性統(tǒng)計分析

    根據(jù)《環(huán)境空氣質(zhì)量標準》(GB3095—2012)中制定的環(huán)境空氣質(zhì)量指數(shù)(AQI)技術規(guī)定,24小時PM2.5的平均值標準值見表2。當24小時的PM2.5濃度超過115 ug/m3時,空氣質(zhì)量處于中度污染以上,對人群心臟、呼吸系統(tǒng)有影響。

    本文以優(yōu)良天氣和中度污染以上天氣做比較。從圖1可以看到,北京2013年優(yōu)良天氣比例為56%,2014年優(yōu)良天氣比例下降為55.6%,但中度污染以上天氣比

    例呈現(xiàn)出下降的趨勢,由2013年的26%下降到2014年的24%。天津優(yōu)良天氣比例有明顯上升的趨勢,2013年優(yōu)良天氣比例為44%,2014年上升為53%;中度污染以上天數(shù)比例有所下降,從2013年的25%下降到23%。石家莊優(yōu)良天氣比例上升幅度是三個城市最明顯的,2013年優(yōu)良天氣比例26.7%,2014年上升到36.7%;中度污染以上天氣比例也有明顯下降,從2013年的53.1%下降到38.6%。

    3.2PM2.5濃度影響因素計量結(jié)果

    排放的大氣污染物會形成聚集,導致PM2.5的濃度累積增加,前期大氣污染物具有聚集形成的滯后效應,前期排放的大氣污染物排放量對當期的PM2.5濃度可能存在正向的顯著影響。本研究運用阿爾蒙(Almon)提出的多項

    δi=γ1+γ2(i-c-)+γ3(i-c-)2……+γn+1(i-c-)n其中,i=1,2,3……m,m>n。

    首先確定模型滯后期長度。本研究采用以下兩種方法決定滯后期:①根據(jù)調(diào)整后的R2決定滯后期數(shù),直到調(diào)整后的R2不再增加;②通過施瓦茨信息準則(SC)、赤池信息準則(AIC)和LR統(tǒng)計量確定滯后期。通過反復檢驗發(fā)現(xiàn)滯后3期最符合大氣污染物擴散效應的實際情況,R2也最大,故選用大氣污染物排放量當期與滯后3期來檢驗對PM2.5的影響,也符合施瓦茨信息準則和赤池信息準則。大氣污染物的滯后系數(shù)選取了服從遠端約束(大氣污染物累積效應距離當期越遠,對當期PM2.5濃度影響越?。┒嘣貧w結(jié)果見表3。

    從表3可知,模型結(jié)果總體上比較良好,調(diào)整R2比較好,達到0.452,即平均相對濕度、平均氣溫、日照時數(shù)、平均相對氣壓、平均風速、降雨量、沙塵天氣、大氣污染物排放量及滯后3期能夠?qū)┙蚣絇M2.5濃度的45.25%做出解釋。模型總體顯著性在1%水平上通過檢驗,F(xiàn)值較大,說明模型中各因素對PM2.5濃度的共同影響顯著。

    在氣象因素中,風力通過稀釋擴散來降低當?shù)豍M2.5的濃度。研究表明,持續(xù)的微風或靜風會抑制大氣污染物的擴散,使近地面層大氣污染物聚集。當風速大于5 m/s,有利于污染物的擴散[17]。我們的統(tǒng)計發(fā)現(xiàn),當京津冀平均風速小于3 m/s時,京津冀PM2.5濃度平均值為111.56 ug/m3;而當京津冀平均風速大于5 m/s,京津冀PM2.5濃度平均值為49.41 ug/m3。

    溫度升高時,有利于上下層空氣產(chǎn)生對流交換,因此有利于大氣污染物的稀釋擴散;而溫度降低,大氣層相對穩(wěn)定不利于污染物的稀釋擴散。平均溫度每增加1℃,PM2.5濃度平均減少3.821 ug/m3。日照時數(shù)對污染物的影響主要體現(xiàn)在太陽輻射強時,加快了地面溫度上升,大氣垂直運動增強,不易產(chǎn)生逆溫層;日照時數(shù)較短時,太陽輻射較弱,氣溫較低,大氣層穩(wěn)定,容易產(chǎn)生逆溫,對污染物的稀釋和擴散不利。因此,日照時數(shù)每增加1 h,PM2.5濃度平均減少2.859 ug/m3。降水對大氣污染物可以起到?jīng)_刷和稀釋清除的作用,在雨水作用下,大氣中的一些大氣污染物能夠溶解在水中,可以降低大氣污染物的濃度,因此,降雨量每增加1mm, PM2.5濃度平均減少2.514 ug/m3。氣壓對污染物濃度主要體現(xiàn)在高氣壓時,空氣在垂直方向的對流較弱,容易形成逆溫,不利于污染物的擴散;低氣壓時,氣流上升較強,大氣污染物容易稀釋和擴散,因此氣壓與PM2.5濃度呈現(xiàn)負向關系。平均氣壓每增加1hpa,PM2.5濃度平均減少2.377 ug/m3。相對濕度與PM2.5有明顯的正向顯著關系。除有降水外,濕度大時,空氣中水氣較多,使得PM2.5附著在水氣中,在空氣中停留,因此平均相對濕度每增加1%,PM2.5濃度平均增加1.399 ug/m3。沙塵暴發(fā)生時,強風把地面沙塵物質(zhì)吹起并卷入空中,造成空氣混濁,水平能見度低。因此,沙塵暴天氣對PM2.5濃度具有顯著的正向影響。沙塵暴天氣發(fā)生時,PM2.5濃度平均增加139.942 ug/m3。而在大氣污染物排放量滯后變量中,當期大氣污染物每增加1t,對當期PM2.5濃度平均增加0.010 2 ug/m3;滯后一期污染物、滯后兩期污染物、滯后三期污染物每增加1t,對當期PM2.5濃度分別增加0.007 6 ug/m3、0.005 1 ug/m3、0.002 6 ug/m3。

    3.3主要污染源影響因素計量分析

    根據(jù)北京、天津、石家莊三個城市2013—2014年《環(huán)境狀況公報》對PM2.5進行源解析,可得知工業(yè)源、機動車、燃煤消耗量對PM2.5貢獻率最大,因此構建了模型(2)。通過檢驗發(fā)現(xiàn)固定效應和隨機效應均優(yōu)于OLS模型(其中固定效應模型的F值為529.29,隨機效應模型chi2(3)檢驗值為21 966.95,因為Hausman檢驗得到chi2(2)為174.4,p值為0,所以應拒絕隨機效應和固定效應的系數(shù)無系統(tǒng)差異的原假設,也就是說,采用固定效應模型更合適)。京津冀大氣污染物排放量的主要污染源影響因素估計結(jié)果見表4。

    由回歸結(jié)果表4可知,固定效應模型的整體顯著性的F值為529.29,P值為0.000,在1%顯著水平上顯著,表明模型(1) 整體是比較顯著的;工業(yè)廢氣重度污染行業(yè)工業(yè)增加值(PII)、燃煤消費量(Coal)、黃標車保有量(Y_car)的t值在1%水平上顯著,說明它們都能很好地解釋大氣污染物排放量的影響因素。模型(2)的整體顯著性的F值為557.7,P值為0,在1%顯著水平上顯著。

    通過固定效應模型可知,工業(yè)廢氣重度污染行業(yè)工業(yè)增加值的對數(shù)估計系數(shù)為0.173,t值檢驗在1%水平上顯著,表明在其他因素不變的情況下,一個地區(qū)工業(yè)廢氣重度污染行業(yè)工業(yè)增加值上升1%,該地區(qū)的大氣污染物排放量將上升0.173%。黃標車保有量的對數(shù)估計系數(shù)為0.199,t值檢驗在1%水平上顯著,表明在其他因素不變

    的情況下,一個地區(qū)黃標車保有量下降1%,該地區(qū)大氣污染物排放量將下降0.199%。煤炭消費量產(chǎn)值的對數(shù)估計系數(shù)為0.469,t值檢驗在1%水平上顯著,表明在其他因素不變的情況下,一個地區(qū)燃煤消費量增加1%,該地區(qū)大氣污染物排放量將上升0.469%,說明目前階段,煤炭消費量是京津冀地區(qū)大氣污染物排放量上升的重要因素。在該模型中經(jīng)過測算發(fā)現(xiàn),煤炭消費量對大氣污染物排放量貢獻率為47.27%,黃標車對大氣污染物排放量貢獻率為18.91%,工業(yè)廢氣重度污染行業(yè)工業(yè)增加值對大氣污染物排放量貢獻率為18.46%。

    4結(jié)論與政策建議

    4.1結(jié)論

    本研究通過使用京津冀PM2.5近兩年的日濃度數(shù)據(jù),提出兩段式的分布滯后模型將自然環(huán)境因素和經(jīng)濟因素結(jié)合起來,對京津冀PM2.5的影響因素進行了綜合分析。在第一段模型中構建了PM2.5和大氣污染物排放量的分布滯后模型,第二段模型中構建了不同的大氣污染源對大氣污染物排放量的影響因素模型。研究得出平均氣溫、平均風速、日照時數(shù)、平均氣壓、降雨量對PM2.5濃度具有顯著的負向作用,平均相對濕度、沙塵暴天氣對PM2.5濃度具有顯著的正向作用,大氣污染物排放量對PM2.5濃度具有聚集的滯后效應,當期大氣污染物排放量、滯后一期、滯后兩期、滯后三期大氣污染物對PM2.5濃度具有顯著的正向作用,且影響依次遞減。最后,構建了大氣污染物排放量的污染源模型,使用固定效應估計得出煤炭消費量、黃標車保有量、工業(yè)廢氣重度污染行業(yè)工業(yè)增加值對大氣污染物排放量具有顯著的正向作用,一個地區(qū)煤炭消費量對數(shù)、黃標車擁有量對數(shù)、工業(yè)廢氣重度污染行業(yè)工業(yè)增加值對數(shù)每增加1%,對該地區(qū)大氣污染物的排放量對數(shù)分別增加0.469%、0.199%和0.173%。

    4.2政策建議

    (1)優(yōu)化能源消費結(jié)構,減少煤炭的消費,大范圍推廣使用集中式供暖,完善清潔能源和PM2.5減排區(qū)域間補貼機制,分別在各空氣污染嚴重地區(qū)的重點行業(yè)進行試點。2014年,煤炭消費占我國能源消費總量的66%,電力、工業(yè)鍋爐、煤化工(煉焦等)以及建材窯爐消費了我國超過90%的煤炭,煤炭(尤其是散煤)的燃燒對PM2.5的貢獻很大。因此一是要完善優(yōu)化能源消費結(jié)構,優(yōu)化燃煤電廠脫硫脫硝的工藝;加大對清潔能源的利用,加強對新能源產(chǎn)業(yè)的補貼,同時建立政府綠色采購和補貼制度。二是要將煤炭為主的能源結(jié)構向清潔能源為主的能源結(jié)構轉(zhuǎn)變,大幅度減少散煤的燃燒。

    (2)加快改造黃標車,大力推廣新能源汽車,控制機動車的尾氣污染。一是加速改造黃標車,要么逐步淘汰黃標車,要么強制性對黃標車的排氣系統(tǒng)進行改造,加裝三元催化轉(zhuǎn)化器,使尾氣排放達國六標準。二是提升我國的油品質(zhì)量標準,降低油品中硫、烯烴及芳烴含量。三是大力推廣新能源汽車,公交、環(huán)衛(wèi)等行業(yè)和政府機關要率先使用新能源汽車、采取財政補貼等措施鼓勵個人購買。

    (3)優(yōu)化產(chǎn)業(yè)結(jié)構,限制工業(yè)廢氣重度污染行業(yè)。京津冀地區(qū)作為我國三大城市群之一,經(jīng)濟發(fā)展長期以來依賴化石能源,資源消耗大、環(huán)境污染嚴重,這也是京津冀大氣環(huán)境質(zhì)量產(chǎn)生問題的經(jīng)濟根源。因此一是調(diào)整優(yōu)化產(chǎn)業(yè)結(jié)構,嚴格限制新建工業(yè)廢氣重度污染行業(yè)的企業(yè),大力發(fā)展質(zhì)量效益型、科技先導型和資源節(jié)約型的綠色產(chǎn)業(yè)。二要將現(xiàn)有的空氣污染密集型企業(yè)向外搬遷,而且搬遷地點不能距離京津冀太近。

    (4)建議全面禁止石油焦的進口,并加強對燃燒高硫石油焦的監(jiān)管。盡管中國政府正在投入巨大的努力治理空氣污染和削減煤炭消費,但是卻出現(xiàn)一個監(jiān)管漏洞: 一些發(fā)電廠、水泥廠、玻璃廠使用石油焦作為替代燃料,在燃煤中故意混燒10-20%的高硫石油焦。這些高硫石油焦含硫量較高,包含多種重金屬污染(如汞、砷、鉻、鎳和鎘),在燃燒時還可能會釋放二惡英、氯化氫和氟化氫等污染物。2013年中國燃燒了3 300萬t石油焦。有關管理部門應出臺政策全面禁止石油焦的進口,并加強燃燒國產(chǎn)高硫石油焦的監(jiān)管,加快重點行業(yè)脫硫、脫硝、除塵改造工程建設,在空氣重度污染地區(qū)不得使用石油焦作為燃料。

    (編輯:劉照勝)

    參考文獻(References)

    [1]吳兌.近十年中國灰霾天氣研究綜述[J].環(huán)境科學學報, 2012,32(2):257-269.[WU Dui. Hazy weather research in China in the last decade: a review[J].Acta scientiae circumstantiae,2012,32(2):257-269.]

    [2]楊洪斌,鄒旭東,汪宏宇,等.大氣環(huán)境中PM2.5的研究進展與展望[J].氣象與環(huán)境學報,2012,28(3):77-82.[YANG Hongbin,ZOU Xudong,WANG Hongyu, et al. Study progress on PM2.5 in atmospheric environment[J].Journal of meteorology and environment,2012,28(3):77-82.]

    [3]梁延剛,WU Manchi,楊敬基.香港能見度、大氣懸浮粒子濃度與氣象條件的關系[J].氣象學報, 2007,66(3):461-469.[LEUNG Yinkong,WU Manchi,YEUNG Kingkay.A study on the relationship among visibility, atmospheric suspended particulate concentration and meteorological conditions in Hong Kong[J].Acta meteorologica sinica,2007,66(3):461-469.]

    [4]CHANG W L,BHAVE P V,BROWN S S, et al. Heterogeneous atmospheric chemistry, ambient measurements,and model calculations of PM2.5: a review[J], Aerosol science and technology,2011(4):665-670.

    [5]鄭玫,張延君,閆才青,等.中國PM2.5來源解析方法綜述[J].北京大學學報(自然科學版),2014,50(6):1141-1154.[ZHENG Mei, ZHANG Yanjun,YAN Caiqing, et al.Review of PM2.5 source apportionment methods in China[J].Acta scientiarum naturalium universitatis pekinensis,2014,50(6): 1141-1154.]

    [6]周靜,劉松華,譚譯,等.蘇州市大氣顆粒物(PM2.5)工業(yè)源排放清單[J].環(huán)境科學學報2016,2(2):1-12.[ZHOU Jing,LIU Songhua,TAN Yi,et al.The emission inventory of atmospheric fine particles from industrial sources in Suzhou City[J].Acta scientiae circumstantiae,2016,2(2):1-12.]

    [7]朱增銀,李冰,趙秋月,等.對國內(nèi)外PM2.5研究及控制對策的回顧與展望[J].環(huán)境科技, 2013,26(1):70-74.[ZHU Zengyin, LI Bing, ZHAO Qiuyue,et al. A review and outlook of the research and pollution control on PM2.5 in domestic and overseas[J].Environmental science and technology,2013,26(1):70-74.]

    [8]郭清彬,程學豐,侯輝,等.冬季大氣中PM10和PM2.5污染特征及形貌分析[J].中國環(huán)境監(jiān)測,2010,26(4):55-58.[GUO Qingbin,CHENG Xuefeng, HOU Hui,et al.Pollution characteristics and morphology analysis of PM10 and PM2.5 in winter[J].Environmental monitoring in China,2010,26(4):55-58.]

    [9]于建華,虞統(tǒng),魏強,等.北京地區(qū)PM10和PM2.5質(zhì)量濃度的變化特征[J].環(huán)境科學研究,2004,17(1):45-47.[YU Jianhua,YU Tong,WEI Qiang,et al.Characteristics of mass concentration variations of PM10 and PM2.5 in Beijing area[J].Research of environmental sciences,2004,17(1):45-47.]

    [10]張淑平,韓立建,周偉奇,等.冬季PM2.5的氣象影響解析[J].生態(tài)學報, 2016,36(24):1-11.[ZHANG Shuping, HAN Lijian,ZHOU Weiqi,e al. Relationships between fine particulate matter and meteorological factors in winter at typical Chinese cities[J].Acta ecologica sinica,2016,36(24):1-11.]

    [11]劉瑞婷,韓志偉,李嘉偉.北京冬季霧霾事件的氣象特征分析[J].氣象與環(huán)境研究, 2014,19(2):164-172.[LIU Ruiting,HAN Zhiwei, LI Jiawei.Analysis of meteorogical characteristics during winter haze events in Beijing[J].Climatic and environmental research,2014,19(2):164-172.]

    [12]孟曉艷,余予,張志富,等.2013年1月京津冀地區(qū)強霧霾頻發(fā)成因初探[J].環(huán)境科學與技術,2014(1):190-194.[MENG Xiaoyan,YU Yu,ZHANG Zhifu,et al. Preliminary study of the dense fog and haze events formation over BeijingTianjinandHebei region in January of 2013[J].Environmental science and technology,2014,37(1):190-194.]

    [13]何楓,馬棟棟.霧霾與工業(yè)化發(fā)展的關聯(lián)研究——中國74個城市的實證研究[J].軟科學,2015(6):110-114.[HE Feng,MA Dongdong. Association study between haze pollution and industrialization[J].Soft science, 2015(6):110-114.]

    [14]馬麗梅,張曉.中國霧霾污染的空間效應及經(jīng)濟、能源結(jié)構影響[J].中國工業(yè)經(jīng)濟, 2014,4(4):19-31.[MA Limei,ZHANG Xiao.The spatial effect of Chinas haze pollution and the impact from economic change and energy structure[J].China industrial economics,2014,4(4):19-31.]

    [15]張家峰.農(nóng)村工業(yè)化與農(nóng)村點源污染關系研究——以江蘇省為例[D].南京:南京農(nóng)業(yè)大學,2013:61-81.[ZHANG Jiafeng,Research on the relationship between rural industrialization and rural industrial point source pollution:taking Jiangsu Province as an example[D].Nanjing:Nanjing Agricultural University,2013:61-81.]

    [16]張洪潮,王澤江,李曉利,等.中國煤炭消費需求波動規(guī)律及成因分析[J].中國人口·資源與環(huán)境,2014,24(1):94-101.[ZHANG Hongchao,WANG Zejiang, LI Xiaoli, et al. The fluctuation pattern of Chinese coal consumption demand and cause analysis[J].China population,resources and environment,2014,24(1):94-101.]

    [17]GROSSMAN G, KRUEGER A. Economic growth and the environmen[J].Quarterly journal of ecomomics,1995,110(2),353-377.

    [18]蔣洪強,張靜,王金南,等.中國快速城鎮(zhèn)化邊際環(huán)境污染效應變化實證分析[J].生態(tài)環(huán)境學報,2012,21(2):293-297.[JIANG Hongqiang,ZHANG Jing,WANG Jinnan,et al.Accounting model and demonstration for marginal environmental effects of the rapid urbanization in China[J].Ecology and environmental sciences,2012,21(2):293-297.]

    [19]楊義彬.成都市大氣污染及氣象條件影響分析[J].四川氣象,2004,89(3):40-43.[YANG Yibin.The air pollution of Chengdu City and the impact analysis of meteorological conditions[J].Journal of Sichuan meteorology,2004,89(3):40-43.]

    [20]SELDEN T, SONG D. Neoclassical growth, the J Curve for abatement and the internet U Curve for pollution[J].Journal of environmental economics and management,1995,29(2),162-168.

    [21]LOPOZ R. The environment as a factor of production:the effects of economic growth and trade liberalization[J].Journal of environmental economics and management,1994,27(2),163-184.

    [22]BROCK W, TAYLOR M. Economics growth and the environment: a review of theory and empirics[J].Handbook of economic growth,2005(1):1749-1821.

    [23]張鳳青.高汽車保有量果真是中國PM2.5污染的主因嗎?[J].經(jīng)濟資料譯叢, 2015,2(4):44-51.[ZHANG Fengqing.Is the high car ownership really the main cause of Chinas PM2.5 pollution?[J].Journal of translation from foreign literature of economics,2015,2(4):44-51.]

    Study on the main influencing factors and their intrinsic relations of PM2.5 in BeijingTianjinHebei

    ZHOU Shudong1,2OUYANG Weiqing1GE Jihong1

    (1.College of Economics and Management, Nanjing Agricultural University, Nanjing Jiangsu 210095, China; 2.China Center for Food Security Studies, Nanjing Agricultural University,Nanjing Jiangsu 210095, China)

    AbstractAir pollutants are the internal causes of haze weather, meteorological conditions are external factors of haze weather. In this study, a twostage distribution hysteresis model for PM2.5 concentration was established, and the influencing factors of PM2.5 were analyzed in combination with natural and economic factors. In the first stage, a distribution hysteresis model for PM2.5 concentration and air pollutants emissions was established. In the second stage, an influencing factors model for air pollutants emission was constructed. Air pollutants sources include industrial emission, life emission, motor vehicles emission, and centralized pollution facilities. In industrial emission sources, the polluting intensive industries for waste gas are the major contributors to air pollutant emissions.In the life emission sources, coal consumption greatly affected air pollutant emissions, which contributed the increase of PM2.5 concentration during winter heating season. In the motor vehicles source, although the proportion of vehicles ownership of wagons with yellow license plates is only about 10%, but they contributed a large part of particulate matters.Based on the daily data of PM2.5 of the representative cities of Beijing, Tianjin and Hebei, the negative and positive impacts on PM2.5 concentration of average temperature, average wind speed, sunshine hours, average air pressure, rainfall, average relative humidity, dust storm and other factors are measured. We find that the air pollutant emission has a lag effect on the accumulation of PM2.5 concentration, e.g. current period air pollutant emissions,lag one, lag two, and lag three have significant positive effects on the PM2.5 concentration, and the influencing effects are decreasing according to the lagging order. A model of influencing factors of air pollutant emission concluded that the consumption of coal, the industrial added value of waste gas polluting intensive industries, and the amount of wagons with yellow license plates do contribute a significant influence on air pollutant emission.Some countermeasures and policy recommendation to optimize energy consumption structure and industrial structure and reduce air pollutant emission are suggested in this paper.

    Key wordsPM2.5; air pollutant emission; Twostage Distribution Hysteresis Model; BeijingTianjinHebei

    猜你喜歡
    灰霾排放量大氣
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    天然氣輸配系統(tǒng)甲烷排放量化方法
    煤氣與熱力(2021年6期)2021-07-28 07:21:40
    黑龍江省碳排放量影響因素研究
    成都城區(qū)冬季一次灰霾期間PM1化學成分特征
    秸稈露天燃燒對北方灰霾天氣的影響分析
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    全國機動車污染物排放量
    ——《2013年中國機動車污染防治年報》(第Ⅱ部分)
    江蘇省火力發(fā)電機組二氧化碳排放量估算
    安徽首座超級灰霾監(jiān)測站將投入使用
    久久久久人妻精品一区果冻| 99久久综合免费| 女性被躁到高潮视频| 丝瓜视频免费看黄片| 亚洲精品一二三| 成人毛片a级毛片在线播放| 欧美人与性动交α欧美精品济南到 | 国产高清国产精品国产三级| 亚洲五月色婷婷综合| 精品国产一区二区三区四区第35| 插逼视频在线观看| 亚洲av男天堂| 中文字幕精品免费在线观看视频 | 欧美精品国产亚洲| av电影中文网址| 婷婷色av中文字幕| 国产又色又爽无遮挡免| 1024视频免费在线观看| 一级,二级,三级黄色视频| 亚洲成人av在线免费| 熟女人妻精品中文字幕| 亚洲欧美中文字幕日韩二区| 久久久久精品性色| √禁漫天堂资源中文www| 久久亚洲国产成人精品v| kizo精华| 久久久久国产网址| 久久久久精品人妻al黑| 亚洲欧美成人精品一区二区| 亚洲精品美女久久久久99蜜臀 | 色94色欧美一区二区| 男女边摸边吃奶| 国产国拍精品亚洲av在线观看| av电影中文网址| av在线app专区| 免费看光身美女| 亚洲精品,欧美精品| 伦理电影免费视频| 超碰97精品在线观看| 两个人免费观看高清视频| 久久鲁丝午夜福利片| 999精品在线视频| 成年av动漫网址| 久久亚洲国产成人精品v| 婷婷成人精品国产| 久久国产精品男人的天堂亚洲 | 欧美日韩精品成人综合77777| 日韩制服骚丝袜av| 色婷婷av一区二区三区视频| 成人18禁高潮啪啪吃奶动态图| 色视频在线一区二区三区| 亚洲欧美清纯卡通| 亚洲国产欧美日韩在线播放| 少妇精品久久久久久久| 两个人免费观看高清视频| 考比视频在线观看| 久久久国产精品麻豆| 制服丝袜香蕉在线| 免费女性裸体啪啪无遮挡网站| 精品人妻一区二区三区麻豆| 久久久久久久亚洲中文字幕| kizo精华| 国产精品久久久久久久电影| 最新的欧美精品一区二区| 日韩中字成人| 国产高清国产精品国产三级| 我的女老师完整版在线观看| 丰满迷人的少妇在线观看| 咕卡用的链子| 久久久久久久国产电影| 亚洲欧美成人精品一区二区| 久久久亚洲精品成人影院| 最近最新中文字幕免费大全7| 人人妻人人添人人爽欧美一区卜| a级毛片黄视频| 高清毛片免费看| 亚洲,欧美精品.| 国产永久视频网站| 少妇的丰满在线观看| 久久韩国三级中文字幕| 啦啦啦啦在线视频资源| 午夜免费观看性视频| 亚洲经典国产精华液单| 亚洲精品自拍成人| 亚洲成国产人片在线观看| 男男h啪啪无遮挡| 亚洲欧美日韩另类电影网站| 久久精品国产鲁丝片午夜精品| 久久这里有精品视频免费| 免费av不卡在线播放| 免费高清在线观看视频在线观看| 爱豆传媒免费全集在线观看| 日本av手机在线免费观看| a 毛片基地| videos熟女内射| 免费看光身美女| 国产永久视频网站| 亚洲综合色惰| 看免费av毛片| 亚洲国产色片| 老女人水多毛片| 蜜桃在线观看..| 免费观看在线日韩| 国产精品无大码| av又黄又爽大尺度在线免费看| 免费观看性生交大片5| 一二三四在线观看免费中文在 | 午夜免费男女啪啪视频观看| 全区人妻精品视频| av免费观看日本| 色吧在线观看| 久久久久久久精品精品| 亚洲成人一二三区av| 午夜视频国产福利| 亚洲精品aⅴ在线观看| xxx大片免费视频| 一本久久精品| 伊人久久国产一区二区| videosex国产| 婷婷色综合大香蕉| 伦理电影免费视频| 亚洲人与动物交配视频| 精品一品国产午夜福利视频| 免费观看av网站的网址| www.色视频.com| 亚洲,一卡二卡三卡| 在线免费观看不下载黄p国产| 久久狼人影院| 久久免费观看电影| 国产精品欧美亚洲77777| 亚洲国产色片| 久久国产精品男人的天堂亚洲 | 日本爱情动作片www.在线观看| 狂野欧美激情性bbbbbb| 亚洲丝袜综合中文字幕| 成人毛片a级毛片在线播放| 高清av免费在线| 亚洲综合精品二区| videos熟女内射| 成人亚洲欧美一区二区av| 街头女战士在线观看网站| 成人漫画全彩无遮挡| 三级国产精品片| 亚洲第一av免费看| 国产男女超爽视频在线观看| 91aial.com中文字幕在线观看| 免费大片18禁| 国产亚洲一区二区精品| 国产高清三级在线| 老司机影院毛片| 国产极品粉嫩免费观看在线| 亚洲精品av麻豆狂野| 国产一级毛片在线| 97在线人人人人妻| 插逼视频在线观看| 亚洲精品乱久久久久久| 久久久久精品性色| 亚洲精华国产精华液的使用体验| 久久99精品国语久久久| 搡老乐熟女国产| 人妻系列 视频| 国产日韩欧美在线精品| 男女免费视频国产| 最近最新中文字幕免费大全7| 美国免费a级毛片| 免费黄频网站在线观看国产| 国产精品.久久久| 精品一区二区三区四区五区乱码 | 亚洲精品日本国产第一区| 极品少妇高潮喷水抽搐| 免费女性裸体啪啪无遮挡网站| 国产激情久久老熟女| 秋霞在线观看毛片| 十八禁网站网址无遮挡| 免费播放大片免费观看视频在线观看| 国产成人a∨麻豆精品| 国产日韩欧美亚洲二区| 永久网站在线| 一级毛片电影观看| 人妻系列 视频| 捣出白浆h1v1| 人体艺术视频欧美日本| 久久精品国产自在天天线| 久久人人爽av亚洲精品天堂| 国产精品久久久久久精品古装| 一级片免费观看大全| 欧美日韩视频精品一区| av天堂久久9| 久久99热6这里只有精品| 精品人妻在线不人妻| 国产永久视频网站| 啦啦啦在线观看免费高清www| 午夜福利影视在线免费观看| 久久精品国产鲁丝片午夜精品| 十八禁网站网址无遮挡| 2022亚洲国产成人精品| √禁漫天堂资源中文www| 国产 精品1| 日韩视频在线欧美| 国产av码专区亚洲av| 欧美人与性动交α欧美精品济南到 | 国产成人午夜福利电影在线观看| 99国产综合亚洲精品| 精品一区二区三区四区五区乱码 | 视频中文字幕在线观看| 秋霞在线观看毛片| 国产精品人妻久久久久久| 日本与韩国留学比较| 久久午夜综合久久蜜桃| 国产无遮挡羞羞视频在线观看| 亚洲人成网站在线观看播放| 蜜臀久久99精品久久宅男| 亚洲情色 制服丝袜| 丝袜在线中文字幕| 熟女av电影| 亚洲av日韩在线播放| 久久久久精品久久久久真实原创| 国产精品三级大全| 丰满饥渴人妻一区二区三| 亚洲国产精品999| 香蕉国产在线看| 在线观看免费高清a一片| 久久精品国产亚洲av天美| 国产有黄有色有爽视频| 免费久久久久久久精品成人欧美视频 | 色视频在线一区二区三区| 激情视频va一区二区三区| 岛国毛片在线播放| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产综合久久久 | 亚洲综合色惰| 国内精品宾馆在线| tube8黄色片| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 香蕉精品网在线| 在线精品无人区一区二区三| 婷婷色麻豆天堂久久| 国产综合精华液| 三级国产精品片| 免费人成在线观看视频色| 最近最新中文字幕大全免费视频 | 青春草亚洲视频在线观看| 国产无遮挡羞羞视频在线观看| 国产精品不卡视频一区二区| 男人操女人黄网站| 欧美性感艳星| 少妇人妻精品综合一区二区| 十八禁高潮呻吟视频| 九九在线视频观看精品| 又大又黄又爽视频免费| 亚洲色图综合在线观看| 国产亚洲一区二区精品| 免费大片黄手机在线观看| 五月伊人婷婷丁香| 久久 成人 亚洲| 亚洲欧美色中文字幕在线| 日本与韩国留学比较| 国语对白做爰xxxⅹ性视频网站| 久久狼人影院| 日本午夜av视频| 精品卡一卡二卡四卡免费| tube8黄色片| 亚洲精品乱码久久久久久按摩| 免费大片18禁| 久久久久久人妻| 国产高清不卡午夜福利| 精品久久国产蜜桃| 欧美国产精品va在线观看不卡| 国产精品国产av在线观看| 成人手机av| 中文字幕av电影在线播放| 亚洲国产精品国产精品| 亚洲av电影在线进入| 在线观看三级黄色| 一二三四在线观看免费中文在 | 高清黄色对白视频在线免费看| 亚洲欧洲精品一区二区精品久久久 | 久热久热在线精品观看| 人人妻人人爽人人添夜夜欢视频| 婷婷色麻豆天堂久久| 多毛熟女@视频| av女优亚洲男人天堂| 精品一区二区三卡| 男女边摸边吃奶| 9色porny在线观看| 国产成人aa在线观看| 欧美变态另类bdsm刘玥| 欧美日韩视频精品一区| 嫩草影院入口| 国产日韩欧美亚洲二区| 满18在线观看网站| 91久久精品国产一区二区三区| 中文字幕亚洲精品专区| 黑丝袜美女国产一区| 色吧在线观看| 少妇人妻精品综合一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美老熟妇乱子伦牲交| 草草在线视频免费看| 人人妻人人添人人爽欧美一区卜| 91成人精品电影| 国产精品偷伦视频观看了| 日韩免费高清中文字幕av| 欧美人与善性xxx| 最近最新中文字幕免费大全7| 国产欧美另类精品又又久久亚洲欧美| 欧美xxxx性猛交bbbb| 亚洲中文av在线| 18禁裸乳无遮挡动漫免费视频| 一级黄片播放器| videosex国产| 日韩精品有码人妻一区| 99久国产av精品国产电影| 不卡视频在线观看欧美| 18禁裸乳无遮挡动漫免费视频| 秋霞在线观看毛片| 看十八女毛片水多多多| 国国产精品蜜臀av免费| 日日啪夜夜爽| 国产一区二区三区av在线| 美国免费a级毛片| 亚洲在久久综合| 久久久久视频综合| 国产成人欧美| 午夜福利网站1000一区二区三区| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| a级毛片在线看网站| 人人澡人人妻人| 精品人妻一区二区三区麻豆| 国产免费一级a男人的天堂| 男女高潮啪啪啪动态图| 99久国产av精品国产电影| 九九爱精品视频在线观看| 99国产综合亚洲精品| 一级毛片电影观看| 亚洲av综合色区一区| 亚洲中文av在线| 日韩中文字幕视频在线看片| 成年人午夜在线观看视频| 欧美精品高潮呻吟av久久| 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久 | 国产成人aa在线观看| 亚洲精品美女久久av网站| 中国三级夫妇交换| 激情视频va一区二区三区| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| a级毛色黄片| 黄色一级大片看看| 五月伊人婷婷丁香| 久久午夜福利片| 在线观看人妻少妇| 永久免费av网站大全| 欧美另类一区| 欧美丝袜亚洲另类| 九草在线视频观看| 国产极品粉嫩免费观看在线| 亚洲国产精品成人久久小说| 国产毛片在线视频| 99香蕉大伊视频| 国产精品久久久久久av不卡| 老女人水多毛片| av网站免费在线观看视频| 91久久精品国产一区二区三区| a级毛色黄片| 十八禁网站网址无遮挡| 欧美精品一区二区免费开放| 毛片一级片免费看久久久久| 在线亚洲精品国产二区图片欧美| 深夜精品福利| 久久久久精品久久久久真实原创| 国语对白做爰xxxⅹ性视频网站| 国产高清三级在线| 国产女主播在线喷水免费视频网站| 亚洲av电影在线观看一区二区三区| 国产成人欧美| 日本-黄色视频高清免费观看| 日韩av免费高清视频| 免费黄色在线免费观看| 国产在线视频一区二区| 免费不卡的大黄色大毛片视频在线观看| 欧美人与性动交α欧美软件 | 精品人妻偷拍中文字幕| www.熟女人妻精品国产 | 1024视频免费在线观看| 国产综合精华液| 久久精品国产鲁丝片午夜精品| 男的添女的下面高潮视频| 国产片内射在线| 精品亚洲成国产av| 伊人亚洲综合成人网| 久久午夜综合久久蜜桃| 亚洲av在线观看美女高潮| 边亲边吃奶的免费视频| 男女国产视频网站| 伊人久久国产一区二区| 大香蕉久久成人网| 高清不卡的av网站| 一级爰片在线观看| av网站免费在线观看视频| 黑人欧美特级aaaaaa片| 极品少妇高潮喷水抽搐| 男女边摸边吃奶| 亚洲欧美日韩卡通动漫| 夜夜爽夜夜爽视频| 亚洲精品乱码久久久久久按摩| 国产精品人妻久久久影院| av卡一久久| 国产免费又黄又爽又色| 亚洲av中文av极速乱| 9色porny在线观看| 欧美精品亚洲一区二区| 人妻一区二区av| 成人国产麻豆网| 亚洲国产精品一区三区| 国产深夜福利视频在线观看| 亚洲欧美一区二区三区黑人 | 日本猛色少妇xxxxx猛交久久| 亚洲国产成人一精品久久久| 男女下面插进去视频免费观看 | 成年人免费黄色播放视频| 午夜久久久在线观看| 乱人伦中国视频| 国产黄频视频在线观看| 婷婷色综合www| 91精品三级在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃国产av成人99| 街头女战士在线观看网站| 亚洲国产精品成人久久小说| 九草在线视频观看| 丰满迷人的少妇在线观看| 亚洲人成网站在线观看播放| 亚洲精品乱久久久久久| 韩国精品一区二区三区 | 99久久综合免费| 精品人妻在线不人妻| 亚洲成国产人片在线观看| 女的被弄到高潮叫床怎么办| 欧美亚洲 丝袜 人妻 在线| 水蜜桃什么品种好| 亚洲精品aⅴ在线观看| 国产片内射在线| 久久精品国产a三级三级三级| 婷婷色综合大香蕉| 黄色 视频免费看| 亚洲天堂av无毛| 久久99热6这里只有精品| 2022亚洲国产成人精品| 王馨瑶露胸无遮挡在线观看| 日韩,欧美,国产一区二区三区| 在线观看美女被高潮喷水网站| 蜜桃在线观看..| 欧美国产精品一级二级三级| 欧美97在线视频| 精品一区二区三区四区五区乱码 | tube8黄色片| 久久久久精品久久久久真实原创| 在线观看www视频免费| 国产视频首页在线观看| 久久精品人人爽人人爽视色| av播播在线观看一区| 亚洲内射少妇av| 亚洲av在线观看美女高潮| 国产成人91sexporn| 久久精品国产鲁丝片午夜精品| 国产一区二区在线观看av| 中文字幕免费在线视频6| 久久精品国产自在天天线| 国产精品蜜桃在线观看| 18在线观看网站| 三级国产精品片| 久久99蜜桃精品久久| 亚洲精品第二区| 日韩av不卡免费在线播放| 国产老妇伦熟女老妇高清| 国产成人精品无人区| 国产片特级美女逼逼视频| 中文字幕人妻丝袜制服| 99久国产av精品国产电影| 永久免费av网站大全| 黄色视频在线播放观看不卡| 男女下面插进去视频免费观看 | 在线观看免费视频网站a站| 亚洲av福利一区| 亚洲第一av免费看| 成人国语在线视频| av线在线观看网站| 69精品国产乱码久久久| a 毛片基地| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频 | 精品卡一卡二卡四卡免费| 色94色欧美一区二区| 美女内射精品一级片tv| 精品人妻偷拍中文字幕| 久久久久久人人人人人| 黑人高潮一二区| 尾随美女入室| 丰满少妇做爰视频| 亚洲一区二区三区欧美精品| 内地一区二区视频在线| 一级a做视频免费观看| 日韩,欧美,国产一区二区三区| 高清欧美精品videossex| 蜜臀久久99精品久久宅男| 乱人伦中国视频| 久久久精品94久久精品| 色94色欧美一区二区| 国产精品一区www在线观看| 欧美日韩精品成人综合77777| 午夜影院在线不卡| 久久综合国产亚洲精品| 美女脱内裤让男人舔精品视频| 亚洲五月色婷婷综合| 最新中文字幕久久久久| 日韩欧美一区视频在线观看| av国产精品久久久久影院| 寂寞人妻少妇视频99o| 成人二区视频| 在线观看三级黄色| 亚洲精品美女久久av网站| 亚洲av免费高清在线观看| 香蕉国产在线看| 久久久久久久亚洲中文字幕| 观看美女的网站| 欧美xxⅹ黑人| 婷婷成人精品国产| 亚洲综合色网址| 99国产精品免费福利视频| 欧美人与性动交α欧美软件 | 精品亚洲成国产av| 欧美性感艳星| 亚洲第一区二区三区不卡| 另类精品久久| 日韩中文字幕视频在线看片| 亚洲精品美女久久久久99蜜臀 | 超色免费av| 极品少妇高潮喷水抽搐| 久久热在线av| 中文字幕人妻丝袜制服| 999精品在线视频| 亚洲av福利一区| 黄色配什么色好看| 国产成人aa在线观看| 亚洲国产av新网站| 日韩一区二区视频免费看| 久久韩国三级中文字幕| 亚洲精品日韩在线中文字幕| 高清在线视频一区二区三区| 9色porny在线观看| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| 高清av免费在线| 性高湖久久久久久久久免费观看| 97精品久久久久久久久久精品| 精品一区二区免费观看| 亚洲,欧美精品.| 日韩免费高清中文字幕av| 香蕉国产在线看| 男人添女人高潮全过程视频| 黄色怎么调成土黄色| 久久99一区二区三区| 天天躁夜夜躁狠狠躁躁| 色哟哟·www| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 制服人妻中文乱码| 国产一区二区三区综合在线观看 | 精品熟女少妇av免费看| 少妇 在线观看| 亚洲天堂av无毛| 国产乱人偷精品视频| 五月伊人婷婷丁香| 日本午夜av视频| 天堂俺去俺来也www色官网| 如日韩欧美国产精品一区二区三区| 狠狠精品人妻久久久久久综合| 国产成人免费观看mmmm| 免费人成在线观看视频色| 青春草亚洲视频在线观看| 伦精品一区二区三区| 久久 成人 亚洲| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区视频在线| 超色免费av| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 国产精品国产av在线观看| 国产av国产精品国产| 国产精品久久久久久久电影| 国产黄色视频一区二区在线观看| 69精品国产乱码久久久| 在线观看人妻少妇| 欧美精品av麻豆av| 女性生殖器流出的白浆| 欧美人与性动交α欧美软件 | 咕卡用的链子| 交换朋友夫妻互换小说| 另类精品久久| 中国三级夫妇交换| 亚洲精品一二三| 日日爽夜夜爽网站| 男女下面插进去视频免费观看 | av福利片在线| 亚洲第一av免费看| 国产精品女同一区二区软件| 人人妻人人澡人人爽人人夜夜| 女性被躁到高潮视频| 在线观看免费日韩欧美大片| 国产在线一区二区三区精| 黑人猛操日本美女一级片| 精品少妇黑人巨大在线播放|