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

    京津冀城市群大氣污染傳輸規(guī)律研究——兩組排放清單的比較分析

    2019-11-28 10:51:06盧亞靈
    中國環(huán)境科學 2019年11期
    關鍵詞:排放量京津冀污染物

    王 媛,李 玥,喬 治*,盧亞靈,2

    京津冀城市群大氣污染傳輸規(guī)律研究——兩組排放清單的比較分析

    王 媛1,李 玥1,喬 治1*,盧亞靈1,2

    (1.天津大學環(huán)境科學與工程學院,天津 300350;2.環(huán)境保護部環(huán)境規(guī)劃院,北京 100012)

    利用WRF/CALPUFF耦合模型,在同樣重污染氣象條件下,選擇了當下模擬應用最廣的兩組排放清單對4種主要污染物(NO,SO2,PM2.5和PM10)進行京津冀城市間區(qū)域傳輸貢獻比較分析.其中一組排放清單來自政府的環(huán)境統(tǒng)計(以下稱:環(huán)統(tǒng)排放清單),另一組排放清單是來自中國多尺度排放清單(以下稱:MEIC排放清單).污染物的濃度空間分布表明,兩種排放清單下污染物濃度均呈現(xiàn)北部以唐山中心,南部以石家莊-邯鄲為中心的分布特征,均是由兩個濃度最高的中心向外逐漸降低.但是環(huán)統(tǒng)排放清單下模擬的污染物濃度高值區(qū)范圍更大,更接近實際監(jiān)測數(shù)據(jù).基于不同的清單輸入,一些城市的傳輸角色存在一些差異.例如,對于4種污染物,兩種排放清單模擬出的滄州與周邊城市的凈傳輸方向完全相反,在MEIC排放清單中,滄州以向外凈傳輸為主,即為源;而在環(huán)統(tǒng)排放清單中,滄州則變成了凈輸入城市,即為匯.這些結論將影響大氣聯(lián)防聯(lián)控中各城市源匯責任的認定,在實際環(huán)境管理中應注重多源數(shù)據(jù)的選擇、驗證和比較.

    CALPUFF模型;京津冀城市群;傳輸規(guī)律;污染物傳輸網絡

    京津冀地區(qū)人口密度大,經濟發(fā)展迅速,加之秋冬季靜穩(wěn)天氣多發(fā),是我國大氣污染防治重點區(qū)域之一[1-2].目前國內外研究者對大氣污染物跨界傳輸問題已經展開了廣泛研究,一些學者從國家的角度研究了污染物在國家間的傳輸規(guī)律[3-4],Itahashi等[3]研究了污染物在日本和中國之間的傳輸規(guī)律,并證明中國東北部的氮氧化物可以傳輸?shù)饺毡揪硟?Lin等[4]將經濟排放和大氣傳輸模型結合,分析了與貿易有關的中國大氣污染物排放對美國大氣環(huán)境的影響.薛文博等[5]將CAMx和PSAT模型相結合,建立了我國31個省市間的PM2.5傳輸矩陣,證實跨區(qū)域傳輸對重點區(qū)域的PM2.5污染具有顯著貢獻,其中,京津冀、長三角、珠三角受外源的貢獻分別達22%、37%和28%.京津冀作為中國大氣污染控制的重點區(qū)域之一,一些學者探討了京津冀城市群空氣污染的傳輸規(guī)律.例如Wang等[6]應用CAMx和PSAT技術來量化在一次重污染事件時,不同地區(qū)的排放對北京PM2.5濃度的貢獻.他們的研究表明北京當?shù)嘏欧诺呢暙I占比83.6%,周邊區(qū)域對它的貢獻占比9.4%.

    大氣重污染三大成因是污染排放、氣象條件和區(qū)域傳輸,其中污染排放是主因和內因,氣象條件是外因,而區(qū)域傳輸受前兩個原因影響.目前,對大氣傳輸?shù)难芯慷鄳每諝赓|量模型進行模擬,但模擬結果常常存在差異.雖然影響因素眾多,但排放數(shù)據(jù)無疑是重要的影響因素之一[7].目前尚未有文章進行過兩組污染物排放清單的比較,探討不同的排放清單對傳輸規(guī)律造成的影響.排放清單選擇不同,對于城市尺度污染責任的界定可能會產生不同的結果.本文就從此角度來探討排放清單對城市尺度的污染責任的界定造成的影響.在同樣氣象條件下,選擇了當下模擬應用最廣泛的兩組排放清單進行傳輸規(guī)律的比較分析,其中一組排放清單來自政府的環(huán)境統(tǒng)計數(shù)據(jù)(以下稱:環(huán)統(tǒng)排放清單),另一組排放清單是來自清華大學開發(fā)的中國多尺度排放清單[8](Multi-resolution Emission Inventory for China,簡稱MEIC排放清單).研究采用WRF/CALPUFF耦合模型進行排放清單的比較分析.選取2012年污染最嚴重的12月為模擬月份,以京津冀地區(qū)為模擬區(qū)域,并選取了SO2、NO、PM2.5和PM10這四種典型的大氣污染物進行傳輸規(guī)律的比較分析.

    本文研究主要從以下兩個方面開展:(1)從濃度空間分布方面對兩個排放清單模擬的結果進行比較;(2)從凈傳輸方向方面對兩個排放清單模擬的結果進行比較,重點分析兩種清單導致的污染物傳輸方向不一致的傳輸流.

    1 研究方法

    1.1 氣象模式

    采用中尺度氣象模式WRF耦合空氣質量模型CALPUFF的氣象模塊CALMET對京津冀的氣象要素進行模擬.選用WRF/CALMET模式系統(tǒng)可以有效地模擬不同氣象條件下地表所表現(xiàn)出的不同變化特點,解決微氣象參數(shù)下的各種效應[8].

    圖1 網格設置示意

    氣象場模擬選取京津冀2012年污染最嚴重的12月進行模擬,采用兩層嵌套的中尺度氣象模型WRF,模型使用蘭伯特正形圓錐(LCC)投影,模擬區(qū)域中心為116.21°E,39.58°N,第一層網格(D1)分辨率為18km,網格數(shù)為60′63,第二層網格(D2)分辨率為6km,網格數(shù)為96′123.第一層網格包括京津冀及其周邊地區(qū),山東、山西、內蒙古以及遼寧的部分地區(qū),第二層網格完整的包括了京津冀地區(qū),兩層網格中心重疊.模擬區(qū)域如圖1所示.

    CALPUFF是大氣環(huán)境影響評價的一種模型[10-11],在空氣質量模擬中具有廣泛的應用[12-15]. CALPUFF模型主要由三部分組成:診斷風場模塊(CALMET),高斯煙團擴散模塊(CALPUFF)以及后處理模塊(CALPOST)[16].在本次研究中,CALPUFF模型的模擬區(qū)域包括京津冀地區(qū)的整個區(qū)域.CALMET/CALPUFF模塊中11個垂直層高度分別設置為0、20、40、80、160、320、640、1200、2000、3000、4000m,網格設置為96′123,東西方向距離為576km,南北方向距離738km.CALMET模塊中的地形高度和土地利用數(shù)據(jù)分別采用SRTM3生成的地形數(shù)據(jù)和美國地質勘探局全球數(shù)據(jù)庫的土地利用數(shù)據(jù)[17-18].本文研究的污染物為SO2、NO、PM2.5以及PM10,基于這4種污染物的化學性質,在CALPUFF模塊中的化學轉換采用MESOPUFFII化學機制[16,19].在該機制中模擬的污染物為SO2、NO、PM2.5、PM10、SO42-、NO3-和HNO3.基于CALPUFF模塊模擬出的污染物濃度場文件,利用CALPOST模塊得到了SO2、NO、PM2.5以及PM10的濃度數(shù)據(jù).

    1.3 城市源匯分區(qū)方法

    對于城市源匯分區(qū)的方法,已經有研究者提出了科學的分類方法.本文應用王燕麗等[20]的研究方法,依據(jù)污染物傳輸矩陣計算傳輸量,并依據(jù)污染物傳輸量對城市群的源匯分區(qū)進行定義:將凈傳輸量(net)定義傳輸輸入量(in)與傳輸輸出量(out)之差.當凈傳輸量大于0時,則該城市受傳輸輸入大于其作為排放源的輸出,稱為凈輸入城市,即為匯;當凈傳輸量小于0時,則該城市作為排放源的輸出大于其受傳輸輸入,稱為凈輸出城市,即為源.公式表示如下:

    net=in-out(1)

    本文依據(jù)此方法,以區(qū)域間污染物的傳輸關系作為邊,以污染物的凈傳輸量作為邊權,應用網絡分析軟件Gephi構造污染物區(qū)域的傳輸網絡,繪制出兩種排放清單下污染物傳輸方向,并對城市源匯分區(qū)的比較分析.網絡分析軟件Gephi是一種交互的可視化平臺,它適用于各種網絡、復雜系統(tǒng)、動態(tài)和層次圖的可視化分析.Gephi主要應用在探索性數(shù)據(jù)分析中,對數(shù)據(jù)的關聯(lián)關系進行實時、直觀化分析[21-22].本文主要應用Gephi軟件來實現(xiàn)對城市間污染物傳輸關系的可視化.

    2 兩種排放清單的比較

    京津冀內有工業(yè)、居民、交通等多種排放源,其中工業(yè)源(包括電廠)占排放總量的絕大部分.根據(jù)中國環(huán)境統(tǒng)計公報,二氧化硫排放總量為2176萬t,其中90.28%來自工業(yè)源.氮氧化物排放總量2337.8萬t,其中70.93%來自工業(yè)源.工業(yè)廢氣煙塵排放量1029.3萬t,其中83.4%來自工業(yè)源[23].因此,我們選擇工業(yè)源進行兩種清單的比較分析,探討不同排放清單對大氣污染物傳輸規(guī)律造成的影響.

    本文所選擇的環(huán)統(tǒng)排放清單來自于2012年環(huán)境統(tǒng)計的12929個工業(yè)點源,主要是從各級環(huán)保局統(tǒng)計收集而來.各點源的空間位置為京津冀城市群中工業(yè)企業(yè)和電廠的實際分布(如圖2(b)所示).另一種是清華大學發(fā)布的2012年MEIC排放清單(工業(yè)源和電力源)[7],排放空間分辨率為0.25′0.25°,它的空間分布較為規(guī)則(如圖2(a)所示).

    圖2 兩種清單下排放源的空間位置分布((a)為MEIC排放清單,(b)為環(huán)統(tǒng)排放清單)

    圖3 兩組排放清單中污染物排放總量對比

    差值指環(huán)統(tǒng)排放清單中各城市的排放量減去MEIC排放清單中各城市的排放量

    從圖3中看出,環(huán)統(tǒng)排放清單中大部分城市的排放量要比MEIC排放清單大8000~60000t左右.尤其是唐山,環(huán)統(tǒng)排放清單的排放量要比MEIC排放清單高出兩倍以上.除此之外,MEIC排放清單中某些城市(北京、滄州、衡水等)的排放量要比環(huán)統(tǒng)排放清單大一萬t左右.從圖3中也可以看出,無論是MEIC排放清單還是環(huán)統(tǒng)排放清單中,4種污染物中唐山均位于13座城市的首位,是污染最為嚴重的城市.衡水、廊坊以及秦皇島是污染排放量較低的城市.

    3 模擬結果比較

    3.1 模型驗證

    2013年前,我國公布的空氣質量數(shù)據(jù)僅包含空氣污染指數(shù)API(Air Pollution Index)、首要污染物、空氣質量級別,缺乏4種污染物的監(jiān)測濃度數(shù)據(jù),無法直接將模擬濃度和監(jiān)測濃度數(shù)據(jù)比對來驗證模型的準確性.我們通過計算,利用API指數(shù)和首要污染物種類來計算得到首要污染物的濃度[24],計算公式如下:

    (其中:CCC1) (2)

    I:第種污染物的污染分指數(shù);C:第種污染物的濃度;I:第i種污染物轉折點的污染分指數(shù);I1:第種污染物1轉折點的污染分指數(shù);C:第種污染物轉折點上的濃度;C1:第種污染物1轉折點上的濃度.各轉折點值參考表1.

    表1 空氣污染指數(shù)API各污染物濃度限值(μg/m3)

    本研究所使用的空氣質量監(jiān)測數(shù)據(jù)主要來自中國生態(tài)環(huán)境保護部數(shù)據(jù)中心(http://datacenter.mep. gov.cn),包含2012年12月京津冀地區(qū)7個重點城市(北京、天津、石家莊、唐山、秦皇島、邯鄲以及保定)的API等污染指標的日監(jiān)測數(shù)據(jù).2012年12月各城市首要污染物為PM10的天數(shù)達到了21d以上,SO2和NO2作為首要污染物僅有1~9d,因而選取PM10來進行模擬濃度與監(jiān)測濃度的對比,一方面來驗證模型的準確性,另一方面來比較哪種排放清單下的模擬濃度與監(jiān)測濃度的準確性更高.

    通過將模擬濃度與同期監(jiān)測濃度進行對比(如圖4所示),可以看出環(huán)統(tǒng)排放清單下污染物的模擬濃度與監(jiān)測濃度的擬合效果更好,其2達到了0.78.而MEIC排放清單的擬合效果稍差,其2為0.60.但是從總體上來看,兩種排放清單下PM10的模擬濃度與監(jiān)測濃度的擬合效果都基本滿足空氣質量模型驗證的相關要求[20,25].

    圖4 模擬值與監(jiān)測值對比((a)為MEIC排放清單,(b)為環(huán)統(tǒng)排放清單)

    3.2 濃度空間分布對比分析

    從空間分布的角度來看(圖5),兩種排放清單下模擬的污染物濃度均呈現(xiàn)北部以唐山中心,南部以石家莊-邯鄲為中心的分布特征,均是由兩個濃度最高的中心向外逐漸降低.在4種污染物的濃度空間分布中,環(huán)統(tǒng)排放清單下模擬的污染物濃度高值區(qū)范圍更大,尤其是針對PM10的模擬結果.從各城市的污染物平均濃度來說(表2),基于環(huán)統(tǒng)排放清單模擬的大部分城市4種污染物的濃度均要高于MEIC排放清單.僅有部分城市(北京等)MEIC排放清單的模擬濃度要高于環(huán)統(tǒng)排放清單.

    圖5 兩種排放清單中各污染物濃度的空間分布((a)~(d)為MEIC排放清單,(e)~(h)為環(huán)統(tǒng)排放清單)

    造成空間分布差異性的原因主要有以下兩點:(1)排放位置的空間分布不同,MEIC排放清單的排放空間分辨率為0.25′0.25°,它的空間分布較規(guī)則,污染物點源排放空間精確位置被放大到了一個網格中.而環(huán)統(tǒng)排放清單的點源空間位置為京津冀城市群中工業(yè)企業(yè)的實際經緯度位置(如圖2所示).(2)排放量不同,通過上一節(jié)的研究,我們發(fā)現(xiàn)環(huán)統(tǒng)排放清單中點源匯總的大部分城市的排放量要高于MEIC排放清單,并且MEIC排放清單中各城市污染物排放量的差異要小于環(huán)統(tǒng)排放清單.

    表2 兩種排放清單下各城市的平均濃度(μg/m3)

    注:表中文字加粗的部分是基于MEIC排放清單模擬的污染物濃度高于環(huán)統(tǒng)排放清單的污染物濃度.

    3.3 污染物傳輸方向的差異分析

    污染物傳輸方向對于污染責任界定具有很大的影響.排放源城市是以污染物凈輸出為主的城市,即為源,而與此對應主要受其他城市影響的城市則為匯.通過對污染物的傳輸方向進行對比分析,可以有效的發(fā)現(xiàn)使用不同排放清單對污染責任方界定造成的影響.根據(jù)式(1),在MEIC排放清單中,邯鄲、石家莊、滄州和唐山是主要的源,以向周邊城市對外輸出為主.其凈傳輸量分別占PM2.5城市間凈傳輸總量的13%、16%、16%和26%.在環(huán)統(tǒng)排放清單中,唐山、石家莊和邯鄲也同樣是13座城市中典型的“源”,向周邊城市凈輸出,其凈傳輸量分別占PM2.5城市間凈傳輸總量的50%、16%和17%.兩個排放清單下計算得到源和匯的結果與王燕麗等[17]研究PM2.5在京津冀城市群傳輸?shù)慕Y果類似.

    雖然在污染物的傳輸網絡圖中,兩種排放清單中城市間污染物傳輸方向一致的凈傳輸量占污染物凈傳輸總量的80%左右,但仍有10%~19%傳輸流方向不一致,從而造成某些城市源和匯的責任認定發(fā)生了變化.我們利用Gephi軟件將兩組排放清單中污染物傳輸方向不一致的傳輸流表示出來,如圖6所示.盡管傳輸方向不一致的凈傳輸量占所有凈傳輸量的比例不大,但仍對于污染責任界定具有重要的影響.

    由圖6可見,對于4種污染物,在兩種排放清單下滄州與周邊城市的凈傳輸方向均存在較大差異.在MEIC排放清單中,滄州以向外凈傳輸為主,即為源;而在環(huán)統(tǒng)排放清單中,滄州則變成了凈輸入城市,即為匯.除此之外,衡水、保定、承德、天津、唐山以及秦皇島這幾座城市在兩組排放清單中某些污染物的凈傳輸方向也不同,例如對于NO而言(如圖6(a)所示),衡水在MEIC排放清單中表現(xiàn)為“源”,而在環(huán)統(tǒng)排放清單中表現(xiàn)為“匯”.衡水—石家莊、邢臺、邯鄲,保定—秦皇島和邯鄲,天津—唐山,廊坊—秦皇島,這幾個城市間的傳輸關系在兩種排放清單下的計算結果完全相反.根據(jù)MEIC排放清單模擬濃度結果顯示,這幾條凈傳輸方向是從“源”到“匯”,而根據(jù)環(huán)統(tǒng)排放清單模擬結果,凈傳輸方向則相反;對于SO2來說(如圖6(b)所示),根據(jù)MEIC排放清單模擬的結果,承德相對于河北省南部城市是凈輸出城市,但根據(jù)環(huán)統(tǒng)排放清單模擬的結果,承德則變成了凈輸入城市;對于PM2.5和PM10來說(如圖6(c)~(d)所示),傳輸方向差異比較類似,秦皇島和天津都是值得注意的城市,根據(jù)MEIC排放清單模擬的結果,秦皇島主要受到北京、保定和邢臺的凈輸入影響,相對這三個城市表現(xiàn)為“匯”,而根據(jù)環(huán)統(tǒng)排放清單,秦皇島則表現(xiàn)為“源”,向北京、邢臺凈輸出;在MEIC排放清單模擬結果中,天津相對于秦皇島、廊坊和邢臺是“源”的角色,但是在環(huán)統(tǒng)排放清單模擬結果中,天津則變成了“匯”.

    造成污染物傳輸方向不一致的原因可能主要有以下兩點:(1)排放量是影響區(qū)域傳輸?shù)闹饕蛩?對于環(huán)統(tǒng)排放清單來說,真實排放源呈現(xiàn)出集中程度高的分布.因而,污染物的傳輸方向均是由排放量高的城市(石家莊等城市)向排放量低的城市(如滄州、衡水、承德等城市)傳輸.例如滄州、衡水、承德屬于排放量低的城市,其排放量均在500~5000t/月,其排放量僅占京津冀地區(qū)總排放量的8%左右,而石家莊、邯鄲這些排放量高的城市,其排放量占總排放量的53%,因而在環(huán)統(tǒng)排放清單中,排放量高的城市(石家莊等城市)向滄州等城市傳輸.(2)對于MEIC排放清單來說,排放源被網格化,集中程度被削弱.污染物的傳輸方向受到地形、風向等因素的影響增強[26-27],其傳輸方向并非均是由排放量高的城市向排放量低的城市傳輸.例如石家莊、邯鄲這些排放量高的城市也只占總排放量的32%,而排放量低的滄州、衡水、承德的排放量均在3000~7000t/月,其排放量也能占京津冀城市群總排放量的16%左右,受到地形、風向等因素的影響造成傳輸方向不一定是從高排放城市向低排放城市傳輸.

    圖6 兩種排放清單下四種污染物(NOx、SO2、PM2.5以及PM10)的方向不同的傳輸流((1)為MEIC排放清單的結果,(2)為環(huán)統(tǒng)排放清單的結果)

    邊和箭頭的不同顏色代表著不同的城市,箭頭方向為污染物的凈輸出方向

    4 討論

    在分析兩種排放清單對傳輸規(guī)律造成的影響基礎上,從排放源的空間分布、時間尺度和可獲性以及行業(yè)和污染物這三個方面對兩種排放清單的優(yōu)缺點進行了比較分析,并且基于我們的研究結果對排放清單的選取進行了推薦.從排放源的空間分布來說,環(huán)統(tǒng)排放清單中各污染點源的空間位置為各企業(yè)的實際經緯度位置;而MEIC排放清單是網格化的排放清單,其空間分辨率為0.25′0.25°,網格化了企業(yè)的實際位置.從時間尺度和可獲性來說,環(huán)統(tǒng)排放清單是一年一更新,但是由于數(shù)據(jù)保密的原因,僅能從各級環(huán)保局申請獲得.而MEIC清單是每兩年更新一次,在清華大學發(fā)布的中國多尺度排放清單的網站(http://www.meicmodel.org)就可以下載,較易獲得.從行業(yè)和污染物來講,環(huán)統(tǒng)排放清單包含30個行業(yè)的三種污染物(NO、SO2和煙粉塵)的排放量.而MEIC清單包含電力、工業(yè)、交通、居民以及農業(yè)五個類別十種大氣污染物(CO、CO2、NH3、NO、OC、PM10、PM2.5、SO2、VOC以及BC),其類別數(shù)雖少,但污染物種類更加豐富.因而,更加適合化學機制復雜的空氣質量模型.本研究僅研究工業(yè)源(包括電廠)污染物的傳輸規(guī)律,兩種清單模擬的濃度空間分布趨勢是類似的,但是依據(jù)本研究所得到的結果顯示:MEIC排放清單下污染物的模擬濃度與實際監(jiān)測濃度的2為0.60;而環(huán)統(tǒng)排放清單下污染物的模擬濃度與實際監(jiān)測濃度的擬合效果反而更好,其R2可達到0.78.因而,在僅研究京津冀區(qū)域點源時選擇環(huán)統(tǒng)排放清單來進行污染物的模擬,更貼近實際情況.

    本研究分析不同排放清單對污染物濃度的模擬、污染物傳輸方向等造成的影響,發(fā)現(xiàn)了排放清單不同會產生不同的結論.其中,污染物的傳輸方向和排放源城市對于污染責任的界定具有重要作用.排放清單選擇不同,對于城市尺度的污染責任的界定可能會產生相反的結果.因而,將空氣質量模擬應用于環(huán)境管理中時,要注意排放清單的選擇.本文選取的兩種排放清單在完全相同的氣象條件下雖然模擬結果總體有較大的一致性,但是對城市級別的精細化管理來說,還是在局部城市間傳輸規(guī)律上具有差異,這種差異可能會影響到大氣聯(lián)防聯(lián)控中各城市源匯責任的認定.在未來的研究中,將選取目前國家正在編制的京津冀及周邊地區(qū)的源清單進行多種排放清單的比較分析.通過比較分析選擇出適合空氣質量模型模擬區(qū)域傳輸?shù)呐欧徘鍐?從而為京津冀地區(qū)環(huán)境管理提供科技支撐.

    5 結語

    從兩種清單的排放量來看,環(huán)統(tǒng)排放清單中點源匯總的大部分城市的排放量要比MEIC排放清單大8000~60000t左右.尤其是唐山,環(huán)統(tǒng)排放清單的排放量要比MEIC排放清單高出兩倍以上.從模擬濃度來看,兩種排放清單下污染物濃度均呈現(xiàn)中心分布特征,均是由兩個濃度最高的中心河北省東北部(唐山市及周邊地區(qū))以及河北省西南部(石家莊、邢臺和邯鄲)向外逐漸降低.在四種污染物的濃度空間分布中,環(huán)統(tǒng)排放清單下污染物由兩個濃度高值中心向周邊城市擴散范圍更大.基于環(huán)統(tǒng)排放清單模擬的大部分城市四種污染物的濃度均要高于MEIC排放清單模擬的結果,僅有部分城市(北京等)MEIC排放清單的模擬濃度要高于環(huán)統(tǒng)排放清單.從污染物的傳輸方向來看,雖然兩種排放清單中城市間污染物傳輸方向一致的凈傳輸量占污染物凈傳輸總量的80%左右,但仍有10%~19%.的傳輸流方向不一致.最典型的是滄州市,全部四種污染物在兩類清單下模擬的傳輸方向都不一致.除此之外,衡水、保定、承德、天津、唐山以及秦皇島這幾座城市在某類污染物的凈傳輸方向也不同.例如對于NO而言,衡水在MEIC排放清單中表現(xiàn)為“源”,而在環(huán)統(tǒng)排放清單中表現(xiàn)為“匯”.造成兩種排放清單污染物傳輸方向不一致的原因主要是排放量的高低以及排放源位置空間分布.對于環(huán)統(tǒng)排放清單來說,其排放源呈現(xiàn)出集中程度高的分布,城市排放量差異顯著.而對于MEIC排放清單來說,排放源被網格化,集中程度被削弱,表現(xiàn)為污染物的傳輸方向受到地形、風向等因素的影響增強.

    [1] Yan D, Lei Y, Shi Y, et al. Evolution of the spatiotemporal pattern of PM2.5concentrations in China-A case study from the Beijing- Tianjin-Hebei region [J]. Atmospheric Environment, 2018,183:225- 233.

    [2] Wu D, Xu Y, Zhang S. Will joint regional air pollution control be more cost-effective? An empirical study of China's Beijing-Tianjin-Hebei region [J]. Journal of Environmental Management, 2015,149:27-36.

    [3] Itahashi S, Uno I, Osada K, et al. Nitrate transboundary heavy pollution over East Asia in winter [J]. Atmospheric Chemistry and Physics, 2017,17(6):3823-3843.

    [4] Lin J, Pan D, Davis S J, et al. China's international trade and air pollution in the United States [J]. Proceedings of the National Academy of Sciences, 2014,111(5):1736-1741.

    [5] 薛文博,付 飛,王金南,等.中國PM2.5跨區(qū)域傳輸特征數(shù)值模擬研究 [J]. 中國環(huán)境科學, 2014,(6):1361-1368. Xue W B, Fu F, Wang J N, et al. Numerical study on the characteristics of regional transport of PM2.5in China [J]. China Environmental Science, 2014,(6):1361-1368.

    [6] Wang Y, Bao S, Wang S, et al. Local and regional contributions to fine particulate matter in Beijing during heavy haze episodes [J]. Science of The Total Environment, 2017,580:283-296.

    [7] 胡亞男,馬曉燕,沙 桐,等.不同排放源對華東地區(qū)PM2.5影響的數(shù)值模擬 [J]. 中國環(huán)境科學, 2018,38(5):1616-1628. Hu Y N, Ma X Y, Sha T, et al. Impact of different emission sources on PM2.5over East China based on numerical study [J]. China Environmental Science, 2018,38(5):1616-1628.

    [8] Li M, Zhang Q, Kurokawa J, et al. MIX: A mosaic Asian anthropogenic emission inventory for the MICS-Asia and the HTAP projects [J]. Atmospheric Chemistry and Physics Discussions, 2015, 15(23):34813-34869.

    [9] Scire J S, Robe F R, Fernau M E, et al. A user's guide for the CALMET meteorological model (Version 5) [J]. Earth Tech, 1998.

    [10] EPA. User's guide for the industrial source complex (ISC2) dispersion models. volume 2: description of model algorithms [Z]. U.S. Environmental Protection Agency Office of Air Quality Planning & Standards Research Triangle, 1992.

    [11] EPA. User's guide for the industrial source complex (ISC2) dispersion models. volume 1: Description of model algorithms [Z]. U.S. Environmental Protection Agency Office of Air Quality Planning & Standards Research Triangle, 1992.

    [12] Visscher A D. Air dispersion modeling: Foundations and applications [M]. A John Wiley & Sons, Inc., Publication, 2014:149-163.

    [13] Abdul-Wahab S, Sappurd A, Al-Damkhi A. Application of California Puff (CALPUFF) model: A case study for Oman [J]. Clean Technologies and Environmental Policy, 2011,13(1):177-189.

    [14] Dresser A L, Huizer R D. CALPUFF and AERMOD model validation study in the near field: Martins Creek revisited [J]. J. Air Waste Manag Assoc, 2011,61(6):647-659.

    [15] Wu H, Zhang Y, Yu Q, et al. Application of an integrated Weather Research and Forecasting (WRF)/CALPUFF modeling tool for source apportionment of atmospheric pollutants for air quality management: A case study in the urban area of Benxi, China [J]. 2018,68(4):347- 368.

    [16] Scire J S, Strimaitis D G, Yamartino R J. A user's guide for the CALPUFF dispersion model [J]. Journal of Non-Newtonian Fluid Mechanics, 2000,223:257.

    [17] 伯 鑫,吳忠祥,王 剛,等. CALPUFF模式的標準化應用技術研究 [J]. 環(huán)境科學與技術, 2014,(S2):530-534. Bo X, Wu Z X, Wang G, et al. Study on the standardized application of CALPUFF model [J]. Environmental Science& Technology, 2014,(S2): 530-534.

    [18] 伯 鑫,王 剛,溫 柔,等.京津冀地區(qū)火電企業(yè)的大氣污染影響 [J]. 中國環(huán)境科學, 2015,35(2):364-373. Bo X, Wang G, Wen R, et al. Air pollution effect of the thermal power plants in Beijing-Tianjin-Hebei region [J]. China Environmental Science, 2015,35(2):364-373.

    [19] Jin X, Fiore A M, Murray L T, et al. Evaluating a space-based indicator of surface ozone-NO-VOC sensitivity over midlatitude source regions and application to decadal trends [J]. Journal of Geophysical Research: Atmospheres, 2017,122(19):10,410-439,461.

    [20] 王燕麗,薛文博,雷 宇,等.京津冀區(qū)域PM2.5污染相互輸送特征 [J]. 環(huán)境科學, 2017,38(12):4897-4904. Wang Y L, Xue W B, Lei Y, et al. Regional transport matrix study of PM2.5in Jingjinji region, 2015 [J]. Environmental Science, 2017, 38(12):4897-4904.

    [21] Cherven K. Network graph analysis and visualization with Gephi [J]. 2013.

    [22] 關迎暉,向 勇,陳 康.基于Gephi的可視分析方法研究與應用 [J]. 電信科學, 2013,29(S1):112-119. Guan Y H, Xiang Y, Chen K. Research and application of visual analysis method based on Gephi [J]. Telecommunications Science, 2013,29(S1):112-119.

    [23] 中華人民共和國環(huán)境保護部.2012中國環(huán)境狀況公報 [R]. 2013. Ministry of Ecology and Environmental Protection of China. Report on the State of the Environment in China: 2012 [R]. 2013.

    [24] 李 薛,龔紹琦,付東洋,等.上海市PM10濃度四季遙感模型研究 [J]. 大氣科學學報, 2016,39(1):126-132. Li X, Gong S Q, Fu D Y, et al. Study on the four—season remote sensing models of PM10concentration in Shanghai [J]. Transactions of Atmospheric Sciences, 2016,39(1):126-132.

    [25] EPA. Clean Air Interstate Rule [Z]. 2013.

    [26] Zhu W, Xu X, Zheng J, et al. The characteristics of abnormal wintertime pollution events in the Jing-Jin-Ji region and its relationships with meteorological factors [J]. Science of The Total Environment, 2018,626:887-898.

    [27] 劉咸德,李 軍,趙 越,等.北京地區(qū)大氣顆粒物污染的風向因素研究 [J]. 中國環(huán)境科學, 2010,30(1):1-6. Liu X D, Li J, Zhao Y, et al. Impact of wind direction on atmospheric particulate matter concentrations in Beijing [J]. China Environmental Science, 2010,30(1):1-6.

    Atmospheric transmission rule on air pollution in Beijing-Tianjin-Hebei urban agglomeration: A comparative analysis of two emission inventories.

    WANG Yuan1, LI Yue1, QIAO Zhi1*, LU Ya-ling1,2

    (1.School of Environmental Science and Engineering, Tianjin University, Tianjin 300350, China;2.Chinese Academy for Environmental Planning, Beijing 100012, China)., 2019,39(11):4561~4569

    In this study, the Weather Research and Forecasting (WRF) model coupled with the California Puff (CALPUFF) air quality model was applied to study the effects of different emission inventories on the regional contribution of atmospheric transmission in the Beijing-Tianjin-Hebei urban agglomeration under heavy pollution weather conditions. The results of four pollutants was compared (NO, SO2, PM2.5and PM10) from two sets of typical emission inventories (the Environmental Statistics Emission Inventory from government and the Multi-resolution Emission Inventory for China (MEIC) from Qinghua University). From the perspective of the spatial distribution of the simulated concentrations, the results based on the two sets of emission inventories in December 2012 were similar. The concentration of pollutants under two emission inventories showed a central distribution, with Tangshan as the centre in north and Shijiazhuang-Handan as centre in the south. However, there were still some significant differences in the transmission roles of some cities based on different emission inventory inputs. For example, for the four pollutants, the transmission direction between Cangzhou and its surrounding cities were completely opposite under the two emission inventories. Cangzhou was the relativer source city among the thirteen cities based on the MEIC emissions inventories. However, based on Environmental Statistics Emission Inventory, Cangzhou was a receptor city. These conclusions will affect the identification of the source and receptors cities in air pollution joint prevention and control. In environmental management, we should pay attention to the verification and comparison of different emissions inventories.

    CALPUFF;Beijing-Tianjin-Hebei urban agglomeration;transmission rule;transfer network of pollutants

    X32,X82

    A

    1000-6923(2019)11-4561-09

    王 媛(1977-),女,天津人,教授,博士,主要研究方向資源與環(huán)境管理.發(fā)表論文20余篇.

    2019-03-31

    國家重點研發(fā)計劃(2018YFC0213600);國家自然科學基金(41871211,41571522);天津市自然科學基金(16JCQNJC08900)

    * 責任作者, 講師, qiaozhi@tju.edu.cn

    猜你喜歡
    排放量京津冀污染物
    菌株出馬讓畜禽污染物變廢為寶
    天然氣輸配系統(tǒng)甲烷排放量化方法
    煤氣與熱力(2021年6期)2021-07-28 07:21:40
    環(huán)境科學研究(2021年6期)2021-06-23 02:39:54
    環(huán)境科學研究(2021年4期)2021-04-25 02:42:02
    你能找出污染物嗎?
    黑龍江省碳排放量影響因素研究
    京津冀大聯(lián)合向縱深突破
    京津冀一化
    養(yǎng)老“京津冀一體化”謹慎樂觀看
    全國機動車污染物排放量
    ——《2013年中國機動車污染防治年報》(第Ⅱ部分)
    下体分泌物呈黄色| 成人18禁高潮啪啪吃奶动态图 | 永久免费av网站大全| 欧美老熟妇乱子伦牲交| 久久亚洲国产成人精品v| 美女中出高潮动态图| 久久精品国产亚洲av天美| 丝袜在线中文字幕| 搡老乐熟女国产| av一本久久久久| 久久久久久伊人网av| 成人国产av品久久久| 久久女婷五月综合色啪小说| 看非洲黑人一级黄片| 国产高清国产精品国产三级| 日韩一区二区三区影片| av专区在线播放| 国产白丝娇喘喷水9色精品| 蜜臀久久99精品久久宅男| 久久99精品国语久久久| 国产伦精品一区二区三区视频9| 亚洲美女视频黄频| 午夜久久久在线观看| 久久ye,这里只有精品| 日韩成人av中文字幕在线观看| 少妇猛男粗大的猛烈进出视频| 精品少妇内射三级| 麻豆精品久久久久久蜜桃| 99热这里只有精品一区| 国产成人午夜福利电影在线观看| 美女脱内裤让男人舔精品视频| 久久热精品热| 日本wwww免费看| 大片电影免费在线观看免费| 国产免费视频播放在线视频| 日韩 亚洲 欧美在线| 晚上一个人看的免费电影| av在线老鸭窝| 国产在线视频一区二区| 这个男人来自地球电影免费观看 | 国产欧美日韩精品一区二区| 日本黄色日本黄色录像| 两个人的视频大全免费| 极品人妻少妇av视频| 有码 亚洲区| 久久午夜综合久久蜜桃| 国内精品宾馆在线| 亚洲国产成人一精品久久久| 国产欧美日韩精品一区二区| 成人亚洲欧美一区二区av| 男女免费视频国产| 精品国产一区二区久久| 中国国产av一级| av国产久精品久网站免费入址| 亚洲欧美一区二区三区国产| 亚洲久久久国产精品| 精品亚洲成a人片在线观看| 久久久欧美国产精品| 人人妻人人澡人人爽人人夜夜| 亚洲欧美精品自产自拍| 日韩精品有码人妻一区| 青春草亚洲视频在线观看| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 秋霞在线观看毛片| 极品少妇高潮喷水抽搐| 久久久久国产网址| 久久99精品国语久久久| 亚洲欧洲日产国产| 一级二级三级毛片免费看| 精品久久久精品久久久| 亚洲av不卡在线观看| 亚洲av成人精品一区久久| 高清视频免费观看一区二区| 天天躁夜夜躁狠狠久久av| 3wmmmm亚洲av在线观看| 午夜免费鲁丝| 97精品久久久久久久久久精品| 欧美bdsm另类| 国产成人91sexporn| √禁漫天堂资源中文www| 成人影院久久| 一级毛片我不卡| 国产 一区精品| 亚洲av电影在线观看一区二区三区| 午夜福利视频精品| 国产精品熟女久久久久浪| 99热这里只有是精品在线观看| 夫妻性生交免费视频一级片| 国产极品天堂在线| 搡老乐熟女国产| 免费看av在线观看网站| 99精国产麻豆久久婷婷| 一边亲一边摸免费视频| 性色av一级| 久热这里只有精品99| 国产永久视频网站| 日本av手机在线免费观看| 亚洲国产欧美日韩在线播放 | 观看美女的网站| 嘟嘟电影网在线观看| 精品少妇黑人巨大在线播放| 美女国产视频在线观看| 国产 精品1| 少妇被粗大的猛进出69影院 | 欧美+日韩+精品| 亚洲天堂av无毛| 亚洲色图综合在线观看| 国产日韩欧美在线精品| 蜜臀久久99精品久久宅男| 久久人人爽av亚洲精品天堂| 一级av片app| 免费大片黄手机在线观看| 欧美bdsm另类| 国产精品久久久久久精品电影小说| 精品一品国产午夜福利视频| 国产 一区精品| 91精品国产国语对白视频| 久久久久久久久久久丰满| 大陆偷拍与自拍| 爱豆传媒免费全集在线观看| 久久精品夜色国产| 女性生殖器流出的白浆| 中国三级夫妇交换| 色婷婷久久久亚洲欧美| 人人妻人人看人人澡| av在线观看视频网站免费| 久久99热这里只频精品6学生| 高清毛片免费看| 精品国产露脸久久av麻豆| 五月玫瑰六月丁香| 另类精品久久| 日本免费在线观看一区| 丰满迷人的少妇在线观看| 欧美xxxx性猛交bbbb| 亚洲久久久国产精品| 国产极品粉嫩免费观看在线 | 欧美精品亚洲一区二区| 黄色配什么色好看| 99热网站在线观看| 在线免费观看不下载黄p国产| 亚洲性久久影院| 美女主播在线视频| 国精品久久久久久国模美| 乱人伦中国视频| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 亚洲av成人精品一区久久| 亚洲精品乱码久久久久久按摩| 亚洲av电影在线观看一区二区三区| 天堂俺去俺来也www色官网| 久久久久国产网址| 有码 亚洲区| 亚洲精品456在线播放app| 色94色欧美一区二区| 中文精品一卡2卡3卡4更新| 少妇人妻精品综合一区二区| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 日韩亚洲欧美综合| 精品少妇黑人巨大在线播放| 一边亲一边摸免费视频| 亚洲精品乱码久久久v下载方式| 国产免费福利视频在线观看| 久久99热这里只频精品6学生| 乱码一卡2卡4卡精品| 中文天堂在线官网| 91精品伊人久久大香线蕉| 国产日韩欧美亚洲二区| 日本与韩国留学比较| 欧美日韩视频高清一区二区三区二| 97在线人人人人妻| 国产精品偷伦视频观看了| 国产免费一区二区三区四区乱码| 免费人成在线观看视频色| 亚洲天堂av无毛| 一本大道久久a久久精品| 国产伦精品一区二区三区四那| 日本午夜av视频| 人妻人人澡人人爽人人| 国产成人午夜福利电影在线观看| 18禁裸乳无遮挡动漫免费视频| 国产美女午夜福利| 久久国产精品大桥未久av | 国产在线免费精品| 国国产精品蜜臀av免费| 乱码一卡2卡4卡精品| 久久久午夜欧美精品| av国产精品久久久久影院| 久久久久久久亚洲中文字幕| 欧美日韩精品成人综合77777| 春色校园在线视频观看| 国产在线免费精品| 黑人猛操日本美女一级片| 女的被弄到高潮叫床怎么办| 久久国产精品大桥未久av | 伦理电影大哥的女人| 亚洲内射少妇av| 美女视频免费永久观看网站| 大陆偷拍与自拍| 精品久久久久久电影网| 国产 精品1| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 国产黄色视频一区二区在线观看| 国产探花极品一区二区| 日本午夜av视频| 中文字幕精品免费在线观看视频 | 乱系列少妇在线播放| 免费看光身美女| 国产成人91sexporn| 亚洲精品成人av观看孕妇| 亚洲经典国产精华液单| 在线观看美女被高潮喷水网站| 97在线人人人人妻| 欧美日韩精品成人综合77777| 色婷婷久久久亚洲欧美| 免费看光身美女| 久久人人爽人人爽人人片va| 极品少妇高潮喷水抽搐| 亚洲熟女精品中文字幕| 美女cb高潮喷水在线观看| 精品一品国产午夜福利视频| av免费在线看不卡| 在线观看人妻少妇| 欧美日韩视频精品一区| 欧美日韩av久久| 亚洲国产日韩一区二区| 亚洲欧洲日产国产| 在线观看人妻少妇| 观看av在线不卡| 大香蕉97超碰在线| 久久鲁丝午夜福利片| 自拍欧美九色日韩亚洲蝌蚪91 | 一区二区三区精品91| 日日爽夜夜爽网站| 97超碰精品成人国产| 亚洲综合色惰| 亚洲国产欧美日韩在线播放 | 99热这里只有是精品在线观看| 在线天堂最新版资源| 久久久精品94久久精品| 中文在线观看免费www的网站| 午夜日本视频在线| 女的被弄到高潮叫床怎么办| 伦理电影大哥的女人| 亚洲精品乱久久久久久| 亚洲成人一二三区av| 大陆偷拍与自拍| 成人特级av手机在线观看| 亚洲美女搞黄在线观看| 欧美日韩在线观看h| 欧美 亚洲 国产 日韩一| 内射极品少妇av片p| 自拍欧美九色日韩亚洲蝌蚪91 | 一本—道久久a久久精品蜜桃钙片| 建设人人有责人人尽责人人享有的| 国产av码专区亚洲av| 欧美少妇被猛烈插入视频| 久久精品国产a三级三级三级| 蜜臀久久99精品久久宅男| 2021少妇久久久久久久久久久| 内射极品少妇av片p| 色哟哟·www| 99热这里只有精品一区| 日本色播在线视频| 国产成人aa在线观看| av天堂中文字幕网| 久久6这里有精品| 国产av一区二区精品久久| 精品久久久久久电影网| 熟妇人妻不卡中文字幕| 黄片无遮挡物在线观看| 日韩 亚洲 欧美在线| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦视频在线资源免费观看| 久久久精品免费免费高清| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网| 午夜福利网站1000一区二区三区| 日韩欧美精品免费久久| 不卡视频在线观看欧美| 久久影院123| 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 国产无遮挡羞羞视频在线观看| 高清黄色对白视频在线免费看 | 天天操日日干夜夜撸| 久久韩国三级中文字幕| 精品人妻熟女毛片av久久网站| 九色成人免费人妻av| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩视频精品一区| 亚洲成人一二三区av| 男女边摸边吃奶| 一本久久精品| 熟女电影av网| av国产久精品久网站免费入址| 国产一区二区三区综合在线观看 | 日韩中文字幕视频在线看片| 免费黄频网站在线观看国产| 一区二区av电影网| 18禁动态无遮挡网站| 亚洲精品,欧美精品| 国产精品99久久99久久久不卡 | 精品久久久久久久久av| 国产一区有黄有色的免费视频| 99久久精品热视频| 三上悠亚av全集在线观看 | 高清在线视频一区二区三区| 日韩精品有码人妻一区| 免费黄网站久久成人精品| 欧美日韩视频精品一区| 建设人人有责人人尽责人人享有的| 国产永久视频网站| 欧美xxxx性猛交bbbb| 成人毛片60女人毛片免费| videossex国产| 国产成人freesex在线| 精品一区在线观看国产| 日韩一区二区视频免费看| 在线观看三级黄色| 大又大粗又爽又黄少妇毛片口| 少妇人妻久久综合中文| 久久精品国产亚洲av天美| 最黄视频免费看| 色视频www国产| 久久免费观看电影| 精品人妻熟女av久视频| 欧美精品一区二区大全| 一区二区三区四区激情视频| 不卡视频在线观看欧美| 大码成人一级视频| 在线观看一区二区三区激情| 久久久久视频综合| 亚洲第一av免费看| 日韩精品有码人妻一区| 久久精品国产a三级三级三级| 少妇被粗大猛烈的视频| 中文字幕制服av| 五月伊人婷婷丁香| 国产精品久久久久久久久免| 日韩伦理黄色片| 另类精品久久| 日韩中字成人| 三上悠亚av全集在线观看 | 精品久久久噜噜| 99热网站在线观看| 又黄又爽又刺激的免费视频.| 日韩制服骚丝袜av| 精品亚洲乱码少妇综合久久| 十分钟在线观看高清视频www | 精品久久久久久久久亚洲| 建设人人有责人人尽责人人享有的| 午夜日本视频在线| 性色av一级| 久久99蜜桃精品久久| 99热6这里只有精品| 亚洲精品亚洲一区二区| 国产亚洲一区二区精品| 亚洲美女黄色视频免费看| 色视频www国产| 久久99蜜桃精品久久| 国产亚洲一区二区精品| 观看av在线不卡| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 亚洲图色成人| 人人妻人人澡人人看| 欧美日韩亚洲高清精品| 国产极品粉嫩免费观看在线 | 日韩精品免费视频一区二区三区 | 一级,二级,三级黄色视频| 亚洲国产色片| 久久久午夜欧美精品| 麻豆精品久久久久久蜜桃| 日本与韩国留学比较| 欧美丝袜亚洲另类| 哪个播放器可以免费观看大片| 欧美3d第一页| 日韩精品有码人妻一区| 久久精品国产亚洲av天美| 精品久久久久久电影网| 一级爰片在线观看| 有码 亚洲区| 欧美变态另类bdsm刘玥| 青青草视频在线视频观看| 精品人妻熟女毛片av久久网站| 日本黄色片子视频| 少妇被粗大猛烈的视频| 中文字幕人妻熟人妻熟丝袜美| 成人国产av品久久久| 久久久国产精品麻豆| 国产免费视频播放在线视频| 美女脱内裤让男人舔精品视频| 亚洲欧美清纯卡通| 国产成人精品无人区| 男人爽女人下面视频在线观看| 欧美精品人与动牲交sv欧美| 99热6这里只有精品| 国产黄频视频在线观看| 香蕉精品网在线| 观看av在线不卡| 成年av动漫网址| 搡女人真爽免费视频火全软件| 又爽又黄a免费视频| 亚洲av男天堂| 日韩精品免费视频一区二区三区 | 国产黄频视频在线观看| 欧美成人精品欧美一级黄| 国内揄拍国产精品人妻在线| 丝袜喷水一区| 精品久久久久久久久av| a级一级毛片免费在线观看| 夫妻午夜视频| 中文资源天堂在线| 久久精品久久精品一区二区三区| 成人亚洲欧美一区二区av| 2021少妇久久久久久久久久久| 男的添女的下面高潮视频| 韩国av在线不卡| 国产 一区精品| 九九爱精品视频在线观看| 国产一区二区在线观看av| 中文字幕人妻熟人妻熟丝袜美| 精品午夜福利在线看| 人妻 亚洲 视频| 亚洲精品国产成人久久av| 亚洲,欧美,日韩| 91久久精品国产一区二区成人| 人体艺术视频欧美日本| 国产精品一区www在线观看| 蜜臀久久99精品久久宅男| 欧美成人精品欧美一级黄| 菩萨蛮人人尽说江南好唐韦庄| 丁香六月天网| 一级毛片电影观看| 欧美激情国产日韩精品一区| 成人无遮挡网站| 精品熟女少妇av免费看| 国产黄色视频一区二区在线观看| 黑人猛操日本美女一级片| 亚洲第一av免费看| 噜噜噜噜噜久久久久久91| 国产成人一区二区在线| 黄色毛片三级朝国网站 | 国产中年淑女户外野战色| 午夜av观看不卡| 久久久久久久国产电影| 五月玫瑰六月丁香| 中文资源天堂在线| 男女国产视频网站| 久久毛片免费看一区二区三区| 观看美女的网站| 免费av不卡在线播放| 一本一本综合久久| 国产黄色免费在线视频| 特大巨黑吊av在线直播| 久久久久视频综合| 黄色欧美视频在线观看| 9色porny在线观看| h日本视频在线播放| 男人和女人高潮做爰伦理| 久久免费观看电影| 中文字幕制服av| 中国三级夫妇交换| a级一级毛片免费在线观看| a级毛片免费高清观看在线播放| 国产成人一区二区在线| 日韩欧美精品免费久久| 交换朋友夫妻互换小说| 妹子高潮喷水视频| 欧美日韩精品成人综合77777| 在线亚洲精品国产二区图片欧美 | 久久午夜福利片| 久久鲁丝午夜福利片| av国产精品久久久久影院| 性高湖久久久久久久久免费观看| 女性被躁到高潮视频| 最新的欧美精品一区二区| 如何舔出高潮| 最新中文字幕久久久久| 日韩视频在线欧美| 国国产精品蜜臀av免费| 国产成人freesex在线| 国产在视频线精品| 日韩av在线免费看完整版不卡| 久久精品久久久久久久性| 日韩精品有码人妻一区| 亚洲丝袜综合中文字幕| 日韩成人av中文字幕在线观看| 美女脱内裤让男人舔精品视频| 狂野欧美白嫩少妇大欣赏| 最黄视频免费看| 日本欧美视频一区| 亚洲情色 制服丝袜| 色94色欧美一区二区| 亚洲国产精品成人久久小说| 边亲边吃奶的免费视频| 国产精品蜜桃在线观看| 欧美精品国产亚洲| 国产成人freesex在线| 中文字幕精品免费在线观看视频 | 91久久精品电影网| 欧美日韩亚洲高清精品| 大片电影免费在线观看免费| 国产一区二区三区综合在线观看 | 王馨瑶露胸无遮挡在线观看| 日韩欧美 国产精品| 香蕉精品网在线| 五月玫瑰六月丁香| 不卡视频在线观看欧美| 老司机影院成人| 街头女战士在线观看网站| videos熟女内射| 丰满乱子伦码专区| 免费av不卡在线播放| 欧美xxxx性猛交bbbb| 人人妻人人澡人人看| 在线观看免费视频网站a站| 午夜免费鲁丝| 男女免费视频国产| 如何舔出高潮| 日韩av免费高清视频| 久久精品国产亚洲av涩爱| 我要看黄色一级片免费的| 国产精品偷伦视频观看了| 涩涩av久久男人的天堂| 成人漫画全彩无遮挡| 美女脱内裤让男人舔精品视频| 欧美日韩av久久| 免费高清在线观看视频在线观看| 三级国产精品片| 亚洲av中文av极速乱| 国产成人a∨麻豆精品| 成年av动漫网址| 亚洲美女视频黄频| 精品一区二区三卡| 18禁在线播放成人免费| 国产一区二区三区综合在线观看 | 高清午夜精品一区二区三区| 一级毛片我不卡| 高清午夜精品一区二区三区| 男人爽女人下面视频在线观看| 婷婷色麻豆天堂久久| 亚洲国产成人一精品久久久| 日韩制服骚丝袜av| 日本91视频免费播放| 国产成人精品婷婷| 成人午夜精彩视频在线观看| 一级a做视频免费观看| 在线观看免费视频网站a站| 欧美性感艳星| 日本免费在线观看一区| 嘟嘟电影网在线观看| 日韩av免费高清视频| 国产白丝娇喘喷水9色精品| 久久久久国产精品人妻一区二区| 国产精品秋霞免费鲁丝片| 蜜桃在线观看..| 国产欧美日韩精品一区二区| 午夜福利视频精品| 国产黄色免费在线视频| 丝袜喷水一区| 久久精品久久精品一区二区三区| av网站免费在线观看视频| av线在线观看网站| 美女xxoo啪啪120秒动态图| 成人亚洲欧美一区二区av| 秋霞在线观看毛片| 久久久久视频综合| 九草在线视频观看| 欧美日韩av久久| 日韩制服骚丝袜av| 麻豆成人av视频| 欧美 亚洲 国产 日韩一| 午夜av观看不卡| 亚洲丝袜综合中文字幕| 国产免费一区二区三区四区乱码| 国产精品偷伦视频观看了| 国产免费又黄又爽又色| 久久97久久精品| 国产成人精品无人区| 午夜精品国产一区二区电影| 五月玫瑰六月丁香| 亚洲在久久综合| 水蜜桃什么品种好| 国产精品秋霞免费鲁丝片| 蜜桃在线观看..| 最近手机中文字幕大全| 日韩中字成人| 午夜免费观看性视频| 人妻人人澡人人爽人人| 国产熟女午夜一区二区三区 | 国产女主播在线喷水免费视频网站| 高清毛片免费看| 国产精品一区二区在线不卡| 久久国产精品大桥未久av | 久久久久久伊人网av| 内射极品少妇av片p| 水蜜桃什么品种好| 日本vs欧美在线观看视频 | 国产综合精华液| 永久网站在线| 久久久久久久国产电影| 久久影院123| 纵有疾风起免费观看全集完整版| 亚洲国产精品国产精品| 偷拍熟女少妇极品色| 亚洲av欧美aⅴ国产| 久久久久久人妻| 久久午夜综合久久蜜桃| 天堂中文最新版在线下载| 女的被弄到高潮叫床怎么办| 成人影院久久|