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

    基于高頻地波雷達(dá)觀測的大亞灣附近海域余流季節(jié)變化特征分析

    2021-05-19 08:53:36袁凱瑞王久娟林佳瑜王凱敏
    關(guān)鍵詞:余流大亞灣風(fēng)場

    袁凱瑞,王久娟,林佳瑜,韓 宇,王凱敏,

    徐宇程3,陳勁松1,2*

    (1.中國科學(xué)院深圳先進(jìn)技術(shù)研究院空間信息研究中心,廣東 深圳 518055;2.深圳海洋環(huán)境大數(shù)據(jù)應(yīng)用與分析工程實驗室,廣東 深圳 518055; 3.深圳市海洋監(jiān)測預(yù)報中心,廣東 深圳 518055)

    大亞灣位于大鵬半島的東部,是南海北部向陸地延伸最深的海灣。大亞灣地理區(qū)位優(yōu)勢突出,交通便利快捷,自然條件優(yōu)越。王友紹等(2004)認(rèn)為近年來大亞灣內(nèi)海岸線的改變、濕地不斷減少以及大量陸源污染物的輸入等,特別是石化區(qū)及核電站的建成,使得大亞灣海域的生態(tài)環(huán)境發(fā)生了較大的變化,海洋生態(tài)環(huán)境惡化趨勢明顯[1]。大亞灣海域水文動力的研究可為該海域的赤潮擴散、水交換能力估算、物質(zhì)運輸、環(huán)境保護(hù)等研究提供重要的數(shù)據(jù)和理論支撐。

    目前對大亞灣水動力的研究主要集中在潮流與余流的特征方面。吳巖等(1998)通過潮流數(shù)值模擬發(fā)現(xiàn),在大亞灣內(nèi)存在著兩個順時針方向的環(huán)流場,且南部環(huán)流強度大于北部環(huán)流[2]。楊國標(biāo)(2001)研究表明,大亞灣的潮流受天文潮和外海潮波的共同作用,灣內(nèi)主要為水平往復(fù)流,灣口及西部的余流較強,灣口余流流向灣外,其余基本沿岸線流動[3]。吳仁豪等(2007)利用HAMSOM模式對大亞灣的潮汐、潮流和余流進(jìn)行了數(shù)值模擬,研究發(fā)現(xiàn)灣內(nèi)海域主要為南北向的往復(fù)流,流速受地形影響強烈,且灣口西側(cè)流速大于東側(cè),夏季灣內(nèi)余流較小,灣口余流大體流向灣外,冬季余流較大,灣口余流大體流向灣內(nèi)[4]。王聰?shù)?2008)采用ECOM-si三維海洋數(shù)值模式模擬了大亞灣海流,認(rèn)為大亞灣潮汐和潮流特性主要受來自西太平洋的潮波制約即受廣東沿岸流控制,在灣口附近海域,余流從大辣甲島東部進(jìn)入,西部流出[5]。Song等(2016) 、武文等(2017)、嚴(yán)聿晗等(2017)借助FVCOM海洋模式建立高時空分辨率的大亞灣三維潮汐潮流數(shù)值模型,分析了大亞灣內(nèi)潮位和潮流“雙峰”現(xiàn)象的原因,發(fā)現(xiàn)了灣內(nèi)不一致的大小潮變化以及灣外所受沿岸流的影響,揭示了大亞灣灣內(nèi)的余流多渦旋結(jié)構(gòu)[6-8]。以上研究多基于現(xiàn)場觀測資料或數(shù)值模擬結(jié)合現(xiàn)場觀測對大亞灣的潮流與余流進(jìn)行分析。

    本研究在大亞灣附近海域基于高頻地波雷達(dá)的長期海流觀測資料,利用Pawlowicz等(2002)發(fā)布的T_tide工具包[9]進(jìn)行潮流調(diào)和分析,實現(xiàn)潮流和余流的分離,進(jìn)而研究該海域余流的季節(jié)分布特征及變化規(guī)律,以期能夠為該海域的海上人員搜救、溢油擴散、污染物漂移等研究工作提供參考。

    1 數(shù)據(jù)

    1.1 數(shù)據(jù)簡介

    1.1.1 高頻地波雷達(dá)數(shù)據(jù) 高頻地波雷達(dá)是利用雷達(dá)發(fā)射的電磁波與海浪作用時產(chǎn)生的Bragg散射和多普勒效應(yīng),探測海表層流場、浪場和風(fēng)場的空間分布,同其他的海洋監(jiān)測設(shè)備相比,它具有監(jiān)測范圍大、探測精度高、受環(huán)境的影響較小、建設(shè)經(jīng)費少以及能夠進(jìn)行連續(xù)無人值守的實時在線監(jiān)測等優(yōu)點,目前已成為了一種重要的海洋觀測手段。OSMAR-S系列便攜式高頻地波雷達(dá)是武漢大學(xué)電波傳播實驗室研制,其長周期工作的穩(wěn)定性和探測精度已被文必洋等(2009)、魏國妹等(2016)與Lai等(2017)的諸多研究認(rèn)可[10-12]。

    該高頻雷達(dá)觀測系統(tǒng)是由兩臺武漢德威斯公司生產(chǎn)的OSMAR-S50(主要技術(shù)指標(biāo)如表1 所示)組成,分別位于深圳大鵬半島的鹿嘴山莊和惠東平海鎮(zhèn)的海龜島(圖1)。本研究采用的表層海流數(shù)據(jù)的起止時間為2015年1月1日至2015 年12月31日,數(shù)據(jù)由深圳市海洋監(jiān)測預(yù)報中心提供。該雷達(dá)系統(tǒng)的波束覆蓋扇角為180°,即每臺雷達(dá)的探測覆蓋范圍為法向方向±90°的半圓區(qū)域。矢量海流為兩臺雷達(dá)的單站探測徑向流合成,空間分辨率為0. 01°×0. 01°,數(shù)據(jù)采樣間隔為20 min。圖1給出了2015年該海域雷達(dá)系統(tǒng)數(shù)據(jù)的空間分布,在兩個雷達(dá)連線附近公共覆蓋區(qū)域的空白區(qū)數(shù)據(jù)系手動剔除,高精度區(qū)主要分布在灣外海域,因此本研究主要關(guān)注灣外海域。

    表1 便攜式高頻地波雷達(dá)(OSMAR-S50)的主要技術(shù)指標(biāo)Tab.1 Main technical indexes of the portable high frequency ground wave radar (OSMAR-S50)

    圖1 大亞灣附近海域雷達(dá)站位及數(shù)據(jù)有效率分布Fig.1 Distribution of the radar station and data availability in adjacent waters of Daya Bay

    1.1.2 風(fēng)場數(shù)據(jù) 風(fēng)場數(shù)據(jù)為美國國家環(huán)境預(yù)測中心和國家大氣研究中心的ERDDAP服務(wù)器中提供的All Metop ASCAT中的月平均數(shù)據(jù)(https://coastwatch.pfeg.noaa.gov/erddap/griddap/erdQMwindmday.html), 空間分辨率為0.25°×0.25°。為了更直觀地分析觀測海域風(fēng)速風(fēng)向的整體變化過程,本研究還下載了單點A(22.25°N,114.75°E)的All Metop ASCAT的日平均數(shù)據(jù),主要用于單點風(fēng)場時間序列分析。

    1.1.3 地形數(shù)據(jù) 地形數(shù)據(jù)為Etopo1的水深數(shù)據(jù),該數(shù)據(jù)是美國的國家海洋和大氣局(NOAA)的國家地球物理數(shù)據(jù)中心(NGDC)公布的數(shù)據(jù)(http://www.ngdc.noaa.gov/mgg/global/),分辨率為1′。

    1.2 數(shù)據(jù)質(zhì)量控制

    高頻地波雷達(dá)在長時間的持續(xù)工作中,不僅會受到天氣、海況的影響,同時雷達(dá)系統(tǒng)本身穩(wěn)定性也會對數(shù)據(jù)質(zhì)量有影響,導(dǎo)致數(shù)據(jù)缺測、數(shù)據(jù)覆蓋率不穩(wěn)定、虛假流速等問題的存在。因此,在對雷達(dá)數(shù)據(jù)進(jìn)行分析前需要對雷達(dá)數(shù)據(jù)進(jìn)行嚴(yán)格的數(shù)據(jù)質(zhì)量控制,這樣才能得到較高質(zhì)量的雷達(dá)數(shù)據(jù),為后期的數(shù)據(jù)應(yīng)用提供可靠保障。本研究依據(jù)朱大勇等(2008)和徐粲等(2014)的經(jīng)驗[13-14]對雷達(dá)數(shù)據(jù)做如下的質(zhì)量控制:①剔除異常值,即剔除觀測中東分量和北分量的流速超過2 m/s的數(shù)據(jù)和由于數(shù)據(jù)傳輸不穩(wěn)定導(dǎo)致的連續(xù)幾組完全相同的明顯異常數(shù)據(jù);②剔除跳躍值,即與平均值之差絕對值大于2倍標(biāo)準(zhǔn)差的數(shù)據(jù)和前后兩次觀測絕對變化值大于2倍標(biāo)準(zhǔn)差的數(shù)據(jù);③剔除孤立值,即剔除前后均無有效采樣點的數(shù)據(jù);④利用時間與空間插值將缺失時刻和缺失的雷達(dá)網(wǎng)格點數(shù)據(jù)補充完整,并進(jìn)行時間序列平滑與空間平滑,最后獲得時空連續(xù)的流場序列。通過以上質(zhì)控手段,認(rèn)為最后獲得時空連續(xù)的雷達(dá)數(shù)據(jù)是可靠的。

    2 結(jié)果與討論

    余流是實際海流濾掉潮流后剩下的部分,包括潮余流、風(fēng)海流和密度流等,余流對海水中物質(zhì)的輸運有重要作用,余流分析是海域泥沙輸運,污染物、赤潮等擴散的重要研究手段。本研究將時空連續(xù)的雷達(dá)表層流數(shù)據(jù)減去調(diào)和分析得到的潮流即得到了余流的時空序列,并對各月余流序列進(jìn)行月平均,得到月均余流場。

    2.1 1—3月余流特征

    如圖2~4所示,1—3月,大亞灣海域整體風(fēng)速為6~8 m/s的東北風(fēng),風(fēng)向逐漸由NE向轉(zhuǎn)向ENE方向(圖2、3)。大亞灣灣口余流總體表現(xiàn)為西進(jìn),且西部流速大于東部流速;灣外東南部外海有一股約28~33 cm/s的強余流,流向為SW向,基本沿著等深線方向流動,流速梯度與水深梯度近乎一致,其中2月余流強度最強,3月開始減弱,這股流在越過114.80°E后隨著等深線的轉(zhuǎn)向也逐漸向陸地側(cè)偏轉(zhuǎn),同時隨著水深的減少,流速也逐漸減弱,流速方向變化趨勢與等深線相同,即余流的整體梯度變化與水深變化一致(圖4)。

    2.2 4—6月余流特征

    如圖5~7所示,4月,季風(fēng)逐漸開始轉(zhuǎn)向,大亞灣海域整體為約5.0 m/s的東風(fēng)(圖5、6),灣外強余流消失,整體余流流速減小到一年之中的最小值(小于5 cm/s),流向指向灣內(nèi),灣口外東側(cè)余流近乎消失;至5月,海域內(nèi)的風(fēng)轉(zhuǎn)為3.5 m/s的南風(fēng),灣口外的強余流逐漸加強,灣口外的流向完全反轉(zhuǎn)為NE向,可能是由于風(fēng)的轉(zhuǎn)向?qū)е?;到?月,海域內(nèi)的風(fēng)轉(zhuǎn)為約3.8 m/s的西南風(fēng),灣外的強余流加強達(dá)到20 cm/s以上,流向為NE向,等流速線與等深線幾乎平行,流速梯度與水深梯度近乎一致,且此時風(fēng)向與流向夾角大約為45°(圖7),與Ekman流的特點一致。此外,在4月流速變化期,余流整體較小,灣外的強余流近乎消失,流速與雷達(dá)本身的測流誤差相當(dāng)。因此,4月的流向、流速可能存在較大誤差,如A點的風(fēng)速矢量顯示4月為風(fēng)速轉(zhuǎn)向期,特別是4月下旬風(fēng)速呈不穩(wěn)定的轉(zhuǎn)向,整體風(fēng)向大體為E向,但余流流向為SW向,與風(fēng)向不一致。

    2.3 7—9月余流特征

    如圖8~10所示,7—8月,海域整體為1.5~2.6 m/s的南風(fēng)和西南風(fēng)(圖8、9),余流分布與6月一致,流向為NW向,8月流速達(dá)到最大,灣外的強余流最大流速可達(dá)23 cm/s,此時灣口外的等流速線與等深線近乎平行;9月海域盛行東風(fēng),灣口東側(cè)有一股約10 cm/s的SW向流,且在越過114.80°E后轉(zhuǎn)向灣內(nèi)(圖10)。同8月相比,余流流向再一次發(fā)生翻轉(zhuǎn),且僅在一個月之內(nèi)完成,比4、5月的翻轉(zhuǎn)時間縮短近一個月,說明該海域的海水對東北季風(fēng)發(fā)展的響應(yīng)迅捷。

    2.4 10—12月余流特征

    如圖11~13所示,10—12月該海域又再次變?yōu)闁|北風(fēng),風(fēng)速為6.5~8.8 m/s(圖11、12),余流整體變化趨勢與1—3月變化一致,流向為SW向,流速可達(dá)30 cm/s以上,且灣外強余流逐漸增強,流速梯度逐漸變大,等流速線逐漸與等深線平行,流速梯度趨向與水深梯度一致,灣外強余流對灣口余流的影響趨于穩(wěn)定(圖13)。

    圖2 2015年1—3月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 2 ASCAT remote sensing wind field in adjacent waters of Daya Bay from January to March, 2015圖中“”是A點(22.25°N,114.75°E)的位置,為All Metop ASCAT日平均數(shù)據(jù)點;箭頭指向表示風(fēng)向;下同。

    圖3 2015年1—3月大亞灣附近海域A點日均風(fēng)矢量變化Fig. 3 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from January to March, 2015

    圖4 2015年1—3月大亞灣附近海域地波雷達(dá)表層余流場Fig. 4 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from January to March, 2015圖中藍(lán)色等值線為水深,白色箭頭指向為余流方向,下同。

    圖5 2015年4—6月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 5 ASCAT remote sensing wind field in adjacent waters of Daya Bay from April to June, 2015

    圖6 2015年 4—6月大亞灣附近海域A點日均風(fēng)矢量變化 Fig. 6 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from April to June, 2015

    圖7 2015年4—6月大亞灣附近海域地波雷達(dá)表層余流場Fig. 7 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from April to June, 2015

    圖8 2015年7—9月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 8 ASCAT remote sensing wind field in adjacent waters of Daya Bay from July to September, 2015

    圖9 2015年7—9月大亞灣附近海域A點日均風(fēng)矢量變化Fig. 9 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from July to September, 2015

    圖10 2015年7—9月大亞灣附近海域地波雷達(dá)表層余流場Fig. 10 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from July to September, 2015

    2.5 討論

    2.5.1 大亞灣附近海域余流季節(jié)變化特征 本研究通過對2015年全年的高頻地波雷達(dá)連續(xù)觀測數(shù)據(jù)分析發(fā)現(xiàn),大亞灣海域表層余流的變化主要呈冬、夏兩種形態(tài),其中,冬季形態(tài)持續(xù)半年以上(從10月到次年3 月),流向為SW向,流速可達(dá)30 cm/s以上;夏季形態(tài)僅在5—8 月間出現(xiàn),流向為NE向,流速可達(dá)23 cm/s。本研究的余流變化特征與朱大勇等[13]的觀測結(jié)論類似,而且這種變化過程與季風(fēng)的轉(zhuǎn)變基本同步。

    2.5.2 大亞灣附近海域余流驅(qū)動機制討論 本研究的余流空間分布數(shù)據(jù)顯示:10月到次年3月,風(fēng)向與強余流的方向基本一致,此時強余流特征與粵東沿岸流基本一致,且余流總在越過114.80°E后逐漸轉(zhuǎn)為向岸流,這與水深等值線在這里發(fā)生向岸轉(zhuǎn)向有關(guān);6—8月,風(fēng)速較小,但風(fēng)向與流向的夾角小于90°,體現(xiàn)了風(fēng)對余流的影響,同時也發(fā)現(xiàn)風(fēng)只能影響表層的流速大小而不能改變其主要流向,這與陳希榮等(2018)在珠江口海域的觀測結(jié)論[15]相同。此外,除4月與9月,余流等值線與水深等值線近乎平行,流速梯度均為NW—SE向,這與研究區(qū)域地形梯度的分布一致,又體現(xiàn)了地形對余流的影響。綜合全年的余流變化特征,不難看出大亞灣附近海域的余流變化主要受東亞季風(fēng)影響,其變化特征與每年10月至次年3月盛行的東北季風(fēng),6—8月盛行的西南季風(fēng)相一致。此外,地形對余流分布也存在重要影響。這與其他學(xué)者的研究相一致,如徐粲等在南黃海輻射沙脊群的研究發(fā)現(xiàn),近岸海域余流分布主要受潮汐與水下地形共同作用的影響,而靠外海域余流格局主要受風(fēng)場的控制[14];史軍強(2014)在長江口海域的研究也同樣證明了風(fēng)場和地形對余流分布的影響[16]。

    此外,粵東沿岸流與珠江沖淡水對大亞灣灣口余流空間分布具有調(diào)控作用。孫湘平(2006)、張志欣(2014)和舒業(yè)強等(2018)的研究結(jié)果認(rèn)為廣東沿岸存在著一股沿岸流即粵東沿岸流,它主要是隨季風(fēng)風(fēng)向變化,冬季為SW向流,流幅窄;夏季為NE向流,流幅寬,且珠江沖淡水在夏季可影響至大亞灣外海域[17-19]。10月至次年3月研究海域主要受到東北風(fēng)影響下的SW向強沿岸流影響,粵東沿岸流在越過114.80°E后,由于受地形的影響,部分流發(fā)生轉(zhuǎn)向作為補償流逐漸轉(zhuǎn)向成為垂直岸線的余流;4月東北風(fēng)逐漸減弱,到5月初已表現(xiàn)為較弱的西南風(fēng),但整體上表現(xiàn)為東北風(fēng),由于東北季風(fēng)的減弱,Ekman效應(yīng)制約沖淡水向外海擴散的能力明顯減弱,沖淡水漸漸影響大亞灣灣口以外的海域; 6—8月隨著珠江徑流量的增加與西南季風(fēng)的建立,在西南季風(fēng)的驅(qū)使及粵東沿岸流的誘導(dǎo)下,表層沖淡水向東擴散至大亞灣附近海域,這與前人的研究結(jié)果[20]是一致的。歐素英(2005)、曾淦寧等(2005)、Ou等(2007)、楊陽等(2014)與Bai等(2015)的研究表明,本研究區(qū)域,大亞灣附近海域的余流受到了珠江沖淡水和粵東沿岸流的共同作用[20-24]。

    2.5.3 高頻地波雷達(dá)在沿岸流觀測中的作用 對于粵東沿岸流的觀測大多基于一次或多次的航次觀測和定點觀測數(shù)據(jù),盡管衛(wèi)星遙感觀測和數(shù)值模擬方法很大程度上豐富了我們對粵東沿岸流的觀測手段,但缺少大范圍連續(xù)觀測數(shù)據(jù),特別是長時間序列的連續(xù)海流觀測數(shù)據(jù),導(dǎo)致我們無法進(jìn)行準(zhǔn)確直觀的觀測。本研究利用的長周期高頻地波雷達(dá)觀測分析余流變化特征與Fang等(1998)對粵東海域沿岸流的研究結(jié)果相似,即沿岸流主要受季節(jié)性反轉(zhuǎn)的風(fēng)場強迫,夏季呈現(xiàn)NE向,冬季為W向流等特征[25],但又略有差異,主要體現(xiàn)在本研究的觀測結(jié)果顯示粵東沿岸流在大亞灣沿岸冬季為SW向,這可能與雷達(dá)觀測結(jié)果比較精細(xì)化有關(guān)。本研究的觀測結(jié)果在一定程度上說明高頻地波雷達(dá)對粵東沿岸流、珠江沖淡水等大尺度流系觀測中的突出作用,對比以往只能通過現(xiàn)場觀測結(jié)合數(shù)值模擬的方式,雷達(dá)有著顯著優(yōu)勢。由于雷達(dá)海流探測精度的限制,在余流較小的4月、5月和9月,余流的分布特征可能存在較大誤差,分析結(jié)論僅供參考。

    圖11 2015年10—12月大亞灣附近海域ASCAT遙感風(fēng)場Fig. 11 ASCAT remote sensing wind field in adjacent waters of Daya Bay from October to December, 2015

    圖12 2015年10—12月大亞灣附近海域A點日均風(fēng)矢量變化 Fig. 12 Change of the average daily wind vector at point A in adjacent waters of Daya Bay from October to December, 2015

    圖13 2015年10—12月大亞灣附近海域地波雷達(dá)表層余流場Fig. 13 Surface residual current field observed by ground wave radar in adjacent waters of Daya Bay from October to December, 2015

    3 結(jié)論

    本研究通過利用T_tide工具包對高頻地波雷達(dá)觀測數(shù)據(jù)進(jìn)行潮流調(diào)和分析,獲得如下結(jié)論:

    (1)大亞灣海域表層余流的變化特征主要呈冬、夏兩種形態(tài),其變化過程與季風(fēng)的轉(zhuǎn)變基本同步。其中冬季形態(tài)持續(xù)半年以上(從10月到次年3月),流向為SW向,流速可達(dá)30 cm/s以上; 夏季形態(tài)僅在5—8月間出現(xiàn),流向為NE向,流速可達(dá)23 cm/s,這與該區(qū)域受到東亞季風(fēng)控制,每年10月至次年3月盛行東北季風(fēng),6—8月盛行西南季風(fēng)相一致。

    (2)大亞灣附近海域的余流分布特征不僅與海表面風(fēng)場變化相關(guān),而且與海底地形也存在顯著相關(guān)關(guān)系。此外,粵東沿岸流與珠江沖淡水對大亞灣灣口余流空間分布具有調(diào)控作用。

    (3)高頻地波雷達(dá)可在粵東沿岸流、珠江沖淡水等大尺度流系的觀測中發(fā)揮重要作用。

    猜你喜歡
    余流大亞灣風(fēng)場
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    大亞灣造礁石珊瑚現(xiàn)狀
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    基于長期觀測的遼東灣口東部海域水動力特征研究
    側(cè)向風(fēng)場中無人機的飛行研究
    夏秋季泉州灣中部海域潮流和余流的變化特征
    基于走航ADCP資料的廈門內(nèi)灣東西口門海域潮流與潮致余流特征分析
    基于SOP規(guī)程的大亞灣堆芯冷卻監(jiān)測系統(tǒng)改造
    基于HJ-1B星熱紅外數(shù)據(jù)的大亞灣海溫反演模型及其應(yīng)用
    大亞灣反應(yīng)堆中微子實驗的快中子本底研究
    av一本久久久久| 国产伦精品一区二区三区四那| 亚洲性久久影院| 亚洲精品aⅴ在线观看| 九九爱精品视频在线观看| 91精品伊人久久大香线蕉| 少妇被粗大猛烈的视频| 中文字幕久久专区| 91成人精品电影| 亚洲高清免费不卡视频| 欧美精品人与动牲交sv欧美| 亚洲成人一二三区av| 爱豆传媒免费全集在线观看| 免费观看在线日韩| 午夜免费观看性视频| 成人美女网站在线观看视频| 成年女人在线观看亚洲视频| 内地一区二区视频在线| 久久久久久久久久久丰满| 中文字幕久久专区| 欧美激情国产日韩精品一区| 亚洲av二区三区四区| 欧美三级亚洲精品| 中文欧美无线码| 亚洲美女黄色视频免费看| 自拍偷自拍亚洲精品老妇| 九色成人免费人妻av| 亚洲va在线va天堂va国产| 亚洲,一卡二卡三卡| 大码成人一级视频| 99热这里只有是精品50| 精品少妇内射三级| 青春草亚洲视频在线观看| 波野结衣二区三区在线| 欧美激情国产日韩精品一区| 久久婷婷青草| 一级毛片aaaaaa免费看小| 人妻 亚洲 视频| 内射极品少妇av片p| 国产男女超爽视频在线观看| 天天操日日干夜夜撸| 91aial.com中文字幕在线观看| 久久久久精品久久久久真实原创| 在线观看三级黄色| 美女xxoo啪啪120秒动态图| 午夜福利视频精品| 丝袜脚勾引网站| 亚洲高清免费不卡视频| 国产 精品1| 精品国产国语对白av| 一区二区三区四区激情视频| 99视频精品全部免费 在线| 国产欧美亚洲国产| 亚洲欧美日韩东京热| 国产探花极品一区二区| 亚洲欧美日韩东京热| 久久久久久久久久人人人人人人| 欧美日韩一区二区视频在线观看视频在线| 免费观看在线日韩| 欧美亚洲 丝袜 人妻 在线| 女人久久www免费人成看片| 国产高清不卡午夜福利| 久久久午夜欧美精品| 国产午夜精品久久久久久一区二区三区| 国产精品一区二区性色av| 亚洲欧美精品自产自拍| 精品国产一区二区三区久久久樱花| 女性生殖器流出的白浆| 尾随美女入室| 亚洲丝袜综合中文字幕| 男女啪啪激烈高潮av片| 亚洲天堂av无毛| av又黄又爽大尺度在线免费看| 青春草亚洲视频在线观看| 日本爱情动作片www.在线观看| www.av在线官网国产| 黄色日韩在线| 插逼视频在线观看| 看十八女毛片水多多多| 亚洲国产日韩一区二区| 亚洲久久久国产精品| 777米奇影视久久| 亚洲av成人精品一二三区| 国产国拍精品亚洲av在线观看| 大片电影免费在线观看免费| 国产亚洲欧美精品永久| 亚洲色图综合在线观看| 亚洲国产色片| 国产91av在线免费观看| 久热这里只有精品99| 亚洲精品视频女| 色婷婷av一区二区三区视频| 欧美性感艳星| 亚洲国产色片| 高清在线视频一区二区三区| 欧美高清成人免费视频www| 人人妻人人爽人人添夜夜欢视频 | 大片电影免费在线观看免费| 色视频www国产| 在线精品无人区一区二区三| 草草在线视频免费看| 99久久综合免费| 日本91视频免费播放| 欧美三级亚洲精品| 亚洲av电影在线观看一区二区三区| 国产精品一二三区在线看| 亚洲,欧美,日韩| 国产 一区精品| 在线观看免费高清a一片| 成年av动漫网址| 全区人妻精品视频| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久精品古装| 亚洲国产日韩一区二区| 国产精品国产三级国产av玫瑰| 曰老女人黄片| 日韩欧美一区视频在线观看 | 一区二区三区四区激情视频| 免费不卡的大黄色大毛片视频在线观看| 视频区图区小说| 精品久久久久久电影网| 视频区图区小说| 中文精品一卡2卡3卡4更新| 七月丁香在线播放| 一级毛片 在线播放| 亚洲伊人久久精品综合| 亚洲怡红院男人天堂| 日本av免费视频播放| 亚洲精品久久久久久婷婷小说| 午夜福利在线观看免费完整高清在| 国产成人精品久久久久久| 国产成人精品福利久久| 精品熟女少妇av免费看| 少妇高潮的动态图| 两个人的视频大全免费| 国产成人精品久久久久久| 夜夜看夜夜爽夜夜摸| 亚洲精品aⅴ在线观看| 国产一区二区在线观看日韩| 亚洲精品乱久久久久久| 久久人人爽av亚洲精品天堂| 中文乱码字字幕精品一区二区三区| 老熟女久久久| 久久久久久久久久成人| 纵有疾风起免费观看全集完整版| 视频区图区小说| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 亚洲国产毛片av蜜桃av| 性高湖久久久久久久久免费观看| 精品一品国产午夜福利视频| 看非洲黑人一级黄片| 欧美日韩国产mv在线观看视频| 99热这里只有是精品50| av播播在线观看一区| 中文字幕av电影在线播放| 免费黄网站久久成人精品| 成年美女黄网站色视频大全免费 | av卡一久久| 欧美一级a爱片免费观看看| 欧美精品一区二区大全| 久久久久久久精品精品| 国产高清三级在线| 欧美精品一区二区免费开放| 毛片一级片免费看久久久久| 欧美老熟妇乱子伦牲交| 亚洲国产av新网站| 久久精品夜色国产| 波野结衣二区三区在线| 9色porny在线观看| 少妇精品久久久久久久| 在线观看免费日韩欧美大片 | 久热久热在线精品观看| 99精国产麻豆久久婷婷| 亚洲欧洲精品一区二区精品久久久 | 久久6这里有精品| 极品少妇高潮喷水抽搐| 亚洲精品日韩在线中文字幕| 最后的刺客免费高清国语| 另类亚洲欧美激情| 一本—道久久a久久精品蜜桃钙片| 夜夜骑夜夜射夜夜干| 香蕉精品网在线| 欧美成人午夜免费资源| 男人狂女人下面高潮的视频| 国产色爽女视频免费观看| 亚洲成人手机| 日韩精品有码人妻一区| 精品久久久久久久久亚洲| 亚洲国产av新网站| 久久久久久久久大av| 夫妻午夜视频| 99久久精品国产国产毛片| 日本91视频免费播放| 在线天堂最新版资源| 五月玫瑰六月丁香| 精品国产一区二区三区久久久樱花| 2022亚洲国产成人精品| 一边亲一边摸免费视频| 99九九在线精品视频 | 不卡视频在线观看欧美| 欧美精品人与动牲交sv欧美| 伦理电影免费视频| 亚洲国产欧美日韩在线播放 | 高清不卡的av网站| 极品少妇高潮喷水抽搐| 国产一区二区三区av在线| 夜夜骑夜夜射夜夜干| 亚洲欧美成人精品一区二区| 五月玫瑰六月丁香| 美女视频免费永久观看网站| 成人毛片a级毛片在线播放| 久久人人爽av亚洲精品天堂| 精品一品国产午夜福利视频| 99热网站在线观看| 18禁动态无遮挡网站| 国产黄片美女视频| 久久综合国产亚洲精品| 观看免费一级毛片| av在线观看视频网站免费| 91精品伊人久久大香线蕉| 国产在线男女| 久久久久精品久久久久真实原创| 国产亚洲一区二区精品| 九九在线视频观看精品| 免费观看无遮挡的男女| 欧美人与善性xxx| 日韩三级伦理在线观看| 亚洲精品一区蜜桃| 汤姆久久久久久久影院中文字幕| 日韩亚洲欧美综合| 人人妻人人澡人人爽人人夜夜| 少妇的逼好多水| 亚洲人成网站在线播| av网站免费在线观看视频| 亚洲成人手机| 亚洲无线观看免费| 夫妻午夜视频| 99热这里只有是精品50| 精品久久久久久久久亚洲| 久久韩国三级中文字幕| 亚洲精品成人av观看孕妇| 日韩不卡一区二区三区视频在线| 亚洲一级一片aⅴ在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲欧美日韩另类电影网站| 成年女人在线观看亚洲视频| 日本免费在线观看一区| 一级毛片 在线播放| 国内精品宾馆在线| 成人综合一区亚洲| 天天操日日干夜夜撸| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 少妇的逼水好多| 久久99精品国语久久久| 青青草视频在线视频观看| 久久久久久久大尺度免费视频| 亚洲欧美日韩东京热| 亚洲欧美成人综合另类久久久| 99久久精品国产国产毛片| 久久这里有精品视频免费| av在线观看视频网站免费| 最近中文字幕高清免费大全6| h视频一区二区三区| 伦理电影免费视频| 午夜91福利影院| 99re6热这里在线精品视频| 一本—道久久a久久精品蜜桃钙片| 欧美丝袜亚洲另类| 欧美精品一区二区大全| 国产精品熟女久久久久浪| www.av在线官网国产| 激情五月婷婷亚洲| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 十八禁高潮呻吟视频 | 久久久久网色| 你懂的网址亚洲精品在线观看| 日本av手机在线免费观看| 国产欧美日韩精品一区二区| 一级毛片 在线播放| 亚洲av不卡在线观看| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区国产| 纵有疾风起免费观看全集完整版| 国产在视频线精品| 看十八女毛片水多多多| 一二三四中文在线观看免费高清| 精品人妻熟女毛片av久久网站| √禁漫天堂资源中文www| 精品亚洲成a人片在线观看| 3wmmmm亚洲av在线观看| 国产日韩欧美视频二区| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 亚洲成人av在线免费| 午夜免费鲁丝| 中文资源天堂在线| 国产亚洲最大av| 久久久a久久爽久久v久久| 亚洲av成人精品一二三区| 妹子高潮喷水视频| 日本av免费视频播放| 亚洲欧美一区二区三区国产| 男女啪啪激烈高潮av片| 乱系列少妇在线播放| 高清午夜精品一区二区三区| 国产免费一区二区三区四区乱码| av福利片在线观看| 免费观看无遮挡的男女| 少妇高潮的动态图| 亚洲欧美日韩卡通动漫| 五月开心婷婷网| 观看美女的网站| av专区在线播放| 在线观看三级黄色| 成人美女网站在线观看视频| 国产男人的电影天堂91| 男人添女人高潮全过程视频| 日韩av免费高清视频| 一个人免费看片子| 我的女老师完整版在线观看| 水蜜桃什么品种好| 国产亚洲av片在线观看秒播厂| 深夜a级毛片| 久久久久国产精品人妻一区二区| 91久久精品电影网| 色婷婷久久久亚洲欧美| 边亲边吃奶的免费视频| 老司机亚洲免费影院| 在线天堂最新版资源| 国产免费视频播放在线视频| 黄色配什么色好看| 日本爱情动作片www.在线观看| 色5月婷婷丁香| a级毛片免费高清观看在线播放| 夜夜爽夜夜爽视频| 亚洲精品色激情综合| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 夜夜骑夜夜射夜夜干| 曰老女人黄片| 日韩亚洲欧美综合| 久久久a久久爽久久v久久| 国产av精品麻豆| 一本色道久久久久久精品综合| 国产精品久久久久成人av| av视频免费观看在线观看| 亚洲欧美精品专区久久| videossex国产| 日韩欧美一区视频在线观看 | 精品少妇内射三级| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 亚洲欧美一区二区三区国产| 精品人妻熟女毛片av久久网站| av又黄又爽大尺度在线免费看| 建设人人有责人人尽责人人享有的| 人体艺术视频欧美日本| av天堂久久9| 亚洲精品中文字幕在线视频 | 在线观看av片永久免费下载| 亚洲精品乱久久久久久| 成人无遮挡网站| 国产91av在线免费观看| 夜夜看夜夜爽夜夜摸| 中文精品一卡2卡3卡4更新| 国产日韩欧美视频二区| 国产欧美亚洲国产| 美女主播在线视频| 国产免费视频播放在线视频| 国产爽快片一区二区三区| 伦理电影大哥的女人| 99热网站在线观看| 国产成人aa在线观看| 精品人妻熟女毛片av久久网站| 国产亚洲91精品色在线| 国产精品伦人一区二区| 久久99热6这里只有精品| 国产色爽女视频免费观看| 国产精品一区二区性色av| 欧美少妇被猛烈插入视频| 欧美xxxx性猛交bbbb| kizo精华| 免费人妻精品一区二区三区视频| 国产综合精华液| 亚洲不卡免费看| 国产免费福利视频在线观看| 成人影院久久| 99热这里只有精品一区| 两个人的视频大全免费| 亚洲精品久久午夜乱码| 久久这里有精品视频免费| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩另类电影网站| 亚洲av二区三区四区| 国产精品欧美亚洲77777| 熟妇人妻不卡中文字幕| 寂寞人妻少妇视频99o| 人人澡人人妻人| 成人美女网站在线观看视频| 少妇人妻一区二区三区视频| 一本大道久久a久久精品| 色5月婷婷丁香| 色哟哟·www| 在线亚洲精品国产二区图片欧美 | av在线老鸭窝| 亚洲电影在线观看av| 十分钟在线观看高清视频www | av在线老鸭窝| 91久久精品电影网| 青春草亚洲视频在线观看| 97在线人人人人妻| 天天躁夜夜躁狠狠久久av| 国产免费一区二区三区四区乱码| 伊人久久精品亚洲午夜| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 免费av不卡在线播放| 99热国产这里只有精品6| 黄色毛片三级朝国网站 | 精品久久久久久久久亚洲| 我要看黄色一级片免费的| 人妻一区二区av| 99久久中文字幕三级久久日本| 青春草国产在线视频| 中文字幕人妻丝袜制服| 久久午夜综合久久蜜桃| 少妇精品久久久久久久| 如日韩欧美国产精品一区二区三区 | 成年美女黄网站色视频大全免费 | 嫩草影院入口| 国产精品人妻久久久久久| 99久久中文字幕三级久久日本| 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 日韩欧美精品免费久久| 久久人妻熟女aⅴ| 我要看黄色一级片免费的| 亚洲欧美一区二区三区黑人 | 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片| 欧美日韩亚洲高清精品| 女的被弄到高潮叫床怎么办| 亚洲四区av| 国产成人一区二区在线| 97超视频在线观看视频| 亚洲欧美一区二区三区国产| 你懂的网址亚洲精品在线观看| 一区二区三区免费毛片| 欧美xxⅹ黑人| 日韩一区二区三区影片| 又大又黄又爽视频免费| 精品久久久噜噜| 夫妻午夜视频| 日韩,欧美,国产一区二区三区| 美女福利国产在线| 亚洲精华国产精华液的使用体验| 有码 亚洲区| 中文在线观看免费www的网站| 欧美日韩在线观看h| 成年人午夜在线观看视频| 亚洲在久久综合| 午夜久久久在线观看| 十八禁高潮呻吟视频 | 国产精品一区www在线观看| 人人妻人人澡人人爽人人夜夜| 人人澡人人妻人| 赤兔流量卡办理| 各种免费的搞黄视频| 午夜精品国产一区二区电影| 国产午夜精品久久久久久一区二区三区| 秋霞在线观看毛片| 亚洲精品亚洲一区二区| 人妻少妇偷人精品九色| 制服丝袜香蕉在线| 在线观看免费高清a一片| 尾随美女入室| 超碰97精品在线观看| 久久久久久久久久久久大奶| 久久6这里有精品| 亚洲久久久国产精品| 少妇裸体淫交视频免费看高清| 精品视频人人做人人爽| 亚洲av成人精品一二三区| 特大巨黑吊av在线直播| 男女啪啪激烈高潮av片| 精品一区二区三区视频在线| 久久国产精品大桥未久av | 精品少妇黑人巨大在线播放| av在线播放精品| 国产有黄有色有爽视频| 麻豆乱淫一区二区| 日韩不卡一区二区三区视频在线| 成年女人在线观看亚洲视频| 日韩成人av中文字幕在线观看| 人体艺术视频欧美日本| 久久久久国产精品人妻一区二区| 最新中文字幕久久久久| 欧美97在线视频| a级毛片免费高清观看在线播放| 国产成人freesex在线| 高清不卡的av网站| 国产白丝娇喘喷水9色精品| 国产成人免费无遮挡视频| 欧美日本中文国产一区发布| 亚洲精品久久午夜乱码| 欧美激情国产日韩精品一区| 日本av手机在线免费观看| 欧美人与善性xxx| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频 | 日韩一区二区视频免费看| 在线观看av片永久免费下载| 人妻人人澡人人爽人人| 黄色配什么色好看| 少妇裸体淫交视频免费看高清| 80岁老熟妇乱子伦牲交| xxx大片免费视频| 涩涩av久久男人的天堂| 最近最新中文字幕免费大全7| 国产亚洲av片在线观看秒播厂| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 18禁裸乳无遮挡动漫免费视频| 少妇精品久久久久久久| xxx大片免费视频| 国产精品麻豆人妻色哟哟久久| av又黄又爽大尺度在线免费看| av线在线观看网站| 一级毛片 在线播放| 久久精品国产鲁丝片午夜精品| 久久精品久久久久久噜噜老黄| 妹子高潮喷水视频| 日日啪夜夜撸| 水蜜桃什么品种好| 我的老师免费观看完整版| 99久久精品国产国产毛片| 插逼视频在线观看| 国模一区二区三区四区视频| 看免费成人av毛片| 国精品久久久久久国模美| 国产探花极品一区二区| 91久久精品国产一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 内地一区二区视频在线| 日韩av不卡免费在线播放| 美女内射精品一级片tv| 国产精品99久久久久久久久| 丰满乱子伦码专区| 欧美老熟妇乱子伦牲交| 在线观看国产h片| 中文资源天堂在线| 少妇裸体淫交视频免费看高清| 18+在线观看网站| 欧美性感艳星| 少妇高潮的动态图| 国产片特级美女逼逼视频| 美女大奶头黄色视频| 51国产日韩欧美| 午夜日本视频在线| 国产精品国产三级国产专区5o| 国产熟女午夜一区二区三区 | 亚洲精品中文字幕在线视频 | 一级二级三级毛片免费看| 少妇被粗大的猛进出69影院 | 午夜免费观看性视频| 亚洲内射少妇av| 亚洲精华国产精华液的使用体验| 男女国产视频网站| 两个人免费观看高清视频 | 国产精品麻豆人妻色哟哟久久| 精品久久国产蜜桃| 美女中出高潮动态图| 一边亲一边摸免费视频| 久久狼人影院| 91久久精品国产一区二区三区| 欧美人与善性xxx| 六月丁香七月| 免费观看在线日韩| 亚洲欧美中文字幕日韩二区| 交换朋友夫妻互换小说| 狂野欧美白嫩少妇大欣赏| 久久久国产一区二区| 一区二区三区乱码不卡18| 久久国产精品大桥未久av | 久久av网站| 久久精品久久精品一区二区三区| 欧美激情极品国产一区二区三区 | 亚洲国产av新网站| 国产亚洲av片在线观看秒播厂| 国产欧美另类精品又又久久亚洲欧美| 大又大粗又爽又黄少妇毛片口| 国产精品三级大全| 亚洲人成网站在线观看播放| 国产成人免费无遮挡视频| 日韩av在线免费看完整版不卡| av黄色大香蕉| 极品人妻少妇av视频| 欧美精品高潮呻吟av久久| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 日日爽夜夜爽网站| 亚洲精品久久久久久婷婷小说| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| 女人久久www免费人成看片| 欧美日韩国产mv在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 国产成人aa在线观看| 午夜av观看不卡| 免费黄网站久久成人精品| 最近最新中文字幕免费大全7| 国产成人精品一,二区|