在线观看国产区-在线观看国产欧美-在线观看国产免费高清不卡-在线观看国产久青草-久久国产精品久久久久久-久久国产精品久久久

美章網 資料文庫 霧天氣高度湍流輸送特征的實驗范文

霧天氣高度湍流輸送特征的實驗范文

本站小編為你精心準備了霧天氣高度湍流輸送特征的實驗參考范文,愿這些范文能點燃您思維的火花,激發您的寫作靈感。歡迎深入閱讀并收藏。

霧天氣高度湍流輸送特征的實驗

《北京大學學報》2016年第二期

摘要

利用天津255米氣象塔層大氣邊界層觀測資料,分析霧日各氣象要素的特征,研究湍流輸送規律。結果表明:霧前,大氣濕度較大,逆溫層高度約為100m,輻射霧發生前的逆溫強于平流霧;霧中,逆溫層持續變強、增厚;霧后,逆溫層出現抬升,大氣呈近中性偏不穩定的層結特征。輻射霧過程的逆溫現象比平流霧明顯;輻射霧霧頂較低,平流霧較高;輻射霧的消散主要受溫度影響,平流霧主要受風速影響。霧天氣過程改變了大氣層結結構,夜間可能呈現不穩定層結狀態;霧前和霧中不同高度湍流垂直輸送微弱,霧過程后期的水平輸送突然增強。輻射霧霧前不同高度的平均動能數值較小,霧中呈增大趨勢,霧消散階段逐漸增大,湍流動能的增大是影響輻射霧消散的重要因素。霧前不同高度平均動能與湍流動能比值的突然增大可能是霧發生的湍流信號之一,比值劇增之后降至霧前水平則為霧消散的信號。

關鍵詞

大氣邊界層;霧天氣;湍流動能;天津

霧是由于空氣中懸浮著大量的微小水滴或冰晶,使水平能見度低于1km的天氣現象。霧是對人類交通活動影響最大的天氣現象之一:霧時的能見度降低會導致很多交通工具(如飛機等)無法正常使用,甚至造成意外事故,還會給人們的日常生活及農作物生長帶來不利影響,也會對人的健康產生很大影響[12]。揭示霧的發生、發展和消散規律,正確認識霧生消過程中的湍流輸送特征,對提高霧的生消預報和人工影響霧實驗方案的設計是非常必要的。由于獲取大氣邊界層觀測資料和快速響應資料較困難以及霧天大氣邊界層結構較復雜等原因,霧的研究多集中在大尺度環流背景[35]、中小尺度地形影響、霧天大氣邊界層特點[610]、霧中微觀物理量演變特征[11]、垂直湍流輸送特征[1213]、霧與環境的相互影響[14]以及不同地區霧過程的特點[1517]等方面,霧中邊界層特征研究和微觀物理量(如霧滴譜等)的研究則集中在輻射霧方面,霧中湍流相似性理論驗證及輸送規律的研究較少涉及。

根據生成條件可以將霧分為輻射霧、平流霧、山地霧和鋒面霧等。霧發生、發展和消散過程與湍流運動密切相關[1819]。湍流運動既可促進霧的發展,也可導致霧的消散。Zhou等[20]認為湍流和輻射及重力沉降之間存在微妙的平衡,由于輻射冷卻而產生的液態水會隨著湍流輸送作用的增強,上升并蒸發消散。Nakanishi[21]研究了湍流輸送對霧的形成和消散機制。吳彬貴[22]指出,霧過程中平均動能較小,但湍流擾動活躍,湍流動量輸送以垂直方向為主;平均動能和湍流動能在霧前出現異常增強,可能是平流霧的啟動信號。輻射霧過程中,低空往往會伴隨著逆溫層的出現,近地面大氣層結常常呈現不穩定特征[23]。Argentini等[24]評估了小風弱切變條件下城市平流霧過程。本文針對霧天氣過程與湍流相互作用研究的相對不足,利用天津255米氣象塔層大氣邊界層探測資料,對比分析城市地區不同霧天氣過程氣象要素廓線演變規律和大氣湍流特征,探討湍流輸送對霧發生、發展和消散的影響機制,試圖找到預測霧發生的湍流信號。

1數據獲取與處理

