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

    基于WRF/UCM的城市氣候高分辨率數(shù)值模擬研究

    2016-10-12 02:10:55
    關(guān)鍵詞:城市熱島熱島氣象站

    郭   飛

    ( 大連理工大學(xué) 建筑與藝術(shù)學(xué)院, 遼寧 大連 116024 )

    ?

    基于WRF/UCM的城市氣候高分辨率數(shù)值模擬研究

    郭 飛*

    ( 大連理工大學(xué) 建筑與藝術(shù)學(xué)院, 遼寧 大連116024 )

    以沿海城市大連2015年7月3日2時(shí)至7月4日2時(shí)的晴朗高溫天氣為研究背景,利用WRF/UCM以1 km的高分辨率進(jìn)行了城市氣候數(shù)值模擬.結(jié)果表明:海風(fēng)、山體對(duì)城市熱島有明顯的緩解作用;中心城區(qū)受山地影響,形成了管道風(fēng)效應(yīng);模擬期間城市熱島中心位于大連開(kāi)發(fā)區(qū)一帶,熱島強(qiáng)度最高達(dá)6.9 ℃,出現(xiàn)在2015年7月3日19時(shí)30分.將模擬結(jié)果與兩個(gè)氣象站和現(xiàn)場(chǎng)觀測(cè)值進(jìn)行了對(duì)比,發(fā)現(xiàn)模擬值與觀測(cè)值變化趨勢(shì)基本一致,地面溫度場(chǎng)模擬結(jié)果優(yōu)于風(fēng)場(chǎng)模擬結(jié)果.溫度場(chǎng)誤差最大-3.10 ℃,最小0.08 ℃;平均誤差最大0.88 ℃,最小0.52 ℃;均方根誤差最大1.23 ℃,最小0.63 ℃.對(duì)于大尺度的城市氣候研究而言誤差在可接受的范圍之內(nèi).對(duì)相應(yīng)的緩解城市熱島的城市規(guī)劃策略也進(jìn)行了初步探討.

    WRF/UCM;城市氣候;數(shù)值模擬;大連

    0 引 言

    中國(guó)目前城市居民已占人口總數(shù)的54%,城市面臨著巨大的環(huán)境壓力.良好的城市氣候是營(yíng)造生態(tài)環(huán)境的重要基礎(chǔ).城市規(guī)劃應(yīng)當(dāng)與城市氣候評(píng)估相結(jié)合,采取科學(xué)措施推動(dòng)城市建設(shè)順應(yīng)氣候規(guī)律.國(guó)外許多城市的成功經(jīng)驗(yàn)表明,這樣能夠有效緩解城市氣候惡化和污染問(wèn)題.

    城市氣候成因復(fù)雜,既受大尺度背景氣候、局地下墊面材質(zhì)影響,也受建筑形態(tài)和能耗情況影響,如何科學(xué)準(zhǔn)確分析城市氣候特征是進(jìn)行城市氣候規(guī)劃的關(guān)鍵環(huán)節(jié)[1].氣象研究模型WRF(weather research and forecasting model),是由美國(guó)環(huán)境預(yù)測(cè)中心(NCEP)、美國(guó)國(guó)家大氣研究中心(NCAR)等科研機(jī)構(gòu)于2000年起共同參與開(kāi)發(fā)的新一代中尺度預(yù)報(bào)系統(tǒng)[2].WRF具有完善的物理方案描述各類復(fù)雜的氣候現(xiàn)象,在全世界范圍內(nèi)廣泛用于天氣預(yù)報(bào)和氣象研究.近年來(lái)利用WRF進(jìn)行城市局地氣候特征分析,并據(jù)此提出城市建設(shè)、規(guī)劃和管理的政策建議及改善策略,逐漸成為研究的熱點(diǎn)之一[3-5].WRF可以耦合Noah陸地下墊面、城市冠層模型(如UCM、BEP和BEM),甚至與微尺度CFD模型耦合來(lái)模擬城市氣候,詳細(xì)描述城市空間形態(tài)的特征、建筑信息、人為產(chǎn)熱的數(shù)據(jù)等,為理解和評(píng)估城市氣候提供了有力的工具[6-7].許多學(xué)者研究了使用WRF和高分辨率的城市信息庫(kù)數(shù)據(jù)模擬城市氣候.例如伍見(jiàn)軍等利用333 m分辨率的GIS數(shù)據(jù)針對(duì)重慶的城市氣候進(jìn)行了研究[8].蒙偉光等利用更新的GIS土地利用資料和城市冠層模型針對(duì)廣州的一次高溫天氣進(jìn)行了數(shù)值實(shí)驗(yàn)[9].Lin等利用MODIS衛(wèi)星1 km分辨率的土地利用數(shù)據(jù)對(duì)中國(guó)臺(tái)灣北部地區(qū)的城市熱島效應(yīng)進(jìn)行了分析[10].Kusaka等將平板模型和UCM的模擬結(jié)果與氣象觀測(cè)數(shù)據(jù)對(duì)比,發(fā)現(xiàn)UCM能夠模擬出熱島效應(yīng)而平板模型沒(méi)有,平板模型比實(shí)際觀測(cè)值最大偏離幅度達(dá)6.2 ℃[11].這些研究表明,WRF可以成功模擬城市熱島效應(yīng)、海風(fēng)環(huán)流等城市氣候現(xiàn)象,其模擬結(jié)果與觀測(cè)數(shù)據(jù)吻合也較好.WRF應(yīng)用于城市總體尺度的氣候評(píng)估具有較好的優(yōu)勢(shì).

    大連是典型的沿海山地城市,自然條件良好,溫度適宜,是北方的旅游熱門城市.

    過(guò)去60余年大連城區(qū)面積逐步擴(kuò)張,從50 km2(1950年)擴(kuò)展為387.7 km2(2009年).由于城市范圍擴(kuò)張、土地高強(qiáng)度開(kāi)發(fā)利用,在中心城區(qū)形成了高層、高密度的形態(tài)特征.過(guò)去60余年城市中心區(qū)年均氣溫逐步升高、風(fēng)速降低、空氣污染逐漸加重.根據(jù)環(huán)境監(jiān)測(cè)數(shù)據(jù),2014年大連空氣質(zhì)量為優(yōu)的天數(shù)僅為69 d[12].利用科學(xué)的工具分析大連城市擴(kuò)張帶來(lái)的氣候影響,據(jù)此提出相應(yīng)的規(guī)劃策略也成為當(dāng)務(wù)之急.

    本文以2015年7月3日為例,利用WRF模擬大連城市溫度場(chǎng)和風(fēng)場(chǎng)特征,并與氣象站和地面現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)進(jìn)行對(duì)比,為了解城市氣候特征并提出規(guī)劃策略提供一種全新思路.

    1 城市下墊面信息處理

    伴隨著城市化的過(guò)程,城市下墊面類型發(fā)生了很大的轉(zhuǎn)變.然而城市土地資料的更新往往跟不上改變的速率,許多研究結(jié)果表明,如果不根據(jù)實(shí)際情況加以修改,就會(huì)造成某些相關(guān)變量產(chǎn)生較大的模擬偏差[13-14].目前WRF提供的土地資料有兩套,一為美國(guó)地質(zhì)調(diào)查局(USGS)1993年全球24種類型的30 s(約50 km)分辨率數(shù)據(jù),二為基于中分辨率成像光譜儀(MODIS)的2001年全球20種類型的15 s(約25 km)分辨率數(shù)據(jù),雖然整體分類精度和質(zhì)量較高,但在許多地區(qū)仍有較為明顯的分類錯(cuò)誤.本文在大連市土地利用現(xiàn)狀圖的基礎(chǔ)之上,利用GIS平臺(tái)對(duì)不符合現(xiàn)狀的地方進(jìn)行了修正.然后將GIS土地利用數(shù)據(jù)處理成WRF的靜態(tài)地理場(chǎng),替換默認(rèn)的USGS和MODIS 土地利用數(shù)據(jù),為WRF數(shù)值試驗(yàn)提供更準(zhǔn)確的城市覆蓋范圍和土地利用邊界條件,修正前后對(duì)比如圖1所示.?dāng)?shù)據(jù)修正和處理流程如圖2所示.

    2 模擬期間背景天氣

    由于靜穩(wěn)天氣有利于城市熱島的形成,適合用來(lái)研究WRF/UCM的模擬效果,選取大連受高氣壓控制的高溫、晴朗的氣象日,即2015年7月2日18時(shí)~2015年7月4日2時(shí),共計(jì)32 h作為模擬的時(shí)間段(圖3).由于前8 h為模型物理過(guò)程的調(diào)整適應(yīng)時(shí)間,本文僅討論后24 h的模擬結(jié)果.為了對(duì)比驗(yàn)證模型設(shè)置和物理方案的有效性,在西尖山和大連理工大學(xué)(以下簡(jiǎn)稱理工)進(jìn)行了現(xiàn)場(chǎng)固定觀測(cè),并獲取了周水子國(guó)際交換站(以下簡(jiǎn)稱周水子)的逐時(shí)觀測(cè)結(jié)果(表1).根據(jù)觀測(cè)結(jié)果,模擬時(shí)段內(nèi)大連區(qū)域風(fēng)向?yàn)楸逼?,風(fēng)速0.2~7.0 ms-1,溫度19.0~25.7 ℃(圖4).

    圖1 大連地區(qū)用地類型分布圖

    圖2土地利用數(shù)據(jù)向WRF靜態(tài)地理場(chǎng)的轉(zhuǎn)換流程

    Fig.2ProcessofconvertinglandusedatatoWRFgeogriddata

    圖32015年7月3日22時(shí)海平面大氣壓等值線圖

    Fig.3Seasurfacepressurecontourat22:00July3rd, 2015

    表1 地面氣象站觀測(cè)數(shù)據(jù)簡(jiǎn)況

    圖4 2015年7月3日12時(shí)背景風(fēng)場(chǎng)及溫度場(chǎng)

    3 模擬參數(shù)設(shè)置

    模擬區(qū)域采用嵌套網(wǎng)格的設(shè)置.嵌套網(wǎng)格的優(yōu)勢(shì)是既可以提高核心關(guān)注區(qū)域的模擬精度,同時(shí)也能充分考慮周邊的氣候影響并節(jié)約計(jì)算資源.本次采用三層嵌套網(wǎng)格,第一層格點(diǎn)數(shù)99×99,網(wǎng)格間距9 km×9 km;第二層格點(diǎn)數(shù)168×168,網(wǎng)格間距3 km×3 km;第三層格點(diǎn)數(shù)156×156,網(wǎng)格間距1 km×1 km.垂直方向設(shè)置54層,總高度16 km.網(wǎng)格中心點(diǎn)坐標(biāo)38.923°N,121.647°E.模擬計(jì)算區(qū)域基本覆蓋了大部分的大連城市建成區(qū)(圖5).

    驅(qū)動(dòng)WRF的氣象數(shù)據(jù)是美國(guó)NCEP的全球3 h 0.25°×0.25°再分析氣象資料,由于時(shí)空分辨率較高,適合用于高分辨率的城市氣候模擬.模擬積分步長(zhǎng)為54 s.本研究采用的物理方案如表2所示.

    圖5 嵌套網(wǎng)格設(shè)置與地形

    表2 WRF的物理方案設(shè)置

    4 模擬結(jié)果分析及驗(yàn)證

    4.1溫度場(chǎng)

    對(duì)城市氣候而言,地面2 m的溫度場(chǎng)對(duì)人體熱舒適影響最重要.地面2 m溫度場(chǎng)變化結(jié)果顯示,模擬期間城市中心區(qū)最高氣溫達(dá)到26.5 ℃,出現(xiàn)在2015年7月3日13時(shí)30分和14時(shí)30分,位于金州區(qū)一帶;7月3日15時(shí)30分起溫度開(kāi)始下降(圖6).夜間最低溫度為13.6 ℃,出現(xiàn)在7月3日4時(shí)30分,位于普蘭店一帶.模擬期間城市熱島最強(qiáng)時(shí)刻出現(xiàn)在2015年7月3日夜間19時(shí)30分(圖7).此時(shí)所有的城市建成區(qū)無(wú)論大小均存在較強(qiáng)的熱島效應(yīng),最高強(qiáng)度達(dá)6.9 ℃,位于大連開(kāi)發(fā)區(qū)一帶.形成的原因主要是風(fēng)受到開(kāi)發(fā)區(qū)北面山地的阻擋,因此城市內(nèi)部蓄積的熱量難以有效地散發(fā)出去.

    從模擬期間溫度場(chǎng)的變化可以看出,白天陸地附近氣溫升高速度快于海面,但是海風(fēng)對(duì)陸地及城市的降溫效果十分明顯.例如2015年7月3日13時(shí)30分,大連北側(cè)渤海海面2 m空氣溫度約為17.7 ℃,北風(fēng)將涼爽空氣吹向陸地,形成較強(qiáng)的降溫作用,渤海岸線附近的氣溫比城市中心區(qū)低5 ℃左右.越往陸地深處氣溫越高,大連城市中心區(qū)、旅順口、開(kāi)發(fā)區(qū)等比普蘭店、金州等地溫度低2~3 ℃.受海面涼爽的風(fēng)影響,城市熱島中心向下風(fēng)向偏移,黃海區(qū)域大連灣附近海面氣溫受上風(fēng)向城市的影響比周邊海面高.

    圖62015年7月3日15時(shí)30分風(fēng)場(chǎng)及溫度場(chǎng)(最內(nèi)層模擬域)

    Fig.6Thewindfieldandtemperaturefieldat15:30July3rd, 2015 (theinnerdomain)

    圖72015年7月3日19時(shí)30分溫度場(chǎng)(最內(nèi)層模擬域)

    Fig.7Thetemperaturefieldat19:30July3rd, 2015 (theinnerdomain)

    為了研究溫度場(chǎng)在高度上的變化,繪制3個(gè)不同剖面,剖切位置如圖8所示.溫度剖面顯示陸地地形及城市下墊面對(duì)大氣產(chǎn)生了較強(qiáng)的影響,地面200 m以下尤為顯著(圖9).A-A剖面熱島中心為大連開(kāi)發(fā)區(qū),模擬期間處于北面山地形成的風(fēng)影區(qū)之內(nèi),存在較強(qiáng)的城市熱島效應(yīng),熱島強(qiáng)度達(dá)3 ℃.B-B剖面熱島中心位于大連甘井子區(qū)、中山區(qū)等處,受海風(fēng)影響熱島強(qiáng)度不高,為1~2 ℃.C-C剖面主要位于高新區(qū)及旅順口區(qū)之間,山體及綠化覆蓋率高,熱島效應(yīng)也在1~2 ℃.

    圖8 最內(nèi)層模擬域地形及溫度場(chǎng)剖面位置

    圖9 2015年7月3日13時(shí)30分溫度場(chǎng)剖面

    4.2風(fēng)場(chǎng)

    對(duì)城市風(fēng)場(chǎng),一般主要評(píng)估地面10 m的風(fēng)場(chǎng)分布情況.風(fēng)速和熱島強(qiáng)度及污染物濃度呈反比關(guān)系,風(fēng)速越大熱島強(qiáng)度和污染物濃度越低[15].模擬期間大連市區(qū)風(fēng)向?yàn)楸逼黠L(fēng),與氣象站觀測(cè)結(jié)果一致.

    模擬結(jié)果顯示,7月3日夜間城市風(fēng)速較低,為1.5~2.5 m·s-1.風(fēng)速最小值1.58 m·s-1,出現(xiàn)在7月4日1時(shí).7月3日9時(shí)至17時(shí)風(fēng)速較高,相應(yīng)地城市中心熱島強(qiáng)度較低.風(fēng)速最大值5.83 m·s-1,出現(xiàn)在7月3日17時(shí).相對(duì)于溫度場(chǎng),風(fēng)場(chǎng)受地形影響更大,例如金州區(qū)和甘井子區(qū)之間兩處山體地勢(shì)較高,在兩者之間形成了管道風(fēng)效應(yīng),風(fēng)速有加強(qiáng)的趨勢(shì)(圖6).受此影響甘井子區(qū)和中山區(qū)的城市熱島強(qiáng)度也有所減弱.

    4.3模擬結(jié)果驗(yàn)證

    為了驗(yàn)證WRF模擬結(jié)果的準(zhǔn)確性,將WRF數(shù)據(jù)與氣象站和現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)加以對(duì)比分析.于2015年7月3日至7月7日進(jìn)行了現(xiàn)場(chǎng)觀測(cè),從8時(shí)到18時(shí)對(duì)西尖山進(jìn)行人工現(xiàn)場(chǎng)連續(xù)觀測(cè),每15 min記錄一組氣象參數(shù),18時(shí)之后的夜間觀測(cè)由放置現(xiàn)場(chǎng)的儀器完成.另外獲取了周水子氣象站、理工氣象站在此期間的逐時(shí)氣象數(shù)據(jù).

    將WRF溫度場(chǎng)模擬結(jié)果分別與現(xiàn)場(chǎng)實(shí)測(cè)和氣象站結(jié)果對(duì)比,可以發(fā)現(xiàn)模擬與觀測(cè)值之間的變化趨勢(shì)基本吻合,且溫度場(chǎng)模擬結(jié)果優(yōu)于風(fēng)場(chǎng)(圖10、11).溫度場(chǎng)觀測(cè)值與模擬值最大誤差-3.10 ℃,最小誤差0.08 ℃,平均誤差0.52~0.88 ℃,均方根誤差0.63~1.23 ℃.風(fēng)速模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)及氣象站觀測(cè)數(shù)據(jù)吻合度低于溫度場(chǎng),最大誤差-5.10 m·s-1,最小誤差0.01 m·s-1,平均誤差1.03~2.19 m·s-1,均方根誤差1.37~2.47 m·s-1(表3).對(duì)于尺度較大的城市而言,模擬結(jié)果尤其是溫度場(chǎng)的準(zhǔn)確性是較為理想的,用于分析城市溫度場(chǎng)、風(fēng)場(chǎng)的基本變化規(guī)律具有較高的可信度.

    圖10地面2 m溫度模擬結(jié)果與氣象站觀測(cè)對(duì)比

    Fig.10Thecomparisonof2msurfacetemperaturesimulationtothestationobserveddata

    圖11地面10 m風(fēng)速模擬結(jié)果與氣象站觀測(cè)對(duì)比

    Fig.11Thecomparisonof10msurfacewindspeedsimulationtothestationobserveddata

    表3 模擬與觀測(cè)數(shù)據(jù)的誤差分析

    5 應(yīng)用前景

    每個(gè)城市的氣候現(xiàn)象都有其獨(dú)特的特征,既受地理?xiàng)l件和背景氣候的影響,也受城市內(nèi)部用地布局和建筑形態(tài)的影響,因此開(kāi)展城市氣候的研究難度較高.而WRF通過(guò)采用豐富完善的物理方案、嵌套網(wǎng)格、真實(shí)的地理和氣象數(shù)據(jù)、建筑和城市形態(tài)參數(shù)等設(shè)置,可以較為準(zhǔn)確地獲取不同時(shí)空分辨率的氣候信息,為研究城市氣候的科學(xué)規(guī)律提供了極大的方便.WRF的應(yīng)用前景較為廣泛,例如應(yīng)用于將規(guī)劃、氣候和地理研究相結(jié)合的城市氣候圖當(dāng)中.德國(guó)魯爾工業(yè)區(qū)20世紀(jì)70年代曾深受大氣污染的困擾,最早開(kāi)始進(jìn)行城市氣候圖研究,采取規(guī)劃措施促進(jìn)弱風(fēng)條件下空氣流動(dòng),緩解大氣污染.通過(guò)編制氣候分析圖,直接指導(dǎo)了許多城市的土地利用規(guī)劃.其內(nèi)容包括:城市氣候分區(qū)(綠地、開(kāi)敞空間、城市功能區(qū)等)、各類污染源(交通、燃?xì)狻⒐I(yè)等)、風(fēng)環(huán)境信息(不同區(qū)域和季節(jié)的風(fēng)玫瑰、地形風(fēng)、管道風(fēng)、風(fēng)屏障等)[16](圖12).WRF獲取的城市氣候數(shù)據(jù)分辨率可從數(shù)百米到數(shù)百千米,將在氣候圖研究中起到重要作用.

    圖12 柏林地區(qū)氣候分析圖局部(2014年,1∶30 000)

    6 結(jié) 語(yǔ)

    本文用WRF對(duì)大連城市進(jìn)行高分辨率的氣候模擬,其結(jié)果與觀測(cè)結(jié)果吻合得較好.模擬期間城市熱島效應(yīng)強(qiáng)度最高時(shí)達(dá)6.9 ℃,但是海風(fēng)對(duì)城市熱島起到了良好的緩解作用.因此大連城市規(guī)劃應(yīng)當(dāng)充分利用黃海和渤海海風(fēng)資源,為城市夏季和污染天氣提供清潔涼爽的空氣.例如保護(hù)海岸線、降低岸線開(kāi)發(fā)強(qiáng)度和密度,規(guī)劃河流沿岸景觀帶、建設(shè)低密度區(qū)并建設(shè)城市通風(fēng)道,將海面上的清潔新鮮空氣引導(dǎo)進(jìn)入城市內(nèi)部,有助于緩解熱島效應(yīng)和大氣污染.

    山體對(duì)城市熱島具有良好的緩解作用.由于一般山體綠化較好,并且氣溫會(huì)隨著高度上升而下降,一定高度的山丘將會(huì)形成冷氣團(tuán)的聚集區(qū),為周邊城市環(huán)境提供山谷風(fēng)和冷空氣.城市規(guī)劃應(yīng)當(dāng)對(duì)此加以保護(hù)和利用,將山體與綠帶、公園相聯(lián)系,有助于分割城市熱島、避免其蔓延擴(kuò)張,對(duì)改善城市環(huán)境具有較大的作用.

    利用WRF模擬結(jié)果對(duì)城市總體溫度場(chǎng)和風(fēng)場(chǎng)變化規(guī)律加以總結(jié),結(jié)合地形分析、土地利用、城市規(guī)劃管理政策等,可以繪制城市氣候圖,并提出改善城市風(fēng)環(huán)境、減緩城市熱島效應(yīng)和空氣污染的策略.

    [1]郭 飛,祝培生,段棟文,等. 高密度城市氣候評(píng)估方法與應(yīng)用[J]. 西部人居環(huán)境學(xué)刊, 2015, 30(6):19-23.

    GUO Fei, ZHU Pei-sheng, DUAN Dong-wen,etal. High-density urban climate evaluation method and application [J]. Journal of Human Settlements in West China, 2015, 30(6):19-23. (in Chinese)

    [2]Michalakes J, Dudhia J, Gill D,etal. Design of a next-generation regional weather research and forecast model [J]. Towards Teracomputing, 1998:117-124.

    [3]GUO Fei, WANG Shi-yuan, ZHU Pei-sheng. Numerical urban climate simulation and preliminary planning strategies of Dalian city based on WRF model [C] // The 9th International Conference of Urban Climate. Toulouse:ICUC, 2015.

    [4]Kong S F, Li L, Li X X,etal. The impacts of firework burning at the Chinese Spring Festival on air quality:Insights of tracers, source evolution and aging processes [J]. Atmospheric Chemistry and Physics, 2015, 15(4):2167-2184.

    [5]Lu X, Pan N, An H. Comparison of traffic noise spatial distribution based on road network - a case study of typical area in Dalian city, China [C] // Proceedings of the 12th Western Pacific Acoustic Conference 2015. Singapore:WESPAC, 2015.

    [6]CHEN Fei, Kusaka H, Bornstein R,etal. The integrated WRF/urban modelling system:development, evaluation, and applications to urban environmental problems [J]. International Journal of Climatology, 2011, 31(2):273-288.

    [7]Tewari M, Kusaka H, CHEN Fei,etal. Impact of coupling a microscale computational fluid dynamics model with a mesoscale model on urban scale contaminant transport and dispersion [J]. Atmospheric Research, 2010, 96(4):656-664.

    [8]伍見(jiàn)軍,王詠薇,朱 彬,等. WRF 模式中城市冠層參數(shù)化方案在重慶氣象環(huán)境模擬中的性能比較[J]. 長(zhǎng)江流域資源與環(huán)境, 2013, 22(12):1627-1634.

    WU Jian-jun, WANG Yong-wei, ZHU Bin,etal. Performance comparison of different urban canopy schemes in WRF model under Chongqing meteorological simulation [J]. Resources and Environment in the Yangtze Basin, 2013, 22(12):1627-1634. (in Chinese)

    [9]蒙偉光,張艷霞,李江南,等. WRF/UCM在廣州高溫天氣及城市熱島模擬研究中的應(yīng)用[J]. 熱帶氣象學(xué)報(bào), 2010, 26(3):273-282.

    MENG Wei-guang, ZHANG Yan-xia, LI Jiang-nan,etal. Application of WRF/UCM in the simulation of a heat wave event and urban heat island around Guangzhou city [J]. Journal of Tropical Meteorology, 2010, 26(3):273-282. (in Chinese)

    [10]LIN Chuan-yao, CHEN Fei, Huang J C,etal. Urban heat island effect and its impact on boundary layer development and land-sea circulation over northern Taiwan [J]. Atmospheric Environment, 2008, 42(22):5635-5649.

    [11]Kusaka H, Chen F, Tewari M,etal. Numerical simulation of urban heat island effect by the WRF model with 4-km grid increment:an inter-comparison study between the urban canopy model and slab model [J]. Journal of the Meteorological Society of Japan, 2012, 90B:33-45.

    [12]中國(guó)環(huán)境監(jiān)測(cè)總站. 2014年大連市區(qū)空氣質(zhì)量概述[EB/OL]. (2015-05-08)[2015-12-15]. http://www.cnemc.cn/publish/totalWebSite/news/news_44673.html.

    China Environmental Monitoring Terminal. Dalian urban area air quality summary[EB/OL]. (2015-05-08)[2015-12-15]. http://www.cnemc.cn/publish/totalWebSite/news/news_44673.html. (in Chinese)

    [13]LIN Chuan-yao, Chen Wan-chin, CHANG Pao-liang,etal. Impact of the urban heat island effect on precipitation over a complex geographic environment in northern Taiwan [J]. Journal of Applied Meteorology and Climatology, 2011, 50(2):339-353.

    [14]Sertel E, Robock A, Ormeci C. Impacts of land cover data quality on regional climate simulations [J]. International Journal of Climatology, 2010, 30(13):1942-1953.

    [15]Morris C J G, Simmonds I, Plummer N. Quantification of the influences of wind and cloud on the nocturnal urban heat island of a large city [J]. Journal of Applied Meteorology, 2001, 40(2):169-182.

    [16]The Association of German Engineers. Environmental Meteorology-Climate and Air Pollution Maps for Cities and Regions:VDI 3787 [S]. Berlin:The Association of German Engineers, 2014.

    Urban climate high resolution numerical simulation research based on WRF/UCM

    GUOFei*

    ( School of Architecture and Fine Art, Dalian University of Technology, Dalian 116024, China )

    Choosing a sunny and high temperature day from 2:00, 3 July to 2:00, 4 July in 2015 as the research background, a 1 km high resolution numerical simulation for the urban climate in the coastal city of Dalian was implemented by WRF/UCM. The experimental results show that the wind from the sea and mountain can significantly moderate the UHI (urban heated island); the wind funnel effect is found in the urban center area caused by the mountain terrain. The UHI center with the maximum intensity reaching 6.9 ℃ is located in the Dalian Development Zone, appeared at 19:30 on 3 July, 2015. The simulation results are compared with the data from two weather stations and the field observation, it is found that the variation tendency of them is consistent. The accuracy of the surface temperature field is better than that of the wind field. The maximum bias of temperature is -3.10 ℃, the minimum is 0.08 ℃;the maximum mean bias of temperature is 0.88 ℃, the minimum is 0.52 ℃; the maximum RSME of temperature is 1.23 ℃, the minimum is 0.63 ℃. For a large-scale urban climate simulation, the deviation is within the acceptable range. The correspondent UHI moderation strategies for urban planning are also analyzed.

    WRF/UCM; urban climate; numerical simulation; Dalian

    1000-8608(2016)05-0502-08

    2016-01-25;

    2016-07-12.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51278078,51308087);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(DUT13RW306).

    郭 飛*(1980-),男,博士,副教授,E-mail:guofei@dlut.edu.cn.

    TU111.19

    A

    10.7511/dllgxb201605010

    猜你喜歡
    城市熱島熱島氣象站
    珠峰上架起世界最高氣象站
    讓城市“發(fā)燒”的熱島效應(yīng)
    基于CiteSpace的城市熱島研究進(jìn)展與熱點(diǎn)分析
    四川建筑(2021年2期)2021-06-04 07:50:20
    基于Modis數(shù)據(jù)的山東省夏季城市熱島效應(yīng)時(shí)空變化分析
    死海不“死”等
    心靈氣象站
    熱熱島真是個(gè)好地方
    探討傳統(tǒng)機(jī)房局部熱島解決方案
    自動(dòng)氣象站應(yīng)該注意的一些防雷問(wèn)題
    自動(dòng)氣象站常見(jiàn)故障判斷與維護(hù)
    河南科技(2014年12期)2014-02-27 14:10:40
    午夜福利一区二区在线看| 一级作爱视频免费观看| av片东京热男人的天堂| 黑人欧美特级aaaaaa片| 一区二区三区激情视频| 日韩精品免费视频一区二区三区| 大陆偷拍与自拍| 久久午夜亚洲精品久久| 久久人妻av系列| 97人妻天天添夜夜摸| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 黄色片一级片一级黄色片| 丰满迷人的少妇在线观看| 97人妻天天添夜夜摸| 精品国产一区二区三区四区第35| 久久中文字幕一级| 中国美女看黄片| 99国产精品一区二区蜜桃av| 国产熟女xx| 一区二区三区国产精品乱码| 精品电影一区二区在线| 露出奶头的视频| 日韩有码中文字幕| 国产激情欧美一区二区| 乱人伦中国视频| 国产精品九九99| 国产av又大| 精品第一国产精品| 国产成人欧美| 露出奶头的视频| 日本一区二区免费在线视频| 咕卡用的链子| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 国产av在哪里看| 国产免费现黄频在线看| 日韩一卡2卡3卡4卡2021年| 看黄色毛片网站| 亚洲全国av大片| 亚洲精品粉嫩美女一区| 亚洲欧美激情综合另类| 两性午夜刺激爽爽歪歪视频在线观看 | 999久久久精品免费观看国产| 午夜a级毛片| 亚洲一码二码三码区别大吗| 国产成人啪精品午夜网站| 欧美国产精品va在线观看不卡| 亚洲欧美日韩无卡精品| av网站免费在线观看视频| 成熟少妇高潮喷水视频| 一级a爱视频在线免费观看| 日韩欧美一区视频在线观看| 亚洲五月婷婷丁香| 黄色毛片三级朝国网站| 最好的美女福利视频网| 午夜亚洲福利在线播放| 99久久综合精品五月天人人| 国产精品一区二区精品视频观看| 一a级毛片在线观看| 欧美成人性av电影在线观看| 高清av免费在线| 亚洲自偷自拍图片 自拍| 成年女人毛片免费观看观看9| 91字幕亚洲| 身体一侧抽搐| 男女午夜视频在线观看| 免费少妇av软件| 成人免费观看视频高清| 天堂中文最新版在线下载| 欧美精品一区二区免费开放| 欧美成人免费av一区二区三区| 人妻丰满熟妇av一区二区三区| 中文字幕精品免费在线观看视频| 精品一区二区三区视频在线观看免费 | av天堂在线播放| 午夜福利在线免费观看网站| 欧美另类亚洲清纯唯美| 色老头精品视频在线观看| 国产亚洲av高清不卡| 亚洲一区二区三区不卡视频| 日韩人妻精品一区2区三区| 男人的好看免费观看在线视频 | 丝袜人妻中文字幕| 热99re8久久精品国产| 美女扒开内裤让男人捅视频| 日韩大尺度精品在线看网址 | 国产99久久九九免费精品| 91在线观看av| 国产精品影院久久| 久久久久久久久久久久大奶| 精品国产乱子伦一区二区三区| 悠悠久久av| 国产亚洲欧美在线一区二区| 十分钟在线观看高清视频www| 欧美人与性动交α欧美软件| 久久久久久久午夜电影 | 一级毛片高清免费大全| 精品人妻1区二区| 午夜福利免费观看在线| 国产亚洲欧美98| 一本大道久久a久久精品| 免费少妇av软件| 亚洲精华国产精华精| 国产高清国产精品国产三级| 伊人久久大香线蕉亚洲五| 国产免费现黄频在线看| 午夜成年电影在线免费观看| 99久久久亚洲精品蜜臀av| 成人免费观看视频高清| 国产精品久久电影中文字幕| 女人被躁到高潮嗷嗷叫费观| 黄频高清免费视频| 国产精品 欧美亚洲| 亚洲 欧美 日韩 在线 免费| 午夜福利在线免费观看网站| 国产精品永久免费网站| 纯流量卡能插随身wifi吗| 一个人免费在线观看的高清视频| 法律面前人人平等表现在哪些方面| 亚洲专区国产一区二区| 国产三级黄色录像| 国产成人一区二区三区免费视频网站| 激情在线观看视频在线高清| 亚洲熟妇熟女久久| 成人三级做爰电影| 女人被躁到高潮嗷嗷叫费观| 免费女性裸体啪啪无遮挡网站| 亚洲精品国产一区二区精华液| 在线观看66精品国产| www.999成人在线观看| 视频区图区小说| 国产成+人综合+亚洲专区| 新久久久久国产一级毛片| 宅男免费午夜| 男女床上黄色一级片免费看| 欧美性长视频在线观看| 色哟哟哟哟哟哟| 亚洲国产看品久久| 亚洲成人国产一区在线观看| 日韩视频一区二区在线观看| 日韩大码丰满熟妇| 桃色一区二区三区在线观看| 久久热在线av| 亚洲午夜理论影院| 妹子高潮喷水视频| 日韩欧美一区视频在线观看| 亚洲性夜色夜夜综合| 亚洲一卡2卡3卡4卡5卡精品中文| √禁漫天堂资源中文www| 国产亚洲精品一区二区www| ponron亚洲| 亚洲,欧美精品.| 久久亚洲真实| 成人手机av| 国产免费现黄频在线看| 香蕉丝袜av| 久久性视频一级片| 制服人妻中文乱码| 欧美日韩瑟瑟在线播放| 亚洲精品粉嫩美女一区| 国产成人欧美在线观看| 欧美日韩亚洲高清精品| 亚洲av第一区精品v没综合| 欧美大码av| 国产欧美日韩一区二区三区在线| 国产伦人伦偷精品视频| 91老司机精品| 人人妻,人人澡人人爽秒播| 大型av网站在线播放| 两个人看的免费小视频| 神马国产精品三级电影在线观看 | 精品卡一卡二卡四卡免费| 久久狼人影院| 黄色片一级片一级黄色片| 岛国在线观看网站| 国产午夜精品久久久久久| 国产97色在线日韩免费| 欧美日韩黄片免| 叶爱在线成人免费视频播放| 国产高清videossex| 色哟哟哟哟哟哟| 国产有黄有色有爽视频| 国产精品一区二区在线不卡| 国产精品影院久久| 久久中文字幕人妻熟女| 亚洲五月天丁香| 黄色女人牲交| 女性生殖器流出的白浆| 午夜免费激情av| 国产在线精品亚洲第一网站| 中文亚洲av片在线观看爽| 人人妻人人爽人人添夜夜欢视频| 国产色视频综合| 国产亚洲欧美精品永久| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩黄片免| 深夜精品福利| 国产精华一区二区三区| 精品第一国产精品| 天堂俺去俺来也www色官网| 亚洲精品国产精品久久久不卡| 别揉我奶头~嗯~啊~动态视频| 久久久久久久久中文| 亚洲精品在线观看二区| 亚洲一码二码三码区别大吗| 操美女的视频在线观看| 9191精品国产免费久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精华国产精华精| 91大片在线观看| 久久久久久人人人人人| 国产伦人伦偷精品视频| 丰满迷人的少妇在线观看| 亚洲伊人色综图| 电影成人av| 一级毛片女人18水好多| 黑人巨大精品欧美一区二区mp4| 国产主播在线观看一区二区| 1024视频免费在线观看| 黑人欧美特级aaaaaa片| 99久久久亚洲精品蜜臀av| 久久久国产精品麻豆| 久久久国产成人精品二区 | 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 免费观看精品视频网站| 国产高清视频在线播放一区| 亚洲av日韩精品久久久久久密| 国产99白浆流出| 啦啦啦免费观看视频1| 亚洲av第一区精品v没综合| 久久久久精品国产欧美久久久| 侵犯人妻中文字幕一二三四区| 丰满的人妻完整版| a在线观看视频网站| 黑人猛操日本美女一级片| 老熟妇乱子伦视频在线观看| 亚洲国产精品sss在线观看 | 国产亚洲欧美98| 日日摸夜夜添夜夜添小说| 亚洲精品在线观看二区| 亚洲精品中文字幕一二三四区| 免费久久久久久久精品成人欧美视频| 亚洲av成人一区二区三| 亚洲精品国产区一区二| videosex国产| 最新在线观看一区二区三区| 久久精品亚洲熟妇少妇任你| 久久人人97超碰香蕉20202| 亚洲中文av在线| 男人操女人黄网站| 91在线观看av| 国产不卡一卡二| 电影成人av| www.精华液| 久久婷婷成人综合色麻豆| 亚洲人成77777在线视频| 一进一出抽搐动态| 国产成人精品久久二区二区免费| 99久久人妻综合| 国产精品久久视频播放| 国产在线精品亚洲第一网站| 999久久久国产精品视频| 久久久久久大精品| 欧美日韩福利视频一区二区| 一本综合久久免费| 久热这里只有精品99| 99riav亚洲国产免费| 黑人欧美特级aaaaaa片| 老熟妇仑乱视频hdxx| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久久毛片| 久久精品影院6| 欧美丝袜亚洲另类 | 岛国视频午夜一区免费看| 一级作爱视频免费观看| 黑人欧美特级aaaaaa片| 99久久久亚洲精品蜜臀av| 黑丝袜美女国产一区| √禁漫天堂资源中文www| 国产精品综合久久久久久久免费 | 午夜福利一区二区在线看| 黄片大片在线免费观看| 成人18禁在线播放| 老熟妇仑乱视频hdxx| 91av网站免费观看| 日本撒尿小便嘘嘘汇集6| 黑人操中国人逼视频| 国产精品久久久人人做人人爽| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看 | 国产精品爽爽va在线观看网站 | 动漫黄色视频在线观看| 免费观看精品视频网站| 国产精品1区2区在线观看.| 丰满的人妻完整版| 男人操女人黄网站| netflix在线观看网站| 国产主播在线观看一区二区| 国产色视频综合| 热re99久久国产66热| 国产高清激情床上av| 一边摸一边抽搐一进一小说| 国产成人免费无遮挡视频| 亚洲国产毛片av蜜桃av| 日本免费一区二区三区高清不卡 | 亚洲,欧美精品.| 中文字幕人妻熟女乱码| 久久精品亚洲精品国产色婷小说| 国产xxxxx性猛交| 欧美精品一区二区免费开放| svipshipincom国产片| 亚洲精品在线美女| 免费观看人在逋| 国产高清国产精品国产三级| 国产一区二区三区视频了| 亚洲av第一区精品v没综合| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 色综合站精品国产| 99久久综合精品五月天人人| 丰满饥渴人妻一区二区三| 黑丝袜美女国产一区| 日本三级黄在线观看| 999久久久精品免费观看国产| 久久久国产精品麻豆| 性欧美人与动物交配| 侵犯人妻中文字幕一二三四区| 亚洲熟妇中文字幕五十中出 | 动漫黄色视频在线观看| videosex国产| 精品久久蜜臀av无| 老司机亚洲免费影院| 免费在线观看完整版高清| www国产在线视频色| 老熟妇仑乱视频hdxx| 亚洲aⅴ乱码一区二区在线播放 | 韩国av一区二区三区四区| 一区二区日韩欧美中文字幕| 大香蕉久久成人网| 窝窝影院91人妻| 欧美在线黄色| 久久精品成人免费网站| 淫妇啪啪啪对白视频| 亚洲av成人不卡在线观看播放网| 国产亚洲精品久久久久5区| 亚洲男人的天堂狠狠| 亚洲国产欧美日韩在线播放| 中文字幕人妻丝袜制服| 中文字幕最新亚洲高清| 很黄的视频免费| 免费在线观看视频国产中文字幕亚洲| 色尼玛亚洲综合影院| 国产精品国产高清国产av| 88av欧美| 久久久国产一区二区| 久久久久亚洲av毛片大全| 亚洲黑人精品在线| 亚洲在线自拍视频| 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 校园春色视频在线观看| 国产精品影院久久| av网站免费在线观看视频| 久久久久久久久中文| 91麻豆精品激情在线观看国产 | 亚洲午夜理论影院| 国产蜜桃级精品一区二区三区| 十分钟在线观看高清视频www| 欧美激情久久久久久爽电影 | 亚洲精品中文字幕一二三四区| 国产成人av激情在线播放| 丁香欧美五月| 一区在线观看完整版| 一本大道久久a久久精品| 日韩欧美三级三区| 自拍欧美九色日韩亚洲蝌蚪91| ponron亚洲| 亚洲一区二区三区色噜噜 | 一边摸一边抽搐一进一小说| 亚洲午夜理论影院| www.999成人在线观看| 欧美另类亚洲清纯唯美| 黑人巨大精品欧美一区二区蜜桃| 久久精品亚洲av国产电影网| 欧美人与性动交α欧美软件| 久久精品影院6| 久久草成人影院| 亚洲成人免费电影在线观看| 国产精品 国内视频| 免费在线观看日本一区| 亚洲国产精品999在线| 午夜福利在线观看吧| www日本在线高清视频| 国产91精品成人一区二区三区| 视频区图区小说| 亚洲三区欧美一区| 午夜亚洲福利在线播放| 日韩大尺度精品在线看网址 | 亚洲欧美日韩无卡精品| 中文欧美无线码| 亚洲午夜理论影院| 亚洲视频免费观看视频| 欧美黄色淫秽网站| 黄色 视频免费看| 一级片免费观看大全| 一区二区三区国产精品乱码| 国产野战对白在线观看| 如日韩欧美国产精品一区二区三区| 老司机亚洲免费影院| 99在线人妻在线中文字幕| 成人亚洲精品av一区二区 | 男人操女人黄网站| 亚洲专区字幕在线| 亚洲久久久国产精品| 女人被狂操c到高潮| 亚洲 欧美 日韩 在线 免费| 天堂中文最新版在线下载| 午夜免费鲁丝| 91成人精品电影| 在线观看午夜福利视频| 欧美性长视频在线观看| 久久久国产精品麻豆| av天堂久久9| 久久狼人影院| 日本wwww免费看| 欧美成人午夜精品| 成人免费观看视频高清| 亚洲情色 制服丝袜| 日韩欧美一区二区三区在线观看| 国产麻豆69| 纯流量卡能插随身wifi吗| 中亚洲国语对白在线视频| 在线国产一区二区在线| 久热这里只有精品99| 99久久综合精品五月天人人| 韩国精品一区二区三区| 大香蕉久久成人网| 老司机深夜福利视频在线观看| 黄片小视频在线播放| 午夜久久久在线观看| 精品电影一区二区在线| 亚洲国产欧美一区二区综合| 一级毛片女人18水好多| 国产一区在线观看成人免费| 国产视频一区二区在线看| 免费看a级黄色片| 如日韩欧美国产精品一区二区三区| 美女国产高潮福利片在线看| 俄罗斯特黄特色一大片| 一级片免费观看大全| 色播在线永久视频| 午夜福利,免费看| 在线播放国产精品三级| 免费少妇av软件| 亚洲人成电影免费在线| 成人特级黄色片久久久久久久| 一夜夜www| 91在线观看av| 日韩三级视频一区二区三区| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看 | 人人妻人人添人人爽欧美一区卜| 国产一卡二卡三卡精品| 级片在线观看| 免费av毛片视频| 人妻久久中文字幕网| 久久久国产欧美日韩av| 极品人妻少妇av视频| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 一级黄色大片毛片| 18禁美女被吸乳视频| 一本大道久久a久久精品| а√天堂www在线а√下载| 免费在线观看视频国产中文字幕亚洲| 大陆偷拍与自拍| 91精品国产国语对白视频| 天堂动漫精品| 国产人伦9x9x在线观看| 欧美黑人精品巨大| 女人被狂操c到高潮| 国产欧美日韩综合在线一区二区| 国产视频一区二区在线看| 夜夜爽天天搞| 侵犯人妻中文字幕一二三四区| 亚洲av熟女| 亚洲自偷自拍图片 自拍| 国产av精品麻豆| 窝窝影院91人妻| 亚洲久久久国产精品| 黑人欧美特级aaaaaa片| 在线免费观看的www视频| 国产精品 欧美亚洲| 久久精品国产亚洲av香蕉五月| 日本三级黄在线观看| 他把我摸到了高潮在线观看| 国产高清国产精品国产三级| 午夜a级毛片| 久久久久国产一级毛片高清牌| 一级a爱片免费观看的视频| 日韩一卡2卡3卡4卡2021年| 三上悠亚av全集在线观看| 国产极品粉嫩免费观看在线| 日韩免费av在线播放| 69精品国产乱码久久久| 一进一出抽搐gif免费好疼 | 亚洲中文字幕日韩| 啪啪无遮挡十八禁网站| 成年版毛片免费区| 老鸭窝网址在线观看| 亚洲精品一卡2卡三卡4卡5卡| 精品福利永久在线观看| 国产精品免费视频内射| 大香蕉久久成人网| 激情视频va一区二区三区| 午夜免费激情av| 久久久久久免费高清国产稀缺| 琪琪午夜伦伦电影理论片6080| 男男h啪啪无遮挡| 久久精品人人爽人人爽视色| 国产99久久九九免费精品| 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久| 国产不卡一卡二| 午夜成年电影在线免费观看| 1024香蕉在线观看| 免费观看精品视频网站| 国产黄a三级三级三级人| 在线播放国产精品三级| 国产av一区二区精品久久| 欧美大码av| 99精品欧美一区二区三区四区| 男男h啪啪无遮挡| 精品免费久久久久久久清纯| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区精品| 国产男靠女视频免费网站| 欧美一级毛片孕妇| 午夜免费观看网址| 中文字幕av电影在线播放| 精品国产乱码久久久久久男人| 777久久人妻少妇嫩草av网站| 国产精品成人在线| 欧美日韩黄片免| 免费在线观看黄色视频的| 日韩欧美国产一区二区入口| 国产97色在线日韩免费| 在线视频色国产色| 精品日产1卡2卡| 国产欧美日韩一区二区三| 国产人伦9x9x在线观看| 日本精品一区二区三区蜜桃| 窝窝影院91人妻| 国产一区二区激情短视频| 亚洲欧美一区二区三区久久| 久久99一区二区三区| 免费在线观看视频国产中文字幕亚洲| 国产有黄有色有爽视频| 国产精品久久久久成人av| 久久天躁狠狠躁夜夜2o2o| 又紧又爽又黄一区二区| 大香蕉久久成人网| 午夜视频精品福利| 97超级碰碰碰精品色视频在线观看| 亚洲第一青青草原| 好看av亚洲va欧美ⅴa在| 黄片大片在线免费观看| 丁香欧美五月| 国产成人免费无遮挡视频| 免费看十八禁软件| 亚洲欧美日韩无卡精品| 国产精品电影一区二区三区| 啪啪无遮挡十八禁网站| 国产av在哪里看| 国产精品一区二区免费欧美| 大型av网站在线播放| 日韩精品青青久久久久久| 国产日韩一区二区三区精品不卡| 亚洲精品中文字幕一二三四区| 日本免费一区二区三区高清不卡 | 一级a爱视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 美女高潮喷水抽搐中文字幕| 欧美国产精品va在线观看不卡| 国产欧美日韩综合在线一区二区| 一进一出好大好爽视频| 亚洲自偷自拍图片 自拍| 久99久视频精品免费| 午夜免费观看网址| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩福利视频一区二区| 国产av一区二区精品久久| 国产伦一二天堂av在线观看| 在线观看日韩欧美| 久久精品国产清高在天天线| 国产伦一二天堂av在线观看| 一级作爱视频免费观看| 欧美午夜高清在线| 久久精品亚洲熟妇少妇任你| 国产成人av教育| 国产精品免费一区二区三区在线| 免费日韩欧美在线观看| 国产精品秋霞免费鲁丝片| 99热只有精品国产| 欧美日韩av久久| 久久久久久亚洲精品国产蜜桃av| 成年人免费黄色播放视频| 国产亚洲av高清不卡| 色综合欧美亚洲国产小说|