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

    南大西洋至幾內(nèi)亞灣冬季海浪傳播特性

    2023-08-02 11:10:34徐福敏牙韓爭(zhēng)周昕偉朱冬琳
    關(guān)鍵詞:模型研究

    徐福敏, 牙韓爭(zhēng), 周昕偉, 朱冬琳

    (1.河海大學(xué) 港口海岸與近海工程學(xué)院,江蘇 南京 210098;2.河海大學(xué) 海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210098;3.華設(shè)設(shè)計(jì)集團(tuán)股份有限公司,江蘇 南京 210014)

    幾內(nèi)亞灣位于赤道附近的大西洋東部,氣候條件特殊,北部屬于赤道低氣壓帶,南部為南半球東南信風(fēng)帶,同時(shí)受南大西洋復(fù)雜氣候的影響[1-2]。獨(dú)特的氣候以及直面大西洋的地理?xiàng)l件,使得幾內(nèi)亞灣近岸頻繁遭受大西洋的海浪襲擊。2007年和2011年2次極端海浪事件造成了科特迪瓦多處房屋被摧毀和約12m岸線的侵蝕[3]。在海浪長(zhǎng)期影響下,幾內(nèi)亞灣海岸功能逐漸退化,高脆弱區(qū)和極高脆弱區(qū)分別可達(dá)整個(gè)岸線的26.98%和11.66%[4]。幾內(nèi)亞灣作為“一帶一路”倡議的重要節(jié)點(diǎn),開(kāi)展南大西洋至幾內(nèi)亞灣海浪傳播特性研究對(duì)幾內(nèi)亞灣沿岸國(guó)家社會(huì)經(jīng)濟(jì)、交通航運(yùn)、海岸工程建設(shè)、海岸減災(zāi)防災(zāi)工作具有重要意義。

    有關(guān)南大西洋至幾內(nèi)亞灣的海浪研究由來(lái)已久。在大洋尺度上,早在20世紀(jì)80年代國(guó)內(nèi)就有南大西洋水文、氣候圖集的報(bào)道[5-6]。劉金芳等[7]根據(jù)1950—1995年間船舶氣象資料,對(duì)南大西洋風(fēng)、浪要素進(jìn)行統(tǒng)計(jì)分析,結(jié)果表明南大西洋高緯度區(qū)(30°S以南)由于西風(fēng)帶的影響常年為暴風(fēng)區(qū),大浪、大涌頻率分別為5%~25%、10%~35%;低緯度區(qū)(30°S以北)海浪和涌浪強(qiáng)度、頻率雖然均低于高緯度區(qū),但其大浪、大涌依然不容忽視(尤其冬季)。Violante-Carvalho等[8]基于浮標(biāo)數(shù)據(jù)研究南大西洋長(zhǎng)周期涌浪對(duì)風(fēng)浪的影響。Liu和Zhao[9]基于SWAN(Simulating WAves Nearshore)海浪模型對(duì)全球海浪分布及傳播規(guī)律進(jìn)行模擬研究,發(fā)現(xiàn)南大西洋涌浪多集中于東部海域。近年來(lái)針對(duì)幾內(nèi)亞灣海浪的研究不斷取得進(jìn)展。Stive等[10]和Almar等[11]研究發(fā)現(xiàn)幾內(nèi)亞灣海域大部分時(shí)間風(fēng)浪和涌浪共存,其中涌浪起主導(dǎo)作用,這些涌浪通常被認(rèn)為來(lái)自于遙遠(yuǎn)的南大西洋咆哮西風(fēng)帶。Olagnon等[12]和Forristall等[13]結(jié)合實(shí)測(cè)數(shù)據(jù)對(duì)西非海域海浪特征進(jìn)行研究,發(fā)現(xiàn)幾內(nèi)亞灣海浪為涌浪主導(dǎo)狀態(tài)。Ondoa等[14]認(rèn)為幾內(nèi)亞灣海浪主要由赤道信風(fēng)區(qū)域風(fēng)浪、咆哮西風(fēng)帶涌浪以及非洲西南海域涌浪組成。王科華和張軍[15]對(duì)幾內(nèi)亞灣若干碼頭項(xiàng)目海浪數(shù)學(xué)模型(MIKE 21)結(jié)果進(jìn)行分析,提出西非的幾內(nèi)亞灣屬長(zhǎng)周期涌浪作用為主的海灣。周昕偉等[2]基于SWAN海浪模型,研究幾內(nèi)亞灣海浪分布特征及傳播規(guī)律,發(fā)現(xiàn)幾內(nèi)亞灣處于支配地位且發(fā)源于南大西洋咆哮西風(fēng)帶。近年來(lái),再分析數(shù)據(jù)也逐漸使用于南大西洋相關(guān)海浪特征研究中,Gramcianinov等[16]使用ERA-5再分析數(shù)據(jù)研究南大西洋西部的極端海浪事件與溫帶氣旋之間的關(guān)系。ERA-5作為于歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasting, ECMWF)最新推出的第五代再分析全球水文氣象數(shù)據(jù)集,其不僅提供了混合浪、風(fēng)浪和涌浪的常規(guī)海浪參數(shù)(波高、波周期、波向),還提供涌浪主成份分析等參數(shù),目前常用于不同尺度的海浪特征研究中[16-19]。

    雖然已有研究結(jié)果表明影響幾內(nèi)亞灣的海浪與南大西洋高緯度區(qū)的咆哮西風(fēng)帶(40°S ~60°S)涌浪有關(guān)[3,10-11],但針對(duì)其影響過(guò)程的研究相對(duì)較少,其影響機(jī)理還尚未深入研究和探討。涌浪從咆哮西風(fēng)帶長(zhǎng)距離傳播至幾內(nèi)亞灣,其傳播機(jī)制仍有待深入研究。本文以南大西洋咆哮西風(fēng)帶至幾內(nèi)亞灣海域?yàn)檠芯繉?duì)象,建立SWAN雙層嵌套海浪模型,分析南大西洋至幾內(nèi)亞灣冬季海浪時(shí)空分布,揭示南大西洋至幾內(nèi)亞灣的海浪傳播特性,在南大西洋至幾內(nèi)亞灣海域混合浪和涌浪間關(guān)系、海浪能量組成和傳播以及幾內(nèi)亞灣涌浪溯源研究方面做進(jìn)一步探討。

    1 海浪數(shù)值模型構(gòu)建

    數(shù)值模擬是目前大洋海浪研究中普遍使用的方法之一,其除了提供與再分析數(shù)據(jù)類似的常規(guī)海浪參數(shù)外,還提供了一維和二維海浪譜等表述海浪組成和能量分布信息,使其廣泛地應(yīng)用于大洋海浪傳播特征的研究。本文采用第三代海浪模型SWAN對(duì)研究區(qū)域進(jìn)行模擬。SWAN模型經(jīng)過(guò)多年的發(fā)展和改進(jìn),目前廣泛應(yīng)用于大、中、小海域尺度的海浪模擬和應(yīng)用研究中[20-23],其模型控制方程為

    1.1 模擬區(qū)域設(shè)置及參數(shù)

    考慮德雷克海峽以及北大西洋海域可能產(chǎn)生的涌浪影響,將模型外層計(jì)算范圍覆蓋整個(gè)大西洋(D1,75°N~ 75°S,83°W~22°E),內(nèi)層范圍為幾內(nèi)亞灣及附近海域(D2,30°N~15°S,20°W~15°E),模型計(jì)算范圍見(jiàn)圖1。水深數(shù)據(jù)采用ETOPO1全球地形模型數(shù)據(jù)。驅(qū)動(dòng)風(fēng)場(chǎng)采用CCMP風(fēng)場(chǎng)數(shù)據(jù),其數(shù)據(jù)空間范圍為78.375°S~78.375°N和0.125°E~359.875°E,空間分辨率為0.25°×0.25°,時(shí)間分辨率為6 h。模型外層空間分辨率為0.25°×0.25°,計(jì)算步長(zhǎng)為15min,內(nèi)層空間分辨率為0.083°×0.083°,計(jì)算步長(zhǎng)為5min。譜頻率范圍為0.04~1.00Hz,相鄰離散頻率(f)間采用對(duì)數(shù)分布(fi+1=γfi,γ取1.106),共劃分33個(gè)。譜空間分布范圍為0°~360°,分辨率為10°。底摩擦采用Hasselmann等[24]的方法,系數(shù)為0.015;海浪破碎采用Battjes等[25]的方法,系數(shù)為0.73;風(fēng)能輸入采用Cavaleri 等[26]的方法;白浪耗散項(xiàng)采用Komen等[27]的方法;三波和四波非線性相互作用分別采用LTA法[28]和DIA算法[29]。

    圖1 嵌套模型計(jì)算范圍、驗(yàn)證點(diǎn)及內(nèi)層Jason-3衛(wèi)星軌跡Fig.1 Nesting model calculation range, verification points, and tracks of validation Jason-3

    1.2 模型驗(yàn)證

    1.2.1 模擬時(shí)段選擇

    南大西洋海浪季節(jié)變化較為明顯,夏季(每年12月至次年3月)為海浪較弱的時(shí)期,冬季(6月至9月)為海浪最活躍的時(shí)期,此時(shí)低緯度區(qū)域的幾內(nèi)亞灣也處于一年中海浪最活躍的時(shí)期[2,7]。Toualy等[3]基于AVISO和NWW3日平均數(shù)據(jù)對(duì)幾內(nèi)亞灣科特迪瓦附近海域有效波高進(jìn)行統(tǒng)計(jì)研究,結(jié)果表明8月為一年之中海浪最大波高出現(xiàn)頻率較大的月份。選取2017年8月為研究時(shí)段,考慮到模型計(jì)算穩(wěn)定和南大西洋海浪傳播時(shí)間,數(shù)值模擬起始時(shí)間在8月1日0時(shí)的基礎(chǔ)上提前9天,實(shí)際模型模擬時(shí)段為(2017年7月22日0時(shí)至2017年11月1日0時(shí))。

    1.2.2 浮標(biāo)觀測(cè)數(shù)據(jù)驗(yàn)證

    大西洋現(xiàn)有的海浪觀測(cè)浮標(biāo)大多位于北美和歐洲海域,赤道附近和南大西洋海浪觀測(cè)公開(kāi)數(shù)據(jù)匱乏。梁小力等[23]對(duì)國(guó)家海洋局全球海浪預(yù)報(bào)模型中大西洋海浪檢驗(yàn)時(shí)使用了美國(guó)國(guó)家海洋和大氣管理局浮標(biāo)數(shù)據(jù)中心(https://www.ndbc.noaa.gov/)包括浮標(biāo)41041(見(jiàn)圖1)的位于美國(guó)東岸的2個(gè)浮標(biāo)數(shù)據(jù)進(jìn)行模型驗(yàn)證。浮標(biāo)41041為與南大西洋和幾內(nèi)亞灣最近且具有公開(kāi)數(shù)據(jù)的觀測(cè)站。基于41041號(hào)浮標(biāo)有效波高和平均波周期,選取2017年8月(冬季、海浪強(qiáng))和2017年10月(冬夏轉(zhuǎn)換)對(duì)模型模擬結(jié)果進(jìn)行驗(yàn)證。圖2、圖3分別為41041號(hào)浮標(biāo)8月、10月有效波高、平均波周期和主波向驗(yàn)證,驗(yàn)證結(jié)果表明模擬結(jié)果和浮標(biāo)觀測(cè)數(shù)據(jù)變化趨勢(shì)基本一致,尤其2017年8月份有效波高、平均波周期和主波向均驗(yàn)證良好。季節(jié)轉(zhuǎn)換期的2017年10月模擬平均波周期和主波向與實(shí)測(cè)值略有偏差,但總體上模型計(jì)算結(jié)果與浮標(biāo)觀測(cè)數(shù)據(jù)吻合良好。

    圖2 41041號(hào)浮標(biāo)與SWAN模擬有效波高、平均波周期和主波向?qū)Ρ龋?017年8月)Fig.2 Comparison of SWAN stimulated significant wave heights, mean wave period, and dominant wave direction with Buoy 41041 in August, 2017

    圖3 41041號(hào)浮標(biāo)與SWAN模擬有效波高、平均波周期和主波向的對(duì)比(2017年10月)Fig.3 Comparison of SWAN stimulated significant wave heights, mean wave period, and dominant wave direction with Buoy 41041 in October, 2017

    1.2.3 衛(wèi)星高度計(jì)數(shù)據(jù)驗(yàn)證

    為彌補(bǔ)浮標(biāo)觀測(cè)數(shù)據(jù)的不足,采用Jason-3衛(wèi)星觀測(cè)數(shù)據(jù)對(duì)模型模擬結(jié)果進(jìn)行補(bǔ)充驗(yàn)證。Jason-3是由歐洲氣象衛(wèi)星組織(EUMESTAT)、美國(guó)國(guó)家航空航天局(NASA)、法國(guó)國(guó)家空間研究中心(CNES)和美國(guó)國(guó)家海洋和大氣管理局(NOAA)聯(lián)合研制的海洋地形衛(wèi)星,數(shù)據(jù)準(zhǔn)確度較高。選取經(jīng)過(guò)幾內(nèi)亞灣海域的6條衛(wèi)星軌跡,如圖1所示。模型模擬時(shí)段為2017年8月、10月,以衛(wèi)星數(shù)據(jù)(有效波高)進(jìn)行海浪模型驗(yàn)證(圖4、圖5)。由于衛(wèi)星經(jīng)過(guò)驗(yàn)證區(qū)域所需的時(shí)間極短(僅約20min),因此,在模型中設(shè)置3個(gè)計(jì)算時(shí)間節(jié)點(diǎn)(2個(gè)計(jì)算時(shí)間步長(zhǎng),即2×15min)的線輸出結(jié)果,對(duì)輸出結(jié)果在時(shí)間和空間上采用線性插值的方法得出所需的驗(yàn)證數(shù)據(jù)。由驗(yàn)證結(jié)果可知SWAN模擬海浪有效波高值與Jason-3觀測(cè)值整體變化趨勢(shì)一致,總體上吻合良好,表明模型可用于幾內(nèi)亞灣海浪數(shù)值模擬研究。

    圖4 Jason-3衛(wèi)星觀測(cè)有效波高與SWAN有效波高模擬值對(duì)比(2017年8月)Fig.4 Comparison of SWAN stimulated significant wave heights with Jason-3 data in August, 2017

    圖5 Jason-3衛(wèi)星觀測(cè)有效波高與SWAN有效波高模擬值對(duì)比(2017年10月)Fig.5 Comparison of SWAN stimulated significant wave heights with Jason-3 data in October, 2017

    2 代表性時(shí)刻大西洋至幾內(nèi)亞灣海浪、涌浪分布特征

    根據(jù)南大西洋東西岸線呈喇叭形狀、幾內(nèi)亞灣橫跨赤道且南大西洋氣旋呈順時(shí)針的特點(diǎn),以幾內(nèi)亞灣典型(最大)有效波高為基點(diǎn)進(jìn)行研究分析。模擬結(jié)果表明,研究時(shí)段內(nèi)幾內(nèi)亞灣海浪最大有效波高為2.8m,出現(xiàn)于8月7日18時(shí)。因此,選取8月7日18時(shí)作為代表性時(shí)刻對(duì)大西洋至幾內(nèi)亞灣海浪特征進(jìn)行分析。

    2.1 風(fēng)場(chǎng)、海浪場(chǎng)分布特征

    圖6、圖7分別為8月7日18時(shí)模型驅(qū)動(dòng)風(fēng)場(chǎng)和大西洋至幾內(nèi)亞灣海浪分布圖。整體上,南大西洋和北大西洋風(fēng)場(chǎng)相對(duì)獨(dú)立,而北大西洋東部直至包括西撒哈拉的北非沿岸海域依然受到南大西洋海浪的影響;南大西洋風(fēng)、浪均強(qiáng)于北大西洋,南大西洋海浪有效波高呈咆哮西風(fēng)帶自高緯度向低緯度海域遞減;幾內(nèi)亞灣海域海浪基本上只和南大西洋海浪相關(guān)。南大西洋至幾內(nèi)亞灣有效波高為8.0~2.8m;咆哮西風(fēng)帶中部有效波高最高,此時(shí)海域有約25m·s-1的強(qiáng)氣旋;由于局部氣旋(最大風(fēng)速約20m·s-1)的影響,南大西洋中部35°S處有約5.0m的有效波高;35°S以北至幾內(nèi)亞灣大部分海域有效波高為2.0~2.8m。

    圖6 大西洋風(fēng)場(chǎng)Fig.6 Wind field of Atlantic

    圖7 大西洋海浪場(chǎng)Fig.7 Wave field of Atlantic

    對(duì)比圖6及圖7中海浪傳播方向和風(fēng)場(chǎng)可見(jiàn),高緯度海域海浪傳播方向基本與風(fēng)場(chǎng)相對(duì)應(yīng),說(shuō)明高緯度海域海浪形成及傳播與當(dāng)?shù)仫L(fēng)場(chǎng)有關(guān)。低緯度海域海浪傳播方向與風(fēng)場(chǎng)差異明顯,海域風(fēng)場(chǎng)總體為SE風(fēng),而海浪基本為S—SW浪,表明低緯度海域除當(dāng)?shù)仫L(fēng)浪外,還有來(lái)源于南大西洋高緯度區(qū)的S—SW涌浪。

    2.2 涌浪場(chǎng)分布特征

    圖8為8月7日18時(shí)大西洋涌浪分布圖。涌浪波高最大值出現(xiàn)于南大西洋咆哮西風(fēng)帶中部,波高超過(guò)6.0m;咆哮西風(fēng)帶以北涌浪波高明顯下降,除南大西洋中部因氣旋產(chǎn)生局部涌浪高值外,大部分海域涌浪波高為2.0~3.0m;幾內(nèi)亞灣此時(shí)有約2.5m高的涌浪。對(duì)比圖7、圖8可發(fā)現(xiàn),南大西洋涌浪與混合浪波高分布高度相關(guān),高緯度區(qū)涌浪波高明顯大于低緯度區(qū)。通過(guò)模擬時(shí)段內(nèi)涌浪波高變化分析得知,低緯度區(qū)的幾內(nèi)亞灣有效波高最大時(shí)刻同樣為涌浪波高最大的時(shí)刻,表明南大西洋涌浪在混合浪中所占比例較大,尤其在風(fēng)速較小的低緯度海域尤為顯著。根據(jù)海浪能量與有效波高平方成正比的關(guān)系[30],此時(shí)幾內(nèi)亞灣涌浪(有效波高約2.5m)能量可達(dá)混合浪(有效波高約2.8m)總能量的約80%,可見(jiàn),涌浪對(duì)低緯度的幾內(nèi)亞灣海浪有顯著的影響。

    圖8 大西洋涌浪場(chǎng)Fig.8 Swell field of Atlantic

    3 南大西洋至幾內(nèi)亞海浪譜特性

    為研究南大西洋至幾內(nèi)亞灣海浪傳播特性,在南大西洋中部至幾內(nèi)亞灣附近海域取3個(gè)研究點(diǎn),其中A點(diǎn)位于南大西洋30°S中部,B點(diǎn)位于南大西洋15°S中部,C點(diǎn)位于幾內(nèi)亞灣附近海域赤道上。繪制各研究點(diǎn)的二維海浪譜,A點(diǎn)、B點(diǎn)海浪譜由模型外層(D1)模擬輸出,C點(diǎn)海浪譜由內(nèi)層(D2)模擬輸出。各研究點(diǎn)位置如圖1及表1所示。

    表1 南大西洋至幾內(nèi)亞灣海浪特征研究點(diǎn)Tab.1 Wave characteristics of research points from South Atlantic to Gulf of Guinea

    圖9為8月7日18時(shí)各研究點(diǎn)的海浪譜。研究點(diǎn)A、B海浪譜呈多波峰的特性,表明該點(diǎn)位海浪由不同海浪組成成分構(gòu)成;C點(diǎn)海浪譜呈單峰狀態(tài)且高峰處NNW—NE有較高頻率能量分布。A、B和C點(diǎn)海浪譜的N—NE方位均存在低頻能量分布,該低頻能量分布集中,在空間和頻率分布范圍均較窄,說(shuō)明南大西洋中部至幾內(nèi)亞灣海域均存在穩(wěn)定的S—SW涌浪系統(tǒng)。從各點(diǎn)位及附近海域風(fēng)場(chǎng)來(lái)看,并沒(méi)有與主涌浪對(duì)應(yīng)的風(fēng)向,顯然該涌浪是由南大西洋高緯度海域傳播而來(lái)。

    圖9 南大西洋至幾內(nèi)亞灣海浪譜特征Fig.9 Wave spectrum characteristics from South Atlantic to Gulf of Guinea

    各研究點(diǎn)高頻區(qū)域均存在部分能量分布,其能量分布方向基本與風(fēng)向一致,此部分能量為當(dāng)?shù)仫L(fēng)場(chǎng)作用產(chǎn)生的風(fēng)浪。A點(diǎn)在海浪譜的ENE—W方位存在寬范圍的高頻能量分布,此為南大西洋中部氣旋引起的風(fēng)浪。B點(diǎn)在海浪譜的W—NNW方位存在高頻能量分布,此為局地風(fēng)生浪成分。從各點(diǎn)低頻、高頻區(qū)域能量及方位分布不難看出,南大西洋至赤道區(qū)域的幾內(nèi)亞灣海浪均包含長(zhǎng)周期涌浪成分,且為S—SW涌浪。南大西洋中部至幾內(nèi)亞灣由于頻繁的強(qiáng)順時(shí)針氣旋和風(fēng)切變,風(fēng)浪能量也不可忽視。

    幾內(nèi)亞灣呈明顯的涌浪主導(dǎo)狀態(tài)。雖然C點(diǎn)海浪譜的高頻區(qū)域有能量分布,但能量相對(duì)低頻區(qū)域明顯較小,此時(shí)幾內(nèi)亞灣海域最大風(fēng)速僅為5.3m·s-1,大部分海域風(fēng)速僅為2 ~3m·s-1。從海浪分布場(chǎng)也可看出,幾內(nèi)亞灣有效波高約2.8m,但涌浪波高就可達(dá)2.5m。對(duì)比幾內(nèi)亞灣和南大西洋海浪譜發(fā)現(xiàn),低頻能量集中方向十分相近(集中在N—NE方位),說(shuō)明幾內(nèi)亞灣海浪主要受南大西洋S—SW涌浪影響。

    4 幾內(nèi)亞灣涌浪溯源

    為了驗(yàn)證前文影響幾內(nèi)亞灣涌浪來(lái)源于南大西洋高緯度區(qū)咆哮西風(fēng)帶的觀點(diǎn),以幾內(nèi)亞灣有效波高較大的時(shí)刻(8月7日18時(shí))為例對(duì)涌浪進(jìn)行溯源分析。圖10為8月1日6時(shí)至8月7日18時(shí)(6.5天)大西洋3個(gè)時(shí)刻(8月1日6時(shí)、8月4日18時(shí)、8月7日18時(shí))的風(fēng)場(chǎng),圖11和圖12為分別為相應(yīng)時(shí)刻的大西洋譜峰周期分布和涌浪有效波高分布。8月1日6時(shí)咆哮西風(fēng)帶存在最大風(fēng)速超過(guò)30m·s-1的順時(shí)針渦旋(圖10a),該渦旋在相應(yīng)海域引起波高6~10m、譜峰周期12~16s的S—SW強(qiáng)浪(圖11a、圖12a)。隨著時(shí)間的推移,此部分海浪不斷傳播并由于速度不同逐漸分離。8月4日18時(shí),南大西洋中部海域已經(jīng)出現(xiàn)譜峰周期14~20s、波高2~3m的涌浪(圖11b、圖12b),譜峰周期等值線沿NW—SE方向規(guī)則排列,表明此部分涌浪均來(lái)源于同一個(gè)涌浪系統(tǒng),即8月1日6時(shí)咆哮西風(fēng)帶生成的強(qiáng)浪傳出強(qiáng)風(fēng)區(qū)后的涌浪。至8月7日18時(shí),該涌浪系統(tǒng)譜峰周期為14~16s的涌浪傳播至幾內(nèi)亞灣,由此產(chǎn)生了幾內(nèi)亞灣8月最強(qiáng)涌浪(圖11c、圖12c)。

    圖10 3個(gè)時(shí)刻大西洋風(fēng)場(chǎng)Fig.10 Wind field of Atlantic in three moments

    圖11 3個(gè)時(shí)刻大西洋譜峰周期分布Fig.11 Peak period distribution of Atlantic in three moments

    圖12 3個(gè)時(shí)刻大西洋涌浪有效波高分布Fig.12 Swell significant wave height distribution of Atlantic in three moments

    8月1日6 時(shí)至8月7日18時(shí)咆哮西風(fēng)帶間斷生成強(qiáng)氣旋并快速向東向移動(dòng)(2~3天內(nèi)即可橫向穿越大西洋),使得南大西洋咆哮西風(fēng)帶內(nèi)產(chǎn)生強(qiáng)海浪。咆哮西風(fēng)帶至20°S局部較強(qiáng)風(fēng)切變(如8月7日18時(shí)南大西洋中部35°S)或較小規(guī)模順時(shí)針氣旋(圖10c)也可能生成較大風(fēng)浪并形成小范圍涌浪。然而,其天氣系統(tǒng)規(guī)模生成的海浪乃至形成的涌浪均無(wú)法和咆哮西風(fēng)帶引起的強(qiáng)風(fēng)、強(qiáng)浪和涌浪相比。與此同時(shí),自8月1日至8月7日18時(shí)6.5天內(nèi)幾內(nèi)亞灣及其鄰近海域均不具有形成幾內(nèi)亞灣高混合浪、涌浪的氣候條件,可以推斷8月7日18時(shí)幾內(nèi)亞灣涌浪高值是由約6.5天前南大西洋咆哮西風(fēng)帶海域的強(qiáng)風(fēng)浪傳播而來(lái)。

    5 結(jié)論

    建立大西洋至幾內(nèi)亞灣的SWAN雙層嵌套海浪模型,以模擬時(shí)段內(nèi)幾內(nèi)亞灣海浪有效波高最大時(shí)刻(2017年8月7日18時(shí))為基點(diǎn),對(duì)冬季代表性時(shí)刻南大西洋至幾內(nèi)亞灣海域海浪特征進(jìn)行研究,得到了風(fēng)、混合浪和涌浪間關(guān)系,明確了海浪、涌浪分布特征;基于模擬得到的南大西洋至幾內(nèi)亞二維海浪譜,表明幾內(nèi)亞灣海浪主要受南大西洋S—SW涌浪影響,揭示了南大西洋至幾內(nèi)亞灣的海浪傳播特性;對(duì)代表性時(shí)刻幾內(nèi)亞灣涌浪進(jìn)行溯源,發(fā)現(xiàn)內(nèi)亞灣涌浪是由約6.5天前南大西洋咆哮西風(fēng)帶海域的強(qiáng)風(fēng)浪傳播而來(lái)。主要結(jié)論如下:

    (1)南大西洋海浪有效波高呈高緯度區(qū)咆哮西風(fēng)帶向低緯度區(qū)遞減的趨勢(shì)。南大西洋咆哮西風(fēng)帶中部有效波高最高可達(dá)約8.0m,除南大西洋中部因氣旋影響產(chǎn)生的局部浪高(約5.0m)外,至幾內(nèi)亞灣海域已降為2.0~2.8m。

    (2)高緯度區(qū)海浪傳播方向與風(fēng)場(chǎng)關(guān)系密切,而低緯度區(qū)則與風(fēng)場(chǎng)差異明顯。低緯度海域風(fēng)場(chǎng)總體為SE風(fēng),但海浪基本為S—SW浪,說(shuō)明低緯度海域不僅有局部風(fēng)浪,更大程度上還受南大西洋高緯度海域涌浪的影響。

    (3)南大西洋涌浪與混合浪波高分布高度相關(guān),且高緯度區(qū)涌浪波高明顯大于低緯度區(qū)。南大西洋咆哮西風(fēng)帶中部涌浪波高超過(guò)6.0m,除南大西洋中部海域因氣旋產(chǎn)生的局部涌浪高值外,低緯度大部分海域涌浪在2.0~3.0m間,幾內(nèi)亞灣涌浪波高已降至約2.5m。根據(jù)海浪能量—有效波高關(guān)系計(jì)算,幾內(nèi)亞灣涌浪能量可達(dá)混合浪總能量的約80%。

    (4)位于南大西洋的A、B和C 3個(gè)海浪研究點(diǎn)二維海浪譜在N—NE方向能量分布集中,且空間和頻率分布范圍均較窄。南大西洋中部至赤道附近的幾內(nèi)亞灣海域均存在穩(wěn)定的S—SW涌浪系統(tǒng),該涌浪由南大西洋高緯度海域傳來(lái)。

    (5)影響幾內(nèi)亞灣的涌浪主要來(lái)源于南大西洋咆哮西風(fēng)帶。通過(guò)對(duì)幾內(nèi)亞灣有效波高最大時(shí)刻的涌浪進(jìn)行溯源分析發(fā)現(xiàn),8月7日18時(shí)幾內(nèi)亞灣涌浪高值為約6.5天前的南大西洋咆哮西風(fēng)帶海域強(qiáng)風(fēng)浪傳播而來(lái)。

    本文選取幾內(nèi)亞灣外海發(fā)生強(qiáng)海浪頻率較大的冬季8月進(jìn)行南大西洋至幾內(nèi)亞灣冬季海浪傳播特性研究,還需對(duì)夏季以及轉(zhuǎn)換期的海浪傳播特性進(jìn)行深入研究,以獲得對(duì)該海域海浪特性的較全面認(rèn)識(shí);南大西洋的公開(kāi)現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)較為匱乏,需要更多和系統(tǒng)觀測(cè)數(shù)據(jù)對(duì)數(shù)值模型進(jìn)行進(jìn)一步驗(yàn)證和機(jī)理研究,并在此基礎(chǔ)上采用不同海浪譜分離機(jī)理的譜海浪模型(比如WWIII等)進(jìn)行對(duì)比模擬研究,以確定模擬南大西洋復(fù)雜海浪組成、特性和傳播的合理模式。

    作者貢獻(xiàn)聲明:

    徐福敏:構(gòu)思框架、論文思路、模型指導(dǎo)、確定文稿。

    牙韓爭(zhēng):處理數(shù)據(jù)、撰寫(xiě)論文、修改論文。

    周昕偉:構(gòu)建模型、審閱論文。

    朱冬琳:數(shù)值模型技術(shù)支持。

    猜你喜歡
    模型研究
    一半模型
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    3D打印中的模型分割與打包
    国内毛片毛片毛片毛片毛片| 国产成人免费无遮挡视频| 久久香蕉激情| 国产精品爽爽va在线观看网站 | av天堂在线播放| 麻豆成人av在线观看| 9热在线视频观看99| 国产精品 欧美亚洲| 99热国产这里只有精品6| 久久久久久久精品吃奶| 午夜福利,免费看| 国产男靠女视频免费网站| 亚洲欧美精品综合一区二区三区| 在线国产一区二区在线| 国内久久婷婷六月综合欲色啪| 国产精品秋霞免费鲁丝片| 大型av网站在线播放| 老汉色av国产亚洲站长工具| 亚洲国产毛片av蜜桃av| 亚洲精品成人av观看孕妇| 人妻丰满熟妇av一区二区三区| 无人区码免费观看不卡| 日本a在线网址| 91成年电影在线观看| 午夜福利一区二区在线看| 婷婷六月久久综合丁香| 欧美 亚洲 国产 日韩一| 欧美大码av| 久99久视频精品免费| 久久中文字幕一级| 丰满饥渴人妻一区二区三| 国内久久婷婷六月综合欲色啪| 91国产中文字幕| 大陆偷拍与自拍| 国产蜜桃级精品一区二区三区| 亚洲专区中文字幕在线| av免费在线观看网站| av免费在线观看网站| 国产精品1区2区在线观看.| 亚洲午夜精品一区,二区,三区| 色综合婷婷激情| 一本综合久久免费| 欧美日韩av久久| 成人av一区二区三区在线看| 亚洲av成人不卡在线观看播放网| 老汉色∧v一级毛片| 午夜两性在线视频| 亚洲专区国产一区二区| 久热爱精品视频在线9| 欧美激情久久久久久爽电影 | 别揉我奶头~嗯~啊~动态视频| 欧美精品一区二区免费开放| 中亚洲国语对白在线视频| 中文字幕人妻丝袜制服| av有码第一页| 19禁男女啪啪无遮挡网站| 免费女性裸体啪啪无遮挡网站| 91国产中文字幕| 日本撒尿小便嘘嘘汇集6| 国产精品亚洲av一区麻豆| 精品久久蜜臀av无| av片东京热男人的天堂| 91麻豆av在线| 操美女的视频在线观看| 99精品久久久久人妻精品| 男人操女人黄网站| 久久性视频一级片| a在线观看视频网站| 中国美女看黄片| 午夜福利在线观看吧| 久久精品91无色码中文字幕| 18禁美女被吸乳视频| 久久精品国产亚洲av高清一级| 免费高清在线观看日韩| 一个人免费在线观看的高清视频| 久久久水蜜桃国产精品网| 欧美黄色片欧美黄色片| 黄色视频不卡| 国产亚洲欧美在线一区二区| 欧美日本中文国产一区发布| 午夜免费观看网址| av在线播放免费不卡| 一级a爱片免费观看的视频| 黑人巨大精品欧美一区二区蜜桃| 女人被躁到高潮嗷嗷叫费观| 国产视频一区二区在线看| 少妇 在线观看| 亚洲午夜精品一区,二区,三区| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区久久| av国产精品久久久久影院| 久久久久精品国产欧美久久久| 美女 人体艺术 gogo| 国产熟女午夜一区二区三区| 如日韩欧美国产精品一区二区三区| 制服诱惑二区| 午夜精品在线福利| 99精品久久久久人妻精品| 制服诱惑二区| 中文欧美无线码| 久久九九热精品免费| 国产黄色免费在线视频| a级毛片在线看网站| 99精品久久久久人妻精品| 欧美乱妇无乱码| 一级片免费观看大全| 在线观看日韩欧美| 在线观看免费午夜福利视频| 在线观看日韩欧美| 美女扒开内裤让男人捅视频| 男女高潮啪啪啪动态图| 熟女少妇亚洲综合色aaa.| 在线av久久热| 成人三级做爰电影| av国产精品久久久久影院| av免费在线观看网站| 最好的美女福利视频网| 18禁国产床啪视频网站| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 级片在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产男靠女视频免费网站| 国产一区二区三区视频了| 亚洲成人久久性| 韩国精品一区二区三区| 欧美+亚洲+日韩+国产| 日本五十路高清| 国产又色又爽无遮挡免费看| 免费av毛片视频| www.999成人在线观看| 91精品三级在线观看| 99国产精品免费福利视频| 国产伦人伦偷精品视频| 美女 人体艺术 gogo| 最新美女视频免费是黄的| 久久影院123| 丰满人妻熟妇乱又伦精品不卡| 在线av久久热| 欧美日韩亚洲国产一区二区在线观看| 大码成人一级视频| 国产亚洲欧美精品永久| 精品国产乱子伦一区二区三区| 国产伦一二天堂av在线观看| 动漫黄色视频在线观看| 嫩草影视91久久| 精品一品国产午夜福利视频| www.999成人在线观看| 亚洲av成人av| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一小说| 99香蕉大伊视频| 视频区图区小说| 久久精品国产99精品国产亚洲性色 | 久久国产精品人妻蜜桃| 在线国产一区二区在线| 欧美精品亚洲一区二区| 免费在线观看影片大全网站| 天天躁狠狠躁夜夜躁狠狠躁| 另类亚洲欧美激情| 亚洲av电影在线进入| 久久人妻熟女aⅴ| 午夜免费激情av| 97超级碰碰碰精品色视频在线观看| 日韩欧美在线二视频| 久久青草综合色| 搡老熟女国产l中国老女人| 亚洲熟女毛片儿| 欧美在线黄色| www.999成人在线观看| 久久精品人人爽人人爽视色| 夜夜爽天天搞| 在线观看日韩欧美| 欧美另类亚洲清纯唯美| 怎么达到女性高潮| 看黄色毛片网站| 亚洲 欧美 日韩 在线 免费| ponron亚洲| 国产精品国产av在线观看| av有码第一页| 国产成人精品久久二区二区91| 两性夫妻黄色片| 精品熟女少妇八av免费久了| 老司机午夜十八禁免费视频| 大码成人一级视频| xxx96com| 不卡av一区二区三区| 神马国产精品三级电影在线观看 | 欧美av亚洲av综合av国产av| 久久精品国产亚洲av香蕉五月| 岛国视频午夜一区免费看| 日韩人妻精品一区2区三区| 中文欧美无线码| 国产国语露脸激情在线看| 精品欧美一区二区三区在线| 久久久久久久久免费视频了| 99精品在免费线老司机午夜| 亚洲中文日韩欧美视频| netflix在线观看网站| 老熟妇乱子伦视频在线观看| 欧美日韩黄片免| 99热只有精品国产| 亚洲aⅴ乱码一区二区在线播放 | 亚洲第一欧美日韩一区二区三区| 国产熟女xx| 精品人妻1区二区| 国产免费现黄频在线看| 午夜老司机福利片| 麻豆一二三区av精品| 悠悠久久av| 国产一卡二卡三卡精品| 欧美在线一区亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 咕卡用的链子| 国产男靠女视频免费网站| 精品久久久久久久毛片微露脸| 视频区图区小说| 在线观看舔阴道视频| 日韩精品青青久久久久久| 亚洲成av片中文字幕在线观看| 亚洲黑人精品在线| 亚洲欧美精品综合久久99| 可以在线观看毛片的网站| 国产精品99久久99久久久不卡| 日韩 欧美 亚洲 中文字幕| cao死你这个sao货| 天堂影院成人在线观看| 亚洲熟女毛片儿| 另类亚洲欧美激情| 18禁黄网站禁片午夜丰满| 超色免费av| 女人精品久久久久毛片| 国产成人av激情在线播放| 亚洲av成人一区二区三| 欧美激情久久久久久爽电影 | 热re99久久精品国产66热6| 欧美午夜高清在线| 欧美成人性av电影在线观看| 如日韩欧美国产精品一区二区三区| 黄色丝袜av网址大全| 色老头精品视频在线观看| ponron亚洲| 国产成年人精品一区二区 | 怎么达到女性高潮| 变态另类成人亚洲欧美熟女 | 国产不卡一卡二| 麻豆av在线久日| 国产99久久九九免费精品| 免费在线观看影片大全网站| 五月开心婷婷网| av福利片在线| 我的亚洲天堂| 黄色成人免费大全| 涩涩av久久男人的天堂| 日韩三级视频一区二区三区| 在线观看午夜福利视频| 久久热在线av| 国产亚洲精品久久久久久毛片| 波多野结衣高清无吗| 1024香蕉在线观看| 亚洲精华国产精华精| 97超级碰碰碰精品色视频在线观看| 成人影院久久| 亚洲va日本ⅴa欧美va伊人久久| 久久人人97超碰香蕉20202| 亚洲视频免费观看视频| 99国产综合亚洲精品| 后天国语完整版免费观看| 精品午夜福利视频在线观看一区| 免费日韩欧美在线观看| 亚洲国产精品sss在线观看 | av超薄肉色丝袜交足视频| 久久精品成人免费网站| 久久草成人影院| 国产精品野战在线观看 | 性色av乱码一区二区三区2| 国产熟女午夜一区二区三区| 午夜成年电影在线免费观看| 在线观看舔阴道视频| 999精品在线视频| 国产一区在线观看成人免费| 久久久久久久午夜电影 | 天天躁夜夜躁狠狠躁躁| 一区二区日韩欧美中文字幕| 欧美日韩亚洲综合一区二区三区_| 亚洲色图 男人天堂 中文字幕| 精品久久久久久成人av| 精品国产国语对白av| 亚洲精品在线美女| 亚洲片人在线观看| 亚洲九九香蕉| 国产av又大| 黄色怎么调成土黄色| 国产高清视频在线播放一区| 黄色视频不卡| 老司机午夜十八禁免费视频| 久久精品国产综合久久久| www.自偷自拍.com| av有码第一页| 亚洲色图 男人天堂 中文字幕| 国产99久久九九免费精品| 日韩免费av在线播放| 国产一区二区三区综合在线观看| 69精品国产乱码久久久| 日日爽夜夜爽网站| 国产av在哪里看| 黄色视频,在线免费观看| 人人妻人人澡人人看| 成人永久免费在线观看视频| 最好的美女福利视频网| 天天添夜夜摸| 最新在线观看一区二区三区| netflix在线观看网站| 久久人人97超碰香蕉20202| 色综合站精品国产| 日本vs欧美在线观看视频| 国产免费男女视频| 亚洲第一青青草原| 嫁个100分男人电影在线观看| 国产成人av激情在线播放| 成人国产一区最新在线观看| 51午夜福利影视在线观看| 人人澡人人妻人| 国产精品影院久久| 1024视频免费在线观看| 久久青草综合色| 日韩成人在线观看一区二区三区| 亚洲五月婷婷丁香| 日韩免费高清中文字幕av| 国产av一区二区精品久久| 欧美人与性动交α欧美精品济南到| 好看av亚洲va欧美ⅴa在| 日韩欧美三级三区| 国产97色在线日韩免费| 久热爱精品视频在线9| 亚洲精品粉嫩美女一区| 操美女的视频在线观看| 精品国产国语对白av| 亚洲自拍偷在线| 一级,二级,三级黄色视频| 人人澡人人妻人| 欧美一区二区精品小视频在线| 久久久水蜜桃国产精品网| 日韩有码中文字幕| 亚洲国产欧美网| 在线观看一区二区三区| 一区在线观看完整版| 国产1区2区3区精品| 91成人精品电影| 亚洲黑人精品在线| 丁香六月欧美| av视频免费观看在线观看| 亚洲五月婷婷丁香| 99国产精品99久久久久| 制服人妻中文乱码| 欧美激情久久久久久爽电影 | 久久精品亚洲熟妇少妇任你| 999久久久国产精品视频| 亚洲av日韩精品久久久久久密| 99国产综合亚洲精品| bbb黄色大片| 免费一级毛片在线播放高清视频 | 国产在线精品亚洲第一网站| 热re99久久精品国产66热6| 午夜福利一区二区在线看| 久久精品成人免费网站| 91大片在线观看| 99久久久亚洲精品蜜臀av| 欧美成人性av电影在线观看| www.熟女人妻精品国产| av免费在线观看网站| 欧美黄色片欧美黄色片| 后天国语完整版免费观看| 久久青草综合色| 男人舔女人的私密视频| videosex国产| av片东京热男人的天堂| 国产熟女xx| 亚洲国产毛片av蜜桃av| 国产精品久久久久久人妻精品电影| 多毛熟女@视频| 亚洲三区欧美一区| 欧美日韩中文字幕国产精品一区二区三区 | 成人18禁高潮啪啪吃奶动态图| 一进一出抽搐gif免费好疼 | 搡老熟女国产l中国老女人| av欧美777| 国产av又大| 美国免费a级毛片| 国产高清激情床上av| 欧美日韩一级在线毛片| 手机成人av网站| 日韩三级视频一区二区三区| 两人在一起打扑克的视频| 国产免费男女视频| 丰满的人妻完整版| 久久精品成人免费网站| 黄色怎么调成土黄色| 12—13女人毛片做爰片一| 国产精品日韩av在线免费观看 | 一级作爱视频免费观看| 高清欧美精品videossex| 午夜福利一区二区在线看| 久久精品国产亚洲av高清一级| 中出人妻视频一区二区| 又黄又爽又免费观看的视频| 他把我摸到了高潮在线观看| 国产av又大| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 久久亚洲真实| 两个人看的免费小视频| 999精品在线视频| 国产成人精品在线电影| 欧美中文日本在线观看视频| 成年版毛片免费区| 欧美激情 高清一区二区三区| 99在线视频只有这里精品首页| 黑人操中国人逼视频| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o| 两个人看的免费小视频| 大型黄色视频在线免费观看| 久久人妻福利社区极品人妻图片| 亚洲男人的天堂狠狠| 精品国产亚洲在线| 国产免费现黄频在线看| 精品午夜福利视频在线观看一区| 侵犯人妻中文字幕一二三四区| 99热只有精品国产| 三级毛片av免费| 色婷婷久久久亚洲欧美| 他把我摸到了高潮在线观看| 久久香蕉激情| 一本大道久久a久久精品| av视频免费观看在线观看| 久久香蕉激情| 高清欧美精品videossex| 日韩欧美一区二区三区在线观看| 99香蕉大伊视频| 亚洲自偷自拍图片 自拍| 丰满的人妻完整版| 日本一区二区免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品99久久久久| 99re在线观看精品视频| 亚洲人成电影观看| 丁香六月欧美| 免费少妇av软件| 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 亚洲av日韩精品久久久久久密| 人人澡人人妻人| 国产精品国产av在线观看| 欧美人与性动交α欧美精品济南到| 丁香六月欧美| 男女下面进入的视频免费午夜 | 黄色丝袜av网址大全| 动漫黄色视频在线观看| 午夜福利免费观看在线| 法律面前人人平等表现在哪些方面| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 欧美日韩黄片免| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 美国免费a级毛片| 88av欧美| 成年女人毛片免费观看观看9| 国产无遮挡羞羞视频在线观看| 欧美激情高清一区二区三区| 搡老乐熟女国产| 97碰自拍视频| 啦啦啦 在线观看视频| 免费在线观看亚洲国产| 亚洲国产精品合色在线| 在线观看一区二区三区| 真人一进一出gif抽搐免费| 老司机亚洲免费影院| 中文字幕人妻丝袜一区二区| 中文字幕色久视频| 免费在线观看影片大全网站| 国产成人精品久久二区二区91| 日日爽夜夜爽网站| 国产蜜桃级精品一区二区三区| 精品无人区乱码1区二区| 国产视频一区二区在线看| 中文字幕人妻丝袜制服| 人妻丰满熟妇av一区二区三区| 一a级毛片在线观看| 国产极品粉嫩免费观看在线| 国产精品98久久久久久宅男小说| 女性生殖器流出的白浆| 丰满人妻熟妇乱又伦精品不卡| 老司机靠b影院| 国产精品 国内视频| 色播在线永久视频| 美女高潮到喷水免费观看| 久久人妻av系列| 一级a爱片免费观看的视频| 如日韩欧美国产精品一区二区三区| 日韩免费高清中文字幕av| 嫩草影院精品99| 美女福利国产在线| 91精品三级在线观看| 脱女人内裤的视频| 99国产精品免费福利视频| 久久久国产一区二区| 亚洲成人免费av在线播放| 在线观看免费午夜福利视频| 精品国产亚洲在线| 久久久久国内视频| 国产精品1区2区在线观看.| 18禁美女被吸乳视频| 免费一级毛片在线播放高清视频 | 80岁老熟妇乱子伦牲交| 成熟少妇高潮喷水视频| 最好的美女福利视频网| 三级毛片av免费| 成人av一区二区三区在线看| 国产成人影院久久av| 欧美激情久久久久久爽电影 | 国产亚洲精品久久久久5区| 久久 成人 亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 国产又色又爽无遮挡免费看| 国产成人欧美在线观看| 精品人妻1区二区| 99香蕉大伊视频| 亚洲成人久久性| 男人的好看免费观看在线视频 | 国产视频一区二区在线看| 在线观看午夜福利视频| 国产精品国产高清国产av| 露出奶头的视频| 免费观看精品视频网站| 精品乱码久久久久久99久播| 免费久久久久久久精品成人欧美视频| 一a级毛片在线观看| 首页视频小说图片口味搜索| 麻豆久久精品国产亚洲av | 精品福利观看| 交换朋友夫妻互换小说| 一本大道久久a久久精品| 精品国产乱码久久久久久男人| 亚洲欧美一区二区三区久久| 久久久国产精品麻豆| 在线国产一区二区在线| 精品久久久久久久毛片微露脸| 中文字幕高清在线视频| 女性生殖器流出的白浆| 国产成人免费无遮挡视频| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图| 国产xxxxx性猛交| 久久精品国产99精品国产亚洲性色 | 久久精品成人免费网站| 一二三四在线观看免费中文在| 一边摸一边抽搐一进一小说| 精品少妇一区二区三区视频日本电影| 久久亚洲真实| a级毛片黄视频| 欧美在线一区亚洲| 午夜激情av网站| 色精品久久人妻99蜜桃| www日本在线高清视频| 99国产综合亚洲精品| 国产亚洲精品综合一区在线观看 | 黑丝袜美女国产一区| 亚洲一区二区三区不卡视频| cao死你这个sao货| 天堂俺去俺来也www色官网| 别揉我奶头~嗯~啊~动态视频| ponron亚洲| 国产精品亚洲av一区麻豆| avwww免费| 午夜久久久在线观看| 国产精品香港三级国产av潘金莲| 老司机午夜福利在线观看视频| 日韩精品青青久久久久久| 操出白浆在线播放| 精品一区二区三区视频在线观看免费 | 亚洲欧美一区二区三区黑人| 亚洲中文日韩欧美视频| 丝袜美腿诱惑在线| 国产成+人综合+亚洲专区| 欧美日韩亚洲综合一区二区三区_| 精品久久久精品久久久| 久久人妻福利社区极品人妻图片| 香蕉久久夜色| 天堂中文最新版在线下载| 热99国产精品久久久久久7| a级毛片在线看网站| 激情视频va一区二区三区| 日韩有码中文字幕| 国产亚洲精品第一综合不卡| 中文字幕色久视频| 久久狼人影院| 搡老乐熟女国产| 免费在线观看完整版高清| 淫妇啪啪啪对白视频| 亚洲精品美女久久av网站| 久久人人爽av亚洲精品天堂| 亚洲av成人av| 精品久久久精品久久久| 国产成人精品无人区| 黄色成人免费大全| 十八禁人妻一区二区| 999久久久精品免费观看国产| 欧美色视频一区免费| 999久久久国产精品视频| 18禁黄网站禁片午夜丰满| 麻豆久久精品国产亚洲av | 97人妻天天添夜夜摸| 久久久久九九精品影院|