天津市位于華北平原東部,鄰近渤海。天津255米氣象鐵塔位于天津南郊(39.06°N,117.1°E),海拔3.3m。該氣象鐵塔周邊500m內沒有高層建筑,300m內建筑物高度控制在3~30m之間。氣象鐵塔東側為商業區(建筑物高10~15m),南側和西側為居民區(建筑物高15~20m),北側約100m處有20~30m高的建筑,基本上代表了城區復雜下墊面條件[25]。氣象鐵塔設置15層觀測平臺,高度分別為5,10,20,30,40,60,80,100,120,140,160,180,200,220和250m。每個觀測平臺在南北兩個方向上各設有伸展臂,長度約4m,進行同步觀測。每層平臺均安裝有風速、風向和溫濕度傳感器。風、溫、濕梯度觀測儀分別采用長春氣象儀器廠生產的風速風向儀和溫度濕度儀。風速、風向、溫度和濕度觀測為全天候、連續和自動觀測,采樣間隔為20s。

采用美國CAMPBELL公司生產的CSAT3型超聲風溫儀,安裝在南伸展臂頂端,開口向東,采樣頻率為10Hz。2006年1月1日,40m高度的超聲風溫儀正式開始自動連續觀測。2007年,在40m高度同方位加裝一臺CAMPBELL公司生產的KH20濕度快速響應儀,采樣頻率為10Hz。能見度觀測采用BelfortModel600能見度監測儀,安裝在距地面2m高度處,采樣間隔為1min。2010年,氣象塔觀測設備全部更新后,塔上15層平臺風向風速、溫度和濕度觀測儀器仍采用長春氣象儀器廠設備;40m平臺有CAMPBELL公司CSAT3型超聲風溫儀和LI-7500二氧化碳水汽分析儀,能見度觀測沿用BelfortModel600能見度監測儀。地面自動站溫濕度傳感器采用芬蘭Väisälä公司生產的HMP45D型傳感器。考慮到霧天氣邊界層氣象要素演變和城市下墊面變化的影響。本文針對2010-2011年秋冬季4次霧天氣過程(表1)分別分析霧發生前、發展中和消散過程的邊界層結構特點。每次霧天氣過程包括霧形成前和霧消散后數小時的資料。所有資料都進行嚴格的質量控制,包括湍流數據預處理、湍流宏觀統計量計算、微觀特征量計算等。其中,大氣湍流數據預處理包括剔除野點、趨勢項去除、數字濾波、坐標旋轉等。湍流數據經過預處理后,去除平均風速小于0.5m/s的資料組,取整點后30min數據滑動進行湍流分析。

2背景場分析

圖1給出2010-2011年10,11和12月的平均溫濕風廓線。整體上看,10月氣溫在0~100m高度表現為下高上低,20m以下高度溫度梯度較大,20~100m高度之間溫度梯度較小,部分層結呈現近中性。100~120m高度和160~220m高度常出現逆溫層。查詢氣溫日變化(圖略)可知:最高氣溫大約出現在北京時15:00左右;最低氣溫一般出現在05:00,較高層的最低氣溫可延遲0.5~1小時出現,高層氣溫隨時間的變化較低層小一些。11和12月與10月的溫度變化規律大體上一致,逆溫層高度和最高/低氣溫出現的時間和量值略有差別。平均風速廓線結果表明:10月風速呈低層小、高層大的趨勢,并且白天(08:00-18:00)高低空風速差較夜間小(白天差值約3m/s,夜間約7m/s,圖略)。低空,風速隨時間變化較小;高空,風速隨時間的變化較大。80~100m高度是上述高低空風速變化趨勢的大致拐點。11和12月風速與10月整體趨勢較一致,主要差別在于高低空風速差值較小,高低空風速變化趨勢的拐點大致在100m高度。平均相對濕度廓線表明:10月平均相對濕度在50%~70%之間。低空相對濕度晝夜差值較大,夜間22:00至次日08:00,相對濕度在65%以上,午時,相對濕度一般在40%~50%;高空相對濕度晝夜差值較小,一般在5%左右(圖略)。11和12月相對濕度變化趨勢與10月基本上一致。

3霧天氣過程溫濕風演變分析

從圖2的能見度和相對濕度變化情況可以確定,2010年10月22日22:00—10月23日06:00和10月24日02:00-06:00期間出現兩次霧過程(個例001和002),最低能見度分別為0.54和0.36km。查詢歷史天氣可知,24日下午12:00-17:00時段內出現小雨天氣。

