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

    基于WRF的城市熱島效應(yīng)高分辨率評(píng)估方法

    2017-03-29 21:50:23郭飛
    土木建筑與環(huán)境工程 2017年1期
    關(guān)鍵詞:城市環(huán)境

    摘要:為解決城市總規(guī)尺度氣候評(píng)估工具缺乏的問題,分析了不同氣候評(píng)估方法的特點(diǎn)和應(yīng)用范圍,研究了氣象學(xué)工具WRF在城市熱島效應(yīng)評(píng)估中的應(yīng)用方法及有效性。建立了針對(duì)城市總規(guī)尺度的高精度城市規(guī)劃信息數(shù)據(jù)庫(kù)、氣象數(shù)據(jù)庫(kù)生成WRF模擬邊界的方法;明確了高精度城市熱島效應(yīng)模擬所需的WRF物理方案。以沿海城市大連夏季典型晴朗高溫天氣為例,利用WRF模型進(jìn)行了300 m的高分辨率模擬。結(jié)果表明,WRF可以較好地模擬出城市熱島效應(yīng)的強(qiáng)度和時(shí)空變化;將模擬結(jié)果與氣象站和現(xiàn)場(chǎng)觀測(cè)值進(jìn)行了對(duì)比,兩者變化趨勢(shì)一致,偏差在可接受的范圍之內(nèi)。

    關(guān)鍵詞:城市環(huán)境;高分辨率模擬;氣候模型;熱島效應(yīng)

    中圖分類號(hào):TU119.4文獻(xiàn)標(biāo)志碼:A文章編號(hào):16744764(2017)01001307

    收稿日期:20160421

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(51308087、51278078)

    作者簡(jiǎn)介:郭飛(1980),男,副教授,博士,主要從事建筑技術(shù)科學(xué)、城市氣候研究,(Email)guofei@dlut.edu.cn。

    Received:20160421

    Foundation item:National Natural Science Foundation of China (No. 51308087, 51278078)

    Author brief:Guo Fei (1980),associate professor,PhD,main research interests:building technology and science,urban climate, (Email) guofei@dlut.edu.cn.Assessment method of urban heat island high resolution based on WRF

    Guo Fei

    (Architecture and Fine Art School, Dalian University of Technology, Dalian 116023, Liao Ning, P. R. China)

    Abstract:Different climate assessment tools features and application scope were analyzed to solve the problem of lack of climate assessment tools in the urban master planning scale. The WRF models application method and validity was studied in urban heat island (UHI) evaluation. According to the demand of urban master planning scale, the method to generate WRF simulation boundary conditions with high resolution urban geographic and land cover data was established. The physics schemes for WRF high resolution simulation were also clarified. A typical summer hot weather day was taken as the research background in the coastal city of Dalian, the WRF model was used to simulate the UHI effect with the resolution of 300 meters. The results showed that WRF model could reproduce the intensity, temporal and spatial variation of UHI effect. The simulation results were compared with the data of weather station and the in situ observation, their change trend was found to be consistent and the deviation was within the acceptable range.

    Keywords:urban environment, high resolution simulation, climate models, urban heat island

    根據(jù)《國(guó)家新型城鎮(zhèn)化規(guī)劃(2014—2020)》,中國(guó)城市人口比例在2020年將達(dá)到60%[1]。城市人口和規(guī)模的擴(kuò)張會(huì)帶來不透水面積增加、人為熱大量排放等問題,出現(xiàn)城市中心區(qū)的溫度高于郊區(qū)的城市熱島現(xiàn)象[2]。在全球氣候變暖的背景下,城市熱島效應(yīng)會(huì)增加環(huán)境的熱壓力,導(dǎo)致高溫天氣的死亡率增加[34],也會(huì)引發(fā)夏季空調(diào)大量使用和電力需求上升,導(dǎo)致更多的污染物產(chǎn)生和各種其他健康威脅[56]。緩解城市熱島效應(yīng)需要綜合采取優(yōu)化城市規(guī)劃布局、促進(jìn)通風(fēng)、建設(shè)綠地和屋頂綠化、優(yōu)化建筑形態(tài)、加強(qiáng)建筑節(jié)能等各類措施。但是由于城市熱島效應(yīng)成因非常復(fù)雜,與風(fēng)、熱、濕、輻射等多種物理因素有關(guān);既受大氣候影響,也受局地地形、下墊面類型、城市和建筑形態(tài)的影響。目前,在城市宏觀層面的科學(xué)評(píng)估工具十分缺乏,難以分析規(guī)劃措施與城市熱島效應(yīng)之間的定量關(guān)系。有鑒于此,許多研究者采用各種研究工具與城市規(guī)劃相對(duì)接,研究了城市熱島的評(píng)估方法。代表性的有現(xiàn)場(chǎng)實(shí)測(cè)、衛(wèi)星遙感影像反演、計(jì)算流體力學(xué)(CFD,Computational Fluid Dynamics)模型和氣象學(xué)模型等。

    例如,王志浩等[7]通過實(shí)測(cè)發(fā)現(xiàn)了重慶城市熱島存在的帶狀特征;齊靜靜等[8]研究了河流對(duì)城市熱島的緩解作用;陳云[9]利用Landsat ETM衛(wèi)星數(shù)據(jù)反演出不同城市熱島效應(yīng)的時(shí)空分布特征和發(fā)展過程,并提出優(yōu)化空間結(jié)構(gòu)、建設(shè)綠化系統(tǒng)、減少人為排放等規(guī)劃建議。現(xiàn)場(chǎng)實(shí)測(cè)對(duì)人力物力要求較高,能夠獲取數(shù)據(jù)的時(shí)間范圍和空間精度受到諸多限制;衛(wèi)星遙感數(shù)據(jù)則限制較少。但這兩種方法都難以對(duì)未來城市熱島的變化趨勢(shì)加以預(yù)測(cè),也無法在方案階段對(duì)不同城市規(guī)劃布局和建筑形態(tài)加以評(píng)估。

    CFD模型的特點(diǎn)是分辨率較高(1~10 m),可以非常精細(xì)地研究城市建筑三維形態(tài)特征對(duì)氣候的影響[10],許多研究者利用Fluent、CFDUrban、StarCD、Phoenics和Envimet等工具進(jìn)行了研究[11]。但是CFD軟件需要針對(duì)城市微氣候復(fù)雜的氣候機(jī)理開發(fā)專門的輻射、蒸發(fā)、人工排熱等模型,限制了其在總規(guī)尺度的使用。

    中尺度氣候模型原本是針對(duì)10~1 000 km尺度天氣預(yù)報(bào)和氣象研究需求產(chǎn)生的工具,主要的有WRF、UKMO、ARPS等[12]。其中WRF氣候模型(Weather Research and Forecast Model)是美國(guó)國(guó)家大氣研究中心等機(jī)構(gòu)開發(fā)的最新一代天氣預(yù)報(bào)和氣象模型,在全世界得到了廣泛應(yīng)用[1316]。城市總規(guī)尺度的氣候現(xiàn)象異常復(fù)雜,WRF模型具有完善的物理方案加以描述。WRF還有優(yōu)秀的跨尺度研究能力,其自身內(nèi)嵌描述城市對(duì)氣候影響的單層城市冠層模型(UCM)、多層城市冠層模型和建筑能耗模型(BEP+BEM)等,可以設(shè)置二十余項(xiàng)參數(shù)描述城市和建筑形態(tài)特征。近年來利用WRF進(jìn)行城市局地氣候特征分析,并據(jù)此提出城市建設(shè)、規(guī)劃和管理的政策建議及改善策略,逐漸成為研究的熱點(diǎn)之一。

    WRF的模擬分辨率可從100 m跨越到100 km。許多學(xué)者研究了將高分辨率的城市信息庫(kù)數(shù)據(jù)應(yīng)用于WRF模擬。例如LIN等利用MODIS(中等分辨率成像光譜儀[17])的1 km分辨率的土地利用數(shù)據(jù)對(duì)臺(tái)北地區(qū)的城市熱島效應(yīng)進(jìn)行了分析[18]。ADACHI等[19]利用NUDAPT(美國(guó)國(guó)家城市數(shù)據(jù)庫(kù)[20])以1 km的分辨率對(duì)芝加哥的模擬結(jié)果研究顯示,WRF耦合的BEP+BEM模型能準(zhǔn)確地再現(xiàn)城市熱島現(xiàn)象,還可以較為準(zhǔn)確地預(yù)測(cè)空調(diào)、采暖的能耗與氣候之間的關(guān)系。PIU[21]利用WRF模擬了香港和珠江三角洲城市的可見度和空氣質(zhì)量,所采用的1 km分辨率城市土地利用信息則來自谷歌衛(wèi)星影像。

    SALAMANCA[22]利用NUDAPT對(duì)美國(guó)休斯頓和德克薩斯城進(jìn)行了研究,發(fā)現(xiàn)詳盡的城市形態(tài)數(shù)據(jù)和信息對(duì)成功模擬城市熱島效應(yīng)十分重要。盡管1 km的分辨率對(duì)于氣象模擬來說已經(jīng)屬于較高的精度,但是D'ALESSANDRO等[23]的研究發(fā)現(xiàn)更高的分辨率(例如250 m)會(huì)提高模擬的準(zhǔn)確性。而且對(duì)于城市控規(guī)、詳規(guī)和城市設(shè)計(jì)而言,較高分辨率的氣候信息便于研究者對(duì)城市形態(tài)與氣候的關(guān)聯(lián)性加以深入分析,其結(jié)果更具有實(shí)踐參考價(jià)值。但是要成功實(shí)現(xiàn)WRF高分辨率模擬需要解決兩個(gè)問題,一是設(shè)置準(zhǔn)確、高精度的邊界條件,包括氣象數(shù)據(jù)、城市規(guī)劃和形態(tài)參數(shù)信息等,以體現(xiàn)城市對(duì)氣候的影響;二是采用合適的物理方案,以準(zhǔn)確捕捉高分辨率條件下城市熱島的特征。

    利用WRF以300 m分辨率模擬了大連夏季典型氣象日的溫度場(chǎng)和風(fēng)場(chǎng)特征,對(duì)WRF高分辨率模擬的邊界條件和物理方案進(jìn)行了分析。并將模擬結(jié)果與氣象站和地面現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)進(jìn)行了對(duì)比。研究對(duì)了解城市熱島效應(yīng)特征并提出相應(yīng)規(guī)劃策略提供了一種思路。1高分辨率評(píng)估方法

    大連東臨黃海、西鄰渤海,屬于海洋性溫暖帶季風(fēng)氣候。冬季主導(dǎo)風(fēng)向?yàn)楸憋L(fēng),平均風(fēng)速5.0 m/s,夏季主導(dǎo)風(fēng)向?yàn)槟掀髂巷L(fēng),平均風(fēng)速為4.0 m/s(圖1),全年溫度在-16.8~37 ℃之間,年平均氣溫為11.19 ℃。大連城市山地多、地形較為復(fù)雜,區(qū)域內(nèi)有長(zhǎng)白山脈的延伸段,主城區(qū)西部山地最高海拔約360 m,東部山地最高海拔約240 m(圖2)。

    圖1大連地區(qū)風(fēng)玫瑰圖

    Fig.1Wind rose diagram of Dalian city圖2大連地區(qū)地形高程圖(90 m分辨率

    DEM數(shù)據(jù),數(shù)據(jù)來源SRTM)

    Fig.2Topographic maps of Dalian city

    (SRTM DEM data of 90m resolution)1.1邊界條件

    研究采用的模擬工具是WRFV3.6版本。驅(qū)動(dòng)WRF模型的數(shù)據(jù)主要包括氣象數(shù)據(jù)、地形數(shù)據(jù)、土地用地類型數(shù)據(jù)3種。WRF的氣象數(shù)據(jù)主要有兩種類型。一是美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NECP)提供的覆蓋全球的再分析氣象資料(FNL)。這個(gè)數(shù)據(jù)來源于全球數(shù)據(jù)同化系統(tǒng)(GDAS),主要由全球數(shù)據(jù)通信網(wǎng)持續(xù)地采集并進(jìn)行實(shí)時(shí)更新,并經(jīng)過數(shù)據(jù)分析和校正后形成,數(shù)據(jù)質(zhì)量最高。二是NECP提供的全球氣象預(yù)報(bào)數(shù)據(jù)(GFS),時(shí)空分辨率比FNL再分析數(shù)據(jù)高,并且能夠提早1 h形成,可用于驅(qū)動(dòng)WRF進(jìn)行中尺度天氣預(yù)報(bào)。由于本研究是針對(duì)歷史已發(fā)生的真實(shí)天氣進(jìn)行模擬,所以采用FNL再分析數(shù)據(jù)(時(shí)間分辨率是6 h、水平空間分辨率是1°×1°約110 km)作為模擬條件。

    地形和土地用地類型數(shù)據(jù)是影響城市熱島特征重要因素,因此,獲取高精度的數(shù)據(jù)十分關(guān)鍵。WRF自帶的數(shù)據(jù)有美國(guó)國(guó)家地質(zhì)調(diào)查局(USGS1993年全球土地分類)和中等分辨率成像光譜儀(MODIS)兩種,盡管質(zhì)量較高,但時(shí)間較老,已難反映最新城市狀況,需要加以更新。目前可獲取的高精度地形數(shù)據(jù)有SRTM 90 m精度(圖2)和ASTERA GDEM 30 m精度的數(shù)據(jù)兩種;土地利用數(shù)據(jù)有中國(guó)國(guó)家基礎(chǔ)地理信息中心開發(fā)的全球土地覆蓋(GLC)30 m精度數(shù)據(jù)和清華大學(xué)開發(fā)的90 m精度土地利用觀測(cè)計(jì)劃(FROM)的數(shù)據(jù)。這些數(shù)據(jù)都是GIS格式,方便了用戶使用、比對(duì)和進(jìn)行數(shù)據(jù)格式的轉(zhuǎn)換。但是WRF能夠使用的數(shù)據(jù)都是專用的二進(jìn)制代碼格式,還需要調(diào)用專門的C程序?qū)IS數(shù)據(jù)加以修改和轉(zhuǎn)換,其流程如圖3所示。修正前后的土地分類圖如圖4所示,可以看出大連城市和海岸線已發(fā)生了明顯的擴(kuò)張和改變。

    圖3地形和用地類型邊界條件生成方法圖示

    Fig.3Conversion procedure of GIS data

    to geogrid static data圖4第3層模擬網(wǎng)格的土地用地分類圖

    Fig.4The land use map of inner grid1.2物理方案

    WRF包括了十幾個(gè)物理方案,較為全面的描述了各種天氣現(xiàn)象,本研究采用的實(shí)驗(yàn)方案如表1。其中對(duì)于高分辨率模擬和城市熱島效應(yīng)最重要的是陸面表面的物理過程。一般低分辨率的研究采用陸面表面物理方案(LSM)就能夠描述城市對(duì)氣候的影響。此時(shí)WRF模型將城市簡(jiǎn)化為一個(gè)均質(zhì)的平板,設(shè)置為均勻的0.8 m粗糙度、表面反照率0.15;假定城市下墊面的材質(zhì)主要是混凝土和瀝青;減少城市部分的綠地面積比例和蒸發(fā)量[14]。但對(duì)于高分辨率的城市熱島模擬而言,需要采用專門的城市物理方案,主要有UCM、BEP和BEM模型3種類型。表1WRF的物理方案設(shè)置

    Table 1WRF physics schemes物理過程物理方案微物理Thompson graupel scheme

    一種包含冰、雪、霰過程的適用于高分辨率模擬的微物理方案長(zhǎng)波輻射Rapid Radiative Transfer Model

    一種高精度的快速輻射傳輸模式短波輻射Goddard Shortwave scheme

    包含氣候態(tài)臭氧和云效應(yīng)的雙速多波段短波方案輻射方案

    計(jì)算間隔10 min近地面

    層物理MoninObukhov (Janjic Eta) Similarity scheme

    基于含有Zilitinkevich熱粗糙長(zhǎng)度和來自查表的標(biāo)準(zhǔn)相似函數(shù)的方案續(xù)表1物理過程物理方案陸面表面Noah LandSurface Model

    含有4個(gè)不同層次土壤溫濕度、積雪覆蓋和凍土物理過程行星邊界層MellorYamadaJanjic (Eta) Turbulence Kinetic Energy (TKE)

    含有局地垂直混合的一維診斷湍流動(dòng)能方案城市物理UCM,單層城市冠層模型

    采用了其中較為常用的UCM城市物理方案,它可以描述城市和建筑的三維特征和人為產(chǎn)熱。UCM假定城市是一個(gè)均勻的、無限長(zhǎng)的街道峽谷,峽谷中陰影、反射、輻射的減少效應(yīng)均被考慮,風(fēng)廓線按照指數(shù)律描述。UCM允許研究者將城市細(xì)分為低密度居住區(qū)、高密度居住區(qū)、工業(yè)和商業(yè)區(qū)3種城市用地類型。采用如建筑平均高度、屋頂和道路平均寬度、人工產(chǎn)熱量、建筑墻面、屋頂和道路的熱容、導(dǎo)熱系數(shù)等參數(shù)分別描述建筑城市對(duì)氣候的影響,本文采用的參數(shù)如表2。

    表2UCM城市物理方案的模擬參數(shù)

    s-1·K-1)低密度住區(qū)58.3 8.314.692.30×1061.25×1061.93×1060.157 50.210.699高密度住區(qū)2512.311.7814.692.30×1061.25×1061.93×1060.157 50.210.699工業(yè)和商業(yè)區(qū)30201414.692.30×1061.25×1061.93×1060.157 50.210.699

    1.3網(wǎng)格設(shè)置

    模擬以大連中心城區(qū)為主要對(duì)象,模擬中心點(diǎn)坐標(biāo)N38.923°,E121.647°。由于采用了300 m高分辨率的網(wǎng)格,需要耗費(fèi)大量的計(jì)算資源。一般來說水平分辨率(網(wǎng)格間距)每提高1倍,模擬的總格點(diǎn)數(shù)將為原來的4倍。此時(shí)積分時(shí)間步長(zhǎng)也需相應(yīng)減小為原來的0.5倍,導(dǎo)致模擬計(jì)算的總積分次數(shù)為原來的2倍。因此一般每提高1倍分辨率,完成模擬所需的CPU總時(shí)間將至少為原來的8倍以上,對(duì)計(jì)算能力的耗費(fèi)非常巨大。針對(duì)這一問題一般可采用嵌套網(wǎng)格的形式,能充分考慮周邊大背景氣候的影響,也保證核心網(wǎng)格的分辨率較高,在一定程度上節(jié)約計(jì)算資源。

    研究采用3層嵌套的方式,第1層格點(diǎn)數(shù)60×60,網(wǎng)格間距2.7 km×2.7 km;第2層格點(diǎn)數(shù)81×81,網(wǎng)格間距0.9 km×0.9 km;第3層格點(diǎn)數(shù)105×105,網(wǎng)格間距0.3 km×0.3 km。計(jì)算區(qū)域基本覆蓋了大部分的城市建成區(qū),其中第3層網(wǎng)格覆蓋的是中山區(qū)、西崗區(qū)等城市中心區(qū)(圖4)。

    參照典型氣象日的標(biāo)準(zhǔn),選取大連高溫、晴朗、南風(fēng)的氣象日,即北京時(shí)間20130706T08:00—20130707T15:00,共計(jì)32 h作為模擬的對(duì)象。

    2模擬結(jié)果

    2.1溫度場(chǎng)

    地面2 m的溫度場(chǎng)結(jié)果顯示,7月6日8時(shí)—7月6日12時(shí),氣溫逐步上升,最高氣溫達(dá)到30.3 ℃;13時(shí)起溫度開始下降,夜間最低是7月7日5時(shí),溫度為22.5 ℃,7月7日6時(shí)起溫度開始上升,至8時(shí)升至25.2 ℃。

    從模擬期溫度場(chǎng)的變化可以看出,海風(fēng)對(duì)陸地及城市的降溫效果十分明顯。南部海面附近空氣溫度低5~10 ℃,南風(fēng)將冷空氣吹向城市,在海洋沿岸形成平行于岸線的等溫線,越往陸地深處氣溫越高。由于海風(fēng)的降溫作用,海岸線附近的氣溫比城市中心區(qū)低5 ℃左右。受南風(fēng)影響,城市熱島中心向下風(fēng)向偏移,位于主城區(qū)北部,大連灣附近海面氣溫受上風(fēng)向城市的影響比周邊海面高(圖5)。

    圖57月6日20時(shí)最內(nèi)層溫度場(chǎng)(AA為剖面位置)

    Fig.5Temperature filed of the inner grid on 20:00

    6th July(AA is the cross section line)地面2 m溫度場(chǎng)剖面變化過程顯示,大連中心區(qū)存在較強(qiáng)的城市熱島效應(yīng)(圖6)。其中剖面上顯示出有2個(gè)熱島中心,分別位于中山區(qū)和甘井子區(qū)。7月6日白天熱島強(qiáng)度較強(qiáng),熱島強(qiáng)度最大為5.6 ℃,出現(xiàn)在7月6日13時(shí)。7月6日15時(shí)之后熱島強(qiáng)度有明顯的回落,夜間城市熱島最小為0.5°,出現(xiàn)在7月7日6時(shí)。夜間山體降溫速度顯著快于城市建成區(qū),例如其中7月6日20時(shí)的溫度場(chǎng)分布顯示,山體氣溫低于城市建成區(qū)約3 ℃,對(duì)城市總體提供了良好的冷源。

    圖6地面2 m溫度剖面隨時(shí)間的變化

    Fig.6Temperature cross section variation with time2.2風(fēng)場(chǎng)

    模擬期內(nèi)風(fēng)向?yàn)槟巷L(fēng),7月6日白天城市中心區(qū)風(fēng)速較低,約為0.5~1.5 m/s之間。風(fēng)速最小值0.5 m/s,出現(xiàn)在7月6日8時(shí)(圖7)。由于7月7日白天風(fēng)速較低,相應(yīng)地城市中心熱島強(qiáng)度較高。7月6日15時(shí)之后風(fēng)速繼續(xù)下降,7月7日凌晨2時(shí)之后風(fēng)速開始加強(qiáng),最大值3.0 m/s,出現(xiàn)在7月7日8時(shí)。

    圖77月6日8時(shí)最內(nèi)層風(fēng)場(chǎng)

    Fig.7The inner wind field at 8:00 a.m. on 6th July2.3模擬結(jié)果驗(yàn)證

    為了驗(yàn)證WRF模擬結(jié)果的準(zhǔn)確性,我們將WRF數(shù)據(jù)與氣象站和現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)加以對(duì)比分析。我們于2013年7月6日—7月7日進(jìn)行了現(xiàn)場(chǎng)觀測(cè),從上午8:00到下午18:00對(duì)大連中山區(qū)進(jìn)行人工現(xiàn)場(chǎng)連續(xù)觀測(cè),每小時(shí)記錄一組氣象參數(shù),18:00之后的夜間觀測(cè)由放置現(xiàn)場(chǎng)的儀器完成(圖8)。人工現(xiàn)場(chǎng)觀測(cè)采用NK4500的手持氣象站。風(fēng)速測(cè)量范圍是0.4~40 m/s,測(cè)量精度是3%;空氣溫度測(cè)量范圍是-29~70 ℃,精度是1 ℃;濕度測(cè)量范圍是0~100%,精度3%。另外在大連理工大學(xué)建筑館屋頂還設(shè)置了一個(gè)HOBO氣象站同步觀測(cè)。風(fēng)速測(cè)量范圍是0~44 m/s,測(cè)量精度是0.5 m/s;空氣溫度測(cè)量范圍是-40~75 ℃,精度是0.2 ℃;濕度測(cè)量范圍是0~100%,精度2.5%。

    圖8現(xiàn)場(chǎng)觀測(cè)及固定氣象站觀測(cè)

    Fig.8In situ observation (left) and weather station (right)將WRF溫度場(chǎng)模擬結(jié)果分別與現(xiàn)場(chǎng)實(shí)測(cè)和氣象站結(jié)果對(duì)比,可以發(fā)現(xiàn)兩者之間的變化趨勢(shì)基本吻合。溫度場(chǎng)觀測(cè)值與模擬值最大偏差2.2 ℃,最小0.1 ℃,平均偏差0.9 ℃(圖9)。風(fēng)速模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)及氣象站觀測(cè)數(shù)據(jù)也基本吻合,最大偏差1.0 m/s,最小0 m/s,平均0.3 m/s(圖10)。對(duì)于尺度較大的城市而言,這些模擬結(jié)果的準(zhǔn)確性是十分理想的,用于分析城市溫度場(chǎng)、風(fēng)場(chǎng)的變化規(guī)律具有較高的可信度。

    圖9溫度場(chǎng)模擬結(jié)果與實(shí)測(cè)對(duì)比

    (7月6日8時(shí)至7月7日8時(shí))

    Fig.9Comparison of temperature field between simulation

    and observations(from 8:00 6th July to 8:00 7th July)圖10風(fēng)場(chǎng)模擬結(jié)果與實(shí)測(cè)對(duì)比

    (7月6日8時(shí)至7月7日8時(shí))

    Fig.10Comparison of wind field between simulation and

    observations(from 8:00 6th July to 8:00 7th July)3結(jié)論

    進(jìn)行高分辨率的氣象模擬可以準(zhǔn)確的研究城市局地氣候的變化特征,為分析城市熱島效應(yīng)、采取相應(yīng)的緩解措施提供了依據(jù)。進(jìn)行WRF高分辨率模擬需要采用高精度城市地形和土地利用數(shù)據(jù)修正模型自帶的數(shù)據(jù)以提供準(zhǔn)確的邊界條件,并選擇UCM等合適的城市物理方案。以大連為例進(jìn)行了300 m高分辨率夏季典型氣象日的模擬,結(jié)果表明模擬期內(nèi)城市存在較強(qiáng)的熱島效應(yīng),最高時(shí)強(qiáng)度達(dá)5.6 ℃,熱島中心位于城市北部即甘井子區(qū)金三角一帶。但是由于陸地和海洋溫差較大,因此海風(fēng)會(huì)對(duì)城市熱島具有較強(qiáng)的緩解作用,風(fēng)速大于2 m/s時(shí)城市熱島效應(yīng)得到了明顯的削弱。因此沿海城市規(guī)劃應(yīng)當(dāng)考慮如何充分利用清潔的海風(fēng)為城市夏季提供致冷源。例如規(guī)劃平行于主導(dǎo)風(fēng)、垂直于岸線的風(fēng)道等措施。

    參考文獻(xiàn):

    [1] 國(guó)務(wù)院.國(guó)家新型城鎮(zhèn)化規(guī)劃(2014—2020年)[EB/OL]. http://www.gov.cn/zhengce/201403/16/content_2640075.htm, 20140316.

    The State Council. National new urbanization plan [EB/OL].http://www.gov.cn/zhengce/201403/16/content_2640075.htm, 20140316. (in Chinese)

    [2] OKE T R. The energetic basis of the urban heat island[J]. Quarterly Journal of the Royal Meteorological Society, 1982, 108(455): 124.

    [3] OLESON K W, MONAGHAN A, WILHELMI O, et al. Interactions between urbanization, heat stress, and climate change[J]. Climatic Change, 2015, 129(3/4): 525541.

    [4] OHASHI Y, IHARA T, KIKEGAWA Y, et al. Numerical simulations of influence of heat island countermeasures on outdoor human heat stress in the 23 wards of Tokyo, Japan[J]. Energy and Buildings, 2015,(5):4044.

    [5] ZHONG S, QIAN Y, ZHAO C, et al. A case study of urbanization impact on summer precipitation in the Greater Beijing Metropolitan Area: Urban heat island versus aerosol effects[J]. Journal of Geophysical Research: Atmospheres, 2015, 120(20):3336.

    [6] FALLMANN J, FORKEL R, EMEIS S. Secondary effects of urban heat island mitigation measures on air quality[J]. Atmospheric Environment, 2016, 125: 199211.

    [7] 王志浩. 山地城鎮(zhèn)熱島特征與測(cè)評(píng)方法研究[D].重慶:重慶大學(xué),2012.

    WANG Z H. Study on characteristics and evaluation methods of heat island in mountainous city[D].Chongqing:Chongqing University,2012. (in Chinese)

    [8] 齊靜靜,劉京,宋曉程,等. 大型城市河流對(duì)城市氣候影響的實(shí)測(cè)研究[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào),2011(10):5659.

    QI J J, LIU J, SONG X C, et al. Field measurement of the influence of large urban river on urban thermal climate [J]. Journal of Harbin Institute of Technology, 2011(10):5659.(in Chinese)

    [9] 陳云. 基于Landsat8的城市熱島效應(yīng)研究初探——以廈門市為例[J]. 測(cè)繪與空間地理信息,2014(2):123128.

    CHEN Y. A preliminary study on the urban heat island effect based on landsat 8:case of xiamen city[J]. Geomatics & Spatial Information Technology,2014(2):123128. (in Chinese)

    [10] 壽亦萱,張大林. 城市熱島效應(yīng)的研究進(jìn)展與展望[J]. 氣象學(xué)報(bào),2012,70(3):338353.

    SHOU Y X, ZHANG D L. Recent advances in understanding urban heat island effects with some future prospects [J]. Acta Meteorologica Sinica, 2012, 70(3): 338353. (in Chinese)

    [11] YANG X, ZHAO L, BRUSE M, et al. An integrated simulation method for building energy performance assessment in urban environments[J]. Energy and Buildings, 2012, 54: 243251.

    [12] 黃菁,張強(qiáng). 中尺度大氣數(shù)值模擬及其進(jìn)展[J]. 干旱區(qū)研究,2012,29(2):273283.

    HUANG J, ZHANG Q. Mesoscale atmospheric numerical simulation and its progress[J]. Arid Zone Research, 2012, 29(2):273283. (in Chinese)

    [13] CHEN F, KUSAKA H, Bornstein R, et al. The integrated WRF/urban modelling system: development, evaluation, and applications to urban environmental problems[J]. International Journal of Climatology, 2011, 31(2): 273288.

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

    GUO F, ZHU P S, DUAN D W, et al. Highdensity urban climate evaluation method and application[J] Journal of Human Settlements in West China, 2015, 30(6):1923. (in Chinese)

    [15] GUO F, WANG S Y, ZHU P S. Numerical urban climate simulation and preliminary planning strategies of dalian city based on WRF model[C]//The 9th International Conference of Urban Climate, Toulouse, France. 2015.

    [16] 郭飛. 基于WRF/UCM的城市氣候高分辨率數(shù)值模擬研究[J]. 大連理工大學(xué)學(xué)報(bào), 2016, 56(5):502509.

    GUO F. Urban climate high resolution numerical simulation research based on WRF/UCM[J]. Journal of Dalian University of Technology,2016,56(5):502509. (in Chinese)

    [17] 李成才,毛節(jié)泰,劉啟漢,等. 利用MODIS研究中國(guó)東部地區(qū)氣溶膠光學(xué)厚度的分布和季節(jié)變化特征[J]. 科學(xué)通報(bào),2003,19:20942100.

    LI C C, MAO J T, LIU Q H et al, Distribution and seasonal variation research of aerosol optical thickness in eastern China using MODIS[J]. Chinese Science Bulletin, 2003,19:20942100. (in Chinese)

    [18] LIN C Y, CHEN F, HUANG J C, et al. Urban heat island effect and its impact on boundary layer development and landsea circulation over northern Taiwan [J]. Atmospheric Environment, 2008, 42(22): 56355649.

    [19] ADACHI S A, KIMURA F, KUSAKA H, et al. Moderation of summertime heat island phenomena via modification of the urban form in the tokyo metropolitan area[J]. Journal of Applied Meteorology & Climatology, 2014, 53(8):18861900.

    [20] CHING J, BROWN M, BURIAN S, et al. National urban database and access portal tool[J]. Bulletin of the American Meteorological Society, 2009, 90(8): 1157.

    [21] PIU N C, HAO S, FAT L Y. Highresolution visibility and air quality forecasting using multilayer urban canopy model for highly urbanized Hong Kong and the Pearl River Delta[C]//EGU General Assembly Conference Abstracts, 2015, 17: 4451.

    [22] SALAMANCA F, MARTILLI A, TEWARI M, et al. A study of the urban boundary layer using different urban parameterizations and highresolution urban canopy parameters with WRF[J]. Journal of Applied Meteorology and Climatology, 2011, 50(5): 11071128.

    [23] D'ALESSANDRO J. Comparisons of anvil cirrus spatial characteristics between airborne observations in DC3 campaign and WRF Simulations[C]//2015 AGU Fall Meeting, 2015.

    猜你喜歡
    城市環(huán)境
    關(guān)于生態(tài)園林與城市環(huán)境保護(hù)的研究
    城市生活污水處理工藝探析
    淺析國(guó)內(nèi)城市環(huán)境中的涂鴉藝術(shù)
    淺析涂鴉藝術(shù)對(duì)城市環(huán)境的影響
    我國(guó)城市環(huán)境污染及其治理措施
    生態(tài)園林與城市環(huán)境發(fā)展之我見
    城市環(huán)境污染問題及對(duì)策
    對(duì)于園林綠化施工設(shè)計(jì)中采用反季節(jié)種植工藝的探討
    簡(jiǎn)述建設(shè)節(jié)約型綠化園林的策略
    生態(tài)理念在園林施工中的融入
    国产免费av片在线观看野外av| 夜夜看夜夜爽夜夜摸| 在现免费观看毛片| 亚洲av熟女| 97人妻精品一区二区三区麻豆| 亚洲av二区三区四区| 男女那种视频在线观看| 国产亚洲欧美98| 女人十人毛片免费观看3o分钟| 国产精品一区二区三区四区久久| 人妻夜夜爽99麻豆av| 日本黄色视频三级网站网址| 黄色日韩在线| 九九在线视频观看精品| 国产免费一级a男人的天堂| 在线天堂最新版资源| 91字幕亚洲| 自拍偷自拍亚洲精品老妇| 亚洲美女黄片视频| 免费看a级黄色片| 午夜福利在线观看免费完整高清在 | 变态另类成人亚洲欧美熟女| 国产精品乱码一区二三区的特点| 亚洲成av人片在线播放无| 国产探花在线观看一区二区| 国产精品野战在线观看| 欧美潮喷喷水| 国产精品不卡视频一区二区 | 他把我摸到了高潮在线观看| 欧美日本亚洲视频在线播放| 99久久精品热视频| 国产成人福利小说| 性色av乱码一区二区三区2| 啪啪无遮挡十八禁网站| 亚洲成人久久爱视频| 国产欧美日韩一区二区精品| 欧美3d第一页| 久久草成人影院| 精品人妻熟女av久视频| 一级作爱视频免费观看| 精品人妻熟女av久视频| 18禁黄网站禁片午夜丰满| 日韩有码中文字幕| 亚洲人成网站在线播| av在线老鸭窝| 黄色一级大片看看| 97碰自拍视频| 亚洲综合色惰| 久久久久久国产a免费观看| 18禁在线播放成人免费| 欧美一区二区亚洲| 亚洲熟妇中文字幕五十中出| 久久中文看片网| 男人舔女人下体高潮全视频| 欧美又色又爽又黄视频| 1000部很黄的大片| 在线天堂最新版资源| 黄色日韩在线| 日本黄大片高清| 免费看a级黄色片| 婷婷精品国产亚洲av在线| 午夜福利在线在线| 男女做爰动态图高潮gif福利片| 日韩欧美三级三区| 国产精品免费一区二区三区在线| 亚洲精品日韩av片在线观看| 国产精品免费一区二区三区在线| 露出奶头的视频| 国产欧美日韩精品亚洲av| 欧美精品国产亚洲| 自拍偷自拍亚洲精品老妇| 欧美最新免费一区二区三区 | 老司机深夜福利视频在线观看| 观看免费一级毛片| 久久久久精品国产欧美久久久| 国产精品永久免费网站| 高潮久久久久久久久久久不卡| 成人精品一区二区免费| 99久久99久久久精品蜜桃| 又黄又爽又免费观看的视频| 自拍偷自拍亚洲精品老妇| 天堂动漫精品| 在线观看美女被高潮喷水网站 | 午夜福利在线观看吧| 女人十人毛片免费观看3o分钟| 欧美成人性av电影在线观看| 狠狠狠狠99中文字幕| 色播亚洲综合网| 国内精品久久久久久久电影| 成年女人毛片免费观看观看9| 国内久久婷婷六月综合欲色啪| 天堂网av新在线| 成人欧美大片| 丁香六月欧美| 亚洲av不卡在线观看| 少妇人妻精品综合一区二区 | 精品乱码久久久久久99久播| 国产精品av视频在线免费观看| 亚洲精品亚洲一区二区| .国产精品久久| 亚洲欧美日韩东京热| 亚洲国产日韩欧美精品在线观看| 亚洲三级黄色毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一卡2卡三卡四卡精品乱码亚洲| 免费在线观看亚洲国产| 精品熟女少妇八av免费久了| 禁无遮挡网站| 亚洲五月天丁香| 国产伦人伦偷精品视频| 国产精品日韩av在线免费观看| 欧美黄色片欧美黄色片| 一进一出好大好爽视频| 少妇被粗大猛烈的视频| 深夜a级毛片| 亚洲乱码一区二区免费版| 欧美+亚洲+日韩+国产| 国产高潮美女av| 色哟哟哟哟哟哟| 动漫黄色视频在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品合色在线| 欧美绝顶高潮抽搐喷水| 国产伦一二天堂av在线观看| 9191精品国产免费久久| 成人精品一区二区免费| 成人欧美大片| 亚洲av中文字字幕乱码综合| 日日夜夜操网爽| 99热精品在线国产| 欧美一级a爱片免费观看看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品国内亚洲2022精品成人| 日韩欧美国产在线观看| 别揉我奶头 嗯啊视频| 亚洲人成网站在线播| 内地一区二区视频在线| 色av中文字幕| 校园春色视频在线观看| 亚洲第一欧美日韩一区二区三区| 久久久国产成人精品二区| 欧美色欧美亚洲另类二区| 桃色一区二区三区在线观看| 精品无人区乱码1区二区| 18禁黄网站禁片午夜丰满| 可以在线观看的亚洲视频| 午夜精品久久久久久毛片777| av女优亚洲男人天堂| 最近在线观看免费完整版| 99精品久久久久人妻精品| 国产淫片久久久久久久久 | 精品人妻一区二区三区麻豆 | 午夜福利在线观看免费完整高清在 | 欧美国产日韩亚洲一区| 日本在线视频免费播放| 午夜久久久久精精品| 内地一区二区视频在线| 日本在线视频免费播放| 日本三级黄在线观看| 国产久久久一区二区三区| 欧美高清性xxxxhd video| 亚洲不卡免费看| 9191精品国产免费久久| 国产乱人伦免费视频| 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 国产欧美日韩精品亚洲av| 欧美成人性av电影在线观看| 日韩亚洲欧美综合| 午夜福利视频1000在线观看| www.www免费av| 一区二区三区免费毛片| 久久天躁狠狠躁夜夜2o2o| 在线观看午夜福利视频| 午夜免费激情av| 在现免费观看毛片| 国产成人aa在线观看| 日韩 亚洲 欧美在线| 很黄的视频免费| 免费在线观看成人毛片| 最新在线观看一区二区三区| 欧美3d第一页| 欧美日韩亚洲国产一区二区在线观看| 精品99又大又爽又粗少妇毛片 | 一进一出抽搐动态| 一本一本综合久久| 老司机午夜十八禁免费视频| 国产在视频线在精品| 亚洲第一电影网av| 久久99热6这里只有精品| 日本 欧美在线| 国产精品日韩av在线免费观看| 少妇丰满av| 国产高清视频在线播放一区| 久久精品国产亚洲av涩爱 | 麻豆成人午夜福利视频| 男女之事视频高清在线观看| 色吧在线观看| 亚洲美女黄片视频| 精品国产亚洲在线| 欧美日韩瑟瑟在线播放| 国产av麻豆久久久久久久| 亚洲avbb在线观看| av女优亚洲男人天堂| 日韩av在线大香蕉| a级毛片免费高清观看在线播放| 精品一区二区免费观看| 亚洲真实伦在线观看| 全区人妻精品视频| 亚洲午夜理论影院| av在线老鸭窝| 99久国产av精品| 午夜激情福利司机影院| 午夜免费激情av| 18禁黄网站禁片午夜丰满| 久久久久久久久中文| 国产精品一区二区性色av| 国产不卡一卡二| 十八禁网站免费在线| 超碰av人人做人人爽久久| 无人区码免费观看不卡| 欧美国产日韩亚洲一区| 午夜免费男女啪啪视频观看 | 亚洲五月天丁香| 18禁黄网站禁片午夜丰满| 国产成人啪精品午夜网站| 免费av观看视频| 国产极品精品免费视频能看的| eeuss影院久久| 99在线人妻在线中文字幕| 男女视频在线观看网站免费| 国产麻豆成人av免费视频| 黄色丝袜av网址大全| 久久久久久久久中文| 在线观看免费视频日本深夜| 老熟妇仑乱视频hdxx| 久久国产乱子伦精品免费另类| 99热只有精品国产| 90打野战视频偷拍视频| 日韩高清综合在线| 日本三级黄在线观看| 天堂av国产一区二区熟女人妻| 成年版毛片免费区| 五月伊人婷婷丁香| 国产成年人精品一区二区| 91在线精品国自产拍蜜月| 久久久久性生活片| 久久久久久久久久成人| 亚洲中文字幕日韩| 亚洲真实伦在线观看| 国产一区二区在线观看日韩| 免费av毛片视频| 性色av乱码一区二区三区2| 国产成人影院久久av| 日本在线视频免费播放| 一个人看视频在线观看www免费| 精品人妻偷拍中文字幕| 久久久久久久久久黄片| 性插视频无遮挡在线免费观看| 很黄的视频免费| 丁香六月欧美| 日本黄色视频三级网站网址| 日本免费一区二区三区高清不卡| 中文字幕人成人乱码亚洲影| 91在线精品国自产拍蜜月| av黄色大香蕉| 在线观看66精品国产| 搡女人真爽免费视频火全软件 | 久久午夜福利片| 高清在线国产一区| 国产老妇女一区| 亚洲专区中文字幕在线| 精品国产三级普通话版| 91久久精品国产一区二区成人| 两人在一起打扑克的视频| 国产aⅴ精品一区二区三区波| 一个人看的www免费观看视频| 我的老师免费观看完整版| 少妇丰满av| 露出奶头的视频| 国产亚洲av嫩草精品影院| 国产v大片淫在线免费观看| 国产欧美日韩精品一区二区| 波野结衣二区三区在线| 久久欧美精品欧美久久欧美| 黄片小视频在线播放| 欧美成人性av电影在线观看| 国产精品美女特级片免费视频播放器| 欧美又色又爽又黄视频| 日本免费一区二区三区高清不卡| 99久久精品热视频| 亚洲av一区综合| 精品人妻熟女av久视频| 中文字幕人妻熟人妻熟丝袜美| 国产免费男女视频| 久久久久久九九精品二区国产| 黄色女人牲交| 深爱激情五月婷婷| 久99久视频精品免费| 午夜久久久久精精品| 欧美日韩国产亚洲二区| 日韩免费av在线播放| 人人妻,人人澡人人爽秒播| 日韩中字成人| 久久精品国产清高在天天线| 久久亚洲真实| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品98久久久久久宅男小说| 村上凉子中文字幕在线| 国产精品日韩av在线免费观看| 日本一二三区视频观看| 欧美高清性xxxxhd video| 国产免费男女视频| 亚洲av成人不卡在线观看播放网| 久久久久久久久久黄片| 亚洲精华国产精华精| 极品教师在线视频| 欧美日本视频| 日韩免费av在线播放| 国产成年人精品一区二区| 最后的刺客免费高清国语| 成熟少妇高潮喷水视频| 亚洲一区二区三区不卡视频| 亚洲精品亚洲一区二区| 免费高清视频大片| 青草久久国产| 级片在线观看| 国产精品98久久久久久宅男小说| 美女免费视频网站| 制服丝袜大香蕉在线| 99久久精品国产亚洲精品| 小蜜桃在线观看免费完整版高清| 色哟哟哟哟哟哟| 日韩欧美三级三区| www.www免费av| 在线播放国产精品三级| 亚洲经典国产精华液单 | 热99re8久久精品国产| 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 高清日韩中文字幕在线| 午夜福利在线观看吧| 亚洲不卡免费看| 白带黄色成豆腐渣| 色噜噜av男人的天堂激情| 亚洲av熟女| 老司机福利观看| 久久久久国产精品人妻aⅴ院| 亚洲欧美日韩无卡精品| 免费看美女性在线毛片视频| 在线看三级毛片| 久久久久久久久中文| 97碰自拍视频| 久久久久久大精品| 日韩 亚洲 欧美在线| 美女被艹到高潮喷水动态| 国产精品一区二区性色av| 欧美黑人欧美精品刺激| 在现免费观看毛片| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 美女黄网站色视频| 亚洲最大成人手机在线| 国产在线精品亚洲第一网站| 日本黄色视频三级网站网址| 免费电影在线观看免费观看| www.www免费av| 观看免费一级毛片| 不卡一级毛片| АⅤ资源中文在线天堂| 成人av在线播放网站| 窝窝影院91人妻| 日本免费一区二区三区高清不卡| 女生性感内裤真人,穿戴方法视频| 别揉我奶头~嗯~啊~动态视频| 99精品在免费线老司机午夜| 日韩欧美 国产精品| 成年女人永久免费观看视频| 一进一出抽搐动态| 小蜜桃在线观看免费完整版高清| 又粗又爽又猛毛片免费看| 午夜福利免费观看在线| 很黄的视频免费| 中文资源天堂在线| 美女免费视频网站| 91久久精品国产一区二区成人| 午夜久久久久精精品| 欧美又色又爽又黄视频| 嫩草影院入口| 少妇的逼水好多| 亚洲欧美日韩无卡精品| 亚洲成av人片免费观看| 国产成+人综合+亚洲专区| 亚洲av美国av| 99热这里只有是精品在线观看 | 久久久精品欧美日韩精品| 免费在线观看成人毛片| 黄色日韩在线| 九色成人免费人妻av| www.www免费av| 99久国产av精品| 少妇的逼水好多| 国内精品美女久久久久久| 一a级毛片在线观看| 久久久精品欧美日韩精品| 国产亚洲精品久久久com| 欧美日韩亚洲国产一区二区在线观看| 少妇的逼水好多| 精品一区二区三区人妻视频| 国产精品av视频在线免费观看| 听说在线观看完整版免费高清| 午夜免费激情av| 日本在线视频免费播放| 老司机深夜福利视频在线观看| 91字幕亚洲| 亚洲无线在线观看| 国产欧美日韩一区二区三| 69人妻影院| 99久久精品国产亚洲精品| av专区在线播放| 国产极品精品免费视频能看的| 男女做爰动态图高潮gif福利片| 男女之事视频高清在线观看| 热99re8久久精品国产| 免费人成在线观看视频色| 午夜福利在线观看免费完整高清在 | 首页视频小说图片口味搜索| 小蜜桃在线观看免费完整版高清| 悠悠久久av| 亚洲自偷自拍三级| 午夜福利视频1000在线观看| 直男gayav资源| av天堂中文字幕网| 欧美性猛交黑人性爽| 国产av一区在线观看免费| 三级国产精品欧美在线观看| 欧美日韩乱码在线| 国产aⅴ精品一区二区三区波| 久久精品人妻少妇| 欧美在线一区亚洲| 看免费av毛片| 亚洲成人中文字幕在线播放| 禁无遮挡网站| 亚洲,欧美,日韩| 日韩欧美三级三区| 99热这里只有是精品在线观看 | 69av精品久久久久久| 亚洲三级黄色毛片| 亚洲精品成人久久久久久| 亚洲最大成人手机在线| 国产伦在线观看视频一区| 国产91精品成人一区二区三区| 中文字幕久久专区| 久久久久国内视频| 国产在视频线在精品| 国产一级毛片七仙女欲春2| 变态另类丝袜制服| aaaaa片日本免费| 亚洲无线观看免费| 男人舔女人下体高潮全视频| 成人无遮挡网站| 看免费av毛片| 久久人人精品亚洲av| 日韩亚洲欧美综合| 免费看光身美女| 国产黄a三级三级三级人| 精品午夜福利视频在线观看一区| 国产精品一区二区性色av| 国产一区二区三区在线臀色熟女| 免费看日本二区| av福利片在线观看| 99在线人妻在线中文字幕| 欧美色欧美亚洲另类二区| 国产毛片a区久久久久| 黄片小视频在线播放| 热99re8久久精品国产| 亚洲欧美日韩东京热| 一进一出好大好爽视频| 精品一区二区三区视频在线| 国产精品国产高清国产av| 亚洲av成人av| 国内久久婷婷六月综合欲色啪| 此物有八面人人有两片| 高清毛片免费观看视频网站| 国产亚洲精品综合一区在线观看| 国产精品,欧美在线| 欧美极品一区二区三区四区| 国内少妇人妻偷人精品xxx网站| 毛片一级片免费看久久久久 | 在现免费观看毛片| 国产精品国产高清国产av| 欧美日韩中文字幕国产精品一区二区三区| 免费黄网站久久成人精品 | 午夜激情福利司机影院| 亚洲人成网站在线播放欧美日韩| 欧美日韩乱码在线| 窝窝影院91人妻| 亚洲综合色惰| 蜜桃亚洲精品一区二区三区| 日本免费一区二区三区高清不卡| 高清毛片免费观看视频网站| 午夜日韩欧美国产| 亚洲国产欧美人成| 日日摸夜夜添夜夜添小说| 国产av不卡久久| 精品国内亚洲2022精品成人| 91麻豆精品激情在线观看国产| 欧美高清性xxxxhd video| 久久精品国产自在天天线| 精品日产1卡2卡| 男人舔女人下体高潮全视频| 精品久久久久久成人av| 少妇被粗大猛烈的视频| 精品久久久久久成人av| 日本与韩国留学比较| 日韩有码中文字幕| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 日韩精品中文字幕看吧| 日韩欧美免费精品| 国产午夜精品论理片| 偷拍熟女少妇极品色| 中文字幕高清在线视频| 女生性感内裤真人,穿戴方法视频| 免费人成在线观看视频色| 97超视频在线观看视频| 一本精品99久久精品77| 精品国产三级普通话版| 日本撒尿小便嘘嘘汇集6| 99久久久亚洲精品蜜臀av| 午夜免费激情av| 婷婷精品国产亚洲av| 18禁在线播放成人免费| 日日夜夜操网爽| 一进一出好大好爽视频| 国产精品亚洲av一区麻豆| 亚洲中文字幕一区二区三区有码在线看| 免费看美女性在线毛片视频| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 成人性生交大片免费视频hd| 亚洲精品影视一区二区三区av| 亚洲欧美日韩东京热| 嫩草影院新地址| 特级一级黄色大片| 一区二区三区四区激情视频 | 国产精品久久久久久亚洲av鲁大| 又爽又黄无遮挡网站| 久久婷婷人人爽人人干人人爱| 久久久精品大字幕| 亚洲av电影在线进入| 日本熟妇午夜| 97热精品久久久久久| 每晚都被弄得嗷嗷叫到高潮| 能在线免费观看的黄片| av中文乱码字幕在线| 欧美极品一区二区三区四区| 国产野战对白在线观看| 麻豆av噜噜一区二区三区| 国产精品综合久久久久久久免费| 18美女黄网站色大片免费观看| 波野结衣二区三区在线| 亚洲精品久久国产高清桃花| 三级男女做爰猛烈吃奶摸视频| 欧美xxxx黑人xx丫x性爽| 黄色日韩在线| 最近最新免费中文字幕在线| 伦理电影大哥的女人| 亚洲黑人精品在线| 老熟妇乱子伦视频在线观看| 免费搜索国产男女视频| 国语自产精品视频在线第100页| 97热精品久久久久久| 99视频精品全部免费 在线| 午夜福利在线观看吧| 女同久久另类99精品国产91| 午夜老司机福利剧场| 亚洲人成网站在线播放欧美日韩| 欧美成人一区二区免费高清观看| 欧美激情久久久久久爽电影| 日本a在线网址| 国产精品亚洲美女久久久| 露出奶头的视频| 丁香六月欧美| 天堂√8在线中文| 亚洲在线观看片| 丰满人妻一区二区三区视频av| 丰满人妻熟妇乱又伦精品不卡| 久久精品久久久久久噜噜老黄 | 国内少妇人妻偷人精品xxx网站| 青草久久国产| 99久久久亚洲精品蜜臀av| 一级作爱视频免费观看| 久久热精品热| 国内精品久久久久久久电影| av天堂在线播放| 国产日本99.免费观看| av天堂在线播放| 又爽又黄无遮挡网站| 欧美在线黄色| 亚洲不卡免费看| 国产私拍福利视频在线观看| av在线老鸭窝| 国产国拍精品亚洲av在线观看| 99久久精品一区二区三区| 欧美性猛交黑人性爽| 91在线观看av| 嫩草影院新地址| 一个人看的www免费观看视频| 最后的刺客免费高清国语| 身体一侧抽搐| 国产中年淑女户外野战色| 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 国产精品免费一区二区三区在线|