• <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网| 免费在线观看影片大全网站| 国产麻豆成人av免费视频| 亚洲av中文av极速乱| 久久国内精品自在自线图片| 国产三级在线视频| 亚洲国产精品久久男人天堂| 久久久久久久久大av| 免费人成在线观看视频色| 最新中文字幕久久久久| 99在线视频只有这里精品首页| 精品熟女少妇av免费看| 国产精品久久久久久久久免| 久久久精品欧美日韩精品| 日本黄大片高清| 男女那种视频在线观看| 男女视频在线观看网站免费| 自拍偷自拍亚洲精品老妇| 啦啦啦韩国在线观看视频| 最近2019中文字幕mv第一页| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 日本在线视频免费播放| 深爱激情五月婷婷| 91av网一区二区| 黄色日韩在线| 国产单亲对白刺激| 一区二区三区四区激情视频 | 午夜免费激情av| 哪里可以看免费的av片| 在线播放无遮挡| 欧美成人一区二区免费高清观看| 国产高潮美女av| 美女被艹到高潮喷水动态| 男女之事视频高清在线观看| 国产淫片久久久久久久久| 色噜噜av男人的天堂激情| 欧美激情在线99| 舔av片在线| 久久精品国产亚洲av涩爱 | 无遮挡黄片免费观看| 成人一区二区视频在线观看| 亚洲婷婷狠狠爱综合网| 色噜噜av男人的天堂激情| 欧美绝顶高潮抽搐喷水| 久久久久性生活片| 国产一区二区三区在线臀色熟女| 日本免费一区二区三区高清不卡| 天堂网av新在线| 亚洲国产色片| 男女边吃奶边做爰视频| 无遮挡黄片免费观看| 精品熟女少妇av免费看| 小蜜桃在线观看免费完整版高清| 男人狂女人下面高潮的视频| 在线观看美女被高潮喷水网站| 日本五十路高清| 国产精品1区2区在线观看.| 大又大粗又爽又黄少妇毛片口| 禁无遮挡网站| 男女边吃奶边做爰视频| 婷婷色综合大香蕉| 麻豆久久精品国产亚洲av| 亚洲最大成人中文| 亚洲欧美日韩高清在线视频| 老司机福利观看| 69av精品久久久久久| 草草在线视频免费看| 国产伦精品一区二区三区四那| 18禁在线无遮挡免费观看视频 | 一边摸一边抽搐一进一小说| 国产精品人妻久久久影院| 在线观看av片永久免费下载| 国产精品人妻久久久久久| 亚洲国产精品久久男人天堂| 深爱激情五月婷婷| eeuss影院久久| 国产视频一区二区在线看| 欧美日韩国产亚洲二区| 国产综合懂色| 成人鲁丝片一二三区免费| 国产成人影院久久av| 国产精品精品国产色婷婷| 国产精品av视频在线免费观看| 神马国产精品三级电影在线观看| h日本视频在线播放| 精品久久久久久久久久免费视频| 三级毛片av免费| 狂野欧美激情性xxxx在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲av中文字字幕乱码综合| 国产成人a∨麻豆精品| 免费看光身美女| 麻豆精品久久久久久蜜桃| 俄罗斯特黄特色一大片| avwww免费| 国产高潮美女av| 国产免费男女视频| 秋霞在线观看毛片| 综合色av麻豆| 91午夜精品亚洲一区二区三区| 看片在线看免费视频| 免费看日本二区| 两个人视频免费观看高清| 色噜噜av男人的天堂激情| 久久久午夜欧美精品| 日本免费一区二区三区高清不卡| 欧美国产日韩亚洲一区| 秋霞在线观看毛片| 国产亚洲av嫩草精品影院| 久久综合国产亚洲精品| 亚洲av免费在线观看| 波多野结衣高清无吗| 亚洲专区国产一区二区| 国产在线精品亚洲第一网站| 成年免费大片在线观看| 一个人观看的视频www高清免费观看| 亚洲熟妇中文字幕五十中出| 亚洲av中文av极速乱| 亚洲精品一卡2卡三卡4卡5卡| 欧美一级a爱片免费观看看| 美女黄网站色视频| 99国产极品粉嫩在线观看| 国产精品一区二区免费欧美| 午夜激情欧美在线| 亚洲五月天丁香| 超碰av人人做人人爽久久| 久久久久久久午夜电影| 丰满人妻一区二区三区视频av| www日本黄色视频网| 成人毛片a级毛片在线播放| 亚洲国产精品合色在线| 人人妻人人看人人澡| 一级毛片aaaaaa免费看小| 淫妇啪啪啪对白视频| 精品99又大又爽又粗少妇毛片| 99热6这里只有精品| 乱系列少妇在线播放| 两个人视频免费观看高清| 国产精品福利在线免费观看| 小说图片视频综合网站| 99热精品在线国产| 国产精品嫩草影院av在线观看| 国产成人福利小说| 我的老师免费观看完整版| 亚洲欧美日韩东京热| 欧美成人精品欧美一级黄| 中国国产av一级| 亚洲真实伦在线观看| 婷婷亚洲欧美| 亚洲久久久久久中文字幕| 国产欧美日韩一区二区精品| 亚洲精品粉嫩美女一区| 欧美+亚洲+日韩+国产| 亚洲精品一卡2卡三卡4卡5卡| 一卡2卡三卡四卡精品乱码亚洲| 丝袜美腿在线中文| 日本 av在线| 日本黄大片高清| 在线播放无遮挡| 亚洲不卡免费看| 美女cb高潮喷水在线观看| 乱人视频在线观看| eeuss影院久久| 亚洲精品一区av在线观看| 国产午夜精品论理片| av在线老鸭窝| 99久国产av精品| 非洲黑人性xxxx精品又粗又长| 久久韩国三级中文字幕| 亚洲性夜色夜夜综合| 成年女人永久免费观看视频| 草草在线视频免费看| 在线播放无遮挡| 国产高清三级在线| 国产黄色视频一区二区在线观看 | 国产91av在线免费观看| 亚洲精品色激情综合| 精品人妻偷拍中文字幕| 国产成人freesex在线 | 晚上一个人看的免费电影| 久久九九热精品免费| 美女xxoo啪啪120秒动态图| 我要搜黄色片| 亚洲最大成人av| 精品免费久久久久久久清纯| 国产成人福利小说| 亚洲18禁久久av| 最近视频中文字幕2019在线8| 美女免费视频网站| 白带黄色成豆腐渣| 国产白丝娇喘喷水9色精品| 亚洲成人av在线免费| 国内精品美女久久久久久| 国产熟女欧美一区二区| 亚洲一区二区三区色噜噜| 不卡一级毛片| 亚洲激情五月婷婷啪啪| 一进一出抽搐gif免费好疼| 亚洲国产日韩欧美精品在线观看| av专区在线播放| 久久久久久久久久成人| 亚洲国产精品合色在线| 日本撒尿小便嘘嘘汇集6| 国产日本99.免费观看| 午夜激情欧美在线| 97超级碰碰碰精品色视频在线观看| 国产v大片淫在线免费观看| 国产日本99.免费观看| 国产精品久久久久久久久免| 欧美另类亚洲清纯唯美| 直男gayav资源| 精品久久久久久久久久久久久| 欧美潮喷喷水| 搡老熟女国产l中国老女人| 欧美日韩国产亚洲二区| 亚洲专区国产一区二区| 国产极品精品免费视频能看的| 亚州av有码| 少妇被粗大猛烈的视频| 国产不卡一卡二| av免费在线看不卡| 国产精品人妻久久久影院| 成人特级av手机在线观看| 亚洲无线在线观看| 在线观看美女被高潮喷水网站| 亚洲aⅴ乱码一区二区在线播放| 给我免费播放毛片高清在线观看| 在线免费十八禁| 国产视频内射| 午夜福利成人在线免费观看| 久久久a久久爽久久v久久| 啦啦啦韩国在线观看视频| 亚洲无线观看免费| 婷婷亚洲欧美| 悠悠久久av| 亚洲精品日韩在线中文字幕 | 搡老熟女国产l中国老女人| 黄色日韩在线| 日本免费一区二区三区高清不卡| 一级a爱片免费观看的视频| 国产精品一及| 色噜噜av男人的天堂激情| 真实男女啪啪啪动态图| 亚洲美女视频黄频| 黄色一级大片看看| 日日摸夜夜添夜夜添小说| 国产v大片淫在线免费观看| 特大巨黑吊av在线直播| 欧美高清性xxxxhd video| 久久久久久久亚洲中文字幕| 国产成人a区在线观看| 赤兔流量卡办理| 看免费成人av毛片| 精品熟女少妇av免费看| 日本五十路高清| 村上凉子中文字幕在线| 亚洲第一区二区三区不卡| 午夜福利高清视频| 麻豆久久精品国产亚洲av| 欧美激情久久久久久爽电影| 又黄又爽又免费观看的视频| 国产亚洲精品综合一区在线观看| 亚洲最大成人中文| 在线天堂最新版资源| 日产精品乱码卡一卡2卡三| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久久国产a免费观看| 91久久精品国产一区二区三区| 色播亚洲综合网| 亚洲成人久久爱视频| 成年免费大片在线观看| 亚洲av电影不卡..在线观看| 精品久久久噜噜| 国产伦一二天堂av在线观看| 免费av不卡在线播放| av在线观看视频网站免费| 精品一区二区免费观看| 亚洲av免费在线观看| 国产精品日韩av在线免费观看| 99九九线精品视频在线观看视频| 国产成人aa在线观看| 免费高清视频大片| 国产精品三级大全| 免费看美女性在线毛片视频| 国产高清有码在线观看视频| 在线观看午夜福利视频| 久久久色成人| 成人鲁丝片一二三区免费| 免费看av在线观看网站| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 大香蕉久久网| 国内精品美女久久久久久| 国内少妇人妻偷人精品xxx网站| 男人舔女人下体高潮全视频| 色噜噜av男人的天堂激情| 亚洲国产欧洲综合997久久,| 偷拍熟女少妇极品色| 亚洲av熟女| 麻豆成人午夜福利视频| 在线看三级毛片| 真实男女啪啪啪动态图| 人人妻人人看人人澡| 日韩精品中文字幕看吧| 白带黄色成豆腐渣| 亚洲电影在线观看av| 2021天堂中文幕一二区在线观| 黄片wwwwww| 成年女人永久免费观看视频| 亚洲精品粉嫩美女一区| 国产成人a∨麻豆精品| 热99re8久久精品国产| 最近最新中文字幕大全电影3| 亚洲无线观看免费| av.在线天堂| 国产真实伦视频高清在线观看| 在线免费观看的www视频| 综合色av麻豆| 最近手机中文字幕大全| 免费看av在线观看网站| 天堂网av新在线| 尾随美女入室| 在线天堂最新版资源| 国产男靠女视频免费网站| 久久久久久久久大av| 亚洲18禁久久av| 搡老妇女老女人老熟妇| 亚洲丝袜综合中文字幕| 桃色一区二区三区在线观看| 日韩欧美国产在线观看| 伦精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲精品一区av在线观看| 久久人人爽人人爽人人片va| 日日干狠狠操夜夜爽| 国产精华一区二区三区| 九九爱精品视频在线观看| 伦精品一区二区三区| 国产爱豆传媒在线观看| 亚洲一级一片aⅴ在线观看| а√天堂www在线а√下载| 亚洲国产色片| 色av中文字幕| 中国美白少妇内射xxxbb| 99久久久亚洲精品蜜臀av| 国产精品久久久久久精品电影| 在线观看美女被高潮喷水网站| 一本久久中文字幕| 亚洲av不卡在线观看| 午夜免费男女啪啪视频观看 | 亚洲av.av天堂| 欧美激情在线99| 别揉我奶头~嗯~啊~动态视频| 女的被弄到高潮叫床怎么办| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片| 91av网一区二区| 成人漫画全彩无遮挡| 久久九九热精品免费| 久久精品国产亚洲av天美| 18+在线观看网站| 12—13女人毛片做爰片一| 草草在线视频免费看| 国产精品爽爽va在线观看网站| 一本一本综合久久| 国产久久久一区二区三区| 免费av观看视频| 免费在线观看影片大全网站| 亚洲图色成人| 麻豆一二三区av精品| 搞女人的毛片| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 久久久国产成人精品二区| 免费一级毛片在线播放高清视频| 日本-黄色视频高清免费观看| 久久久久久久久大av| 日韩av在线大香蕉| 成人鲁丝片一二三区免费| 亚洲人成网站在线观看播放| 亚洲欧美精品自产自拍| 国产蜜桃级精品一区二区三区| 天美传媒精品一区二区| a级毛色黄片| 久久久欧美国产精品| 又粗又爽又猛毛片免费看| 亚洲av美国av| 亚洲精品国产av成人精品 | 久久精品人妻少妇| 亚洲熟妇熟女久久| 最近在线观看免费完整版| 国产精品野战在线观看| 日本黄色视频三级网站网址| 免费av毛片视频| 99riav亚洲国产免费| 午夜久久久久精精品| h日本视频在线播放| 色噜噜av男人的天堂激情| 亚洲av熟女| 国产日本99.免费观看| 国产在线男女| 国产精品综合久久久久久久免费| 国产欧美日韩精品亚洲av| 18禁在线无遮挡免费观看视频 | 亚洲精品一区av在线观看| 国产一区二区三区在线臀色熟女| 99热全是精品| 日韩成人伦理影院| 自拍偷自拍亚洲精品老妇| 99九九线精品视频在线观看视频| 久久人人爽人人片av| 看片在线看免费视频| 国产精品久久久久久亚洲av鲁大| 亚洲av成人av| 禁无遮挡网站| 国产av麻豆久久久久久久| 国产高潮美女av| 日韩欧美三级三区| 国产亚洲精品久久久com| 国产av麻豆久久久久久久| 国产高清视频在线播放一区| 精品久久国产蜜桃| 国产精品美女特级片免费视频播放器| 亚洲欧美日韩高清在线视频| 国产成人精品久久久久久| 日韩欧美三级三区| 亚洲欧美中文字幕日韩二区| 国产成人a区在线观看| 高清午夜精品一区二区三区 | 成人av一区二区三区在线看| 狠狠狠狠99中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产午夜福利久久久久久| 日韩精品有码人妻一区| 少妇裸体淫交视频免费看高清| 成人亚洲精品av一区二区| 人人妻人人澡人人爽人人夜夜 | 六月丁香七月| 精品日产1卡2卡| 亚洲专区国产一区二区| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 观看免费一级毛片| 亚洲中文字幕一区二区三区有码在线看| 国产精品人妻久久久影院| 亚洲欧美日韩高清专用| 久久久欧美国产精品| 日韩强制内射视频| 男人的好看免费观看在线视频| 成年av动漫网址| 亚洲人成网站在线播放欧美日韩| 哪里可以看免费的av片| 在线看三级毛片| 国产成人freesex在线 | 啦啦啦观看免费观看视频高清| 少妇丰满av| 夜夜夜夜夜久久久久| 美女高潮的动态| 久久亚洲国产成人精品v| 麻豆国产97在线/欧美| 欧美一级a爱片免费观看看| 亚洲美女视频黄频| 中出人妻视频一区二区| 2021天堂中文幕一二区在线观| 99热这里只有是精品在线观看| 久久精品人妻少妇| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av| 久久韩国三级中文字幕| 免费看光身美女| or卡值多少钱| 看片在线看免费视频| 久久久精品94久久精品| 青春草视频在线免费观看| 亚洲av美国av| 又黄又爽又刺激的免费视频.| 精品乱码久久久久久99久播| 99久久成人亚洲精品观看| 久久精品国产自在天天线| 亚洲欧美精品综合久久99| 日本五十路高清| 97超级碰碰碰精品色视频在线观看| 嫩草影视91久久| 在现免费观看毛片| 国产极品精品免费视频能看的| 亚洲自偷自拍三级| 老司机午夜福利在线观看视频| 久久人人精品亚洲av| 亚洲国产精品成人久久小说 | 欧美日本视频| 97在线视频观看| 看片在线看免费视频| 亚洲av成人av| 日本一本二区三区精品| 国产蜜桃级精品一区二区三区| 又爽又黄无遮挡网站| 亚洲激情五月婷婷啪啪| 日本五十路高清| 日韩欧美免费精品| 内地一区二区视频在线| 99久久九九国产精品国产免费| 久久精品影院6| 欧美丝袜亚洲另类| 国内精品一区二区在线观看| 亚洲一级一片aⅴ在线观看| 能在线免费观看的黄片| 少妇熟女aⅴ在线视频| 亚洲国产欧洲综合997久久,| 久久草成人影院| 亚洲中文字幕一区二区三区有码在线看| 成人漫画全彩无遮挡| 天堂影院成人在线观看| 嫩草影视91久久| 国产精品久久久久久av不卡| a级一级毛片免费在线观看| 男人舔女人下体高潮全视频| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 日本黄色视频三级网站网址| 中文字幕熟女人妻在线| 桃色一区二区三区在线观看| 成人一区二区视频在线观看| 亚洲成人精品中文字幕电影| 成人漫画全彩无遮挡| 桃色一区二区三区在线观看| 国产av在哪里看| 最近2019中文字幕mv第一页| 又黄又爽又免费观看的视频| 久久精品影院6| 欧美国产日韩亚洲一区| 国产男靠女视频免费网站| 天天躁日日操中文字幕| 搡老熟女国产l中国老女人| 亚洲综合色惰| 麻豆av噜噜一区二区三区| 亚洲成人av在线免费| 美女xxoo啪啪120秒动态图| 黄色欧美视频在线观看| 人人妻人人澡欧美一区二区| 亚洲婷婷狠狠爱综合网| 在线国产一区二区在线| 亚洲av成人av| 别揉我奶头 嗯啊视频| 人人妻,人人澡人人爽秒播| 九九热线精品视视频播放| 搡老妇女老女人老熟妇| av卡一久久| 久久6这里有精品| 国产精品爽爽va在线观看网站| 久久久精品94久久精品| 国产亚洲91精品色在线| 日本黄色片子视频| 国产三级中文精品| 亚洲最大成人手机在线| 人人妻人人澡人人爽人人夜夜 | 成人亚洲精品av一区二区| 又粗又爽又猛毛片免费看| 午夜福利视频1000在线观看| 亚洲色图av天堂| 久久久精品欧美日韩精品| 大又大粗又爽又黄少妇毛片口| 日韩精品有码人妻一区| 成人二区视频| 亚洲精品日韩在线中文字幕 | 99九九线精品视频在线观看视频| 成人一区二区视频在线观看| 九九久久精品国产亚洲av麻豆| 亚洲一区高清亚洲精品| 男人舔奶头视频| 好男人在线观看高清免费视频| 亚洲成人av在线免费| av卡一久久| 色综合色国产| 免费大片18禁| 久久亚洲国产成人精品v| 麻豆国产97在线/欧美| 久久久精品大字幕| eeuss影院久久| 免费看av在线观看网站| 国产免费一级a男人的天堂| 性色avwww在线观看| 免费无遮挡裸体视频| 欧美一区二区亚洲| 日韩在线高清观看一区二区三区| 天堂√8在线中文| 欧美一级a爱片免费观看看| 少妇丰满av| 欧美成人免费av一区二区三区| 日产精品乱码卡一卡2卡三| 少妇人妻精品综合一区二区 | 亚洲av电影不卡..在线观看| 18禁黄网站禁片免费观看直播| 国产精品一二三区在线看| 中文字幕av成人在线电影| 国产真实乱freesex| 成年女人看的毛片在线观看| 一本久久中文字幕| 插逼视频在线观看| 俺也久久电影网| 三级男女做爰猛烈吃奶摸视频| 久久久精品欧美日韩精品| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 男女那种视频在线观看| 久久久成人免费电影| 日韩欧美国产在线观看| 丰满人妻一区二区三区视频av| 久久久久久久久久久丰满| 九九久久精品国产亚洲av麻豆| 亚洲久久久久久中文字幕| 一区福利在线观看| 欧美中文日本在线观看视频| 网址你懂的国产日韩在线|