3.1個例001分析綜合圖3中相對濕度、比濕、氣溫和風速的時間‒高度分布情況,分析22日22:00—23日06:00的霧過程的邊界層特征以及霧體的特點,可以得到以下結果。霧發生前數小時,整層大氣比濕保持在一個較大的數值,約9g/kg;相對濕度同時出現明顯增大,達到85%以上,低層大氣甚至達到90%。150m左右高度出現逆溫層,并逐漸發展、變厚、下壓。低空風速持續低于1m/s,有利于霧的形成,并且霧發生前十幾個小時風速都較小。風向由霧前的北風逐漸轉為東偏南風(圖略)。由于天津地區東臨渤海,東風帶來的暖濕氣流是霧形成的水汽條件,也是此次霧形成的可能原因之一。

在此次霧發生發展過程中,大氣的相對濕度幾乎處于飽和狀態,比濕則比霧前有減小的趨勢。在整個霧過程中,比濕與相對濕度存在一定的反對應關系,即相對濕度明顯變大的時刻和高度,比濕一般都呈現減小的趨勢,但兩者都比平時的水平高,這與前人關于輻射霧的結論[25]一致。霧發生前中期,逆溫層持續變強、增厚,并下壓逐漸接地;而霧過程的后期,逆溫層底出現抬升的情況,大氣呈現弱不穩定層結特征。整個霧發生過程中,霧區低層風速較小(約1m/s),直到霧臨近消散,才有緩慢增加的趨勢。根據相對濕度,低于80%可以判定為非霧區,不難得出結論:此次霧過程的大部分時間段,霧頂高度在100~150m之間,整個霧層較淺,只在霧發展的后期,霧頂才略有抬升。霧消散前,隨著逆溫層底的抬升,逆溫強度也在減弱,直至霧消散前1~2小時,大氣呈現弱不穩定、近中性的特征,對流增強,加快了霧的消散速度。06:00以后,隨著太陽輻射逐漸出現,大氣緩慢升溫,湍流交換變強,霧加速消散。并且,從23日06:00開始,高層的極大風速逐漸下壓,低層的風速也在緩慢增加,空氣流動加強,促使霧的消散。可以認為,溫度和風速是影響本次霧消散的重要因素。

3.2個例002分析2010年10月24日02:00-06:00的霧過程(個例002)與22日22:00-23日06:00的霧過程(個例001)時間間隔比較短,為了更好地對比不同霧日的特點,兩次霧過程均呈現在圖3中。對個例002的相對濕度、比濕、氣溫和風速的變化進行分析,可以看出002號霧過程具有以下特點。霧前3~4個小時,大氣的相對濕度也明顯增大,達到80%以上,低層大氣的相對濕度霧前接近90%。霧前逆溫層開始出現并逐漸發展、增厚、下壓。低空風速持續較低,有些時段風速低于1m/s,甚至無風。霧前風向雖然持續變化,但是都有偏東的風向,為此次霧的形成提供了水汽。此次霧發生發展過程中,大氣的相對濕度幾乎處于飽和狀態,相比于霧前比濕則略微減小。在霧發生一小時后,逆溫層底出現抬升的情況,大氣呈現弱不穩定層結特征,逆溫層減弱。與個例001不同的是,此次霧過程逆溫層較弱,逆溫層底的下壓、抬升也不明顯,不易通過逆溫層判斷霧頂的高度。但是,霧前以及霧中逆溫層的出現和發展依然是霧過程的重要特點。從02:00開始,霧區風速持續增大,達到6m/s左右;06:00以后的霧消散期間,在200~250m高空出現風的極大值。本次霧過程期間,風速一直較大,霧中一段時間有東偏南風向,暖濕氣流提供了水汽輸送。可以判定此次霧為平流霧。此次霧過程的大部分時間段,霧頂高度在250m以上。霧發展的初期高度較低,之后逐漸升高,直到大于250m。霧消散前,大氣呈現弱不穩定,近中性,對流增強,利于霧的消散。06:00以后,隨著太陽輻射逐漸增強,大氣緩慢升溫,湍流交換變強,加快了霧的消散速度。同時風速在高層達到一個極大值,空氣流動加強,促使霧消散。可以認為風速是影響此次霧消散的重要因素。

