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

    聲屏障對(duì)城市道路環(huán)境空氣質(zhì)量的影響規(guī)律研究

    2019-08-28 02:45:18秦曉春
    中國(guó)環(huán)境科學(xué) 2019年8期
    關(guān)鍵詞:風(fēng)向屏障風(fēng)速

    秦曉春,劉 睿,韓 瑩

    聲屏障對(duì)城市道路環(huán)境空氣質(zhì)量的影響規(guī)律研究

    秦曉春*,劉 睿,韓 瑩

    (北京交通大學(xué)土木建筑工程學(xué)院,北京 100044)

    為研究聲屏障對(duì)道路交通污染擴(kuò)散和空氣質(zhì)量的影響規(guī)律,基于三維計(jì)算流體動(dòng)力學(xué)理論,運(yùn)用Fluidyn-Panache仿真及CFD求解器進(jìn)行數(shù)值模型計(jì)算分析,對(duì)比有、無(wú)聲屏障情況下的風(fēng)速矢量圖與PM2.5、CO、NO3種主要污染物濃度分布圖,并對(duì)仿真區(qū)域污染物濃度進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè),得出聲屏障對(duì)城市道路環(huán)境空氣質(zhì)量影響規(guī)律.研究結(jié)果顯示,有聲屏障時(shí)道路內(nèi)側(cè)風(fēng)速由2m/s降低至0.02m/s;CO、NO濃度上升40%~50%,道路內(nèi)側(cè)PM濃度下降50%;道路外側(cè)CO、NO濃度下降20%~50%,PM濃度下降20%~38%.聲屏障的存在會(huì)阻擋風(fēng)的流動(dòng),減弱風(fēng)的強(qiáng)度,使風(fēng)向從垂直于道路上升轉(zhuǎn)變?yōu)檠刂缆贩较蛄鲃?dòng),風(fēng)向改變使污染物聚集在道路內(nèi),增加道路內(nèi)的氣體污染物濃度,減弱風(fēng)往道路外側(cè)流動(dòng),同時(shí)也降低了氣體污染物向道路外側(cè)擴(kuò)散的程度.聲屏障對(duì)PM2.5的影響顯著,設(shè)置聲屏障的道路內(nèi)側(cè)和道路外側(cè)建筑群內(nèi),PM的濃度相較于無(wú)聲屏障均得到大幅度降低.

    道路;聲屏障;環(huán)境空氣質(zhì)量;污染物擴(kuò)散

    我國(guó)城市道路在飛速建設(shè)發(fā)展的同時(shí)帶來(lái)了嚴(yán)重的環(huán)境污染問(wèn)題,車(chē)流量的增加造成道路嚴(yán)重?fù)矶?汽車(chē)排放的尾氣也給環(huán)境帶來(lái)了嚴(yán)重負(fù)荷,甚至已經(jīng)成為很多城市的首要污染物[1].機(jī)動(dòng)車(chē)排放顯著增加了PM2.5和CO、NO等污染物濃度,據(jù)調(diào)查,公路路側(cè)污染物濃度明顯高于幾百米外城市區(qū)域[2].近公路區(qū)域高水平的空氣污染還帶來(lái)了健康影響,包括出生缺陷、肺部疾病和心血管疾病[3].因此,對(duì)道路交通污染擴(kuò)散的治理和對(duì)空氣質(zhì)量的改善是十分必要的.

    國(guó)外在交通污染防治方面的研究,主要集中在汽車(chē)尾氣污染物排放對(duì)于城市道路環(huán)境的影響,以及道路擁擠引起的環(huán)境影響方面,國(guó)外在城市道路交通規(guī)劃對(duì)環(huán)境影響的研究與實(shí)踐較多,在進(jìn)行城市道路交通規(guī)劃的時(shí)候,美國(guó)更加注重土地的利用率和經(jīng)濟(jì)情況分析,對(duì)道路交通環(huán)境污染因素及城市道路環(huán)境空氣質(zhì)量的研究較少[4-6].

    路側(cè)植被屏障可以有效降低路側(cè)空氣污染水平,目前國(guó)內(nèi)外在城市道路路側(cè)植被屏障方面開(kāi)展了大量研究[7-8].通過(guò)實(shí)體構(gòu)造物減緩城市道路交通污染影響的研究尚不多見(jiàn).在使用示蹤氣體的現(xiàn)場(chǎng)研究中,發(fā)現(xiàn)路側(cè)聲屏障將順風(fēng)方向污染物濃度降低20%~50%,聲屏障的結(jié)構(gòu)改變了附近的風(fēng)動(dòng)氣流場(chǎng),將車(chē)輛排放的空氣污染物由地面帶入障礙物上方區(qū)域[9].因此,聲屏障有助于降低空氣污染物的額外分散度,使空氣污染物水平相較于附近道路開(kāi)放區(qū)域降低約15%~80%[10].

    在城市道路交通污染物控制方面,我國(guó)相對(duì)于國(guó)外研究較晚.目前我國(guó)在對(duì)于汽車(chē)尾氣擴(kuò)散方面的研究中,主要運(yùn)用的模型大部分從國(guó)外引進(jìn).將國(guó)外的模型引入到我國(guó),對(duì)該模型進(jìn)行分析,并分析其對(duì)于我國(guó)污染物擴(kuò)散情況的適用性[11].劉永紅等[12]建立了一個(gè)有關(guān)于機(jī)動(dòng)車(chē)尾氣擴(kuò)散的模型,該模型在高斯擴(kuò)散理論的基礎(chǔ)上,針對(duì)于線(xiàn)源擴(kuò)散的特點(diǎn),主要用于計(jì)算城市復(fù)雜道路、城市建筑群、高架橋上的汽車(chē)尾氣污染物濃度.張文杰等[13]人也指出了幾個(gè)在高斯模型的理論基礎(chǔ)上建立的有關(guān)污染物擴(kuò)散的模型,并且根據(jù)各自模型的特點(diǎn),分析了我國(guó)將來(lái)建立有關(guān)污染物擴(kuò)散新模型亟待解決的難題.國(guó)內(nèi)學(xué)者研究城市街道峽谷內(nèi)機(jī)動(dòng)車(chē)排放污染物擴(kuò)散規(guī)律也是在國(guó)外學(xué)者研究的基礎(chǔ)上發(fā)展起來(lái)的.劉寶章[14]基于澳門(mén)復(fù)雜街區(qū)進(jìn)行了風(fēng)洞試驗(yàn),分析了交通堵塞和交通順暢時(shí)汽車(chē)排放的污染物擴(kuò)散規(guī)律.周洪昌[15]通過(guò)風(fēng)洞實(shí)驗(yàn)分析了不同風(fēng)向來(lái)流以及不同形狀街道峽谷對(duì)街道峽谷內(nèi)機(jī)動(dòng)車(chē)排放污染物濃度場(chǎng)的影響.何澤能等[16]應(yīng)用雷諾應(yīng)力瑞流模型模擬了不同幾何結(jié)構(gòu)的城市街道峽谷的氣體流場(chǎng).

    汽車(chē)尾氣對(duì)城市大氣環(huán)境造成污染,嚴(yán)重威脅,人類(lèi)健康,其中CO、NO、HC與NO等危害最大,此類(lèi)污染物易引發(fā)生理障礙、中毒等多種疾病,威脅著人類(lèi)健康;近年來(lái),PM造成的霧霾問(wèn)題也不斷引起重視.

    盡管?chē)?guó)內(nèi)目前已經(jīng)開(kāi)展交通尾氣影響的大量研究,但在汽車(chē)尾氣對(duì)于城市道路環(huán)境空氣污染的減緩方法和措施上還較為局限.在城市交通環(huán)境空氣污染防治主要集中在政策制定、燃料改革、減少排放量、增加綠化面積或使用功能性路面等手段來(lái)解決,有關(guān)實(shí)體設(shè)施對(duì)城市道路交通污染控制的研究基本還是空白.因此本研究以北京城市道路典型聲屏障為例,運(yùn)用Fluidyn分析軟件,研究聲屏障對(duì)道路交通污染擴(kuò)散和空氣質(zhì)量的影響規(guī)律.

    (1)考慮污染物源強(qiáng)、氣象條件和道路條件等,選取典型路段,建立地形數(shù)值模型和氣象數(shù)據(jù)庫(kù).運(yùn)用Fluidyn軟件對(duì)典型路段主要交通污染物(CO、NO、PM)擴(kuò)散進(jìn)行仿真模擬,提出基于仿真分析的聲屏障對(duì)于道路交通污染擴(kuò)散和空氣質(zhì)量的影響規(guī)律.

    (2)對(duì)典型路段聲屏障區(qū)域交通污染狀況進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè),提出基于現(xiàn)場(chǎng)實(shí)驗(yàn)監(jiān)測(cè)的聲屏障對(duì)于道路交通污染擴(kuò)散和空氣質(zhì)量的影響規(guī)律.

    (3)結(jié)合仿真結(jié)果與現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù),考慮交通量、車(chē)速、車(chē)型等影響因素,確定聲屏障對(duì)城市道路環(huán)境空氣質(zhì)量影響規(guī)律.

    1 計(jì)算原理與方法

    本研究采用法國(guó)Fluidyn軟件對(duì)典型路段交通污染擴(kuò)散進(jìn)行仿真模擬分析,采用三維計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)模型計(jì)算復(fù)雜地形的風(fēng)速矢量場(chǎng)和污染物擴(kuò)散情況由于三維計(jì)算流體力學(xué)(CFD)模型能夠考慮障礙物(城市中建筑)、復(fù)雜地形、底和高風(fēng)速、近場(chǎng)和運(yùn)場(chǎng),因此在計(jì)算復(fù)雜地形的風(fēng)速矢量場(chǎng)和污染物擴(kuò)散等情況都是基于CFD模型.

    1.1 大氣邊界層(PBL)模型

    PBL位于地表附近的大氣區(qū)域.大氣與外界相互作用是以湍流運(yùn)動(dòng)體現(xiàn)的,是小尺度、次網(wǎng)格的,有必要進(jìn)行參數(shù)化建模來(lái)分析.

    PBL用于創(chuàng)立CFD求解器里需要的氣象觀(guān)測(cè)數(shù)據(jù)和邊界條件之間的接口.它由兩部分組成:(1)微氣象模型,是由觀(guān)測(cè)數(shù)據(jù)計(jì)算PBL的基本物理參數(shù)模型;(2)邊界層模型,是一個(gè)規(guī)定湍流、風(fēng)速和溫度的垂直剖面模型.

    1.2 邊界條件

    在求解時(shí),主區(qū)域邊界、地面以及障礙物上需要給定邊界條件,邊界條件分為3類(lèi),如表1所示.

    表1 邊界條件

    注:——速度矢量;y——物種m的質(zhì)量分?jǐn)?shù);——溫度;——濃度;——耗散率;——速率.

    頂部界限通常認(rèn)為是出流邊界,地域的側(cè)邊界通常憑借風(fēng)向認(rèn)為是入口或出口界限,如圖1所示.

    1.2.1 入口 在入口的邊界處,需給出溫度、物種的濃度、湍流變量、速率,壓力由區(qū)域內(nèi)部計(jì)算得到,用物種的背景濃度來(lái)確定物種的濃度.入口邊界的速度設(shè)為:

    1.2.2 出口 如果出口的邊界由一個(gè)界限所確定,但因?yàn)槭艿搅斯趯踊蛘系K物等因素的影響,那么風(fēng)仍然有能力由這里流動(dòng)進(jìn)入?yún)^(qū)域中.因此, PANACHE軟件不僅會(huì)確定出來(lái)一些變量,還會(huì)從域內(nèi)推導(dǎo)出一些變量.表2列出了幾類(lèi)出口邊界條件.

    環(huán)境的壓力一般是由出口邊界的靜壓力來(lái)確定.

    1.2.3 壁面 PANACHE軟件是運(yùn)用壁面函數(shù)來(lái)計(jì)算位于湍流邊界層中固體壁面的拖曳力,而壁面函數(shù)則是在平衡條件下,由湍流邊界層的N-S方程中計(jì)算來(lái)獲得.

    表2 出口邊界條件

    1.3 CFD求解

    PANACHE的CFD求解器在三維空間和時(shí)間上求解控制方程,而這些控制方程則可以用以下一般形式的對(duì)流擴(kuò)散方程來(lái)表示:

    式中:為待求變量(動(dòng)量、溫度、濃度等);G為的交換系數(shù)(粘度,熱擴(kuò)散系數(shù)等);為的源項(xiàng)(污染物排放等);其中,(Ⅰ)為時(shí)間微分項(xiàng),(Ⅱ)為對(duì)流項(xiàng),(Ⅲ)為耗散項(xiàng),(Ⅳ)為源項(xiàng).

    1.3.2 空間差分 空間離散的實(shí)現(xiàn)通過(guò)由單元(控制體)構(gòu)成的立體三維網(wǎng)格來(lái)完成的.控制體或整體平衡方法()用于構(gòu)造每個(gè)單元的有限差分近似(控制體),以保證差分變量的局部守恒.PANACHE求解速度向量笛卡爾坐標(biāo)系下分量的NS方程.所有變量(壓力、速度分量、溫度、濃度、湍流等)在同一控制量()內(nèi)求解得到.

    二階精度方法用于對(duì)耗散項(xiàng)的離散.通常情況下,對(duì)流項(xiàng)的離散對(duì)于求解的精度有著決定性的作用.數(shù)值格式是指對(duì)對(duì)流項(xiàng)的離散.當(dāng)使用有限體積法離散時(shí)對(duì)流項(xiàng)表示如下:

    式中:f為控制體的面;A為面積;N為面上的法向量;為流過(guò)面的質(zhì)量通量;由連續(xù)性方程求解得到.因此ff的計(jì)算決定了所采用計(jì)算格式的精度.變量f在面f附近的變化如圖2所示:

    C為上風(fēng)單元,D為下風(fēng)單元,U為遠(yuǎn)處上風(fēng)單元.可由的函數(shù)表示

    2 交通污染物擴(kuò)散仿真模型建立

    2.1 確定仿真地段

    通過(guò)google earth在北京市域內(nèi)找到有聲屏障的2~3條典型城市道路,通弄過(guò)分析,確定高粱橋斜街聲屏障路段是研究聲屏障對(duì)城市道路環(huán)境空氣質(zhì)量的典型路段,如圖3(a)所示.

    綜上所述選取高粱橋斜街路段聲屏障作為本次課題的研究對(duì)象.結(jié)合google earth與實(shí)地測(cè)量,得到建筑本身的長(zhǎng)與寬、建筑之間的距離、道路的寬度等數(shù)據(jù),如圖3(b).

    將數(shù)據(jù)標(biāo)記在地形圖上,并且建立坐標(biāo)系,得出各個(gè)建筑物的坐標(biāo),確定各建筑物的空間位置.

    圖3 典型路段信息

    2.2 建立地形模型

    圖4 地形模型

    考慮到CFD算法劃分網(wǎng)格需設(shè)定一定范圍的計(jì)算域.Panache軟件計(jì)算域范圍的設(shè)定需要通過(guò)外部文件來(lái)獲得.設(shè)定600m*600m的數(shù)值作為計(jì)算域的平面范圍,并設(shè)置計(jì)算域的高度為最高建筑高度(70m)的3倍,即為210m,從而開(kāi)展CFD模擬大氣擴(kuò)散的仿真研究.

    通過(guò)Terrain Objects繪制各個(gè)建筑物,并輸入已測(cè)得的對(duì)應(yīng)地形數(shù)據(jù).通過(guò)繪制建筑物和不規(guī)則建筑物(圖中紅圈所標(biāo)),并繪制監(jiān)測(cè)點(diǎn)(圖中藍(lán)圈所標(biāo)),以便于后期處理中記錄數(shù)據(jù),得出所有建筑和監(jiān)測(cè)點(diǎn),如圖4所示.

    2.3 設(shè)定污染物參數(shù)

    研究表明,汽車(chē)尾氣排放主要污染物一般為NO、CO和PM,監(jiān)測(cè)道路車(chē)輛以小型汽車(chē)為主,經(jīng)查閱文獻(xiàn),得到小型汽車(chē)污染物排放量,NO濃度為0.18g/km,CO濃度為1.0g/km,PM濃度為0.0045g/km[8]:

    于2018年4月10日對(duì)研究道路進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè),得出平峰11:00~13:00和高峰18:00~20:00 2個(gè)時(shí)期的4條道路車(chē)流量,如表3所示.

    表3 平、高峰時(shí)期車(chē)流量

    計(jì)算得出排放源污染物排放量[g/(km×h)]如表4所示:

    表4 排放源污染物排放量[g/(km×h)]

    圖5 地形模型

    在確定污染物種類(lèi)和污染物排放數(shù)據(jù)后,對(duì)污染物排放源進(jìn)行設(shè)置,將城市道路設(shè)為污染物排放源.在計(jì)算域中將道路作為排放源,輸入排放源在平峰和高峰時(shí)的污染物排放量數(shù)據(jù),根據(jù)實(shí)際距離在道路上通過(guò)Terrain繪制聲屏障,得到最終地形模型,如圖5所示.

    2.4 設(shè)定氣象參數(shù)

    用AirVisual軟件每30min對(duì)高粱橋斜街無(wú)聲屏障路段的氣象數(shù)據(jù)進(jìn)行實(shí)時(shí)監(jiān)測(cè),如表5所示,根據(jù)表所列數(shù)據(jù)分別設(shè)定平、高峰時(shí)的氣象參數(shù).

    表5 平峰、高峰時(shí)期氣象數(shù)據(jù)

    2.5 劃分網(wǎng)格

    采用Panache內(nèi)嵌的網(wǎng)格劃分工具劃分網(wǎng)格,通過(guò)不同劃分方式對(duì)比確定網(wǎng)格參數(shù).在無(wú)需計(jì)算結(jié)果的地方將網(wǎng)格設(shè)置較為稀疏;而在排放源附近網(wǎng)格則設(shè)置較為緊密,網(wǎng)格可以準(zhǔn)確捕捉到聲屏障,從而得出劃分后的網(wǎng)格圖,如圖6所示.

    3 基于計(jì)算流體力學(xué)CFD的仿真分析

    表6 監(jiān)測(cè)點(diǎn)平、高峰時(shí)期污染物濃度(mg/g)

    CFD求解器不同于很多其它常用傳統(tǒng)工具,它可以考慮城市中復(fù)雜建筑、復(fù)雜地形、低和高風(fēng)速、近場(chǎng)和遠(yuǎn)場(chǎng),進(jìn)行風(fēng)場(chǎng)計(jì)算和污染物濃度擴(kuò)散的計(jì)算.針對(duì)本次研究對(duì)象,有兩組數(shù)據(jù)共四種情況需要進(jìn)行仿真模擬.分別為平峰和高峰時(shí)期有無(wú)聲屏障時(shí)的污染物擴(kuò)散情況對(duì)比.

    各監(jiān)測(cè)點(diǎn)污染物濃度數(shù)據(jù)如表6所示:

    對(duì)所得結(jié)果進(jìn)行分析處理,需進(jìn)行分析的結(jié)果主要有四種,包括:(1)某時(shí)刻研究域中,風(fēng)場(chǎng)、污染物濃度、溫度及其它流動(dòng)特性的空間分布;(2)監(jiān)測(cè)點(diǎn)的風(fēng)場(chǎng)、污染物濃度、溫度及其它流動(dòng)特性的時(shí)間分布;(3)給定時(shí)間間隔內(nèi),地面濃度的平均時(shí)間與空間分布;(4)數(shù)值求解的收斂性的結(jié)果.

    對(duì)于本次研究對(duì)象,分別針對(duì)本分分別平峰和高峰兩部分,進(jìn)行后處理分析.

    3.1 平峰時(shí)期

    3.1.1 殘差圖 繪制殘差圖,用于在模擬過(guò)程中檢查求解的收斂情況,如圖7所示.

    殘差繪圖沿縱坐標(biāo)采用對(duì)數(shù)模式,可以看出,所有圖形最后的縱坐標(biāo)數(shù)值均小于-3,所以污染物濃度的觀(guān)察值與擬合值之差較小,擬合程度很高,數(shù)據(jù)均較為精確,誤差較小,可以進(jìn)行進(jìn)一步處理分析.

    3.1.2 矢量圖 速度矢量顯示了風(fēng)場(chǎng),矢量用箭頭表示,箭頭的方向指風(fēng)向,大小和顏色代表風(fēng)速.分別繪制有聲屏障時(shí)和無(wú)聲屏障的時(shí)水平地面平面風(fēng)速矢量圖,如圖8.

    圖8 有、無(wú)聲屏障時(shí)速度矢量圖(km/h)

    對(duì)比可以發(fā)現(xiàn),有、無(wú)聲屏障對(duì)城市道路內(nèi)風(fēng)速大小有較大影響.有聲屏障時(shí)城市道路水平面上的風(fēng)速基本集中在0.03~3m/s之間,而無(wú)聲屏障區(qū)域則大多集中在3~6m/s之間,聲屏障的存在減少了道路內(nèi)的水平面風(fēng)速.而對(duì)于處于道路兩側(cè)、建筑群內(nèi)的空間,同樣可由圖看出,相同的位置有聲屏障算例里風(fēng)速仍然小于無(wú)聲屏障算例風(fēng)速.

    在風(fēng)速矢量圖上,除了速度大小外,風(fēng)速的方向也是重要的指標(biāo)之一.以下通過(guò)局部對(duì)比觀(guān)察聲屏障的存在對(duì)于風(fēng)速方向的影響,如圖9所示.

    圖9 有、無(wú)聲屏障時(shí)局部風(fēng)向

    可以得到聲屏障的存在使得風(fēng)向有了很大的改變.無(wú)聲屏障時(shí),風(fēng)向幾乎垂直于道路流動(dòng),沒(méi)有受到任何阻擋;當(dāng)聲屏障存在時(shí),會(huì)阻擋風(fēng)的流動(dòng),不僅減小了風(fēng)的流動(dòng)速度,還改變了風(fēng)的流動(dòng)方向,使得原本垂直道路的風(fēng)向改為沿著道路的方向流動(dòng).

    3.1.3 濃度分布圖 濃度分布圖可以顯示變量的空間分布.在濃度分布圖中,相同顏色的線(xiàn)或條帶上的所繪變量為常數(shù).().繪制平峰時(shí)有、無(wú)聲屏障情況的3種污染物濃度濃度分布圖,來(lái)比較聲屏障對(duì)于各種污染物濃度擴(kuò)散的影響.

    如圖10(a)(b)所示,有聲屏障時(shí),道路內(nèi)CO的濃度明顯高于無(wú)聲屏障道路內(nèi)相同位置的CO濃度.而在道路外側(cè)的建筑群內(nèi)內(nèi),在道路的順風(fēng)邊,即道路東西邊的建筑領(lǐng)域,有聲屏障時(shí)相同位置CO的濃度低于無(wú)聲屏障時(shí)CO的濃度.在道路的逆風(fēng)方向,即西南向,2種情況CO的濃度都接近于0.

    如圖10(c)(d)所示,NO在有、無(wú)聲屏障情況下的擴(kuò)散情況與CO類(lèi)似.在有聲屏障的時(shí)候,道路內(nèi)NO的濃度同樣高于無(wú)聲屏障道路內(nèi)相同位置的NO濃度.在道路外側(cè)的建筑群中,在道路的順風(fēng)邊,即道路東北邊的建筑領(lǐng)內(nèi),有聲屏障時(shí)相同位置NO的濃度低于無(wú)聲屏障時(shí)NO的濃度.在道路的逆風(fēng)方向,即西南向,2種情況的CO的濃度都接近于0.

    如圖10(e)(f)所示,聲屏障對(duì)于PM的擴(kuò)散影響程度遠(yuǎn)大于CO和NO.在有聲屏障時(shí),不論在道路內(nèi),還是在道路外側(cè)順風(fēng)邊的建筑群內(nèi),相同位置的PM濃度均遠(yuǎn)低于無(wú)聲屏障的情況.而且,在無(wú)聲屏障的情況下,由于缺少聲屏障對(duì)PM的阻擋作用,PM在較遠(yuǎn)處的建筑群區(qū)域仍有非常大的濃度.但在道路的逆風(fēng)方向,即西南向,兩種情況的PM濃度都接近于0.

    綜上所述,在平峰的時(shí)間段里,聲屏障對(duì)于污染物的擴(kuò)散有較大的影響.有聲屏障時(shí),聲屏障改變了風(fēng)速的大小和方向,使得道路內(nèi)的氣體污染物(CO、NO)無(wú)法被吹散,聚集在道路內(nèi)部,相對(duì)于無(wú)聲屏障的情況下,道路內(nèi)的污染物濃度顯著增加.但同時(shí)由于聲屏障的存在,改變了風(fēng)的方向,降低了污染物向道路兩側(cè)擴(kuò)散的程度,使得道路兩側(cè)的污染物濃度相比于無(wú)聲屏障的情況下得到了一定程度上的降低.由于PM為顆粒狀污染物,相對(duì)更易擴(kuò)散,因此聲屏障對(duì)其擴(kuò)散的抑制作用更加明顯,極大地減少了PM的擴(kuò)散,PM擴(kuò)散濃度的降低高達(dá)50%以上.

    3.2 高峰時(shí)期

    3.2.1 殘差圖 繪制殘差圖,如圖11所示.

    同平峰時(shí)間段情況一樣,圖中最后的縱坐標(biāo)數(shù)值均小于-3,所以污染物濃度的觀(guān)察值與擬合值之差較小,擬合程度很高,數(shù)據(jù)均較為精確,可進(jìn)行進(jìn)一步處理分析.

    圖11 無(wú)聲屏障時(shí)3種污染物濃度殘差圖Fig.11 Residual history of three pollutants concentration with and without noise barriers

    3.2.2 矢量圖 分別繪制有、無(wú)聲屏障時(shí)水平地面的風(fēng)速矢量圖如圖12所示.

    如圖12可以發(fā)現(xiàn),現(xiàn)有、無(wú)聲屏障對(duì)城市道路內(nèi)部區(qū)域的風(fēng)速大小有一些影響,但影響程度無(wú)平峰時(shí)大.聲屏障的存在仍然減少了道路內(nèi)部水平面的風(fēng)速.

    圖13為同一條道路的同一位置有、無(wú)聲屏障時(shí)風(fēng)的方向圖.

    圖12 有、無(wú)聲屏障時(shí)速度矢量圖(km/h)

    圖13 有、無(wú)聲屏障時(shí)局部風(fēng)向

    圖14 道路兩側(cè)建筑群

    聲屏障在高峰時(shí)期和平峰時(shí)期對(duì)風(fēng)向的改變相同,無(wú)聲屏障時(shí),風(fēng)向幾乎垂直于道路流動(dòng),未受到任何阻擋;當(dāng)設(shè)置聲屏障時(shí),會(huì)阻擋風(fēng)的流動(dòng),不僅減小了風(fēng)的流動(dòng)速度,還改變了風(fēng)的流動(dòng)方向,使得原本垂直道路的風(fēng)向改為沿著道路的方向流動(dòng).

    通過(guò)對(duì)高峰時(shí)期聲屏障對(duì)風(fēng)速影響的分析可得,聲屏障對(duì)于風(fēng)速大小和方向的影響與平峰時(shí)期類(lèi)似,均是減小風(fēng)速,改變風(fēng)向.但影響程度從圖中可以看出,聲屏障對(duì)風(fēng)速在高峰時(shí)期的影響程度要比平峰時(shí)期要小很多,根據(jù)對(duì)地形條件和氣象條件的分析,主要有兩點(diǎn)原因:(1)高峰時(shí)期的風(fēng)速為3.6km/h,僅為平峰時(shí)期的一半,風(fēng)速減小,影響程度也減少.(2)平峰時(shí)期風(fēng)向?yàn)槲髂巷L(fēng),而高峰時(shí)期為東北風(fēng).由地形圖可以看出,位于道路西南方向建筑群建筑寬度均較小,且建筑之間縫隙較多,如圖14所示,風(fēng)從西南方向吹來(lái)時(shí),建筑對(duì)于風(fēng)的阻擋較弱,所以平峰時(shí)期道路附近風(fēng)的強(qiáng)度較大,聲屏障對(duì)風(fēng)的影響也較大.而位于道路東北方的建筑群建筑寬度較大,且建筑之間縫隙較少,風(fēng)從東北方吹來(lái)的時(shí)候,建筑對(duì)于風(fēng)的阻擋較強(qiáng),所以高峰時(shí)期道路附近風(fēng)較小,聲屏障對(duì)風(fēng)的影響也較小.所以?xún)煞N情況之所以會(huì)有較大差別主要是由平、高峰時(shí)期風(fēng)速、風(fēng)向的不同以及道路兩側(cè)建筑構(gòu)造的不同這兩種因素共同導(dǎo)致的結(jié)果.

    3.2.3 濃度分布圖 繪制高峰時(shí)有、無(wú)聲屏障時(shí)3種污染物濃度濃度分布圖,比較聲屏障對(duì)于各種污染物濃度擴(kuò)散的影響.

    如圖15(a)(b)所示,有聲屏障時(shí),道路內(nèi)部有些位置CO的濃度高于無(wú)聲屏障道路內(nèi)部相同位置的CO濃度,也有一小部分位置CO的濃度是低于無(wú)聲屏障時(shí)的.而在道路外側(cè)的建筑群內(nèi),在道路的順風(fēng)邊,即道路西南邊的建筑領(lǐng)域,有聲屏障時(shí),相同位置CO的濃度低于無(wú)聲屏障時(shí)CO的濃度.

    如圖15(c)(d)所示,在有聲屏障時(shí),道路內(nèi)NO的濃度同樣高于無(wú)聲屏障道路內(nèi)相同位置的NO濃度.對(duì)于道路外側(cè)的建筑群,在道路的順風(fēng)邊,即道路西南邊的建筑領(lǐng)域,有聲屏障時(shí),相同位置NO濃度低于無(wú)聲屏障時(shí)NO濃度.在道路的逆風(fēng)方向,即東北向,兩種情況NO濃度都接近于0.

    如圖15(e)(f)所示,在高峰期,聲屏障對(duì)于PM擴(kuò)散的影響程度相對(duì)于CO和NO相對(duì)較大.在有聲屏障時(shí),不管是在道路內(nèi)部,還是在道路外側(cè)順風(fēng)邊的建筑群內(nèi),相同位置PM濃度均遠(yuǎn)低于無(wú)聲屏障情況下PM濃度.同樣,在道路的逆風(fēng)方向,即東北向,兩種情況PM的濃度都接近于0.

    綜上所示,高峰時(shí)期聲屏障對(duì)于CO和NO擴(kuò)散的影響作用與平峰時(shí)期相近,同樣是聲屏障使得道路內(nèi)部的大氣污染物(CO、NO)無(wú)法被吹散,聚集在道路內(nèi)部,相對(duì)于無(wú)聲屏障的情況下,增加了道路內(nèi)部的污染物濃度,但同時(shí)由于聲屏障的存在也改變了風(fēng)的方向,降低了污染物往道路兩側(cè)擴(kuò)散的程度,使得道路兩側(cè)污染物濃度相比于無(wú)聲屏障的情況得到了一定程度上的降低.而對(duì)于PM擴(kuò)散影響,聲屏障對(duì)于PM擴(kuò)散的抑制作用相比于CO、NO更加明顯,極大地減少了PM的擴(kuò)散,對(duì)PM擴(kuò)散濃度的降低高達(dá)50%以上.

    聲屏障對(duì)污染物擴(kuò)散的影響程度在高峰時(shí)期沒(méi)有平峰時(shí)期那么顯著,由于氣象條件和地形條件的不同,使得風(fēng)的強(qiáng)度在高峰時(shí)期相對(duì)于平峰時(shí)期有了大幅度的下降,很大程度上的減弱了污染物的擴(kuò)散、流動(dòng)能力,所以也使得聲屏障在高峰時(shí)期對(duì)于污染物濃度的擴(kuò)散影響遠(yuǎn)不如平峰時(shí)期那么顯著.

    4 交通污染物現(xiàn)場(chǎng)監(jiān)測(cè)

    采用YT-HPC3000A空氣凈化檢測(cè)儀于2018年5月7日對(duì)高梁橋斜街聲屏障區(qū)域進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè),采樣時(shí)間為60s.

    分別在平峰和高峰時(shí)期選擇相同10個(gè)監(jiān)測(cè)點(diǎn),在所選監(jiān)測(cè)點(diǎn)用儀器進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè),獲取PM和PM10在各個(gè)監(jiān)測(cè)點(diǎn)的濃度,如表7所示.

    表7 平、高峰時(shí)期監(jiān)測(cè)點(diǎn)污染物濃度(μg/m3)

    由于受到儀器設(shè)備的限制,只能現(xiàn)場(chǎng)監(jiān)測(cè)PM2.5和PM10的濃度擴(kuò)散值.由表7的數(shù)據(jù)所知,聲屏障對(duì)PM2.5和PM10污染物濃度的擴(kuò)散有著非常大的影響.根據(jù)比較,監(jiān)測(cè)點(diǎn)1、3、5、7、9中監(jiān)測(cè)點(diǎn)3的PM2.5和PM10濃度最低,9和7處的污染物濃度最高.由監(jiān)測(cè)結(jié)果計(jì)算出,平峰時(shí)期,有聲屏障道路附近的PM2.5濃度降低10%左右,PM10濃度降低50%左右;高峰時(shí)期,有聲屏障道路附近的PM2.5、PM10濃度均降低50%左右.同時(shí),監(jiān)測(cè)點(diǎn)2、4、6的污染物濃度均比監(jiān)測(cè)點(diǎn)1、3、5處的濃度要低,由于聲屏障的存在改變了風(fēng)向,抑制了污染物的擴(kuò)散,使得遠(yuǎn)離道路的位置PM2.5和PM10濃度降低.而監(jiān)測(cè)點(diǎn)8和10均比7、9處的污染物濃度高,是由于這兩處缺少了聲屏障對(duì)于污染物的阻擋作用,經(jīng)過(guò)風(fēng)的擴(kuò)散作用,離道路稍遠(yuǎn)的地方污染物的濃度反而增大,這也表明聲屏障對(duì)于PM顆粒污染物的擴(kuò)散有著很大的抑制作用,產(chǎn)生較大影響.

    5 結(jié)論

    5.1 聲屏障對(duì)風(fēng)速及風(fēng)向有很大影響.針對(duì)監(jiān)測(cè)道路,平峰期風(fēng)向主要為西向,垂直道路方向?yàn)橹饕L(fēng)向,高峰期風(fēng)向主要為東北向,沿道路方向?yàn)橹饕L(fēng)向.通過(guò)矢量圖分析,聲屏障的存在會(huì)阻擋風(fēng)的流動(dòng),減弱了風(fēng)的強(qiáng)度.平峰無(wú)聲屏障時(shí),道路內(nèi)的風(fēng)速大小相對(duì)于有聲屏障時(shí)從0.03m/s提高到了3m/s,甚至6m/s,風(fēng)速增加100到200倍.無(wú)聲屏障時(shí)道路內(nèi)部的風(fēng)速相對(duì)于有聲屏障時(shí)從0.02m/s上升到近2m/s,風(fēng)速增加100倍.除了改變風(fēng)的強(qiáng)度,聲屏障的存在還使得風(fēng)向從原本垂直于道路上升轉(zhuǎn)變?yōu)檠刂缆贩较蛄鲃?dòng).

    5.2 有聲屏障時(shí),道路內(nèi)部風(fēng)的強(qiáng)度和方向發(fā)生很大改變,因此會(huì)影響污染物的擴(kuò)散情況.風(fēng)向改變使得汽車(chē)排放的氣體污染物聚集在道路內(nèi),增加了道路內(nèi)部的氣體污染物濃度.平峰有聲屏障時(shí),道路內(nèi)部CO的濃度相對(duì)于無(wú)聲屏障時(shí)約從0.025′10-6上升到0.05′10-6,NO從0.005′10-6上升到0.01′10-6,上升比均為50%.高峰時(shí)期由于風(fēng)向沿道路方向,影響程度較小,高峰有聲屏障時(shí),CO濃度從0.03′10-6上升到0.05′10-6左右,NO則從0.005′10-6上升到0.008′10-6左右,上升比約為40%.

    5.3 聲屏障的存在不僅改變了道路內(nèi)部污染物排放濃度,對(duì)道路外側(cè)污染物擴(kuò)散濃度也有較大的影響.由于聲屏障改變了風(fēng)向,減弱了風(fēng)往道路外側(cè)的流動(dòng),降低了氣體污染物向道路外側(cè)擴(kuò)散的程度,從而使得道路外側(cè)氣體污染物濃度相對(duì)于無(wú)聲屏障時(shí)得到了一定程度上的降低.平峰有聲屏障時(shí),道路外側(cè)CO的濃度相對(duì)于無(wú)聲屏障的時(shí)候從0.025′10-6下降到約0.0125′10-6,NO濃度則從0.005′10-6下降到0.0025′10-6,下降比均為50%.從濃度分布圖中可知,聲屏障對(duì)道路外側(cè)污染物擴(kuò)散的影響在距道路50~80米內(nèi)較為明顯.高峰時(shí)期由于受地形和氣象條件的響,聲屏障對(duì)于道路外側(cè)污染物擴(kuò)散影響沒(méi)有平峰時(shí)明顯.聲屏障存在的時(shí)候,使得CO的濃度從0.04′10-6降到約0.03′10-6;NO則從0.0075′10-6降到0.006′10-6,下降比約為20%.這種影響在距離道路外側(cè)50m到80m內(nèi)效果明顯,聲屏障在靠近道路處側(cè)對(duì)污染物擴(kuò)散的影響效果并不明顯.

    5.4 結(jié)合仿真結(jié)果和現(xiàn)場(chǎng)監(jiān)測(cè),聲屏障對(duì)PM(顆粒污染物)的影響顯著.有聲屏障時(shí),道路內(nèi)和道路外側(cè)的建筑群內(nèi),PM的濃度均遠(yuǎn)低于無(wú)聲屏障時(shí).平峰有聲屏障時(shí),道路內(nèi)側(cè)部分區(qū)域的PM濃度為0.00015′10-6以下,無(wú)聲屏障的時(shí)候PM濃度為0.0003′10-6以上,聲屏障的存在使得PM濃度下降了超過(guò)一倍.高峰時(shí)期,道路內(nèi)側(cè)部分區(qū)域的PM濃度為0.00025′10-6左右,而無(wú)聲屏障時(shí)PM濃度為0.000375′10-6左右,影響程度依然很高.平峰時(shí)期,道路外側(cè)靠近道路的地方,PM2.5和PM10的濃度在有屏障時(shí)比無(wú)聲屏障時(shí)降低了將近一倍.同樣在有聲屏障的地方,距離道路50m處PM2.5和PM10的濃度相對(duì)于道路附近的污染物濃度分別下降了將近37%和20%.而在無(wú)聲屏障的地方,距離道路50m處PM2.5和PM10的濃度相對(duì)于道路附近處的污染物濃度分別上升了38%和25%.可見(jiàn)聲屏障的存在極大抑制了PM的擴(kuò)散作用,不僅使得道路附近的PM濃度下降,還阻擋了PM的擴(kuò)散.

    [1] 王 宇.北京市近10年空氣環(huán)境質(zhì)量變化特征及影響因素分析 [C]//中國(guó)環(huán)境科學(xué)學(xué)會(huì).2013中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第四卷).中國(guó)環(huán)境科學(xué)學(xué)會(huì), 2013:8. Wang Y. Analysis of the characteristics and influencing factors of air quality change in Beijing in recent 10years [C]//Chinese Society of Environmental Sciences. Papers of the academic annual meeting of the Chinese Society of Environmental Sciences in 2013 (Volume 4). Chinese Society of Environmental Sciences: August 2013.

    [2] Kim J J, Smorodinsky S, Ostro B, et al. Traffic-related air pollution and respiratory health: the east bay Children's respiratory health study [J]. Epidemiology, 2002,13:S100.

    [3] Beelen R, Hoek G, vanden Brandt P A, et al. long-term effects of traffic-related airpollution on mortality in a Dutch cohort (NLCS-AIR study [J]). Environ. Health Perspect., 2008,116:196-202.

    [4] Real-world fuel efficiency and exhaust emissions of light-duty diesel vehicles and their correlation with road conditions [J]. Journal of Environmental Sciences, 2012,5:865-874.

    [5] Cantarella G E. Macroscopic vs. Mesoscopic traffic flow models in Signal Setting Design [C]//Proceedings of 17th International IEEE Conference on Intelligent Transportation Systems Ⅲ.IEEE ITSS Society, IEEE, Chinese Association of Automation, the Institute of Automation, Chinese Academy of Sciences, Xi'an Jiaotong University, Qingdao Academy of Intelligent Industries:,2014:6.

    [6] Wang Guihua, Bai Song, Ogden Joan M. Identifying contributions of on-road motor vehicles to urban air pollution using travel demand model data [J]. Transportation Research Part D: Transport and Environment, 2009,14(3):168-179.

    [7] Tong Z, Baldauf R W, Isakov V, et al. Roadside vegetation barriers designs to mitigate near-road air pollution impacts [J]. Sci. Total Environ, 2016,541:920-927

    [8] 王雪艷.基于道路綠化帶影響的街道峽谷內(nèi)尾氣擴(kuò)散的數(shù)值模擬 [D]. 濟(jì)南:山東大學(xué), 2015. Wang X Y. Numerical simulation of exhaust gas diffusion in street canyons based on the influence of road greenbelt [D]. Ji¢nan: Shandong University, 2015.

    [9] Richard Baldauf. Roadside vegetation design characteristics that can improvelocal, near-road air quality [J]. Transportation Research Part D 52, 2017,5:354–361.

    [10] S. Amini, F. E. Ahangar, N. Schulte, A. Venkatram using models to interpret the impact of roadside barriers on near-road air quality [J]. Atmos. Environ., 2016,138:55-64.

    [11] 劉春玲.淺談空氣污染監(jiān)測(cè)和建模 [J]. 黑龍江科技信息, 2015, (9):16. Liu C L. Discuss on air pollution monitoring and modeling [J]. Heilongjiang Science and Technology Information, 2015,(9):16.

    [12] 劉永紅,畢索陽(yáng),周 兵,等.佛山市中心城區(qū)機(jī)動(dòng)車(chē)限行對(duì)污染物削減效果的分析 [J]. 中國(guó)環(huán)境科學(xué), 2010,30(11):1563-1567. Liu Y H, Bi S Y, Zhou B, et al. Analysis of the effect of motor vehicle restriction on pollutant reduction in downtown Foshan [J]. China Environmental Science, 2010,30(11):1563-1567.

    [13] 張文杰.道路交通尾氣污染對(duì)我國(guó)城市居民生活環(huán)境影響研究 [J].市場(chǎng)研究, 2016,(11):6-8. Zhang W J. Study on the impact of road traffic exhaust pollution on the living environment of urban residents in China [J]. Market Research, 2016,(11):6-8.

    [14] 劉寶章,康 凌,林官明.澳門(mén)荷蘭園街汽車(chē)尾氣擴(kuò)散規(guī)律實(shí)驗(yàn)研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2000,20(S1):27-33. Liu B Z, Kang L, Lin G M. Experimental study on vehicle exhaust diffusion in Dutch Garden Street, Macau [J]. Journal of Environmental Sciences, 2000,(S1):27-33.

    [15] 周洪昌,高延令,吳曉琰.街道峽谷地面源污染物擴(kuò)散規(guī)律的風(fēng)洞試驗(yàn)研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 1994,(4):389-396.Zhou H C, Gao Y L, Wu X Y. Wind tunnel experimental study on the diffusion law of ground source pollutants in street canyon [J]. Journal of Environmental Science, 1994,(4):389-396.

    [16] 何澤能,高陽(yáng)華,李永華,等.城市街道峽谷對(duì)稱(chēng)性對(duì)內(nèi)部氣流場(chǎng)的影響研究 [J]. 氣象與環(huán)境學(xué)報(bào), 2008,(2):62-67.He Z N, Gao Y H, Li Y H, et al. Study on the influence of symmetry of urban street canyons on internal airflow field [J]. Journal of Meteorology and Environment, 2008,(2):62-67.

    [17] 王 樂(lè),張?jiān)苽?顧兆林.動(dòng)態(tài)風(fēng)場(chǎng)及交通流量下街道峽谷內(nèi)污染物擴(kuò)散模擬 [J]. 中國(guó)環(huán)境科學(xué), 2012,32(12):2161-2167. Wang L, Zhang Y W, Gu Z L. Simulation of pollutant diffusion in street canyons under dynamic wind field and traffic flow [J]. China Environmental Science, 2012,32(12):2161-2167.

    [18] 周 禮.汽車(chē)尾氣污染防治綜述 [J]. 工程設(shè)計(jì)與研究, 2011,(2):42- 44. Zhou Li. Overview of prevention and control of automobile exhaust pollution [J]. Engineering Design and Research, 2011,(2):42-44.

    [19] Lin M Y, Hagler G, Baldauf R, et al. A. Khlystov The effects of vegetation barriersbarriers on near-road ultrafine particle number and carbon monoxide concentrations [J]. Sci. Total Environ., 2016,553: 372-379.

    The impact of sound barriers on air quality of urban road environment.

    QIN Xiao-chun*, LIU Rui, HAN Ying

    (School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China)., 2019,39(8):3558~3571

    To precisely quantify the impact of sound barriers on air quality of urban road, we applied the fluidyn-Panache simulation and CFD solver to launch numerical models and analysis based on the three-dimensional computational fluid dynamics theory. We compared the wind velocity vector maps and contaminant concentration nephograms of PM2.5, CO and NOnear the road under both presence and absence of noise barriers, monitored the contaminant concentration in the simulation area on the spot, and evaluated the regularity of impacts of barriersbarriers on air quality of urban road environment. Our results showed that 1) the wind speed inside the road decreased from 2m/s to 0.02m/s; 2) the concentration of CO and NOinside the road increased by 40% to 50%;3) the concentration of PM inside the road dropped by 50%;4) the concentration of CO and NO3outside the road decreased by 20% to 50%; and 5) the concentration of PM outside the road reduced by 20% to 38%. We found that 1) the noise barriers blocked the flow of wind, weakened the intensity of wind, and changed the direction of wind from vertical to upward to along the road; 2) by changing the direction of wind pollutants gather in the road, the concentration of gas pollutants rises; and 3) the flow of wind to the outside, gas pollutants out of the road, as well as the extent of lateral diffusion drop. Noise barriers has a significant effect on thePM2.5concentration in the buildings on the inside and outside of the road that PM2.5decreases significantly compared with that of the silent barriers.

    road;sound barriers;ambient air quality;pollutant diffusion

    X51

    A

    1000-6923(2019)08-3558-14

    秦曉春(1982-),女,內(nèi)蒙古包頭人,副教授,博士,主要從事交通環(huán)保與景觀(guān)、綠色可持續(xù)交通等方面的研究工作.發(fā)表論文40篇.

    2019-02-25

    國(guó)家自然科學(xué)基金資助面上項(xiàng)目(51878039)

    * 責(zé)任作者, 副教授, xcqin@bjtu.edu.cn

    猜你喜歡
    風(fēng)向屏障風(fēng)速
    咬緊百日攻堅(jiān) 筑牢安全屏障
    屏障修護(hù)TOP10
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
    基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
    一道屏障
    維護(hù)網(wǎng)絡(luò)安全 筑牢網(wǎng)絡(luò)強(qiáng)省屏障
    自然與風(fēng)Feeling Nature
    基于GARCH的短時(shí)風(fēng)速預(yù)測(cè)方法
    行業(yè)統(tǒng)計(jì)帶來(lái)哪些風(fēng)向?
    考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
    超碰成人久久| 99国产精品99久久久久| 国产高清激情床上av| 亚洲欧美一区二区三区久久| 母亲3免费完整高清在线观看| 亚洲三区欧美一区| 午夜久久久在线观看| 日本黄色日本黄色录像| 手机成人av网站| 操美女的视频在线观看| 欧美成狂野欧美在线观看| 80岁老熟妇乱子伦牲交| 亚洲精品中文字幕在线视频| 两个人看的免费小视频| 亚洲色图 男人天堂 中文字幕| 日本精品一区二区三区蜜桃| 日本黄色视频三级网站网址 | 另类精品久久| 成人18禁在线播放| 久久天躁狠狠躁夜夜2o2o| 国产在线免费精品| 国产欧美日韩精品亚洲av| 在线观看一区二区三区激情| 国产日韩欧美亚洲二区| 激情在线观看视频在线高清 | 啦啦啦在线免费观看视频4| 亚洲精品久久午夜乱码| 中文字幕最新亚洲高清| 搡老乐熟女国产| 国产高清videossex| 多毛熟女@视频| 多毛熟女@视频| 最近最新中文字幕大全免费视频| 另类亚洲欧美激情| 国产精品久久久久久精品古装| 色播在线永久视频| av天堂在线播放| 夜夜夜夜夜久久久久| 男女之事视频高清在线观看| 人妻一区二区av| 80岁老熟妇乱子伦牲交| 国产一区二区三区视频了| 99热网站在线观看| 精品久久久久久电影网| 亚洲成a人片在线一区二区| tocl精华| 亚洲国产看品久久| 欧美乱码精品一区二区三区| 欧美成狂野欧美在线观看| 啦啦啦在线免费观看视频4| 日韩 欧美 亚洲 中文字幕| 高清av免费在线| 美女主播在线视频| 亚洲性夜色夜夜综合| 国产精品久久久人人做人人爽| 午夜两性在线视频| 中文字幕av电影在线播放| 人成视频在线观看免费观看| 国产精品免费一区二区三区在线 | 男男h啪啪无遮挡| 国产亚洲精品第一综合不卡| 视频在线观看一区二区三区| 一级黄色大片毛片| 狂野欧美激情性xxxx| 肉色欧美久久久久久久蜜桃| 丝袜喷水一区| 国产成人欧美| 亚洲人成伊人成综合网2020| 国产淫语在线视频| 久久中文字幕一级| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频| 一级片'在线观看视频| 91麻豆精品激情在线观看国产 | 99久久人妻综合| 国产高清视频在线播放一区| xxxhd国产人妻xxx| av又黄又爽大尺度在线免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 大型av网站在线播放| 老司机影院毛片| av片东京热男人的天堂| 久久久国产精品麻豆| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 亚洲全国av大片| 久久久国产一区二区| 制服诱惑二区| 国产精品影院久久| 国产aⅴ精品一区二区三区波| 蜜桃在线观看..| 日韩欧美国产一区二区入口| 日韩免费av在线播放| 日韩一卡2卡3卡4卡2021年| 少妇被粗大的猛进出69影院| 亚洲国产欧美网| 久久免费观看电影| 欧美人与性动交α欧美精品济南到| 久久精品亚洲精品国产色婷小说| 国产片内射在线| 久久香蕉激情| 日本黄色日本黄色录像| 亚洲男人天堂网一区| 亚洲精品久久成人aⅴ小说| 老司机在亚洲福利影院| 高清欧美精品videossex| 国产福利在线免费观看视频| 精品国产亚洲在线| 91麻豆av在线| 亚洲精品av麻豆狂野| 99re在线观看精品视频| aaaaa片日本免费| 免费看十八禁软件| 欧美黑人欧美精品刺激| av视频免费观看在线观看| www.自偷自拍.com| 成人影院久久| 亚洲av国产av综合av卡| 超碰97精品在线观看| 日本一区二区免费在线视频| 国产日韩欧美视频二区| 一本色道久久久久久精品综合| 亚洲av第一区精品v没综合| 久久久水蜜桃国产精品网| 80岁老熟妇乱子伦牲交| 久久精品国产a三级三级三级| 亚洲欧美一区二区三区久久| 91成年电影在线观看| 亚洲人成电影免费在线| 亚洲avbb在线观看| 亚洲,欧美精品.| 在线观看www视频免费| 亚洲自偷自拍图片 自拍| 考比视频在线观看| 黄色视频,在线免费观看| 国产免费视频播放在线视频| 国产精品美女特级片免费视频播放器 | 国产精品美女特级片免费视频播放器 | 精品免费久久久久久久清纯 | 国产成人啪精品午夜网站| 日本av手机在线免费观看| www.精华液| 精品少妇久久久久久888优播| 色婷婷久久久亚洲欧美| 在线观看人妻少妇| 国产不卡一卡二| kizo精华| 久久久欧美国产精品| 国产在线视频一区二区| 一二三四社区在线视频社区8| 黄色视频不卡| 大香蕉久久网| 岛国毛片在线播放| 女性生殖器流出的白浆| 精品国产乱码久久久久久男人| 精品国产一区二区三区久久久樱花| 久久久久国内视频| 成人手机av| 99香蕉大伊视频| 老司机福利观看| 黄色毛片三级朝国网站| 亚洲精品国产区一区二| 免费少妇av软件| 欧美日韩中文字幕国产精品一区二区三区 | 精品第一国产精品| 免费一级毛片在线播放高清视频 | 菩萨蛮人人尽说江南好唐韦庄| 欧美成人免费av一区二区三区 | 黄色片一级片一级黄色片| 宅男免费午夜| 精品人妻在线不人妻| 午夜福利一区二区在线看| 日本精品一区二区三区蜜桃| 精品国产乱码久久久久久男人| 一进一出抽搐动态| 大香蕉久久成人网| 国产一区二区在线观看av| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人啪精品午夜网站| 国产精品欧美亚洲77777| 久久午夜综合久久蜜桃| 欧美人与性动交α欧美软件| 国产高清videossex| 波多野结衣一区麻豆| 亚洲欧美一区二区三区黑人| 91精品三级在线观看| 免费少妇av软件| 在线观看一区二区三区激情| 桃花免费在线播放| 夜夜夜夜夜久久久久| 欧美一级毛片孕妇| 久久久久国产一级毛片高清牌| 久久ye,这里只有精品| 日韩免费高清中文字幕av| 自线自在国产av| 日韩欧美一区二区三区在线观看 | 狠狠狠狠99中文字幕| 久久国产精品大桥未久av| 日本精品一区二区三区蜜桃| 少妇裸体淫交视频免费看高清 | 久久99热这里只频精品6学生| 91精品国产国语对白视频| 色婷婷久久久亚洲欧美| 免费在线观看视频国产中文字幕亚洲| 怎么达到女性高潮| 国产精品国产高清国产av | 夜夜骑夜夜射夜夜干| 国产精品免费视频内射| 妹子高潮喷水视频| 日韩免费av在线播放| 女同久久另类99精品国产91| 亚洲专区中文字幕在线| 亚洲av成人不卡在线观看播放网| av线在线观看网站| 久久久国产欧美日韩av| 精品人妻熟女毛片av久久网站| 免费女性裸体啪啪无遮挡网站| 欧美精品人与动牲交sv欧美| 亚洲免费av在线视频| 最新在线观看一区二区三区| 国产精品 欧美亚洲| 精品国内亚洲2022精品成人 | 男人舔女人的私密视频| 老熟女久久久| 两个人免费观看高清视频| 国产有黄有色有爽视频| av在线播放免费不卡| 久久久久精品人妻al黑| 狠狠狠狠99中文字幕| 一区在线观看完整版| 法律面前人人平等表现在哪些方面| 精品国产乱子伦一区二区三区| 亚洲自偷自拍图片 自拍| 真人做人爱边吃奶动态| 国产1区2区3区精品| 无限看片的www在线观看| 国产精品久久久久久精品电影小说| 美女国产高潮福利片在线看| 91精品国产国语对白视频| 少妇 在线观看| 欧美日韩av久久| 欧美激情久久久久久爽电影 | 中文亚洲av片在线观看爽 | 成年女人毛片免费观看观看9 | 亚洲精品中文字幕在线视频| 亚洲美女黄片视频| 欧美日韩亚洲高清精品| 亚洲色图av天堂| av电影中文网址| 美国免费a级毛片| 99九九在线精品视频| 黄色毛片三级朝国网站| 91大片在线观看| 欧美成人免费av一区二区三区 | 国产99久久九九免费精品| 桃花免费在线播放| 中文欧美无线码| a在线观看视频网站| 露出奶头的视频| 国产av又大| 自拍欧美九色日韩亚洲蝌蚪91| 欧美精品av麻豆av| 黑人操中国人逼视频| 国产精品一区二区精品视频观看| 精品亚洲乱码少妇综合久久| 国产91精品成人一区二区三区 | 国产片内射在线| 久久精品国产综合久久久| 天堂俺去俺来也www色官网| 久久精品熟女亚洲av麻豆精品| 日韩欧美免费精品| 日本vs欧美在线观看视频| 天天躁日日躁夜夜躁夜夜| 丰满饥渴人妻一区二区三| 国产高清视频在线播放一区| 天堂俺去俺来也www色官网| 天天躁狠狠躁夜夜躁狠狠躁| 成年人免费黄色播放视频| 亚洲精品自拍成人| 91麻豆av在线| 色播在线永久视频| 久久精品亚洲精品国产色婷小说| 美女视频免费永久观看网站| a级片在线免费高清观看视频| 国产精品 国内视频| 狂野欧美激情性xxxx| www日本在线高清视频| 婷婷成人精品国产| 少妇被粗大的猛进出69影院| 久久久久久久精品吃奶| 视频区图区小说| 黑人欧美特级aaaaaa片| 国产免费视频播放在线视频| 国产一区二区三区视频了| 亚洲视频免费观看视频| 一边摸一边抽搐一进一小说 | 91麻豆精品激情在线观看国产 | 久久久欧美国产精品| 午夜福利一区二区在线看| 久久久精品免费免费高清| 免费在线观看完整版高清| 日韩视频在线欧美| 国产精品免费一区二区三区在线 | 最近最新中文字幕大全免费视频| 极品教师在线免费播放| 99精品欧美一区二区三区四区| 极品少妇高潮喷水抽搐| 精品人妻1区二区| 精品国产超薄肉色丝袜足j| 精品亚洲成国产av| h视频一区二区三区| 丁香欧美五月| 亚洲国产看品久久| 国产男靠女视频免费网站| 国产深夜福利视频在线观看| 成人国产一区最新在线观看| 亚洲欧美色中文字幕在线| 国产精品国产av在线观看| 亚洲 国产 在线| 精品亚洲乱码少妇综合久久| 国产一卡二卡三卡精品| 国内毛片毛片毛片毛片毛片| 女人精品久久久久毛片| 国产1区2区3区精品| 亚洲精品在线观看二区| 91成年电影在线观看| 国产成人欧美在线观看 | 亚洲国产av新网站| 丝袜美足系列| 欧美+亚洲+日韩+国产| 深夜精品福利| 色精品久久人妻99蜜桃| e午夜精品久久久久久久| 亚洲午夜理论影院| 女人高潮潮喷娇喘18禁视频| 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 久久 成人 亚洲| 国产精品久久久久久精品古装| 亚洲中文av在线| 精品一区二区三区四区五区乱码| 一级毛片女人18水好多| 伊人久久大香线蕉亚洲五| 国产精品国产av在线观看| 视频在线观看一区二区三区| 久久国产精品人妻蜜桃| 悠悠久久av| 亚洲人成77777在线视频| 操出白浆在线播放| 午夜福利影视在线免费观看| 丝袜在线中文字幕| 啦啦啦 在线观看视频| 亚洲国产av新网站| 在线十欧美十亚洲十日本专区| 男女无遮挡免费网站观看| 久久精品91无色码中文字幕| 国产精品欧美亚洲77777| 热re99久久精品国产66热6| 在线天堂中文资源库| 男人操女人黄网站| 精品一区二区三区av网在线观看 | 久久婷婷成人综合色麻豆| 啦啦啦 在线观看视频| 在线观看免费日韩欧美大片| 我的亚洲天堂| 欧美 日韩 精品 国产| 高清视频免费观看一区二区| 欧美在线一区亚洲| 高清在线国产一区| 黄网站色视频无遮挡免费观看| 91麻豆精品激情在线观看国产 | 欧美日韩亚洲高清精品| 叶爱在线成人免费视频播放| 人妻久久中文字幕网| 亚洲三区欧美一区| 久久天堂一区二区三区四区| 欧美日韩福利视频一区二区| 99国产精品一区二区三区| h视频一区二区三区| 欧美精品一区二区免费开放| 99九九在线精品视频| 精品一区二区三卡| 成年人黄色毛片网站| 亚洲成人手机| 丁香欧美五月| 国产av又大| 欧美人与性动交α欧美软件| 亚洲男人天堂网一区| 国产97色在线日韩免费| 国产激情久久老熟女| 国产精品免费视频内射| 欧美日本中文国产一区发布| 最新在线观看一区二区三区| 一二三四在线观看免费中文在| 欧美午夜高清在线| 美女主播在线视频| a级毛片黄视频| 捣出白浆h1v1| 久久久久视频综合| 国产亚洲av高清不卡| 国产亚洲欧美在线一区二区| 精品久久蜜臀av无| 建设人人有责人人尽责人人享有的| 日韩三级视频一区二区三区| 国产成人免费观看mmmm| 亚洲五月色婷婷综合| 亚洲人成77777在线视频| 成人三级做爰电影| 亚洲国产成人一精品久久久| 大型av网站在线播放| 色视频在线一区二区三区| 99国产极品粉嫩在线观看| 最近最新免费中文字幕在线| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 无限看片的www在线观看| 欧美老熟妇乱子伦牲交| 日韩一卡2卡3卡4卡2021年| 国产精品99久久99久久久不卡| 日韩三级视频一区二区三区| √禁漫天堂资源中文www| 三上悠亚av全集在线观看| 国产成人精品无人区| 久久青草综合色| 精品亚洲成国产av| 久久久国产成人免费| 曰老女人黄片| 99精品欧美一区二区三区四区| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频 | 日本av免费视频播放| 中文字幕高清在线视频| 成人精品一区二区免费| 精品国产亚洲在线| 欧美日韩亚洲综合一区二区三区_| 欧美另类亚洲清纯唯美| 亚洲精品国产区一区二| 免费在线观看完整版高清| 这个男人来自地球电影免费观看| 美女扒开内裤让男人捅视频| 老司机午夜十八禁免费视频| 丁香六月欧美| 悠悠久久av| 欧美在线黄色| 在线观看免费午夜福利视频| av不卡在线播放| videos熟女内射| 日韩欧美一区二区三区在线观看 | 免费久久久久久久精品成人欧美视频| 手机成人av网站| 极品少妇高潮喷水抽搐| 亚洲,欧美精品.| 国产欧美日韩精品亚洲av| 91av网站免费观看| 久久精品亚洲av国产电影网| 俄罗斯特黄特色一大片| 汤姆久久久久久久影院中文字幕| 美女主播在线视频| 亚洲成a人片在线一区二区| av片东京热男人的天堂| 精品一区二区三区四区五区乱码| 91国产中文字幕| 热99re8久久精品国产| 2018国产大陆天天弄谢| 老司机午夜福利在线观看视频 | 精品一区二区三卡| 曰老女人黄片| 欧美 亚洲 国产 日韩一| 免费在线观看日本一区| 天天躁日日躁夜夜躁夜夜| 一边摸一边抽搐一进一出视频| 999精品在线视频| 麻豆乱淫一区二区| 热99久久久久精品小说推荐| 久久免费观看电影| 狂野欧美激情性xxxx| 高清在线国产一区| 老熟妇仑乱视频hdxx| 天堂中文最新版在线下载| 操美女的视频在线观看| 午夜福利视频精品| 亚洲精品av麻豆狂野| 亚洲五月色婷婷综合| 国产精品免费一区二区三区在线 | 久久九九热精品免费| 久久 成人 亚洲| 天天躁夜夜躁狠狠躁躁| 十八禁网站免费在线| 亚洲人成伊人成综合网2020| 久久精品亚洲精品国产色婷小说| 他把我摸到了高潮在线观看 | 女警被强在线播放| 热99国产精品久久久久久7| 好男人电影高清在线观看| avwww免费| 中文欧美无线码| 91麻豆av在线| 丁香六月欧美| 久久国产精品影院| 精品人妻1区二区| 交换朋友夫妻互换小说| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利影视在线免费观看| 美女视频免费永久观看网站| 丝袜在线中文字幕| av视频免费观看在线观看| 天堂中文最新版在线下载| 18禁黄网站禁片午夜丰满| 高潮久久久久久久久久久不卡| 中文亚洲av片在线观看爽 | 人妻 亚洲 视频| 韩国精品一区二区三区| 成人国产一区最新在线观看| 在线永久观看黄色视频| 狠狠精品人妻久久久久久综合| 精品乱码久久久久久99久播| 午夜免费成人在线视频| 国产精品香港三级国产av潘金莲| 精品一品国产午夜福利视频| av有码第一页| 成人亚洲精品一区在线观看| 国产有黄有色有爽视频| 亚洲欧洲日产国产| 无限看片的www在线观看| 久久久久网色| 狠狠精品人妻久久久久久综合| 人人妻人人爽人人添夜夜欢视频| 中文字幕另类日韩欧美亚洲嫩草| 最近最新免费中文字幕在线| 久久国产精品人妻蜜桃| 水蜜桃什么品种好| 丝袜在线中文字幕| a级毛片黄视频| 国产在线一区二区三区精| 999久久久国产精品视频| 在线观看一区二区三区激情| 欧美日韩中文字幕国产精品一区二区三区 | 人人妻,人人澡人人爽秒播| 国产成人精品久久二区二区91| 亚洲五月婷婷丁香| 91老司机精品| 757午夜福利合集在线观看| 80岁老熟妇乱子伦牲交| 超色免费av| 99九九在线精品视频| 国产av一区二区精品久久| 国产欧美日韩一区二区精品| 国产精品熟女久久久久浪| 性高湖久久久久久久久免费观看| 制服诱惑二区| 国产一卡二卡三卡精品| 精品国产一区二区三区四区第35| 国产黄频视频在线观看| 国产精品免费视频内射| 中文字幕另类日韩欧美亚洲嫩草| 国产av又大| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| 国产在线视频一区二区| 十八禁网站网址无遮挡| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 天堂动漫精品| 高清在线国产一区| 18禁观看日本| 中亚洲国语对白在线视频| 美女福利国产在线| 国产在线观看jvid| 捣出白浆h1v1| 国产日韩一区二区三区精品不卡| 日本欧美视频一区| 一级毛片电影观看| 乱人伦中国视频| 日韩人妻精品一区2区三区| 国产精品一区二区免费欧美| 一区二区三区激情视频| 国产在线免费精品| 国产单亲对白刺激| 亚洲欧美激情在线| 久久亚洲精品不卡| 又大又爽又粗| 法律面前人人平等表现在哪些方面| 在线十欧美十亚洲十日本专区| 亚洲精品国产区一区二| 80岁老熟妇乱子伦牲交| 亚洲精华国产精华精| 91精品国产国语对白视频| 99精品欧美一区二区三区四区| 久久久久视频综合| 亚洲av片天天在线观看| 97在线人人人人妻| 国产亚洲精品久久久久5区| av国产精品久久久久影院| 黄色成人免费大全| 丁香六月天网| 精品久久久久久久毛片微露脸| 免费不卡黄色视频| 一本综合久久免费| 久久狼人影院| 欧美精品一区二区免费开放| 日本av手机在线免费观看| 999久久久精品免费观看国产| 精品乱码久久久久久99久播| av超薄肉色丝袜交足视频| 久久久久久久久免费视频了| 另类精品久久| 精品国产一区二区三区久久久樱花| 午夜免费鲁丝| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久久久大奶| 777米奇影视久久| 热99国产精品久久久久久7| 国产一区二区三区在线臀色熟女 | 国产精品熟女久久久久浪| 久久亚洲精品不卡|