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

    北京大氣熱力和動(dòng)力結(jié)構(gòu)對(duì)污染物輸送和擴(kuò)散條件的影響

    2017-05-23 11:04:18孫兆彬李梓銘廖曉農(nóng)唐宜西趙秀娟中國(guó)氣象局北京城市氣象研究所北京00089京津冀環(huán)境氣象預(yù)報(bào)預(yù)警中心北京00089
    中國(guó)環(huán)境科學(xué) 2017年5期
    關(guān)鍵詞:逆溫北京地區(qū)熱力

    孫兆彬,李梓銘,廖曉農(nóng),唐宜西,趙秀娟(.中國(guó)氣象局北京城市氣象研究所,北京 00089;.京津冀環(huán)境氣象預(yù)報(bào)預(yù)警中心,北京 00089)

    北京大氣熱力和動(dòng)力結(jié)構(gòu)對(duì)污染物輸送和擴(kuò)散條件的影響

    孫兆彬1,2*,李梓銘2,廖曉農(nóng)2,唐宜西2,趙秀娟1(1.中國(guó)氣象局北京城市氣象研究所,北京 100089;2.京津冀環(huán)境氣象預(yù)報(bào)預(yù)警中心,北京 100089)

    利用北京、天津、河北地區(qū)常規(guī)氣象數(shù)據(jù)、PM2.5濃度、激光雷達(dá)數(shù)據(jù),鐵塔數(shù)據(jù),結(jié)合風(fēng)廓線數(shù)據(jù)對(duì)2015年10月13~17日發(fā)生在北京的空氣重污染過(guò)程中大氣熱力和動(dòng)力結(jié)構(gòu)進(jìn)行分析.結(jié)果表明,白天熱力條件轉(zhuǎn)好后,邊界層低層垂直方向氣流由上-下運(yùn)動(dòng)轉(zhuǎn)變?yōu)閮A斜運(yùn)動(dòng),同時(shí)水平方向偏南風(fēng)增大,高濃度氣團(tuán)向北京移動(dòng),污染物輸送增強(qiáng),造成這種空間結(jié)構(gòu)的原因主要是熱力條件加強(qiáng)后等溫面由準(zhǔn)東-西向轉(zhuǎn)變?yōu)槲鞅?東南向,由于山地加熱激發(fā)的重力波振幅在白天增強(qiáng),高空水平方向輸送至北京的污染物在垂直方向出現(xiàn)向上和向下的傳輸.16日白天地面PM2.5濃度的快速升高,這是在PM2.5存在一定垂直交換的情況下,近地面偏南風(fēng)輸送的作用所致,以上出現(xiàn)的大氣污染物輸送和擴(kuò)散機(jī)制可以在一定程度上解釋區(qū)域重污染背景下,白天混合層高度升高而近地層PM2.5濃度卻隨之升高的現(xiàn)象.

    熱力和動(dòng)力結(jié)構(gòu);污染物輸送和擴(kuò)散;PM2.5;激光雷達(dá);重力波

    近些年北京地區(qū)空氣重污染事件頻發(fā),給居民身體健康帶來(lái)極大危害,因此研究重污染天氣形成機(jī)制對(duì)于預(yù)報(bào)、預(yù)警的及時(shí)發(fā)布具有重要意義.大氣混合層是人類生活和生產(chǎn)的主要空間,許多重要天氣和重污染事件均發(fā)生在這一空間中[1-8],因此混合層中大氣動(dòng)力和熱力過(guò)程對(duì)污染物濃度具有重要影響.此外北京地區(qū)三面環(huán)山,大氣的結(jié)構(gòu)除受天氣尺度系統(tǒng)影響外,還受局地山谷風(fēng)影響,天氣尺度系統(tǒng)與山谷風(fēng)耦合在一起所形成的三維大氣熱力和動(dòng)力結(jié)構(gòu)在很大程度上決定著污染物的輸送和擴(kuò)散.

    諸多研究結(jié)果表明,對(duì)混合層高度估算的方法較多,例如利用觀測(cè)到的氣象數(shù)據(jù)進(jìn)行估算、基于數(shù)值模式預(yù)報(bào)結(jié)果進(jìn)行估算、利用激光雷達(dá)和微波輻射計(jì)的監(jiān)測(cè)資料進(jìn)行估算[9-14].對(duì)于混合層高度與污染物之間關(guān)系的研究已有很多[15-19],研究表明較高的最大混合層高度出現(xiàn)頻率較大的月,空氣質(zhì)量為優(yōu)-良出現(xiàn)的頻率較高,即混合層高度高,大氣環(huán)境容量加大,利于污染物的擴(kuò)散,污染物濃度低,反之混合層高度降低,大氣環(huán)境容量減小,污染物堆積壓縮不易擴(kuò)散,污染物濃度升高.總結(jié)起來(lái),混合層高度與大氣污染物濃度整體呈負(fù)相關(guān),這一結(jié)論已經(jīng)長(zhǎng)期指導(dǎo)我國(guó)環(huán)境氣象及空氣質(zhì)量預(yù)報(bào)業(yè)務(wù).本文的研究在經(jīng)過(guò)長(zhǎng)時(shí)間大量數(shù)據(jù)統(tǒng)計(jì)計(jì)算后發(fā)現(xiàn)這一結(jié)論仍然成立,但是在對(duì)北京地區(qū)PM2.5濃度數(shù)據(jù)和混合層數(shù)據(jù)進(jìn)行分析后發(fā)現(xiàn),在區(qū)域空氣重污染背景下的某些時(shí)段(例如白天,以小時(shí)為時(shí)間分辨率),隨著混合層升高,逆溫減弱,客觀上大氣對(duì)污染物局地容納能力和擴(kuò)散條件好轉(zhuǎn),而監(jiān)測(cè)到的污染物濃度卻呈現(xiàn)升高的趨勢(shì).因此本研究選擇 2015年 10月13~17日北京地區(qū)已經(jīng)出現(xiàn)了空氣重污染的時(shí)段開(kāi)展研究,希望在一定程度上揭示混合層高度升高而PM2.5濃度也隨之升高的大氣動(dòng)力和熱力機(jī)制.

    1 材料與方法

    1.1 PM2.5數(shù)據(jù)、氣象數(shù)據(jù)和激光雷達(dá)數(shù)據(jù)

    PM2.5數(shù)據(jù)來(lái)自北京市環(huán)保局(http://www. bjmemc.com.cn/g68.aspx)、河北省環(huán)保廳(http://121. 28.49.85:8080/)、天津市環(huán)保局(http://air.tjemc.org. cn/)所公布的數(shù)據(jù).氣象資料來(lái)自北京市氣象局,鐵塔資料來(lái)自于中國(guó)科學(xué)院大氣物理研究所.激光雷達(dá)數(shù)據(jù)為北京地區(qū)EV-LINAR型微脈沖激光雷達(dá)在垂直方向的觀測(cè)數(shù)據(jù).

    1.2 混合層高度的計(jì)算方法

    可用地面常規(guī)氣象參數(shù)來(lái)估算平均混合層高度,本文根據(jù)羅氏法計(jì)算混合層厚度[13-14],計(jì)算公式如下:

    式中: T為地面氣溫,℃;Td為地面露點(diǎn)溫度,℃;uz為z高度處的平均風(fēng)速,m/s;z0為地面粗糙度,m;f為柯氏參數(shù),f= 2?sinΦ,? 為地轉(zhuǎn)角速度,Φ 為地理緯度;P為帕斯圭爾穩(wěn)定度級(jí)別(大氣穩(wěn)定度級(jí)別為A~F時(shí),P值依次為1~6).

    風(fēng)廓線數(shù)據(jù)目前常被用來(lái)研究混合層內(nèi)各種天氣現(xiàn)象發(fā)生的規(guī)律,在混合層內(nèi)常有較強(qiáng)的湍流,而混合層外湍流較弱,會(huì)在混合層頂形成湍流的極大值,本文依據(jù)風(fēng)廓線雷達(dá)重要湍流參數(shù)-折射率結(jié)構(gòu)常數(shù)(CN2)來(lái)計(jì)算混合層高度,CN2(zi)為 zi高度處的折射率結(jié)構(gòu)常數(shù),湍流能量法的計(jì)算公式如下:

    2 污染物時(shí)空分布特征

    2.1 氣象要素時(shí)間序列變化

    10月13~17日,北京地區(qū)近地面受偏南風(fēng)和偏東風(fēng)影響.偏南風(fēng)主要出現(xiàn)在每日的午后至傍晚,將河北中南部污染物向北京、天津、河北東北部一帶輸送,導(dǎo)致污染物大量積累;其他時(shí)段為偏東風(fēng)影響,對(duì)北京來(lái)說(shuō),偏東風(fēng)主要來(lái)向?yàn)槲廴疚镆呀?jīng)大量積累的天津、河北東北部地區(qū),因此偏南風(fēng)、偏東風(fēng)兩只氣流交替將污染物向北京輸送,導(dǎo)致北京平原地區(qū)污染物快速積累和維持.垂直方向上,北京中高空為弱的下沉氣流,將偏南風(fēng)壓制在較低層次,使其所輸送的污染物無(wú)法翻越太行山、燕山山脈,只能在山前積累,同時(shí)在高濕的條件下,顆粒物吸濕增長(zhǎng),導(dǎo)致能見(jiàn)度降低(圖1).

    圖1 2015年10月13~17日北京逐小時(shí)風(fēng)向、風(fēng)速、相對(duì)濕度Fig.1 Beijing hourly wind direction, wind speed, relative humidity October 2015 13~17

    2.2 PM2.5濃度時(shí)間序列變化

    如圖2所示,13~17日京津冀地區(qū)出現(xiàn)一次空氣污染事件.能見(jiàn)度從13日下午的30km下降到14日的 5km, PM2.5的濃度逐漸上升,14日 20:00達(dá)到150μg/m3左右.14日夜間至 15日受弱冷空氣影響,PM2.5濃度有所下降,能見(jiàn)度略有好轉(zhuǎn).16日白天空氣質(zhì)量急劇惡化,其中北京地區(qū) PM2.5濃度連續(xù)24h超過(guò)250μg/m3,17日下午超過(guò)300μg/m3,能見(jiàn)度一直維持在2km以下,局地小于1km,17日后半夜冷空氣南下,此次污染過(guò)程自北向南逐漸結(jié)束.

    圖2 2015年10月13~17日京津冀地區(qū)PM2.5濃度(a)和能見(jiàn)度(b)時(shí)間演變特征Fig.2 Beijing, Tianjin and Hebei PM2.5concentration (a) and visibility (b) variation during October 2015 13~17

    2.3 京津冀地區(qū)地面風(fēng)場(chǎng)和PM2.5空間分布

    由圖3可見(jiàn)16日08:00~10:00時(shí),北京東南部、河北廊坊、天津3地交界處有一個(gè)PM2.5濃度為300~500μg/m3的高濃度氣團(tuán),空間尺度較小(約為幾十km),午后伴隨著近地層偏南風(fēng)明顯增大,此高濃度氣團(tuán)向北移動(dòng),于13:00前后到達(dá)北京城區(qū),這一時(shí)間節(jié)點(diǎn)與圖2中PM2.5濃度明顯升高的時(shí)間相一致,受排放源分布及偏南風(fēng)等不利于污染物擴(kuò)散條件的影響,京津冀平原沿山地區(qū)形成了常見(jiàn)的 “廠”字形污染帶.

    圖3 2015年10月16日京津冀地區(qū)08:00~19:00的地面自動(dòng)站10m風(fēng)場(chǎng)和PM2.5濃度(μg/m3)Fig.3 10m wind field and PM2.5concentration distribution in Beijing Tianjin and Hebei region (08:00~19:00) (μg/m3)

    2.4 近地面 PM2.5濃度急劇升高時(shí)段北京上空的消光系數(shù)

    為了研究偏南風(fēng)下污染物水平和垂直擴(kuò)散條件,本研究選取PM2.5濃度急劇增長(zhǎng)的16日白天,分析激光雷達(dá)消光系數(shù)的時(shí)空演變趨勢(shì).1km以上北京上空無(wú)明顯的消光系數(shù)增大現(xiàn)象,由北京地區(qū) 1km 內(nèi)的消光系數(shù)演變(圖 4)可見(jiàn), 13:00~15:00出現(xiàn)了3個(gè)消光系數(shù)較高的層次,這3個(gè)層次分別處于200~600m、100m左右、近地層,消光系數(shù)數(shù)值最大區(qū)位于 200~600m,大值中心的出現(xiàn)時(shí)間為13:00前后,這與圖3中PM2.5高濃度氣團(tuán)到達(dá)北京城區(qū)的時(shí)間相一致,說(shuō)明16日白天外源輸送很可能是導(dǎo)致北京地區(qū)PM2.5濃度快速上升的主要原因,而且由于高層和低層消光系數(shù)大值中心同時(shí)出現(xiàn),所以外源輸送很可能在高層和低層同時(shí)進(jìn)行.15:00左右原本位于 200~600m的消光系數(shù)大值區(qū)強(qiáng)度減弱,而且有向上移動(dòng)的趨勢(shì),消光系數(shù)大值區(qū)與高層大氣最大梯度由600m上移至900m,整層大氣消光系數(shù)有趨于均一的趨勢(shì),伴隨著局地?zé)崃l件好轉(zhuǎn),1000m大氣混合作用增強(qiáng),15:00 ~16:00整層大氣消光系數(shù)達(dá)到了 16日白天最為均一的狀態(tài),16:00~17:00在800m左右的高度再次出現(xiàn)消光系數(shù)大值中心,而且 17:00~19:00近地層的消光系數(shù)在明顯增大,更為重要的是消光系數(shù)大值區(qū)隨著時(shí)間出現(xiàn)了高值區(qū)向下伸展的現(xiàn)象,由 800m向下伸展至 200m,甚至更低的高度,對(duì)于這一現(xiàn)象出現(xiàn)的原因?qū)⒃诒疚?.2章節(jié)進(jìn)行分析.

    圖4 北京地區(qū)2015年10月16日11:00~20:00激光雷達(dá)消光系數(shù)Fig.4 Beijing area lidar extinction coefficient from 11:00~20:00 on the 16th

    3 大氣動(dòng)力過(guò)程影響

    3.1 混合層高度變化

    混合層高度的計(jì)算具有很大不確定性,因此本文采用羅氏法和湍流能量法計(jì)算了混合層高度.如圖5(a)和(b)所示,混合層高度均具有明顯的日變化特征,其中在14~16日空氣重污染期間,兩種算法的混合層高度均在2000m以下波動(dòng),其中15日受弱冷空氣影響,混合層高度最高,16日夜間至白天混合層高度均呈升高趨勢(shì).在14~16日的空氣重污染期間,兩種算法的混合層高度均在2000m以下波動(dòng),兩者的區(qū)別是羅氏法的計(jì)算結(jié)果對(duì)于近地層大氣的動(dòng)力、熱力條件及大氣穩(wěn)定度更加敏感,湍流能量法對(duì)于低層大氣垂直方向大氣動(dòng)力結(jié)構(gòu)的改變更加敏感,兩種算法顯示16日夜間至白天混合層高度均呈升高趨勢(shì).

    圖5 北京地區(qū)混合層高度Fig.5 Mixed layer height

    3.2 區(qū)域風(fēng)場(chǎng)日變化分析

    圖6 2015年10月16日京津冀地區(qū)10m風(fēng)場(chǎng)分布Fig.6 10m wind field distribution in Beijing Tianjin Hebei region in October 16, 2015

    BJ-RUC v2.0系統(tǒng)的空間區(qū)域共有2個(gè),分辨率為9km(D1區(qū)域)和3km(D2區(qū)域),水平網(wǎng)格點(diǎn)數(shù)分別為400×649和424×550,垂直方向分為38層,BJ-RUC v2.0系統(tǒng)主要物理過(guò)程具體如下:WSM6顯式微物理方案;Kain-Fritsch對(duì)流參數(shù)化方案(3km區(qū)域無(wú)對(duì)流參數(shù)化方案);YSU邊界層方案;RRTM長(zhǎng)波輻射方案;Goddard短波輻射方案;輻射方案每15min計(jì)算一次;Noah LSM陸面模式.

    圖6(a)~ (f)反映了北京及周邊地區(qū)10月16日地面風(fēng)場(chǎng)變化,通過(guò)數(shù)值模式的預(yù)報(bào)結(jié)果能夠看到近地層風(fēng)場(chǎng)精細(xì)特征,經(jīng)過(guò)綜合分析我們可以得到如下規(guī)律:

    (1)京津冀地區(qū)夜間近地層受輻射冷卻效應(yīng)影響,冷空氣下泄,平原沿山地區(qū)夜間多為山風(fēng)控制,但風(fēng)速較小,在北京東北部的密云區(qū)、平谷區(qū)由于均三面環(huán)山,在其平原地區(qū)山風(fēng)對(duì)應(yīng)的冷空氣明顯下泄,山風(fēng)風(fēng)速較大,常形成異于北京平原地區(qū)的氣象條件.

    (2)早晨日照輻射作用加強(qiáng)后,近地層熱條件轉(zhuǎn)好,谷風(fēng)率先出現(xiàn)山區(qū)(圖6b),北京平原地區(qū)風(fēng)場(chǎng)由山風(fēng)控制轉(zhuǎn)變?yōu)闊o(wú)固定的風(fēng)場(chǎng),京津冀沿山地區(qū)為風(fēng)速較小的弱谷風(fēng)控制,而此弱谷風(fēng)控制區(qū)的東側(cè)則為較明顯的偏南風(fēng)影響,自沿山地區(qū)向東偏南風(fēng)的分量逐漸增大,會(huì)促進(jìn)京津冀平原地區(qū)正渦度的形成,有利于沿山地區(qū)動(dòng)力輻合作用的增強(qiáng),這是京津冀沿山地區(qū)大氣熱力作用先變化,而后促使動(dòng)力條件發(fā)生改變,從而有利于沿山地區(qū)污染物積累的有力證據(jù).

    (3)午后京津冀沿山地區(qū)受山地和平原上空空氣熱力梯度加強(qiáng)的影響,平原沿山地區(qū)的山風(fēng)明顯加強(qiáng)為較強(qiáng)的偏南風(fēng),在熱力和動(dòng)力條件雙重加強(qiáng)的條件下,平原沿山地區(qū)混合層高度也隨之升高(圖5),大氣對(duì)污染物的輸送作用加強(qiáng).

    3.3 風(fēng)廓線數(shù)據(jù)資料

    如圖 7(a)所示,0~2000m 高度上 16日10:00~20:00為一致的偏南風(fēng),高層風(fēng)速較大,低層風(fēng)速較小,大于 12m/s的低空急流在 16日04:00~10:00出現(xiàn)在 1000m 以上的高空,在10:00~11:00、13:00前后、16:00~17:00急流的強(qiáng)度減弱,急流在時(shí)間序列演變中發(fā)生了斷裂現(xiàn)象,這段時(shí)間也正是北京本地?zé)崃l件逐漸轉(zhuǎn)好的時(shí)段.16日00:00~06:00北京近地層為一致的上升氣流,高層為下沉氣流,而在 06:00-10:00熱力條件開(kāi)始轉(zhuǎn)好的時(shí)段,2000m以下的大氣出現(xiàn)了明顯的垂直上升和下沉的氣流,但在 10:00~17:00熱力條件繼續(xù)加強(qiáng)的階段,垂直氣流開(kāi)始向遠(yuǎn)離山地的一側(cè)傾斜,出現(xiàn)了垂直方向的斜升氣流,總體看來(lái)上升運(yùn)動(dòng)的強(qiáng)度在2000m內(nèi)最大可大于 0.75m/s,下沉運(yùn)動(dòng)則可大于 0.45m/s,上升運(yùn)動(dòng)主要為熱力強(qiáng)迫引起,而下沉運(yùn)動(dòng)更多為動(dòng)力強(qiáng)迫造成,所以其強(qiáng)度小于上升運(yùn)動(dòng),上升和下沉運(yùn)動(dòng)影響著北京地區(qū)污染物擴(kuò)散和輸送條件.

    圖7 16日00:00~20:00的北京風(fēng)廓線Fig.7 Beijing wind profile from 00:00~20:00 on the 16th

    4 大氣熱力過(guò)程影響

    4.1 大氣所320m鐵塔數(shù)據(jù)分析

    由圖8可見(jiàn),早晨08:00 250~320m高度間仍然存在逆溫,至10:00受熱力條件轉(zhuǎn)好影響,逆溫層頂向上抬升至320m以上的高度,320m以下仍存在較弱的逆溫,14:00~16:00期間北京上空320m內(nèi)逆溫基本消失,由圖4可見(jiàn)消光系數(shù)大值區(qū)出現(xiàn)也正好是320m內(nèi)逆溫逐漸減弱消失的時(shí)段.伴隨著午后熱力條件顯著轉(zhuǎn)好(圖 9c),相對(duì)濕度的降低(圖 9d),偏南風(fēng)明顯增大(圖 9a和 9b),高層的偏南風(fēng)在 13:00~14:00明顯增大,與消光系數(shù)(圖4)有較好的對(duì)應(yīng)關(guān)系,再次從風(fēng)場(chǎng)角度印證了偏南風(fēng)對(duì)污染物有明顯的輸送作用.

    圖8 320m鐵塔16日溫度廓線Fig.8 320m tower temperature profile on the 16th

    圖9 320m鐵塔16日11:00~20:00氣象條件變化過(guò)程Fig.9 320m tower meteorological conditions from 11:00to 20:00 on the 16th

    4.2 垂直運(yùn)動(dòng)轉(zhuǎn)變?yōu)閮A斜運(yùn)動(dòng)的機(jī)制以及對(duì)污染物輸送和擴(kuò)散的影響

    由圖10所示16日08:00北京上空有逆溫存在,平原上空等溫面呈準(zhǔn)東-西向分布,16日 10:00逆溫層頂抬高,逆溫的強(qiáng)度減弱,北京西部地區(qū)率先增溫,西部地區(qū)等溫面呈西北-東南向傾斜,這說(shuō)明山區(qū)在無(wú)明顯天氣系統(tǒng)影響時(shí),日出以后太陽(yáng)輻射加熱使山谷壁附近的空氣變暖,其變暖速度要快于相同海拔高度上平原地區(qū)上空的空氣,16日12:00和14:00北京上空的逆溫已經(jīng)消失,溫度場(chǎng)受西部、北部山區(qū)熱力加熱的作用變?yōu)橐恢碌奈鞅?東南向.結(jié)合圖10及圖11,當(dāng)北京上空存在逆溫期間(08:00~10:00)垂直方向的運(yùn)動(dòng)方向以上-下為主,當(dāng)逆溫消失期間(11:00~17:00),伴隨著等溫面由準(zhǔn)東西向轉(zhuǎn)變?yōu)槲鞅?東南向,垂直方向的運(yùn)動(dòng)向東傾斜.即垂直運(yùn)動(dòng)的方向總是傾向于盡可能與等溫面正交的方向,這種垂直運(yùn)動(dòng)的機(jī)制,保證了近地層被加熱的空氣以最大的效率穿過(guò)等溫面,在分析天氣尺度系統(tǒng)對(duì)擴(kuò)散條件影響時(shí),垂直運(yùn)動(dòng)影響比水平運(yùn)動(dòng)相對(duì)較弱,可以忽略不計(jì),但是在本研究中空間尺度相當(dāng)于中尺度運(yùn)動(dòng)的情況來(lái)說(shuō),垂直運(yùn)動(dòng)對(duì)擴(kuò)散條件的影響則更為重要.

    圖10 BJ-RUC數(shù)值模式16日08:00~14:00溫度場(chǎng)Fig.10 Temperature field of BJ-ruc numerical simulation from 08:00to 14:00 on the 16th

    在15日夜間至16日08:00北京處于穩(wěn)定層結(jié)的影響,由于多層逆溫(圖 10中 08:00)的存在且風(fēng)速較小,大氣對(duì)污染物的輸送較弱,隨著上文提到的等溫面向東側(cè)傾斜,以及垂直方向斜升氣流的出現(xiàn),水平方面上谷風(fēng)加強(qiáng).嚴(yán)格來(lái)講,該谷風(fēng)是中尺度系統(tǒng)對(duì)應(yīng)的風(fēng),與天氣尺度高空槽前對(duì)應(yīng)的偏南風(fēng)疊加后,就是本研究在鐵塔和風(fēng)廓線數(shù)據(jù)所觀測(cè)的偏南風(fēng).由圖9及圖11中的V風(fēng)分量可見(jiàn),11:00~17:00 2000m 內(nèi)偏南風(fēng)增大,這也正是激光雷達(dá)顯示的高消光系數(shù)所出現(xiàn)的時(shí)間,伴隨著偏南風(fēng)的增大,垂直運(yùn)動(dòng)的振幅也隨之增大(圖 12),污染物由河北中南部輸送至北京地區(qū),水平方向出現(xiàn)消光系數(shù)大值區(qū),而且消光系數(shù)大梯度區(qū)向高空和近地層伸展,說(shuō)明這種垂直運(yùn)動(dòng)振幅隨著熱力條件逐漸增大,而且垂直方向傾斜運(yùn)動(dòng)的氣流會(huì)促進(jìn)來(lái)自水平方向高濃度氣團(tuán)向高空和近地層擴(kuò)散,近地層監(jiān)測(cè)站所觀測(cè)到的濃度也正是在這種極為復(fù)雜多尺度系統(tǒng)、多維度的運(yùn)動(dòng)的結(jié)果.

    圖11 BJ-ruc數(shù)值模式15日20:00~16日20:00U風(fēng)(上)和V風(fēng)(下)Fig.11 Wind field of BJ-ruc numerical simulation from 20:00 on the 15th to 20:00 on the 16th

    圖12 BJ-ruc數(shù)值模式16日垂直速度Fig.12 Vertical Speed of BJ-ruc on the 16th

    圖 13為本個(gè)例中北京地區(qū)在無(wú)明顯天氣系統(tǒng)影響下垂直運(yùn)動(dòng)的概念模型,夜間北京大氣層結(jié)穩(wěn)定,有逆溫存在,早晨太陽(yáng)輻射加熱近地層及山區(qū),但熱力條件仍較弱,相同海拔高度下,山區(qū)表面空氣的加熱速度快于同海拔高度的平原地區(qū),上升運(yùn)動(dòng)首先在山區(qū)及沿山地區(qū)出現(xiàn),上升運(yùn)動(dòng)則會(huì)將近地層較為溫暖的空氣和污染物向高空輸送,直到遇到逆溫層,此時(shí)被地面帶入高空氣塊的溫度小于上層的逆溫層溫度,因此受到負(fù)浮力的作用,此氣塊向下運(yùn)動(dòng),被帶入高層的污染物則在逆溫層以下均勻混合,而此氣塊向下運(yùn)動(dòng)的同時(shí)也將溫度較高逆溫層中的暖空氣向下夾卷,逆溫層以下大氣受到地表加熱及高層夾卷共同作用,溫度上升,逆溫的強(qiáng)度隨之減弱,隨著太陽(yáng)輻射作用的加強(qiáng),北京地區(qū)由于地形的作用,等溫面準(zhǔn)東西向轉(zhuǎn)變?yōu)槲鞅?東南向,并且垂直氣流的振幅增大(圖 12),逆溫強(qiáng)度減弱,逆溫層以下大氣污染物混合更加充分和均勻,垂直方向氣流由上-下運(yùn)動(dòng)轉(zhuǎn)變?yōu)橄驏|側(cè)傾斜的斜升和斜降運(yùn)動(dòng),以使高低層大氣之間的熱交換效率達(dá)到最大.對(duì)于污染物的輸送,在夜間穩(wěn)定層結(jié)下,逆溫層以下水平風(fēng)速較小,對(duì)污染物的輸送能力較弱,即便逆溫層以上由于夜間地表對(duì)動(dòng)量消耗減少而形成的低空急流,其能否明顯將河北中南部沿山地區(qū)的污染物輸送至北京仍存在著極大的不確定性,以往鮮有研究,即便果真出現(xiàn)此夜間的低空急流對(duì)污染物的輸送,由于逆溫層的阻隔,逆溫層以上低空急流所帶來(lái)的污染物較難向底層輸送,但是早晨 08:00~10:00北京上空逆溫在減弱的同時(shí),逆溫層以下的偏南風(fēng)都在增大,逆溫的存在雖然會(huì)阻隔高低空污染物的混合和交換,但是依然不會(huì)明顯影響水平方向污染物的傳輸,這就是 16日白天無(wú)論逆溫強(qiáng)弱和有無(wú),近地面監(jiān)測(cè)站污染物濃度一直持續(xù)上升的原因,當(dāng)然這其中仍然存在熱力條件激發(fā)的振幅不斷放大的重力波在逆溫消失前后所進(jìn)行垂直方向大氣污染物濃度混合的作用.

    此外還值注意的是,在等溫面由于山體加熱作用發(fā)生傾斜后,垂直方向氣流發(fā)生傾斜運(yùn)動(dòng),山體右側(cè)溫度梯度加大,水平方向南風(fēng)分量加大,總體的結(jié)果就是,在太行山沿山一帶排放源密、污染最為嚴(yán)重地區(qū)的近地層的高濃度氣團(tuán)被這種傾斜上升而且向北運(yùn)動(dòng)的力量驅(qū)使,流向北京.

    圖14是Hooke[20]所提出的經(jīng)典的重力波能量和位相傳播的示意圖,t=t1時(shí)刻由于北京西部和北部山地受太陽(yáng)輻射熱力作用強(qiáng)迫,開(kāi)始激發(fā)重力波在波包中的波動(dòng)產(chǎn)生了切變,到達(dá)t2時(shí)刻(t2>t1),波包背離原地向右上方移動(dòng).在 t3時(shí)刻(t3>t2),波包繼續(xù)向右上方移動(dòng),波的能量繼續(xù)向這一方向傳播,而波的位相則向右下方傳播,在實(shí)際大氣中,重力波的位相和能量的傳播受諸多因素控制,如風(fēng)切變、溫度層結(jié)等因素,情況較為復(fù)雜.關(guān)于重力波的研究主要出現(xiàn)在對(duì)降水天氣的動(dòng)力學(xué)診斷和分析中,在大氣環(huán)境領(lǐng)域所引起的關(guān)注較少,可重力波又是對(duì)大氣污染物擴(kuò)散和傳輸具有重要影響的一類波動(dòng),應(yīng)當(dāng)引起重視.對(duì)于短周期波動(dòng)k12<<k22,k1為垂直方向的波數(shù),k2為水平方向的波數(shù),位相傳播接近于水平方向,量能傳播接近于垂直方向,若大氣不可壓,則質(zhì)點(diǎn)的運(yùn)動(dòng)與波的移動(dòng)方向垂直,在本研究中由于地面可假設(shè)為剛體,所以波動(dòng)不可能向地下傳播,因此風(fēng)廓線觀測(cè)到的垂直方向的運(yùn)動(dòng)為上下或傾斜運(yùn)動(dòng),位相則為自西向東.

    圖13 16日08:00~14:00北京地區(qū)邊界層垂直氣流結(jié)構(gòu)特征(+為上升運(yùn)動(dòng),-為下沉運(yùn)動(dòng))Fig.13 Boundary layer vertical airflow structure in Beijing Region (+: rising movement, -: sinking movement)

    圖14 重力波能和位相的傳播Fig.14 Propagation of gravity wave energy and phase

    5 不確定性討論

    由于本研究主要研究區(qū)域?yàn)楸本┑貐^(qū),未詳細(xì)討論河北地區(qū)污染物的積累和傳輸,而北京地區(qū)污染物的外源性輸送主要來(lái)自于河北地區(qū),所以河北地區(qū)污染物如何從近地層被輸送至高空而向北京運(yùn)動(dòng)文中未詳細(xì)分析,對(duì)此展開(kāi)討論需要結(jié)合空氣質(zhì)量模式及多地的風(fēng)廓線、激光雷達(dá)等特種資料加以分析.

    本研究雖討論了伴隨偏南風(fēng)增大,熱力擾動(dòng)所造成污染物垂直方向傳輸?shù)木唧w機(jī)制,但是其在垂直方向的傳輸數(shù)量仍然沒(méi)有量化的結(jié)果,這需要利用數(shù)值模式進(jìn)一步模擬驗(yàn)證.

    目前對(duì)城市空氣質(zhì)量的監(jiān)測(cè)站點(diǎn)主要集中于近地層,所監(jiān)測(cè)到的污染物濃度是水平和垂直傳輸?shù)墓餐Y(jié)果,針對(duì)這一方面開(kāi)展的詳細(xì)討論的研究成果極少,本研究也僅是利用各種常規(guī)和特種資料進(jìn)行了一次嘗試,仍然需要更多的資料和研究結(jié)果加以佐證.

    6 結(jié)論

    6.1 2015年10月13~17日北京地區(qū)出現(xiàn)了一次空氣重污染過(guò)程,重污染過(guò)程期間,北京地區(qū)白天熱力條件轉(zhuǎn)好后,逆溫減弱,垂直方向由上-下運(yùn)動(dòng)轉(zhuǎn)變?yōu)閮A斜運(yùn)動(dòng),同時(shí)偏南風(fēng)增大,污染物輸送增強(qiáng),造成這種動(dòng)力結(jié)構(gòu)的原因主要是熱力條件加強(qiáng)后等溫面由準(zhǔn)東-西向轉(zhuǎn)變?yōu)槲鞅?東南向,使氣流穿越等溫面對(duì)近地面熱量向上輸送效率最大造成的.

    6.2 空氣重污染過(guò)程期間,由山地加熱激發(fā)的重力波振幅在白天增強(qiáng),造成高空輸送至北京的污染物在垂直方向出現(xiàn)向上和向下的傳輸, 16日白天地面 PM2.5濃度的快速升高,是在PM2.5存在一定垂直交換的情況下,近地面偏南風(fēng)輸送所形成的大氣污染物輸送和擴(kuò)散機(jī)制可以在一定程度上解釋區(qū)域重污染背景下,白天混合層高度升高,逆溫減弱而近地層 PM2.5濃度隨之升高的現(xiàn)象.

    [1] 丁一匯,馬淑芬.1991年江淮暴雨期地敢通量與混合層結(jié)構(gòu)的研究 [J]. 氣象學(xué)報(bào), 1991,55(3):257-270.

    [2] 張福青,杜華武,蔣全榮.成熟臺(tái)風(fēng)邊界層作用的數(shù)值研究 [J].熱帶氣象學(xué)報(bào), 1994,10(2):107-114.

    [3] 劉曉慧,馬翠平,趙玉廣,等.2014年河北中南部?jī)纱沃伥蔡鞖獬梢蚍治?[J]. 中國(guó)環(huán)境科學(xué), 2016,36(6):1621-1629.

    [4] 王曉玲,王海燕,王珊珊,等.邊界層準(zhǔn)靜止干線觸發(fā)的中尺度暴雨機(jī)理分析 [J]. 高原氣象, 2015,34(5):1310-1322.

    [5] 徐祥德,王寅鈞,趙天良,等.中國(guó)大地形東側(cè)霾空間分布“避風(fēng)港”效應(yīng)及其“氣候調(diào)節(jié)”影響下的年代際變異 [J]. 科學(xué)通報(bào), 2015(12):1132-1143.

    [6] 吳 兌,廖碧婷,吳 蒙,等.環(huán)首都圈霾和霧的長(zhǎng)期變化特征與典型個(gè)例的近地層輸送條件 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(1): 1-11.

    [7] 劉 建,范紹佳,吳 兌,等.珠江三角洲典型灰霾過(guò)程的邊界層特征 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(6):1664-1674.

    [8] 蘇福慶,高慶先,張志剛,等.北京邊界層外來(lái)污染物輸送通道[J]. 環(huán)境科學(xué)研究, 2004,17(1):26-29.

    [9] 賀千山,毛節(jié)泰.北京城市大氣混合層與氣溶膠垂直分布觀測(cè)研究 [J]. 氣象學(xué)報(bào), 2005,63(3):374-384.

    [10] 林常青,楊東偉,李成才,等.北京地區(qū)大氣氣溶膠的激光雷達(dá)觀測(cè)及反演算法研究 [J]. 北京大學(xué)學(xué)報(bào):自然科學(xué)版, 2013,49(3): 426-434.

    [11] 劉思波,何文英,劉紅燕,等.地基微波輻射計(jì)探測(cè)大氣邊界層高度方法 [J]. 應(yīng)用氣象學(xué)報(bào), 2015,26(5):626-635.

    [12] 劉夢(mèng)娟,陳 敏.BJ-RUC系統(tǒng)對(duì)北京夏季邊界層的預(yù)報(bào)性能評(píng)估 [J]. 應(yīng)用氣象學(xué)報(bào), 2014,25(2):212-221.

    [13] 程水源,張寶寧.大氣混合層高度的確定與計(jì)算方法研究 [J].中國(guó)環(huán)境科學(xué), 1997,17(6):512-516.

    [14] 馬 金,鄭向東.混合層厚度的經(jīng)驗(yàn)計(jì)算及與探空觀測(cè)對(duì)比分析[J]. 應(yīng)用氣象學(xué)報(bào), 2011,22(5):567-576.

    [15] 孫 揚(yáng),王躍思,劉廣仁,等.北京地區(qū)一次大氣環(huán)境持續(xù)嚴(yán)重污染過(guò)程中SO2的垂直分布分析 [J]. 環(huán)境科學(xué), 2006,27(3):408-414.

    [16] 尤煥苓,劉偉東,譚江瑞.北京地區(qū)平均最大混合層厚度的時(shí)間變化特征 [J]. 氣象, 2010,36(5):51-55.

    [17] 劉 梅,嚴(yán)文蓮,張 備,等.2013年1月江蘇霧霾天氣持續(xù)和增強(qiáng)機(jī)制分析 [J]. 氣象, 2014(7):835-843.

    [18] 郭 利,張艷昆,劉樹華,等.北京地區(qū)PM10質(zhì)量濃度與邊界層氣象要素相關(guān)性分析 [J]. 北京大學(xué)學(xué)報(bào):(自然科學(xué)版), 2011, 47(4):607-612.

    [19] 鄭 飛,張 鐳,朱 江.復(fù)雜地形城市冬季邊界層對(duì)氣溶膠輻射效應(yīng)的響應(yīng) [J]. 大氣科學(xué), 2006,30(1):171-179.

    [20] Hooke W H. Gravity Waves [M]// Mesoscale Meteorology and Forecasting. American Meteorological Society, 1986.

    致謝:本研究得到中國(guó)科學(xué)院大氣物理研究所鐵塔數(shù)據(jù)的支持,在此表示感謝!

    Influence of Beijing atmospheric thermal and dynamics structure on pollutants transport and diffusion conditions.

    SUN Zhao-bin1,2*, LI Zi-ming2, LIAO Xiao-nong2, TANG Yi-xi2, ZHAO Xiu-juan1(1.Institute of Urban Meteorology, Beijing, China Meteorological Administration, Beijing 100089, China;2.Environment Meteorology Forecast Center of Beijing-Tianjin-Hebei, China Meteorological Administration, Beijing 100089, China). China Environmental Science, 2017,37(5):1693~1705

    Based on the conventional meteorological data, PM2.5concentration, lidar data, tower data combined with wind profile data on 13~17 October 2015 in Beijing, Tianjin and Hebei, thermal and dynamic structure of heavy air pollution process have been researched. The results showed that the thermal conditions turn better during the day, vertical direction in the lower boundary layer changes from top-down motion into tilting motion, while southerly winds strengthened, pollution was transported enhanced, the main reasons for this spatial structure was that quasi east-west isotherm changed into the northwest-southeast direction after the thermal conditions strengthened and gravity wave caused by mountain heating excited amplitude increased during the day, resulted in air pollutants transmitted to Beijing at high-altitude appeared up and down transmission in the vertical direction, PM2.5concentration increased rapidly, there was vertical exchange under the condition of near ground southerly transports on November 16th, above transport and dispersion mechanism may explain the phenomenon, which mixing height increased with an increased PM2.5concentration in the daytime under the background of region heavy pollution.

    thermal and dynamic structure;pollutant transport and dispersion;PM2.5;laser radar;gravity waves

    X511

    A

    1000-6923(2017)05-1693-13

    孫兆彬(1985-),男,黑龍江省黑河市人,碩士,工程師,主要從事空氣重污染形成機(jī)制和大氣污染對(duì)人體健康影響評(píng)估的研究.發(fā)表論文20余篇.

    2016-10-06

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41305130,41575010);國(guó)家科技支撐計(jì)劃資助項(xiàng)目(2014BAC23B03,2014BAC16B04);北京市自然科學(xué)基金資助項(xiàng)目(8161004);北京市科技計(jì)劃課題(Z151100002115045)

    * 責(zé)任作者, 工程師, szb850804@163.com

    猜你喜歡
    逆溫北京地區(qū)熱力
    小諸葛漫“話”天氣之逆溫
    濟(jì)南市低空逆溫特征及其對(duì)空氣質(zhì)量影響分析
    熱力工程造價(jià)控制的影響因素及解決
    熱力站設(shè)備評(píng)測(cè)分析
    酒泉市大氣邊界層逆溫特征及其與沙塵天氣的關(guān)系
    周六福520愛(ài)跑節(jié)1000人登陸西安城墻 熱力開(kāi)跑
    北京地區(qū)居民鎘攝入量評(píng)估
    2011—2013年長(zhǎng)春市逆溫特征分析
    1949—1966年北京地區(qū)貞操觀的變革——解放與進(jìn)步
    北京地區(qū)地下水人工回灌簡(jiǎn)介
    欧美另类亚洲清纯唯美| 欧美在线黄色| 国产高清视频在线观看网站| 国产精品 国内视频| 欧美黄色淫秽网站| 免费av观看视频| 国产一区二区三区视频了| eeuss影院久久| 18禁在线播放成人免费| e午夜精品久久久久久久| 国产精品电影一区二区三区| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 观看美女的网站| 精品无人区乱码1区二区| 中文在线观看免费www的网站| 免费无遮挡裸体视频| 亚洲精品日韩av片在线观看 | 99精品久久久久人妻精品| 国产在视频线在精品| 成人高潮视频无遮挡免费网站| 日本五十路高清| 亚洲国产欧美人成| 国产成人av教育| 国产精品电影一区二区三区| 波多野结衣高清无吗| 久久精品影院6| 亚洲五月婷婷丁香| xxx96com| 成人国产综合亚洲| 国产精品99久久久久久久久| 免费看a级黄色片| 亚洲狠狠婷婷综合久久图片| 欧美黄色淫秽网站| 国产精品一区二区三区四区免费观看 | 一区福利在线观看| 久久6这里有精品| 色视频www国产| 人妻久久中文字幕网| 国产在视频线在精品| 好看av亚洲va欧美ⅴa在| 国产免费一级a男人的天堂| 国产私拍福利视频在线观看| 久久性视频一级片| 亚洲在线自拍视频| 91九色精品人成在线观看| 久久人妻av系列| 欧美国产日韩亚洲一区| 制服人妻中文乱码| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 中文字幕高清在线视频| 人人妻人人澡欧美一区二区| 国产精品女同一区二区软件 | 真人一进一出gif抽搐免费| 国产aⅴ精品一区二区三区波| 国产精品99久久99久久久不卡| www日本黄色视频网| 久久久国产精品麻豆| 亚洲熟妇中文字幕五十中出| 小说图片视频综合网站| 亚洲av二区三区四区| 老司机深夜福利视频在线观看| 757午夜福利合集在线观看| 老汉色av国产亚洲站长工具| 免费大片18禁| 日本撒尿小便嘘嘘汇集6| 日韩欧美在线乱码| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| 制服人妻中文乱码| 国产高潮美女av| 少妇的逼好多水| 亚洲国产色片| 中文资源天堂在线| 丁香欧美五月| 亚洲av熟女| 99久久综合精品五月天人人| 黄色日韩在线| 亚洲国产欧洲综合997久久,| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 国产99白浆流出| 岛国视频午夜一区免费看| 久久性视频一级片| 一区二区三区国产精品乱码| 亚洲欧美日韩无卡精品| 亚洲欧美日韩卡通动漫| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 色综合站精品国产| 亚洲精品456在线播放app | 国产精品野战在线观看| 日韩欧美精品v在线| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 麻豆成人av在线观看| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 国产探花在线观看一区二区| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 亚洲精品456在线播放app | 91在线观看av| 午夜福利欧美成人| 久久性视频一级片| 在线观看舔阴道视频| 亚洲自拍偷在线| 黄色视频,在线免费观看| a在线观看视频网站| 免费人成视频x8x8入口观看| 午夜福利视频1000在线观看| 蜜桃亚洲精品一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲在线自拍视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 在线观看舔阴道视频| 操出白浆在线播放| 在线观看午夜福利视频| 99热6这里只有精品| 国内久久婷婷六月综合欲色啪| 天堂av国产一区二区熟女人妻| 精品福利观看| 嫩草影视91久久| 久久99热这里只有精品18| 国产一区二区三区在线臀色熟女| 青草久久国产| 欧美日韩亚洲国产一区二区在线观看| 99久久精品一区二区三区| 少妇高潮的动态图| 亚洲美女黄片视频| 国产三级黄色录像| 久久精品91无色码中文字幕| 国产色爽女视频免费观看| 国产精品综合久久久久久久免费| 18禁美女被吸乳视频| 亚洲男人的天堂狠狠| 亚洲av二区三区四区| 亚洲av美国av| 日韩有码中文字幕| 国内精品一区二区在线观看| 亚洲av美国av| 国产单亲对白刺激| 欧美三级亚洲精品| 偷拍熟女少妇极品色| 亚洲欧美精品综合久久99| 国产视频一区二区在线看| 岛国视频午夜一区免费看| 午夜激情欧美在线| 在线观看66精品国产| 88av欧美| 非洲黑人性xxxx精品又粗又长| 内射极品少妇av片p| 男女做爰动态图高潮gif福利片| 久久久色成人| 亚洲欧美日韩东京热| 他把我摸到了高潮在线观看| 两人在一起打扑克的视频| 757午夜福利合集在线观看| 色吧在线观看| e午夜精品久久久久久久| 波野结衣二区三区在线 | 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 深爱激情五月婷婷| 男女那种视频在线观看| 久久国产乱子伦精品免费另类| 久久久久久久精品吃奶| 看片在线看免费视频| 人人妻人人澡欧美一区二区| 免费电影在线观看免费观看| 免费搜索国产男女视频| 亚洲男人的天堂狠狠| 亚洲国产高清在线一区二区三| 亚洲无线在线观看| 在线播放国产精品三级| 日本五十路高清| 99久久无色码亚洲精品果冻| 国产色婷婷99| 男人舔女人下体高潮全视频| 美女 人体艺术 gogo| 757午夜福利合集在线观看| 宅男免费午夜| 色综合欧美亚洲国产小说| 日韩中文字幕欧美一区二区| 舔av片在线| 国产av一区在线观看免费| 日韩免费av在线播放| 一夜夜www| 制服人妻中文乱码| 美女免费视频网站| 在线观看日韩欧美| 久久香蕉精品热| 亚洲狠狠婷婷综合久久图片| 国产aⅴ精品一区二区三区波| 人人妻,人人澡人人爽秒播| 两个人视频免费观看高清| 欧美日韩福利视频一区二区| 免费看光身美女| 88av欧美| 最新在线观看一区二区三区| 人妻夜夜爽99麻豆av| а√天堂www在线а√下载| xxx96com| 在线视频色国产色| 亚洲成av人片免费观看| 婷婷亚洲欧美| 欧美性感艳星| 真人做人爱边吃奶动态| 久久中文看片网| 色视频www国产| 欧美三级亚洲精品| а√天堂www在线а√下载| 天堂网av新在线| 五月伊人婷婷丁香| 女警被强在线播放| 成人特级黄色片久久久久久久| 噜噜噜噜噜久久久久久91| 高清在线国产一区| 欧美av亚洲av综合av国产av| 亚洲国产色片| 夜夜看夜夜爽夜夜摸| 99国产综合亚洲精品| 精品人妻一区二区三区麻豆 | 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 国产三级中文精品| 亚洲成av人片免费观看| 老熟妇乱子伦视频在线观看| 一个人免费在线观看电影| 国产色婷婷99| 亚洲,欧美精品.| 欧美丝袜亚洲另类 | 成人永久免费在线观看视频| 九色成人免费人妻av| 欧美极品一区二区三区四区| 此物有八面人人有两片| 国产精品久久视频播放| 日韩精品中文字幕看吧| 亚洲狠狠婷婷综合久久图片| av专区在线播放| 国产伦精品一区二区三区视频9 | 久久久久久九九精品二区国产| 成人午夜高清在线视频| 欧美乱码精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 欧美日韩亚洲国产一区二区在线观看| 国产精品影院久久| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 国产黄a三级三级三级人| 亚洲国产中文字幕在线视频| 午夜福利18| 一本一本综合久久| 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 久久久久久久久中文| 亚洲,欧美精品.| 亚洲av不卡在线观看| 国产高清视频在线播放一区| 国产亚洲精品久久久com| 国产免费男女视频| 亚洲狠狠婷婷综合久久图片| 免费一级毛片在线播放高清视频| 中文字幕人妻丝袜一区二区| 99热这里只有精品一区| eeuss影院久久| 国产中年淑女户外野战色| 欧美三级亚洲精品| 俺也久久电影网| 免费看光身美女| 午夜福利免费观看在线| 日日夜夜操网爽| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 成人午夜高清在线视频| 伊人久久精品亚洲午夜| 极品教师在线免费播放| 小蜜桃在线观看免费完整版高清| 国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 亚洲片人在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品日韩av片在线观看 | 欧美日韩黄片免| 最近最新免费中文字幕在线| 色老头精品视频在线观看| 黄色视频,在线免费观看| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久,| 国产成人av激情在线播放| 久久99热这里只有精品18| 婷婷丁香在线五月| 日韩欧美精品v在线| 三级毛片av免费| 身体一侧抽搐| 十八禁人妻一区二区| 欧美精品啪啪一区二区三区| 亚洲欧美一区二区三区黑人| 伊人久久精品亚洲午夜| 少妇裸体淫交视频免费看高清| 成人欧美大片| 99在线人妻在线中文字幕| www.熟女人妻精品国产| 黄色女人牲交| 蜜桃久久精品国产亚洲av| 欧美日本亚洲视频在线播放| aaaaa片日本免费| 人妻夜夜爽99麻豆av| 一本久久中文字幕| 在线观看免费视频日本深夜| 一区二区三区免费毛片| 欧美日韩综合久久久久久 | 日韩欧美一区二区三区在线观看| 校园春色视频在线观看| 一区福利在线观看| 精品99又大又爽又粗少妇毛片 | 日韩免费av在线播放| 给我免费播放毛片高清在线观看| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 亚洲无线观看免费| 免费看美女性在线毛片视频| 久久精品亚洲精品国产色婷小说| 国产精品女同一区二区软件 | 窝窝影院91人妻| 亚洲18禁久久av| 天堂影院成人在线观看| 美女黄网站色视频| 最新美女视频免费是黄的| 婷婷六月久久综合丁香| 尤物成人国产欧美一区二区三区| 国产蜜桃级精品一区二区三区| 国产v大片淫在线免费观看| 人人妻人人澡欧美一区二区| 首页视频小说图片口味搜索| 97人妻精品一区二区三区麻豆| 老司机午夜十八禁免费视频| 久久久久亚洲av毛片大全| 黄色片一级片一级黄色片| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频| 午夜两性在线视频| 国产成人影院久久av| 在线观看66精品国产| 国产黄色小视频在线观看| 91在线观看av| 国产黄a三级三级三级人| 一本精品99久久精品77| 在线免费观看的www视频| 精品久久久久久久久久久久久| 最新中文字幕久久久久| 国产久久久一区二区三区| 色噜噜av男人的天堂激情| 亚洲国产色片| 波野结衣二区三区在线 | 天堂影院成人在线观看| 国产av在哪里看| 国产成人av教育| 嫩草影视91久久| 亚洲国产精品sss在线观看| 日韩有码中文字幕| 午夜福利在线在线| 两人在一起打扑克的视频| 亚洲av成人不卡在线观看播放网| 可以在线观看的亚洲视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| 日韩精品青青久久久久久| 国产美女午夜福利| 亚洲精品一区av在线观看| 免费在线观看亚洲国产| 亚洲欧美日韩无卡精品| 亚洲人成网站在线播| 中文字幕人成人乱码亚洲影| 国产精品久久电影中文字幕| 国产精品一及| 精品一区二区三区av网在线观看| 床上黄色一级片| 看片在线看免费视频| 九九热线精品视视频播放| 18+在线观看网站| 极品教师在线免费播放| 欧美日韩瑟瑟在线播放| 国模一区二区三区四区视频| 亚洲中文字幕一区二区三区有码在线看| 午夜精品在线福利| 好男人电影高清在线观看| 一进一出抽搐动态| 亚洲美女黄片视频| 亚洲avbb在线观看| 99久久精品热视频| 欧美大码av| 欧美黄色片欧美黄色片| 啦啦啦观看免费观看视频高清| 亚洲熟妇熟女久久| 中文字幕高清在线视频| 国产伦人伦偷精品视频| 国产精品一区二区三区四区免费观看 | 看免费av毛片| 国产一区二区亚洲精品在线观看| 91九色精品人成在线观看| 超碰av人人做人人爽久久 | 欧美黄色片欧美黄色片| 天堂av国产一区二区熟女人妻| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 亚洲aⅴ乱码一区二区在线播放| 国产精品1区2区在线观看.| 国模一区二区三区四区视频| 精品一区二区三区人妻视频| 国产成人av教育| 人人妻人人澡欧美一区二区| 欧美另类亚洲清纯唯美| 日韩欧美三级三区| 午夜视频国产福利| 亚洲狠狠婷婷综合久久图片| 欧美日韩瑟瑟在线播放| 成人午夜高清在线视频| 一级黄色大片毛片| 午夜福利视频1000在线观看| 又粗又爽又猛毛片免费看| 色视频www国产| 3wmmmm亚洲av在线观看| 尤物成人国产欧美一区二区三区| 国产av不卡久久| 色综合婷婷激情| 性色av乱码一区二区三区2| 亚洲国产精品合色在线| 在线观看午夜福利视频| 老熟妇仑乱视频hdxx| 在线免费观看不下载黄p国产 | 91久久精品电影网| 成年人黄色毛片网站| 日韩欧美 国产精品| 国产成人av教育| 日韩成人在线观看一区二区三区| 日本在线视频免费播放| 久久久久久久久久黄片| 亚洲国产欧美网| 我要搜黄色片| 国语自产精品视频在线第100页| 亚洲欧美日韩卡通动漫| 精品国产亚洲在线| 免费在线观看亚洲国产| 久99久视频精品免费| 久久精品国产综合久久久| 亚洲欧美一区二区三区黑人| 亚洲av日韩精品久久久久久密| 成人精品一区二区免费| 狠狠狠狠99中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 给我免费播放毛片高清在线观看| 国产乱人视频| 精品日产1卡2卡| 3wmmmm亚洲av在线观看| 国产在线精品亚洲第一网站| 亚洲内射少妇av| 亚洲一区二区三区不卡视频| 日本 欧美在线| 国产精品98久久久久久宅男小说| 午夜精品在线福利| 久久欧美精品欧美久久欧美| 欧美成人免费av一区二区三区| 国产蜜桃级精品一区二区三区| 日韩大尺度精品在线看网址| 制服丝袜大香蕉在线| 黄色日韩在线| 可以在线观看的亚洲视频| 少妇人妻一区二区三区视频| 最近最新中文字幕大全电影3| 床上黄色一级片| 特级一级黄色大片| 桃红色精品国产亚洲av| 婷婷精品国产亚洲av在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av美国av| 国产亚洲精品一区二区www| 日本 欧美在线| 亚洲成人久久性| 首页视频小说图片口味搜索| 宅男免费午夜| netflix在线观看网站| 两个人视频免费观看高清| 日韩中文字幕欧美一区二区| xxxwww97欧美| 村上凉子中文字幕在线| 在线播放国产精品三级| 亚洲无线在线观看| 网址你懂的国产日韩在线| 久9热在线精品视频| 欧美av亚洲av综合av国产av| 日本一本二区三区精品| 超碰av人人做人人爽久久 | 我要搜黄色片| 夜夜看夜夜爽夜夜摸| 99热6这里只有精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲在线观看片| 国产精品自产拍在线观看55亚洲| 一级黄色大片毛片| x7x7x7水蜜桃| 搡女人真爽免费视频火全软件 | 久久亚洲精品不卡| 欧美最黄视频在线播放免费| 岛国在线免费视频观看| 两个人的视频大全免费| 成熟少妇高潮喷水视频| 国产精品久久久久久精品电影| 亚洲不卡免费看| 草草在线视频免费看| 亚洲久久久久久中文字幕| 日本五十路高清| 99久久精品国产亚洲精品| 亚洲av成人不卡在线观看播放网| 亚洲成av人片免费观看| 99视频精品全部免费 在线| 欧美在线一区亚洲| 欧美日本视频| 一个人免费在线观看的高清视频| 国产av麻豆久久久久久久| 一二三四社区在线视频社区8| 淫妇啪啪啪对白视频| 男女床上黄色一级片免费看| 怎么达到女性高潮| 国产国拍精品亚洲av在线观看 | 91字幕亚洲| 在线看三级毛片| 亚洲美女黄片视频| 老司机午夜福利在线观看视频| www日本黄色视频网| 欧美成人免费av一区二区三区| 啦啦啦免费观看视频1| 桃色一区二区三区在线观看| 精品久久久久久,| 久久久国产成人精品二区| 亚洲成人中文字幕在线播放| 国产精品久久电影中文字幕| 日韩欧美一区二区三区在线观看| 在线播放国产精品三级| 久久久国产精品麻豆| 国产精品1区2区在线观看.| 少妇人妻精品综合一区二区 | 国产精品三级大全| 又黄又爽又免费观看的视频| 久9热在线精品视频| 国产不卡一卡二| 国产激情偷乱视频一区二区| 亚洲黑人精品在线| 男人舔奶头视频| 中文在线观看免费www的网站| 18禁在线播放成人免费| 欧美日韩黄片免| 色综合欧美亚洲国产小说| 久久香蕉国产精品| 白带黄色成豆腐渣| 亚洲人成伊人成综合网2020| 丰满人妻一区二区三区视频av | 久久精品综合一区二区三区| 久久国产精品影院| 国产97色在线日韩免费| 国产高清视频在线观看网站| 国产毛片a区久久久久| 国产真实乱freesex| 男人舔女人下体高潮全视频| 人人妻人人看人人澡| 18美女黄网站色大片免费观看| 国产男靠女视频免费网站| 日本一本二区三区精品| 成人午夜高清在线视频| 成人国产一区最新在线观看| 制服人妻中文乱码| 51午夜福利影视在线观看| 麻豆久久精品国产亚洲av| 亚洲男人的天堂狠狠| 中文字幕熟女人妻在线| 免费av不卡在线播放| 亚洲人成伊人成综合网2020| 亚洲精品色激情综合| 天天一区二区日本电影三级| 一进一出抽搐动态| 高清日韩中文字幕在线| 成人永久免费在线观看视频| 国产淫片久久久久久久久 | 麻豆国产97在线/欧美| 精品福利观看| 成年免费大片在线观看| 国产精品一及| 成人高潮视频无遮挡免费网站| 狂野欧美白嫩少妇大欣赏| 久久精品国产清高在天天线| 成人永久免费在线观看视频| 在线观看日韩欧美| 亚洲激情在线av| 天天一区二区日本电影三级| 无人区码免费观看不卡| 在线免费观看的www视频| 色老头精品视频在线观看| x7x7x7水蜜桃| www.999成人在线观看| 婷婷亚洲欧美| 两人在一起打扑克的视频| АⅤ资源中文在线天堂| 一级黄色大片毛片| 久久草成人影院| 九九热线精品视视频播放| 国产极品精品免费视频能看的| 内地一区二区视频在线| 乱人视频在线观看| 99久久久亚洲精品蜜臀av| 日韩欧美国产在线观看| 成人午夜高清在线视频|