3.3普適性分析為了驗證上述結果的普適性,我們另外選取兩次霧天氣個例進行分析。分別是2010年11月19日02:00-08:00的輻射霧過程(記為003)和2011年11月28日04:00-10:00的平流霧過程(記為004)。將003與001兩次輻射霧過程進行對比(圖略),相同點是:1)霧前比濕、相對濕度較平均狀態大,風速較平均狀態小,利于霧的形成;霧前出現逆溫層,并且持續發展、增厚;2)霧中相對濕度持續增大,但比濕略有減小,保持較小風速,利于霧的維持,霧初期逆溫層變強、增厚、下壓,霧后半段逆溫層底抬升,逆溫層減弱;3)溫度變化與001號霧過程比較一致,溫度是導致霧消散的主要原因。不同點是:1)003號霧過程的霧前比濕小于001號霧過程的9g/kg,相對濕度也較小,003號霧過程能見度只略小于1km;2)霧中相對濕度略小于90%,與001號霧過程中95%以上的高濕狀態有一定區別,霧頂高度在200m左右,略低于001號霧過程;3)霧消散階段風速增大現象不如001號霧過程明顯。將004與002兩次平流霧過程進行對比(圖略),相同點是:1)霧前比濕、相對濕度較平均狀態大,逆溫層和小風速均利于霧的形成;2)霧中相對濕度持續增大,比濕略有減小,風速有增大趨勢,利于平流霧的發生,霧初期逆溫層變強、增厚,并下壓逐漸接地,霧后半段逆溫層底抬升,逆溫層減弱;3)004號霧過程霧消散時的溫度和風速變化與002號霧過程趨勢一致,風速是導致霧消散的主要原因。不同點是:1)004號霧過程霧前比濕和風速略有差異,2)霧中風速尤其是低層的風速比002號霧過程中小,平流不強,初期霧頂高度在150m左右,很快霧頂就上移到250m以上。

4秋冬季霧天的湍流輸送特征研究

4.1個例001分析圖4給出三層高度(40,120和220m)水平熱量輸送和垂直熱量輸送。可以看出,不論是霧前還是霧持續過程中,垂直輸送都非常弱,水平輸送則一直存在,尤其是霧過程的后半段,40和120m高度還出現較強的水平熱量輸送,有利于霧的消散。檢查三層高度穩定度參數(圖略),2010年10月22日白天,大氣層結都以不穩定層結為主。霧前一小時直至霧開始消散時間段,三層穩定度都呈現明顯的穩定層結與不穩定間層結振蕩的趨勢,前半時間段以穩定層結為主,后半段時間以不穩定間層結為主。從霧前后層結穩定度參數的演變看,本次霧天氣的存在改變了大氣原有的層結結構,降低了夜間的層結穩定度,使夜間部分時段呈現不穩定層結狀態。這一點與前人研究結果相似:輻射霧過程中由穩定層結演變為不穩定層結[2627]。然而,水汽輸送結果(圖5)顯示,不論是霧前還是霧持續過程中,垂直水汽輸送都非常弱,以水平輸送為主,尤其是霧過程的前段和后段,120m高度和220m高度出現較強的水平水汽輸送,說明本次霧的水汽來源以平流輸送為主。即使在霧結束階段,水平輸送量也能達到垂直輸送的3倍,甚至更多。值得注意的是,23日02:00至霧開始消散的這段時間,三層水汽水平輸送都有增大的趨勢,但是40m高度是正方向,而120m高度和220m高度是負方向。圖6(a)是三層高度平均動能隨時間的變化。整體來看,16:00以前,三層高度平均動能都不大。

16:00以后,平均動能有一個快速的升高過程,并很快達到峰值,40m高度在4m2/s2左右,120m高度在7m2/s2左右,220m高度達到8m2/s2以上。之后又有降低的趨勢,其中40m高度和120m高度平均動能的降低一直持續到霧前的22:00,基本上達到最小值。但是,220m高度的平均動能卻在19:00開始再次上升,并且上升趨勢持續到霧前,并在霧后1小時左右達到極大值。霧的形成需要在較低的平均動能背景下,此次霧過程由于120m及以上高度的平均動能過大,不利于霧的形成和保持。我們可以推測霧頂的高度應該在100m左右,只在中后期短時內霧頂高度可能有所提升。霧消散前,平均動能較霧中期的平均值小,在開始消散后的2小時左右才達到較大的值。可以認為,平均動能對此次的霧消散有影響,但不是決定性因素。圖6(b)是三層高度湍流動能隨時間的變化。可以看出,22日白天三層湍流動能的波動都較大,波動趨勢基本上一致。霧前的幾個小時,三層湍流動能均開始降低;霧中期三層湍流動能均在波動中穩步上升,在霧臨近消散前達到一個較大的數值。湍流動能減小到一定數值是霧發生的條件之一,可以認為,湍流動能的增大是影響霧消散的重要因素。圖7是40,120和220m三層高度平均動能與湍流動能的比值演變情況。不難發現,霧前數小時,比值突然增大到30倍左右,并保持一定的時間,直到霧發生。霧持續的中間時段,比值有一個劇增,可達100倍以上。之后,比值降低到霧前的水平,霧開始消散,消散期間,比值較小。在霧頂高度以下的40m和120m高度,這一特征更為明顯。吳彬貴等[28]針對平流霧的研究指出:霧前的湍流動能和平均動能異常,且前者大于后者,可能是平流霧的啟動信號。根據對三層高度平均動能與湍流動能比值的分析,可將二者比值的突然增大看做輻射霧發生的一個信號,而霧過程中比值在劇增之后降低到霧前的水平,可以作為輻射霧即將消散的信號(通過下文的分析可以看出,這一比值也可以作為平流霧發生和消散的信號)。

