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

    P-RSHMU蒸發(fā)波導(dǎo)預(yù)測與湍流通量相關(guān)性分析

    2023-11-10 06:35:50張利軍張超王紅光郭相明韓杰張守寶
    海洋預(yù)報 2023年5期
    關(guān)鍵詞:模型

    張利軍,張超,王紅光,郭相明,韓杰,張守寶

    (中國電波傳播研究所,山東青島 266107)

    0 引言

    近海面水汽會隨著高度迅速下降而形成蒸發(fā)波導(dǎo),該類型大氣波導(dǎo)發(fā)生概率高,因此受到雷達(dá)以及電波傳播學(xué)界的極大關(guān)注。這種反常的傳播環(huán)境會極大地影響海洋環(huán)境中的雷達(dá)、通信、電子對抗設(shè)備的性能,尤其會使得雷達(dá)實現(xiàn)超視距傳播或者形成雷達(dá)盲區(qū)等。從20 世紀(jì)60 年代開始,伴隨著微氣象學(xué)的發(fā)展,國內(nèi)外的科研工作者基于海上氣象觀測數(shù)據(jù),以相似理論為基礎(chǔ)開始了蒸發(fā)波導(dǎo)模型的研究分析,相繼提出了JESKES 模型[1]、ROTHERAM 模 型[2]、FAIRALL 模 型[3]、PAULUSJESKES 模型[4]、RSHMU(Russian State Hydrometeorological University,俄羅斯國立水文氣象大學(xué))模型[5]、MGB(Musson-Gauthier-Bruth)模 型[6]、BYC(Babin-Young-Carton)模 型[7]、NPS(Navy Postgraduate School,海軍研究生院)模型[8]以及偽折射率模型[9]等。

    RSHMU 模型公開發(fā)表時間較晚,針對此模型的適用性分析文獻(xiàn)較少。由于文獻(xiàn)[5]中的溫度濕度剖面穩(wěn)定度函數(shù)存在編輯錯誤,導(dǎo)致張永剛等[10-11]在錯誤地引用了這些剖面穩(wěn)定度函數(shù)后,得到RSHMU 模型預(yù)測風(fēng)速特征參數(shù)的結(jié)果較差以及溫度、濕度剖面函數(shù)在南海海區(qū)不適用的結(jié)論。張利軍等[12]更正了溫度、濕度剖面穩(wěn)定度函數(shù)的編輯錯誤,并基于正確的剖面穩(wěn)定度函數(shù)給出了該模型在不穩(wěn)定條件下的高度預(yù)測,結(jié)果優(yōu)于美國NPS 模型。進(jìn)一步地,張利軍等[13]借鑒了PAULUSJESKES 的訂正方法,給出了改進(jìn)的RSHMU 模型(即P-RSHMU 模型),該模型在不穩(wěn)定條件下的預(yù)測性能保持不變,穩(wěn)定條件下的預(yù)測性能獲得極大提高,因此基于P-RSHMU 蒸發(fā)波導(dǎo)模型,能夠給出適合穩(wěn)定條件以及不穩(wěn)定條件下的合理蒸發(fā)波導(dǎo)高度。

    郭相明等[14]基于海上梯度塔觀測數(shù)據(jù),分析了多種蒸發(fā)波導(dǎo)預(yù)測模型的適用性。郭相明等[15]還應(yīng)用湍流通量分析了蒸發(fā)波導(dǎo)的形成機(jī)理,并利用TOGA-COARE 實驗期間測試的1992 年11 月25—26 日的一組不穩(wěn)定大氣時的實驗數(shù)據(jù),定性分析了蒸發(fā)波導(dǎo)的生成和變化與代表湍流運動的湍流通量的關(guān)系。眾所周知,海氣界面通量交換是海氣相互作用的重要環(huán)節(jié),海氣界面通量的確定在海氣相互作用研究方面具有重要意義,受到海洋學(xué)界和氣象學(xué)界的廣泛重視。張超等[16]從基本的氣象觀測資料出發(fā),利用經(jīng)海上觀測數(shù)據(jù)驗證有效的耦合海洋-大氣響應(yīng)試驗(Coupled Ocean-Atmosphere Response Experiment,COARE)通量算法[17-18]以及改進(jìn)的NPS蒸發(fā)波導(dǎo)預(yù)測模型,計算獲得了南海北部海域海面湍流通量以及蒸發(fā)波導(dǎo)高度,并定量分析了二者之間的相關(guān)性。由于穩(wěn)定條件下低風(fēng)速時改進(jìn)的NPS 模型計算可能出現(xiàn)蒸發(fā)波導(dǎo)高度達(dá)40 m 的情形,進(jìn)而可能影響穩(wěn)定條件下蒸發(fā)波導(dǎo)高度與湍流通量的相關(guān)性計算,因此本文引入預(yù)測精度較高的P-RSHMU 蒸發(fā)波導(dǎo)模型,基于海上氣象水文數(shù)據(jù)計算蒸發(fā)波導(dǎo)高度,有效抑制了穩(wěn)定條件下蒸發(fā)波導(dǎo)高度預(yù)測不合理的現(xiàn)象。進(jìn)一步地,討論了PRSHMU 模型預(yù)測高度與張超等[16]給出的湍流通量結(jié)果的相關(guān)性。這些結(jié)果可為深入理解蒸發(fā)波導(dǎo)的形成機(jī)理和內(nèi)涵提供參考。

    1 湍流通量、蒸發(fā)波導(dǎo)高度計算方法

    1.1 湍流通量計算方法

    海氣界面通量的計算方法通常分為直接計算方法和整體參數(shù)化(Bulk Parameter)方法。直接計算方法是根據(jù)通量的定義直接測量氣象要素的湍流脈動量并計算其二階矩,該方法是相對最準(zhǔn)確的方法,但由于實際中會受到人力、物力以及觀測手段的影響,因此通常采用整體參數(shù)化方法來估計海氣界面的湍流通量[17]。整體參數(shù)化方法依賴于風(fēng)速、溫度、濕度等基本氣象參數(shù),目前普遍使用的是計算湍流通量的整體公式,較為先進(jìn)且普遍認(rèn)可的通量算法是COARE 算法,該算法可以計算風(fēng)速、溫度、比濕的特征參數(shù),進(jìn)而計算海氣湍流通量。海氣湍流通量關(guān)系式定義如下[18]:

    式中:M為動量通量,單位:N/m2;HS為感熱通量,又稱顯熱通量,單位:W/m2;HE為潛熱通量,單位:W/m2;ρ為空氣密度,單位:kg/m3;cp為空氣定壓熱容,單位:m2/(s2·K);LE是水的蒸發(fā)潛熱,單位:J/kg;w'、u'、θ'、q'分別表示垂直風(fēng)速、水平風(fēng)速、位溫、比濕的湍流脈動值;u*、θ*、q*分別為莫寧奧布霍夫相似(Monin-Obukhov Similarity,MOS)理論中的風(fēng)速、位溫、比濕特征量。

    1.2 P-RSHMU蒸發(fā)波導(dǎo)高度計算方法

    蒸發(fā)波導(dǎo)高度采用P-RSHMU 模型[13]計算。近地層的溫度剖面以及比濕剖面計算公式如下[19]:

    式中:Tsea、qs分別為海表溫度(單位:K)以及海面比濕(單位:kg/kg);θ*、q*分別為位溫特征量、比濕特征量;κ為卡曼常數(shù),取0.4;z0θ、z0q分別為溫度粗糙度和比濕粗糙度分別為溫度普適函數(shù)和比濕普適函數(shù);Γd為干絕熱遞率,取0.009 76 K/m;z為高度,單位:m。

    需要特別指出的是,溫度和濕度剖面穩(wěn)定度函數(shù)的正確表達(dá)式為[12]:

    式中:β1= 8;y1=(1 -β2ξ)13;β2= 35;γ2= 6。

    大氣壓力剖面由流體靜力學(xué)公式以及理想氣體定律聯(lián)合確定,積分得到如下壓力剖面公式[19]:

    式中:R為干空氣理想氣體率常數(shù)287.04 J/(kg·K);為層高度z1和z2處平均虛位溫。

    由比濕以及壓力剖面計算水汽壓e剖面[19]:

    式中:ε為干空氣氣體常數(shù)以及水汽常數(shù)之比,取0.621 7)。一旦確定了溫度、水汽壓、壓強(qiáng)剖面,修正折射率[19]剖面可以按照下式計算得到:

    通過尋找M 剖面中的最小值即可確定對應(yīng)的蒸發(fā)波導(dǎo)高度。進(jìn)一步地,當(dāng)氣海溫差小于-1 ℃時,采用RSHMU 模型計算方法;當(dāng)氣海溫差大于-1 ℃時,保持海溫不變,利用RSHMU 模型分別計算氣海溫差等于-1 ℃以及等于0 ℃的蒸發(fā)波導(dǎo)高度h-1和h0,如果h0>h-1,則蒸發(fā)波導(dǎo)高度為h-1,否則依照原來的氣溫與海溫來計算蒸發(fā)波導(dǎo)高度[13]。該方法即為蒸發(fā)波導(dǎo)P-RSHMU模型的計算方法。

    2 數(shù)據(jù)資料與計算結(jié)果

    2.1 數(shù)據(jù)資料

    依托廣東茂名博賀海洋氣象科學(xué)試驗基地的海上試驗平臺[20],采用要素集成的氣象傳感器以及海溫傳感器,以分鐘級為單位,采集了2017 年11 月26 日—2018 年5 月31 日的氣象數(shù)據(jù),包含氣溫、海溫、相對濕度、氣壓、風(fēng)速、風(fēng)向、降雨量等數(shù)據(jù),共計262 769 組。剔除湍動分量的影響,對分鐘級的氣象要素進(jìn)行小時平均,獲得小時的氣象要素分布,共計4 395組。

    基于P-RSHMU 模型計算獲得小時蒸發(fā)波導(dǎo)高度,結(jié)果見圖1。由圖可見,大部分蒸發(fā)波導(dǎo)高度位于0~20 m,沒有出現(xiàn)蒸發(fā)波導(dǎo)高度為40 m的現(xiàn)象,波導(dǎo)高度曲線呈現(xiàn)連續(xù)變化。而張超等[16]基于改進(jìn)的NPS模型計算出的蒸發(fā)波導(dǎo)高度在5月中下旬風(fēng)向穩(wěn)定的時段,可能出現(xiàn)極其穩(wěn)定的條件,此時蒸發(fā)波導(dǎo)高度為40 m,且呈現(xiàn)跳躍變化。為對比起見,圖2 給出該模型計算結(jié)果與改進(jìn)的NPS 模型的結(jié)果對比圖,由圖可見,大部分情形下二者一致性較好,P-RSHMU 模型有效抑制了蒸發(fā)波導(dǎo)高度較高尤其預(yù)測的蒸發(fā)波導(dǎo)高度為40 m的情形。

    圖1 蒸發(fā)波導(dǎo)高度計算結(jié)果Fig.1 Calculation result of evaporation duct height

    圖2 不同模型計算結(jié)果對比Fig.2 Comparison of calculation result between different models

    2.2 相關(guān)性分析

    將計算的蒸發(fā)波導(dǎo)高度與潛熱通量、感熱通量、動量通量作散點圖,由圖3 可見,蒸發(fā)波導(dǎo)高度與潛熱通量以及動量通量呈現(xiàn)某種線性關(guān)系,而與感熱通量的線性關(guān)系不明顯。進(jìn)一步地,計算蒸發(fā)波導(dǎo)高度與潛熱通量、感熱通量、動量通量的相關(guān)系數(shù)分別為0.92、-0.12、-0.41。

    圖3 湍流通量與蒸發(fā)波導(dǎo)高度散點圖Fig.3 Scatter plot of turbulent flux and evaporation duct height

    以大氣溫度與海水溫度之差為標(biāo)準(zhǔn)區(qū)分穩(wěn)定條件,并分別給出不穩(wěn)定條件(氣海溫差<0 ℃)、中性條件(氣海溫差=0 ℃)、穩(wěn)定條件(氣海溫差>0 ℃)下的蒸發(fā)波導(dǎo)高度與湍流通量的散點圖,結(jié)果如圖4—6。不穩(wěn)定條件(見圖4)以及中性條件下(見圖5)的散點圖與張超等[16]研究結(jié)果類似,而穩(wěn)定條件下(見圖6)的散點圖顯著不同于文獻(xiàn)[16],圖7給出了穩(wěn)定條件下文獻(xiàn)[16]的散點圖進(jìn)行對比。

    圖4 不穩(wěn)定條件下湍流通量與蒸發(fā)波導(dǎo)高度散點圖Fig.4 Scatter plot of turbulent flux and evaporation duct height under unstable condition

    圖5 中性條件下湍流通量與蒸發(fā)波導(dǎo)高度散點圖Fig.5 Scatter plot of turbulent flux and evaporation duct height under neutral condition

    圖6 穩(wěn)定條件下湍流通量與蒸發(fā)波導(dǎo)高度散點圖Fig.6 Scatter plot of turbulent flux and evaporation duct height under stable condition

    圖7 穩(wěn)定條件下的計算結(jié)果(引自文獻(xiàn)[16])Fig.7 Calculation result under stable condition(cited from literature[16])

    張超等[16]指出,40 m 的蒸發(fā)波導(dǎo)高度存在于0 N/m2動量通量、-3~0 W/m2感熱通量以及0~40 W/m2潛熱通量附近。動量通量較小,對應(yīng)風(fēng)速較小,感熱通量為負(fù),對應(yīng)感熱由上面的大氣向下面的海洋傳輸能量,即氣海溫差>0 ℃,因此在穩(wěn)定條件下,較低的風(fēng)速導(dǎo)致蒸發(fā)波導(dǎo)高度的預(yù)測效果較差。這可能與該條件下湍流輸運不夠旺盛進(jìn)而導(dǎo)致蒸發(fā)波導(dǎo)預(yù)測理論MOS 中假定近海面空間為常通量不成立有關(guān)。

    文中引入的P-RSHMU 模型則有效地避免了穩(wěn)定條件下預(yù)測蒸發(fā)波導(dǎo)高度不合理的現(xiàn)象,在穩(wěn)定條件下同樣可以獲得與不穩(wěn)定條件、中性條件一致的預(yù)測結(jié)果。區(qū)分不同穩(wěn)定性條件下蒸發(fā)波導(dǎo)高度與湍流通量的相關(guān)系數(shù)見表1,其中穩(wěn)定條件計算樣本數(shù)為3 969,中性條件為16,不穩(wěn)定條件為410。由表中數(shù)據(jù)可以看到,在不區(qū)分穩(wěn)定性條件下,潛熱通量與蒸發(fā)波導(dǎo)的高度最相關(guān),動量通量次之;不穩(wěn)定條件下,兩者具有一致的相關(guān)性特征;中性條件下,感熱通量的相關(guān)性次之;穩(wěn)定條件下,動量通量的相關(guān)性次之。整體來說,潛熱通量對應(yīng)著水汽相變直接影響蒸發(fā)波導(dǎo)的生成,動量通量又或快或慢地影響海面水汽蒸發(fā)進(jìn)而間接影響蒸發(fā)波導(dǎo)的生成,感熱通量僅僅代表由于溫度差異導(dǎo)致海洋向大氣或者大氣向海洋輸送的熱量,對整體的湍流通量貢獻(xiàn)較小,因此潛熱通量、動量通量與蒸發(fā)波導(dǎo)高度的相關(guān)系數(shù)大于感熱通量。

    2.3 對比分析

    作為對比,表1還列出了文獻(xiàn)[16]的相關(guān)系數(shù)結(jié)果。由表中數(shù)據(jù)可以看到,引入的P-RSHMU 模型顯著改變了穩(wěn)定條件下即氣海溫差>0 ℃時湍流通量與P-RSHMU 蒸發(fā)波導(dǎo)高度的相關(guān)系數(shù),計算結(jié)果較為合理。

    需要特別指出的是,穩(wěn)定條件下改進(jìn)的NPS 模型在氣海溫差>0 ℃時以及動量通量較小(或者風(fēng)速較?。r會出現(xiàn)蒸發(fā)波導(dǎo)高度為40 m 或者計算波導(dǎo)高度偏高的現(xiàn)象,導(dǎo)致該條件下計算的蒸發(fā)波導(dǎo)高度與潛熱通量的相關(guān)系數(shù)為0.73,而P-RSHMU模型則抑制了這種現(xiàn)象的出現(xiàn),進(jìn)而二者的相關(guān)系數(shù)提高到0.92,且保持了與不穩(wěn)定條件以及中性條件一致的計算結(jié)果,即蒸發(fā)波導(dǎo)與潛熱通量顯著相關(guān)。

    通常情況下蒸發(fā)波導(dǎo)被認(rèn)為是由于近海面水汽蒸發(fā)從而引起大氣濕度隨高度銳減而形成的一種大氣層結(jié)。近海面水汽可認(rèn)為是飽和的,隨著高度的增加,相對濕度迅速下降進(jìn)而形成蒸發(fā)波導(dǎo)。而潛熱通量是由于水汽相變向大氣傳輸?shù)臒崃客浚蚨闹杏嬎愕恼舭l(fā)波導(dǎo)高度與潛熱通量的相關(guān)系數(shù)為0.92,有力支撐了蒸發(fā)波導(dǎo)的形成原因。

    3 結(jié)論

    本文基于海上平臺觀測的水文氣象數(shù)據(jù),利用不同穩(wěn)定性條件下預(yù)測效果都較好的P-RSHMU 模型,計算獲得不同時刻的蒸發(fā)波導(dǎo)高度,有效抑制了穩(wěn)定條件下蒸發(fā)波導(dǎo)預(yù)測高度為40 m 或者預(yù)測偏高的不合理現(xiàn)象。進(jìn)一步地,分析了該模型預(yù)測蒸發(fā)波導(dǎo)高度與動量通量、感熱通量、潛熱通量等湍流通量參數(shù)的相關(guān)性,結(jié)果表明蒸發(fā)波導(dǎo)高度與潛熱通量高度顯著相關(guān),相關(guān)系數(shù)可達(dá)0.92,與動量通量相關(guān)性次之,相關(guān)系數(shù)為-0.41,與感熱通量的相關(guān)系數(shù)最小,為-0.12。這些結(jié)果進(jìn)一步為豐富蒸發(fā)波導(dǎo)的內(nèi)涵以及形成機(jī)理提供重要數(shù)據(jù)支撐。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    国产精品一区二区三区四区免费观看 | 夜夜看夜夜爽夜夜摸| 午夜免费成人在线视频| 十八禁网站免费在线| 亚洲av熟女| 午夜老司机福利片| 看黄色毛片网站| 亚洲电影在线观看av| 亚洲美女视频黄频| 在线观看午夜福利视频| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 2021天堂中文幕一二区在线观| 亚洲中文av在线| 真人一进一出gif抽搐免费| www.熟女人妻精品国产| 国语自产精品视频在线第100页| 欧美av亚洲av综合av国产av| 亚洲av第一区精品v没综合| 97超级碰碰碰精品色视频在线观看| 亚洲 欧美一区二区三区| 欧美三级亚洲精品| 免费高清视频大片| 国产欧美日韩一区二区精品| 丁香六月欧美| 久久久久亚洲av毛片大全| 亚洲精品国产精品久久久不卡| 青草久久国产| 国产精品电影一区二区三区| www.自偷自拍.com| 露出奶头的视频| 国产精品美女特级片免费视频播放器 | 夜夜躁狠狠躁天天躁| 国产亚洲精品久久久久久毛片| av在线播放免费不卡| 日日干狠狠操夜夜爽| 首页视频小说图片口味搜索| 欧美激情久久久久久爽电影| 久久精品国产亚洲av香蕉五月| av视频在线观看入口| 老司机午夜福利在线观看视频| 好看av亚洲va欧美ⅴa在| 又粗又爽又猛毛片免费看| 亚洲国产欧洲综合997久久,| 操出白浆在线播放| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| av在线天堂中文字幕| 一进一出抽搐动态| 在线国产一区二区在线| 最好的美女福利视频网| 成人午夜高清在线视频| 熟女少妇亚洲综合色aaa.| 日本熟妇午夜| 在线永久观看黄色视频| 一本一本综合久久| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 久久久精品国产亚洲av高清涩受| 又黄又粗又硬又大视频| 亚洲精品久久成人aⅴ小说| 69av精品久久久久久| 国产97色在线日韩免费| 天天一区二区日本电影三级| 男女做爰动态图高潮gif福利片| 久久久久久久精品吃奶| 91大片在线观看| 国产精品亚洲一级av第二区| 五月玫瑰六月丁香| 成人av在线播放网站| 国产av麻豆久久久久久久| 欧美一区二区精品小视频在线| 床上黄色一级片| 在线观看免费午夜福利视频| 国产成+人综合+亚洲专区| 在线视频色国产色| 欧美zozozo另类| 精品久久久久久成人av| 精品国产亚洲在线| 91成年电影在线观看| 很黄的视频免费| 黑人操中国人逼视频| 人人妻人人澡欧美一区二区| 搡老妇女老女人老熟妇| 国产精品永久免费网站| 动漫黄色视频在线观看| 久久香蕉国产精品| 国产视频内射| 国产乱人伦免费视频| 久久99热这里只有精品18| 日韩高清综合在线| 看片在线看免费视频| www.自偷自拍.com| 国产精品久久久久久人妻精品电影| 国产精品久久电影中文字幕| 老司机靠b影院| 亚洲中文av在线| 欧美色欧美亚洲另类二区| 久久精品亚洲精品国产色婷小说| 日韩欧美一区二区三区在线观看| 午夜精品久久久久久毛片777| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看 | 免费在线观看视频国产中文字幕亚洲| 91麻豆精品激情在线观看国产| 桃红色精品国产亚洲av| 91在线观看av| 亚洲片人在线观看| 麻豆成人午夜福利视频| 国产亚洲精品一区二区www| 村上凉子中文字幕在线| 国产精品久久视频播放| 精品高清国产在线一区| 美女免费视频网站| 制服诱惑二区| 又黄又粗又硬又大视频| 久久精品夜夜夜夜夜久久蜜豆 | www日本在线高清视频| 男女下面进入的视频免费午夜| 午夜福利在线观看吧| 黄色a级毛片大全视频| 久久久久久人人人人人| 国产成人精品无人区| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| 久久精品aⅴ一区二区三区四区| 嫁个100分男人电影在线观看| 看免费av毛片| 美女午夜性视频免费| 18禁裸乳无遮挡免费网站照片| 高潮久久久久久久久久久不卡| 97碰自拍视频| 在线观看午夜福利视频| 中文资源天堂在线| 国产私拍福利视频在线观看| 国产视频内射| 日本黄色视频三级网站网址| 免费看美女性在线毛片视频| 久久热在线av| 欧美av亚洲av综合av国产av| 好男人在线观看高清免费视频| 男男h啪啪无遮挡| 日韩成人在线观看一区二区三区| 成人手机av| 色综合婷婷激情| 一级a爱片免费观看的视频| a级毛片在线看网站| 国产成人av激情在线播放| 神马国产精品三级电影在线观看 | 亚洲精品中文字幕在线视频| av视频在线观看入口| 久久午夜综合久久蜜桃| 国产野战对白在线观看| 久久性视频一级片| 99国产精品一区二区三区| 国产精品久久久久久精品电影| 99在线视频只有这里精品首页| 久久久久免费精品人妻一区二区| 欧美黑人欧美精品刺激| 欧美极品一区二区三区四区| 久久婷婷成人综合色麻豆| 欧美日韩亚洲综合一区二区三区_| 亚洲avbb在线观看| 18禁观看日本| 高清毛片免费观看视频网站| 搞女人的毛片| 国产精品九九99| 精品久久蜜臀av无| 久久久久久人人人人人| 中国美女看黄片| 天堂av国产一区二区熟女人妻 | 三级毛片av免费| 在线十欧美十亚洲十日本专区| 香蕉av资源在线| 18禁黄网站禁片午夜丰满| 天天一区二区日本电影三级| 国产一区二区激情短视频| 亚洲第一欧美日韩一区二区三区| 999久久久国产精品视频| 两个人视频免费观看高清| 欧美久久黑人一区二区| 五月伊人婷婷丁香| bbb黄色大片| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线观看免费| 狠狠狠狠99中文字幕| 九色国产91popny在线| 国产男靠女视频免费网站| 久久精品91无色码中文字幕| 国产成+人综合+亚洲专区| 宅男免费午夜| 欧美另类亚洲清纯唯美| 久久久久亚洲av毛片大全| 午夜福利在线观看吧| 精品高清国产在线一区| 亚洲精品色激情综合| 亚洲成a人片在线一区二区| 色老头精品视频在线观看| 国产乱人伦免费视频| 日韩 欧美 亚洲 中文字幕| 99在线人妻在线中文字幕| 午夜精品久久久久久毛片777| 99热这里只有是精品50| 精品一区二区三区视频在线观看免费| 亚洲精品国产精品久久久不卡| 美女 人体艺术 gogo| 一级毛片高清免费大全| 在线视频色国产色| 亚洲国产欧美一区二区综合| 观看免费一级毛片| 一区二区三区激情视频| 他把我摸到了高潮在线观看| 大型av网站在线播放| 欧美一级毛片孕妇| 欧美高清成人免费视频www| 久久久久国产精品人妻aⅴ院| 欧美又色又爽又黄视频| 亚洲国产欧洲综合997久久,| 丝袜人妻中文字幕| 欧美另类亚洲清纯唯美| 欧美日韩瑟瑟在线播放| 亚洲真实伦在线观看| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 免费一级毛片在线播放高清视频| 欧美+亚洲+日韩+国产| 国产一区二区三区在线臀色熟女| 狂野欧美白嫩少妇大欣赏| 50天的宝宝边吃奶边哭怎么回事| 777久久人妻少妇嫩草av网站| 日本 av在线| 国产精品美女特级片免费视频播放器 | 波多野结衣巨乳人妻| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 正在播放国产对白刺激| 一a级毛片在线观看| 国产成人aa在线观看| 亚洲欧美精品综合一区二区三区| 制服丝袜大香蕉在线| 黄色a级毛片大全视频| 国产又色又爽无遮挡免费看| 欧美3d第一页| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 手机成人av网站| 九九热线精品视视频播放| 欧美大码av| 美女免费视频网站| 夜夜夜夜夜久久久久| 色哟哟哟哟哟哟| 国产精品久久久久久人妻精品电影| 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 日本三级黄在线观看| av中文乱码字幕在线| 亚洲成av人片免费观看| 黑人巨大精品欧美一区二区mp4| 免费看美女性在线毛片视频| 亚洲国产欧美网| 深夜精品福利| 国产人伦9x9x在线观看| 日本 av在线| 国产97色在线日韩免费| 日本免费一区二区三区高清不卡| 两性夫妻黄色片| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| 香蕉久久夜色| 国产一区二区三区在线臀色熟女| www.熟女人妻精品国产| 99国产精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 老熟妇乱子伦视频在线观看| 国产免费av片在线观看野外av| 亚洲色图 男人天堂 中文字幕| 欧美成人免费av一区二区三区| 黄色成人免费大全| 男女床上黄色一级片免费看| 久久久久国内视频| 亚洲avbb在线观看| 777久久人妻少妇嫩草av网站| 亚洲国产精品999在线| 亚洲国产欧美一区二区综合| 校园春色视频在线观看| 女同久久另类99精品国产91| 久久精品国产清高在天天线| 亚洲成人久久爱视频| 搞女人的毛片| 国产精品久久久人人做人人爽| 91国产中文字幕| 欧美高清成人免费视频www| 亚洲美女视频黄频| 国产真实乱freesex| 亚洲国产欧美网| 国产亚洲精品第一综合不卡| 免费观看人在逋| 久9热在线精品视频| 国产1区2区3区精品| 日本免费一区二区三区高清不卡| 久久精品综合一区二区三区| 无人区码免费观看不卡| 亚洲午夜理论影院| 人人妻人人澡欧美一区二区| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 岛国在线观看网站| 国产欧美日韩一区二区三| 给我免费播放毛片高清在线观看| 国产99久久九九免费精品| 亚洲 国产 在线| 国产三级黄色录像| 日本撒尿小便嘘嘘汇集6| 婷婷亚洲欧美| 久久久久久久久久黄片| 91大片在线观看| 国产精品av视频在线免费观看| 久久久国产成人精品二区| 男人舔奶头视频| 最近在线观看免费完整版| 国产亚洲精品久久久久久毛片| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产成人免费| 中文亚洲av片在线观看爽| 久久精品91无色码中文字幕| 精品少妇一区二区三区视频日本电影| 精品久久久久久久末码| 在线观看免费视频日本深夜| 十八禁人妻一区二区| 国产精品1区2区在线观看.| 国产成人精品无人区| 又黄又粗又硬又大视频| 精品免费久久久久久久清纯| www日本在线高清视频| 91av网站免费观看| 日本熟妇午夜| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| 在线观看美女被高潮喷水网站 | 床上黄色一级片| 夜夜躁狠狠躁天天躁| 99热只有精品国产| 91老司机精品| 长腿黑丝高跟| 丁香六月欧美| 亚洲人成77777在线视频| 国产1区2区3区精品| 日日夜夜操网爽| 婷婷六月久久综合丁香| 久久伊人香网站| 亚洲片人在线观看| 高清在线国产一区| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 久久久久久久久久黄片| 国产成人av激情在线播放| 99热只有精品国产| 日韩大尺度精品在线看网址| 黄色视频不卡| 久久久久久久久中文| videosex国产| 91九色精品人成在线观看| 国产精品自产拍在线观看55亚洲| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 欧美精品啪啪一区二区三区| 一级片免费观看大全| 欧美精品亚洲一区二区| 又大又爽又粗| 国产一区二区三区视频了| 午夜精品在线福利| 男人舔女人的私密视频| 国产日本99.免费观看| av天堂在线播放| 在线看三级毛片| 久久久久亚洲av毛片大全| 99久久精品国产亚洲精品| 亚洲精品一区av在线观看| 国产视频内射| 日韩 欧美 亚洲 中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| √禁漫天堂资源中文www| 久久99热这里只有精品18| 国产人伦9x9x在线观看| www.999成人在线观看| 亚洲 欧美一区二区三区| 久久久久久国产a免费观看| 久久久久久久久中文| 亚洲狠狠婷婷综合久久图片| 免费在线观看视频国产中文字幕亚洲| 免费人成视频x8x8入口观看| 一级黄色大片毛片| 亚洲欧美日韩高清专用| 久久欧美精品欧美久久欧美| 久久草成人影院| 美女午夜性视频免费| 亚洲国产中文字幕在线视频| 国产精品一及| 日韩大尺度精品在线看网址| 亚洲av电影在线进入| 亚洲一区二区三区不卡视频| 亚洲一码二码三码区别大吗| 精品不卡国产一区二区三区| 久9热在线精品视频| 亚洲全国av大片| 精品无人区乱码1区二区| 村上凉子中文字幕在线| 法律面前人人平等表现在哪些方面| 免费av毛片视频| 久久久久久久久免费视频了| av有码第一页| 国产免费男女视频| 久久亚洲真实| 精品国产超薄肉色丝袜足j| 男插女下体视频免费在线播放| 亚洲精品中文字幕在线视频| 免费在线观看成人毛片| 亚洲国产精品sss在线观看| 99热这里只有精品一区 | 女生性感内裤真人,穿戴方法视频| 精品久久久久久久人妻蜜臀av| АⅤ资源中文在线天堂| 窝窝影院91人妻| 舔av片在线| 欧美 亚洲 国产 日韩一| 国产成人精品久久二区二区91| 黄色成人免费大全| 99在线视频只有这里精品首页| 97碰自拍视频| 午夜福利成人在线免费观看| 国产91精品成人一区二区三区| 午夜精品一区二区三区免费看| 欧美日韩中文字幕国产精品一区二区三区| 国产免费男女视频| 久久久久国内视频| 日韩免费av在线播放| 国产区一区二久久| 日本 av在线| 男女那种视频在线观看| 午夜福利免费观看在线| 岛国在线观看网站| 国产亚洲av嫩草精品影院| 亚洲全国av大片| 又爽又黄无遮挡网站| 国产高清激情床上av| 久久人妻av系列| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 欧美成人性av电影在线观看| 成人特级黄色片久久久久久久| 日韩有码中文字幕| 美女免费视频网站| 九色国产91popny在线| 国产精品日韩av在线免费观看| 又黄又粗又硬又大视频| 日韩大码丰满熟妇| 男人舔女人的私密视频| 色综合欧美亚洲国产小说| 国产免费av片在线观看野外av| 午夜精品一区二区三区免费看| a级毛片a级免费在线| 亚洲最大成人中文| 久久性视频一级片| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 欧美黄色片欧美黄色片| 中文资源天堂在线| 18美女黄网站色大片免费观看| 我的老师免费观看完整版| 桃色一区二区三区在线观看| 国产精品久久久人人做人人爽| 国产99久久九九免费精品| 叶爱在线成人免费视频播放| 在线a可以看的网站| 亚洲欧美日韩高清专用| 久久国产乱子伦精品免费另类| 色综合亚洲欧美另类图片| 观看免费一级毛片| 久久久精品欧美日韩精品| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av高清一级| 日本成人三级电影网站| 99在线视频只有这里精品首页| 免费看a级黄色片| 午夜成年电影在线免费观看| 亚洲精品国产精品久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看 | netflix在线观看网站| 亚洲国产精品成人综合色| 三级国产精品欧美在线观看 | 一本大道久久a久久精品| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站 | 免费在线观看完整版高清| 午夜日韩欧美国产| 99久久无色码亚洲精品果冻| 老熟妇仑乱视频hdxx| 成年女人毛片免费观看观看9| 免费电影在线观看免费观看| 亚洲av成人一区二区三| 午夜福利在线观看吧| 大型黄色视频在线免费观看| 亚洲男人的天堂狠狠| 哪里可以看免费的av片| 岛国在线观看网站| 欧美色视频一区免费| 亚洲国产精品sss在线观看| 最近最新中文字幕大全电影3| 欧美一区二区国产精品久久精品 | 99国产综合亚洲精品| av福利片在线观看| 非洲黑人性xxxx精品又粗又长| 国产主播在线观看一区二区| 国产1区2区3区精品| 亚洲国产精品成人综合色| 日韩欧美 国产精品| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 无限看片的www在线观看| 桃红色精品国产亚洲av| 色在线成人网| 亚洲av成人精品一区久久| 丝袜美腿诱惑在线| 在线观看www视频免费| 桃色一区二区三区在线观看| 国产精品,欧美在线| 9191精品国产免费久久| 中国美女看黄片| 久久欧美精品欧美久久欧美| 精品久久久久久久久久久久久| 国产精品爽爽va在线观看网站| 亚洲专区国产一区二区| 久久久久性生活片| 国产精品免费一区二区三区在线| 欧美乱码精品一区二区三区| 国产精华一区二区三区| 日韩成人在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 狠狠狠狠99中文字幕| 悠悠久久av| 老司机福利观看| 国产欧美日韩精品亚洲av| 少妇熟女aⅴ在线视频| 麻豆国产av国片精品| 欧美在线黄色| 窝窝影院91人妻| 老鸭窝网址在线观看| 他把我摸到了高潮在线观看| 亚洲av成人一区二区三| 嫁个100分男人电影在线观看| 久久久国产欧美日韩av| 人妻丰满熟妇av一区二区三区| av视频在线观看入口| 国内少妇人妻偷人精品xxx网站 | 一边摸一边做爽爽视频免费| 久久久久久亚洲精品国产蜜桃av| 五月玫瑰六月丁香| 国内精品一区二区在线观看| 亚洲人成电影免费在线| 99久久精品热视频| 欧美日本视频| 成年人黄色毛片网站| 美女免费视频网站| 国产精品野战在线观看| 免费看十八禁软件| 亚洲美女黄片视频| 中文字幕久久专区| 久久亚洲真实| av中文乱码字幕在线| 亚洲av成人精品一区久久| av视频在线观看入口| 亚洲 国产 在线| 一本精品99久久精品77| 一夜夜www| 欧美一区二区精品小视频在线| 国产精品 欧美亚洲| 中文字幕熟女人妻在线| 欧美三级亚洲精品| 亚洲av电影不卡..在线观看| 淫秽高清视频在线观看| 一区二区三区激情视频| 日本撒尿小便嘘嘘汇集6| 在线观看一区二区三区| 精品一区二区三区av网在线观看| 精品久久久久久久毛片微露脸| 精品午夜福利视频在线观看一区| 9191精品国产免费久久| 亚洲av电影在线进入| 欧美激情久久久久久爽电影| 久久亚洲真实| 人妻夜夜爽99麻豆av| av片东京热男人的天堂| 久久精品国产亚洲av香蕉五月| 中文在线观看免费www的网站 | 亚洲成a人片在线一区二区| 男人舔女人下体高潮全视频| 久久人妻av系列| 波多野结衣高清作品| 看片在线看免费视频| a级毛片a级免费在线| 国产精品乱码一区二三区的特点| 久久精品91无色码中文字幕| 一区二区三区高清视频在线| 桃色一区二区三区在线观看| 久久亚洲精品不卡| 男女那种视频在线观看| 视频区欧美日本亚洲| www日本在线高清视频| 给我免费播放毛片高清在线观看|