4.2個例002分析參照圖5~7后半時段的結果,我們對002個例中平流霧過程的湍流輸送特征進行分析。分析圖4中三層熱通量垂直輸送的方向和強度,可以看出:04:00-08:00,40m高度向上輸送,120m和220m高度向下輸送,120m高度熱量垂直輸送強度大于220m和40m高度;08:00-12:00,三層高度都向上輸送,40m高度的輸送強度大于120m和220m高度的輸送強度。觀察溫度的梯度情況(圖3),可以看出120m高度的溫度梯度較大,而220m和40m高度較弱,這可能是導致輸送強度不同的重要原因。08:00-12:00,在太陽輻射的加熱作用下,整層大氣呈現弱不穩定層結,熱量向上輸送,但強度較弱。分析圖4中水平方向三層高度熱量輸送的情況可知,23日08:00-14:00,40m高度以負方向輸送為主,平均值在0.1K•m/s左右,120m和220m高度則以正方向輸送為主,輸送強度比40m高度大,其中220m高度的峰值達到0.4K•m/s,平均值在0.2K•m/s左右。14:00以后一直到霧發生的前兩個小時,從整體上看,水平輸送很微弱,大部分時段幾乎為零。04:00至霧消散階段,三層的水平熱量輸送相對之前都增強,其中40m高度基本上保持負方向輸送,數值持續增大;120m和220m高度則表現出正負方向交替,峰值和平均值都較大。

檢查三層高度穩定度參數(圖略),2010年10月23日20:00-24日02:00(霧前),三層高度都呈現明顯的穩定層結與不穩定間層結振蕩的趨勢,前半時段以穩定層結為主,后半時段以不穩定層結為主。霧持續期間,40m高度由不穩定逐漸過渡到近中性并持續到霧消散,而120m和220m高度則基本上保持穩定層結。這與持續的較大風速有關,也受到溫度廓線的影響。直到08:00以后,在更大的風速影響下,三層高度都調整到近中性并保持到12:00以后開始降水。從三層高度水汽輸送對比(圖5)可以看出,霧前數小時直至霧前半時段,垂直輸送一直都非常弱,以水平輸送為主。霧后半時段,各層都出現較大的垂直輸送和水平輸送,但這可能與霧滴附著在儀器表面有關。整體上以水平輸送為主,平均值約為垂直輸送的3倍,甚至更多。通過比較3個高度層的平均動能(圖6(a))可以看出,三層高度的平均動能在霧前都有一個相當長時間的低值區,霧的出現也伴隨著平均動能的迅速增大。當平均動能增大到一定數值之后,霧開始消散。由此也可以判定,霧的形成是由比較強的平流引起的,而霧的消散同樣是大風引起的,這非常符合平流霧的定義和特征。從圖6(b)后半段002個例三層高度湍流動能來看,霧前較長時間湍流動能很小,隨著霧的發生迅速增大。當湍流動能增大到一定數值后,霧開始消散。這也是典型的平流霧的特征。值得注意的是,在霧的消散階段,可以發現湍流動能呈現上低下高的狀態,即40m高度湍流動能>120m高度湍流動能>220m高度湍流動能。從圖7后半段三層高度平均動能與湍流動能的比值演變情況可以看到,23日傍晚一直到霧前,比值都有變大的趨勢,三層高度的比值平均增大3倍左右。霧中期,40m高度和220m高度的比值明顯有一個突增,分別可達40和150,之后又迅速降低到霧前的水平,很快霧開始消散,消散期間的比值也基本上保持在霧前較大的水平。120m高度變化不明顯,但也可以看到這個趨勢。比值的突然增大可以看做霧發生的一個信號;比值在劇增之后降低到霧前的水平,可以作為霧即將消散的信號。

4.3普適性分析同樣將003號與001號霧過程、004號與002號霧過程進行對比分析(圖略),得到以下一致性的結論。1)綜合三層熱量輸送情況來看,不論是霧前還是霧持續過程中,垂直輸送都非常弱,水平輸送則一直存在,可見霧中的熱量輸送以水平輸送為主。2)不論是霧前還是霧持續過程中,垂直水汽輸送都非常弱,以水平水汽輸送為主。3)霧過程中,三層高度的平均動能都有增大趨勢;霧消散階段,三層的平均動能先略有減小然后逐漸增大。可以看出,霧的形成需要在較低的平均動能背景下。霧消散前,平均動能數值增大到較大的數值,對此次的霧消散有很大的影響。4)霧前幾小時,三層高度的湍流動能一致降低;霧持續過程中,雖然各層湍流動能先在短時間內降低,但之后穩步上升,在達到一定的數值后,霧開始消散。湍流動能減小到一定的數值是霧發生的條件之一,湍流動能的增大是影響霧消散的重要因素。5)霧頂以下高度的平均動能與湍流動能比值在霧消散前均增大。

5結論

本文采用天津市255米氣象鐵塔不同高度的常規氣象和湍流觀測資料,對秋冬季節4次霧過程的大氣邊界層結構和湍流輸送特征進行分析。得到以下結論。1)霧前數小時,大氣邊界層整層的比濕和相對濕度數值都較大;逆溫層發生并逐漸發展,其中輻射霧霧前的逆溫強度大于平流霧霧前的逆溫強度;霧前風速較小,平流霧發生前1~2小時呈增大趨勢。霧中逆溫層厚度和強度持續發展,大氣呈弱不穩定層結特征,輻射霧的逆溫強度高于平流霧。輻射霧霧區的低風速現象持續到霧消散時才略有增加;平流霧后期則往往出現風速的增大。輻射霧霧頂高度(150~200m)略低于平流霧(大于250m)。輻射霧的消散主要受溫度的影響,而平流霧消散的主要影響因素是風速。2)霧前和霧過程中,熱量和水汽的垂直輸送較弱,水平輸送相對較強。逆溫層的發展導致霧前期的大氣層結狀態以穩定層結為主。霧天氣的存在降低了夜間的大氣層結穩定度,部分時段呈現弱不穩定層結狀態。霧過程中,不同高度的平均動能有增大趨勢;霧消散前,平均動能數值有明顯增加。湍流動能的增加是霧消散的重要因素之一,平均動能與湍流動能比值的突然增大可以視為霧發生的湍流信號,而霧過程中出現的比值數值的迅速減小可以作為霧即將消散的湍流信號。

作者:葉鑫欣 魏偉 李航 張宏升 單位:北京大學物理學院大氣與海洋科學系 氣候與海氣實驗室

主站蜘蛛池模板: 中文国产成人精品少久久 | 亚洲精品福利在线观看 | 国产日本在线播放 | 亚洲免费观看网站 | 久久国产精品99久久小说 | 亚洲国产成人久久一区久久 | 亚洲视频中文字幕在线观看 | 亚洲fuli在线观看 | 久久精品国产精品2020 | 亚洲精品视 | 四虎在线最新永久免费播放 | 婷婷爱五月天 | 国色天香社区视频免费版 | 亚洲欧美久久婷婷爱综合一区天堂 | 久热这里只有精品视频6 | 欧美视频亚洲视频 | 在线a亚洲视频播放在线观看 | 在线亚洲免费 | 91精品综合久久久久久五月天 | 在线黄色小视频 | 六月婷婷七月丁香 | 精品国产综合区久久久久99 | 成人羞羞视频国产 | 亚洲精品蜜桃久久久久久 | 日韩欧美视频二区 | 亚洲免费一区 | 日韩成人免费在线视频 | 亚洲人成人网毛片在线播放 | 亚洲网站免费看 | 国色天香社区视频免费版 | 欧美亚洲一区二区三区四 | 九色亚洲 | 久久综合一区二区三区 | 国产噜噜噜视频在线观看 | 久久精品资源 | 美女被爆羞羞视频网站视频 | 亚洲午夜视频在线 | 男女网站免费 | 尤物精品视频一区二区三区 | 亚洲国产成人在人网站天堂 | 精品一区二区三区水蜜桃 |