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

    珊瑚礁破碎帶附近波浪演化和波生流實(shí)驗(yàn)研究

    2017-05-24 08:38:42唐政江杜睿超蔣昌波
    海洋科學(xué) 2017年2期
    關(guān)鍵詞:坪上波浪諧波

    姚 宇, 唐政江, 杜睿超, 蔣昌波

    ?

    珊瑚礁破碎帶附近波浪演化和波生流實(shí)驗(yàn)研究

    姚 宇1, 2, 唐政江1, 杜睿超1, 蔣昌波1, 2

    (1. 長沙理工大學(xué)水利工程學(xué)院, 湖南長沙 410004; 2. 水沙科學(xué)與水災(zāi)害防治湖南省重點(diǎn)實(shí)驗(yàn)室, 湖南長沙 410004)

    為了研究珊瑚岸礁破碎帶附近波浪演化和波生流特性, 通過水槽實(shí)驗(yàn)對規(guī)則波浪作用下珊瑚岸礁上沿礁分布的水位和流速進(jìn)行了詳細(xì)的測量。在典型卷破波條件下, 測試了礁冠存在與不存在的兩種情況。實(shí)驗(yàn)結(jié)果表明, 多重波浪反射作用引起岸礁上形成不完全駐波, 而破碎帶附近的淺化作用則產(chǎn)生高次諧波, 波浪破碎所耗散的波能主要來源于主頻波, 礁坪上透射波成分中二次諧波與主頻波的能量相當(dāng); 礁冠的存在引起破碎帶寬度減小、礁坪上增水變大以及礁坪上各次諧波變小; 礁冠不存在時, 岸礁上波生流的沿礁分布與平直海岸相似, 而礁冠的存在一定程度上阻礙了礁坪上水體向外海的回流。研究成果將豐富和發(fā)展珊瑚礁水動力學(xué)理論, 并為島礁工程的建設(shè)和維護(hù)提供一定的理論參考。

    高次諧波; 波生流; 破碎帶; 珊瑚岸礁

    珊瑚礁(coral reef)是由碳酸鈣組成的珊瑚蟲骨骼沉積形成的巖體, 廣泛分布于熱帶和亞熱帶淺海地區(qū)。我國擁有豐富的珊瑚礁資源, 主要分布于西沙、南沙、臺灣島和海南島沿岸。近年來, 我國在南海地區(qū)珊瑚島礁周圍進(jìn)行填礁造陸的活動日益增多, 建設(shè)了諸如機(jī)場、碼頭、燈塔、通信、氣象等相關(guān)設(shè)施。這些設(shè)施在復(fù)雜的遠(yuǎn)海海洋動力環(huán)境下的穩(wěn)定性問題日益受到學(xué)者們的關(guān)注, 其中由于臺風(fēng)暴潮引起的波浪作用是最重要的海洋動力因素之一, 因此研究珊瑚礁附近與波浪相關(guān)的水動力問題對于珊瑚島礁工程設(shè)施的建設(shè)和維護(hù)具有指導(dǎo)意義。

    岸礁(fringing reef)是珊瑚礁的常見類型, 主要由連接深海海床的礁前斜坡和延伸向海岸的水平礁坪組成, 在礁前斜坡和礁坪相接的礁緣處還有可能存在一個隆起的礁冠。與平直海岸類似, 波浪由深海傳至岸礁, 由于礁前斜坡上水深急劇變淺, 波浪因淺化作用而變陡, 進(jìn)而在礁緣附近發(fā)生破碎, 破碎過程如圖1所示, 破碎波在礁坪上傳播一段距離消耗大量的能量后重新生成行進(jìn)波, 因此礁坪上并不存在類似于平直海岸的沖泄區(qū)。同時波浪破碎也會引起礁坪上產(chǎn)生波浪增水和波生流[1-2], 并且部分波生流以海底回流的形式返回到外?!,F(xiàn)場觀測也發(fā)現(xiàn)伴隨著礁前斜坡和礁緣附近的波浪淺化和破碎作用, 礁坪上波譜顯著變寬, 產(chǎn)成高頻波[3-4]和低頻波[5]。

    國內(nèi)外文獻(xiàn)中對于岸礁上波浪演化的物理模型實(shí)驗(yàn)研究主要集中于波浪傳播變形、破碎以及增水方面的問題, 例如: 張慶河等[6]對規(guī)則波在臺階地形上的破碎規(guī)律進(jìn)行了實(shí)驗(yàn)研究。Yao等[7]用實(shí)驗(yàn)和數(shù)值方法研究了岸礁對于波浪傳播的影響,驗(yàn)證了礁頂水深是控制波浪破碎和增水產(chǎn)生的決定性參數(shù)。Yao等[8]采用物理模型系統(tǒng)地研究了破波類型、破波位置和破碎帶寬度等破碎波的特征隨礁盤水深和礁前斜坡坡度的變化規(guī)律, 并給出了相關(guān)的經(jīng)驗(yàn)公式。梅弢和高峰[9]實(shí)驗(yàn)?zāi)M了常年平均波浪和重現(xiàn)期為50 a的兩種波浪條件, 研究了波浪在珊瑚礁上的傳播規(guī)律。柳淑學(xué)等[10]通過物理模型實(shí)驗(yàn)對規(guī)則波和不規(guī)則波在珊瑚礁作用下的波浪破碎及波高變化進(jìn)行了研究, 并給出波浪破碎指標(biāo)。姚宇等[11-12]通過物理模型實(shí)驗(yàn)分別研究了珊瑚礁破碎帶附近的波浪反射、透射以及能量衰減規(guī)律和礁冠的存在對于珊瑚礁上波浪傳播變形以及增水產(chǎn)生的影響?,F(xiàn)有的實(shí)驗(yàn)研究幾乎沒有涉及破碎帶附近波浪的詳細(xì)演化過程, 對波生流問題也研究甚少。Gourlay[13]通過一系列物理模型實(shí)驗(yàn)首次研究了規(guī)則波作用下礁坪上的波生流問題, 然而其僅在礁坪上單點(diǎn)位置對岸礁上的波生流進(jìn)行測量, 無法分析流速的垂向分布和沿礁變化。

    因此, 本文通過波浪水槽物理模型實(shí)驗(yàn), 對珊瑚岸礁破碎帶附近沿礁分布的波和流在空間上進(jìn)行詳細(xì)的測量, 以此研究波浪演化和波生流特性。具體來說, 以某一典型的卷破波工況為例, 分別探討兩種礁形(存在和不存在礁冠)破碎帶附近波高的變化、高頻波的產(chǎn)生和波生流的分布。研究成果將豐富和發(fā)展珊瑚礁水動力學(xué)理論, 為島礁工程的建設(shè)和維護(hù)提供一定的理論參考。

    1 實(shí)驗(yàn)設(shè)置

    物理模型實(shí)驗(yàn)在新加坡南洋理工大學(xué)的長36 m, 寬0.55 m, 高0.6 m的波浪水槽中進(jìn)行, 實(shí)驗(yàn)設(shè)置如圖2a所示, 水槽左端為一臺具備主動吸收功能的推板式造波機(jī), 另一端距造波機(jī)約32 m處設(shè)置坡度為1︰8的斜面模擬礁后岸灘。物理模型按照弗洛伊德相似定律, 根據(jù)現(xiàn)場觀測的原型礁冠[14]按 1︰20制作, 對應(yīng)的時間比尺為1︰4.5。在距造波機(jī)16.35 m處設(shè)置坡度為1︰6的斜坡模擬礁前斜面, 斜面后接長度為7 m的水平平臺模擬礁坪。礁冠寬50 cm, 高5 cm, 長55 cm, 長度與水槽寬度相匹配, 礁冠固定于礁坪平臺上, 與礁坪邊緣對齊。實(shí)驗(yàn)中沿礁布置12個浪高儀(G1~G12)來測量自由液面的變化。對于礁冠存在和不存在的情況, 實(shí)驗(yàn)均測試了礁坪水深范圍為0.05~0.10 m、波浪周期范圍為1~2 s、深水波高范圍為0.04~0.12 m的組合規(guī)則波工況, 據(jù)上述模型比尺對應(yīng)的原型波浪條件為礁坪水深1~2 m、波浪周期4.5~9 s、深水波高0.8~2.4 m, 符合Hench等[14]的現(xiàn)場觀測情況, 詳細(xì)的實(shí)驗(yàn)設(shè)置見姚宇等[12]。

    1.1 波浪的測量

    實(shí)驗(yàn)中觀察到, 當(dāng)水深增大時, 破碎點(diǎn)由礁前斜坡移向礁緣, 大部分波況下波浪以卷破的形式破碎。破碎帶寬度定義為從破碎帶起點(diǎn)(破碎點(diǎn))到破碎帶終點(diǎn)(破碎停止并重新生成行進(jìn)波)之間的水平距離[11], 由內(nèi)破碎帶和外破碎帶組成, 如圖1所示。為研究波流的沿礁分布, 本文選擇一個代表性波況(遠(yuǎn)海靜水深0=0.45 m, 深水波高0=0.095 m, 波周期=1.25 s)沿礁進(jìn)行詳細(xì)的測量, 此工況下波浪破碎為典型的卷破波, 所觀測到的破碎帶的特征位置如表1所示。在此工況下選擇破碎帶附近60個測點(diǎn)位置對礁冠存在和不存在的情況均進(jìn)行詳細(xì)的波浪測量, 如圖2b所示: 測量區(qū)域始于礁前斜坡上的淺化帶(距離礁緣外海側(cè)0.75 m), 終于破碎帶后再生波帶(距離礁緣向岸側(cè)3.675 m), 相鄰測點(diǎn)間距75 mm。實(shí)驗(yàn)中通過在同次運(yùn)行中沿水槽移動浪高儀G4對60個位置的自由液面依次進(jìn)行測量, 每個測點(diǎn)的采樣時間為3 min, 浪高儀采樣頻率為50 Hz。由于卷破波沖擊點(diǎn)附近存在強(qiáng)烈的摻氣作用, 測量誤差較大, 因此沖擊點(diǎn)附近的測點(diǎn)在本文中不作分析。

    1.2 水流的測量

    本文采用二維電磁流速儀測量波浪作用下沿礁水流的分布, 流速儀直徑為5 mm, L形探頭正對波浪傳播方向放置以測量瞬時水平流速和垂向流速。

    測量時須保證流速儀探頭始終位于水面以下, 故波谷以上的流速未能測量。如圖2b所示, 實(shí)驗(yàn)中沿水槽中軸線選取13條垂向測線(1~3位于淺化帶,4~11位于破碎帶,12~13位于破碎帶后再生波帶, 具體位置見表2), 詳細(xì)測量每條垂線上槽底上方到波谷下方之間每隔5 mm的流速分布。破碎帶內(nèi)測線的具體位置根據(jù)有無礁冠兩種情形時破碎點(diǎn)、沖擊點(diǎn)、濺波點(diǎn)位置以及破碎帶寬度的差異(見表1)略有不同。由于礁緣或礁冠附近的流速變化較大, 對此區(qū)域測線進(jìn)行了局部加密, 同樣由于卷破波的沖擊點(diǎn)和濺波點(diǎn)之間的區(qū)域存在強(qiáng)烈的摻氣作用, 故此區(qū)域不進(jìn)行測量。

    表1 觀測的破碎帶特征位置(規(guī)則波工況h0=0.45 m, H0=0.095 m, T=1.25 s, 圖3、圖4同)

    a: 原點(diǎn)為礁緣位置, 礁緣向岸側(cè)為正向(表2同)

    實(shí)驗(yàn)中電磁流速儀自造波機(jī)啟動10 min后開始采樣, 以保證波浪和水流達(dá)到穩(wěn)定狀態(tài), 采樣頻率為30 Hz。每條垂線上的流速分布在同一次運(yùn)行中進(jìn)行測量, 而后在新的運(yùn)行中測量下一條垂線, 各采樣點(diǎn)的采樣時間設(shè)為5 min, 相鄰兩次運(yùn)行至少間隔30 min以消除可能存在的余流影響。

    2 結(jié)果分析

    2.1 波浪演化

    類似于臺階式地形(如潛堤、沙壩等), 波浪作用于岸礁的礁前斜坡和礁冠時會產(chǎn)生高次諧波, 以往的物理實(shí)驗(yàn)和數(shù)值模擬研究(如Brossard等[15])表明, 超過60%的入射波能可通過高頻波傳遞到結(jié)構(gòu)物下游。本文擬通過傅里葉變換對測得的自由液面時間序列進(jìn)行諧波分析來研究能量從主頻波到高頻波的轉(zhuǎn)移。圖3和圖4分別展示了代表性工況0=0.45 m,0=0.095 m,=1.25 s下礁冠不存在與存在時總波幅(t)、各次(精確到五次)諧波波幅(t~5)以及平均水位(MWL)的沿礁變化, 圖中MWL通過對所測各點(diǎn)自由液面的時間序列求平均獲得。

    表2 流速測量位置

    圖3a和圖4a中波幅t的沿礁變化表明: 波浪破碎點(diǎn)均位于礁緣離岸一側(cè)(此時波幅達(dá)到最大); 在礁前斜坡上, 由于波浪反射, 波幅在空間上存在起伏(形成了不完全駐波), 當(dāng)礁冠存在時波浪反射更強(qiáng), 故駐波現(xiàn)象更為明顯; 相較于礁冠不存在的情況, 礁冠的存在使局部水深變淺, 破碎帶變窄, 隨之造成波幅沿礁衰減加快; 在礁坪上破碎帶后的再生波帶, 雖然兩種礁形下礁坪靜水深相同(h=0.1 m), 但礁冠存在時再生波的波幅仍略小; 由于礁坪兩端和礁后岸灘造成的多重波浪反射的影響, 礁坪上的波幅亦出現(xiàn)空間上的起伏。理論上駐波波長為相同條件下行進(jìn)波波長的一半, 對于上述工況, 礁坪上淺水透射波波長, 則礁坪上不完全駐波的波長估算為0.63 m, 與圖中觀測到的駐波波長(礁坪上起伏的波幅相鄰峰值間的距離)比較接近。

    圖3a和圖4a中各次諧波振幅(1~5)的沿礁分布表明: 在淺化帶, 由于波浪的淺化作用, 波能由主頻波向高頻波轉(zhuǎn)移, 高頻波波幅沿礁逐漸增大, 在破碎點(diǎn)附近達(dá)到峰值, 此時二次諧波波幅(2≈0.02 m)達(dá)到入射波幅(t≈0.06 m)的33%, 三次諧波波幅(3≈0.01 m)大約為t的17%, 而四到五次諧波則相對較小, 由于斜坡的強(qiáng)反射作用, 主頻波1在礁前斜坡上產(chǎn)生駐波現(xiàn)象(尤其是礁冠存在時); 在破碎帶內(nèi), 由于波浪破碎耗能, 主頻波1迅速減小, 而高次諧波2~5僅略微減小; 在破碎帶后礁坪上的再生波帶, 透射波能量中高次諧波的能量所占據(jù)了相當(dāng)?shù)牟糠? 其中二次諧波與主頻波能量接近, 而三次到五次諧波則幾乎可以忽略, 特別是第五次諧波僅表現(xiàn)為入射波幅的2%, 所以本文不再分析更高次的諧波。同時各次諧波的幅值仍和上述波高一樣出現(xiàn)空間上的波動(即駐波現(xiàn)象)。礁冠的存在對于諧波沿礁變化有一定的影響, 由于礁冠一定程度上增強(qiáng)了破碎波能量的消耗, 因此礁冠存在時礁坪上各次諧波的能量略有減小。

    圖3b和圖4b表明, 波浪在礁前斜坡由于變淺作用造成MWL下降產(chǎn)生減水(即MWL小于靜水位), 其直至波浪破碎點(diǎn)達(dá)到最低值(此為破碎帶起點(diǎn)), 此時減水值分別是–0.7 cm(無礁冠)和–2.3 cm(有礁冠), 隨后MWL開始沿破碎帶上升, 直到破碎停止重新生成行進(jìn)波時(此為破碎帶終點(diǎn)), MWL表現(xiàn)為增水(MWL大于靜水位0), 這與上述波幅的變化趨勢恰好相反, 符合Longuet-Higgins和Stewart[16]的輻射應(yīng)力理論的描述; 隨后由于破碎過程結(jié)束, 在礁坪上的再生波帶增水趨于常值, 其中無礁冠時波浪破碎導(dǎo)致的增水為0.6 cm, 對應(yīng)的相對增水值(礁坪增水最大值與礁坪靜水深的比值)為6%, 有礁冠時的增水為1.7 cm, 對應(yīng)的相對增水值為17%; 值得注意的是, MWL的沿礁分布以及增水的產(chǎn)生并未受到駐波的影響而出現(xiàn)起伏。對比有無礁冠兩種礁形下MWL的變化特征, 發(fā)現(xiàn)當(dāng)礁冠存在時破碎點(diǎn)附近平均水位的變化同樣更為劇烈。根據(jù)圖中MWL變化特征, 估算有無礁冠時的破碎帶寬度(即圖3b和圖4b中箭頭所示的破碎點(diǎn)和破碎帶終點(diǎn)之間的距離)分別為1.2 m和2.7 m, 與實(shí)驗(yàn)中的觀測值1.38 m和2.92 m (由表1計(jì)算)基本一致。

    2.2 波生流

    珊瑚岸礁波浪破碎帶內(nèi)垂向水流結(jié)構(gòu)特征對于珊瑚沙輸運(yùn)等海岸動力過程具有重要影響。本文通過對沿波浪傳播方向的水流在垂向上的分布進(jìn)行詳細(xì)的測量, 來分析岸礁上波生流的變化規(guī)律。以選取的工況為例, 圖5對比了礁冠存在與不存在時若干有代表性的垂線上中間水深處水平沿礁流速(調(diào)整到相位一致)的時間序列。在淺化帶(1), 實(shí)測流速呈現(xiàn)明顯的周期性變化, 且振幅在有無礁冠兩種情況下基本相同; 在外破碎帶的破碎波沖擊點(diǎn)附近(5), 瞬時流速呈鋸齒狀不對稱分布, 具有典型的非線性波特征; 在內(nèi)破碎帶的濺波點(diǎn)附近(9), 流速則不甚規(guī)則, 一方面是由于波浪破碎引起的紊流和摻氣作用對流速測量產(chǎn)生了影響, 另一方面是由于高頻波的產(chǎn)生; 相較于1和5, 由于礁冠的阻水作用,9位置在礁冠存在時的流速幅值較礁冠不存在時小; 在破碎帶后的再生波帶(13), 流速變化傾向有序, 由于高頻波的影響存在次波峰, 而流速幅值幾乎與9處相同。

    圖6和圖7分別展示了礁冠不存在和存在時沿礁平均水平流速()隨水深的變化規(guī)律。本文通過集合平均法求得平均流速, 即先對每個采樣點(diǎn)所測的流速在幾十個波周期內(nèi)取相位平均, 然后求一個波周期內(nèi)相位平均流速的平均值得到總時均流速, 圖6和圖7中采用誤差條表示各測點(diǎn)時均流速的測量誤差, 誤差通過相位平均流速的標(biāo)準(zhǔn)差計(jì)算。測點(diǎn)的位置由離礁緣的距離(礁緣向岸側(cè)為正向)和相對高程/確定(為測點(diǎn)高程,=0為靜水面,為當(dāng)?shù)仂o水深), 時均流速由當(dāng)?shù)販\水波波速進(jìn)行無量綱化。由于破碎波的摻氣作用, 在濺波點(diǎn)附近流速測量誤差較大(如圖6和圖7中的6~9)。

    圖6表明, 當(dāng)礁冠不存在時, 淺化帶(1~3)平均水流流向離岸方向, 即海底回流(波谷以上未能測到區(qū)域的平均水流為向岸流), 且隨著水深增大呈線性增大, 這是非破碎波的典型特征; 在外破碎帶的沖擊點(diǎn)離岸側(cè)(4~6), 流速呈非線性分布, 床面附近接近于0, 在破碎波波谷附近最大, 其中在5位置, 波谷附近的最大流速達(dá)到0.15; 在內(nèi)破碎帶(濺波區(qū)向岸側(cè))(6~11), 破碎波完全發(fā)展為涌波(bore), 大部分測點(diǎn)的流速為負(fù)值, 即流向離岸方向, 僅在7和8處的波谷附近流速為正值, 即流向海岸方向。實(shí)驗(yàn)中觀察到, 從6到11波浪破碎強(qiáng)度逐漸減小, 涌波尺度逐漸減小, 由于向岸的質(zhì)量輸移主要由行進(jìn)的涌波產(chǎn)生, 故觀察到的平均流的強(qiáng)度也在沿礁減小; 在破碎帶后的再生波帶(12和13), 垂向平均流速基本為零, 表明破碎帶后近岸區(qū)的波生流幾乎可以忽略不計(jì)。

    圖7表明, 當(dāng)礁冠存在時,1和2處垂向水流分布與礁冠不存在時相似; 但在3處, 波浪臨近破碎, 礁冠附近的垂向水流結(jié)構(gòu)明顯改變; 在礁冠頂部(4~6), 同樣存在離岸方向的平均流(海底回流), 但強(qiáng)度約為無礁冠時的兩倍(0.3), 最大海底回流出現(xiàn)在礁緣附近, 約為0.45; 在濺波點(diǎn)附近(7), 流速分布與6處相似, 但是幅值明顯較小; 在內(nèi)破碎帶(8~11), 波生流在底床附近流向離岸方向, 波谷附近則流向海岸方向, 同時強(qiáng)度逐漸減小; 在破碎帶后的再生波帶(12和13), 海底回流的強(qiáng)度幾乎可以忽略。綜上所述, 相較于無礁冠的情況, 當(dāng)礁冠存在時海底回流在礁冠附近明顯較大, 而在礁冠后隨著破碎帶的變窄而出現(xiàn)海底回流的范圍變小。這表明, 類似于明渠流中寬頂堰的阻水作用, 礁冠能夠阻礙近岸側(cè)水體向外海側(cè)的回流, 從而影響礁坪和遠(yuǎn)海之間長期的水體交換, 這對礁坪上珊瑚沙的運(yùn)動會產(chǎn)生重要的影響。

    值得注意的是, 本文測得的礁冠不存在時典型卷破波作用下沿礁的波生流垂向分布規(guī)律與已有文獻(xiàn)[17]中報(bào)道的在平直海岸情況下的測量結(jié)果比較一致, 表明典型的一維珊瑚岸礁和平直海岸上波浪破碎、增水和波生流的產(chǎn)生機(jī)理具有一致性。由于實(shí)驗(yàn)中的波浪水槽為封閉容器, 向岸方向的質(zhì)量輸移流必須由離岸方向的海底回流補(bǔ)償, 以保證水槽任一截面凈質(zhì)量輸移為零。實(shí)際的珊瑚礁系統(tǒng)可能在沿岸方向分布并不均勻, 往往存在著某些裂口, 近似于沙質(zhì)海岸的裂流系統(tǒng), 從而形成近岸水平環(huán)流。當(dāng)此類環(huán)流存在時, 向岸流主要是通過裂流回到外海而海底回流會顯著減小, 甚至有實(shí)驗(yàn)室測量中發(fā)現(xiàn), 垂向水流僅表現(xiàn)為向岸流而不存在海底回流的情況[18], 因此本文的結(jié)論適用于沿岸方向分布較均勻的大片岸礁(不存在裂口)或遠(yuǎn)離裂口的區(qū)域。

    3 結(jié)論

    本文通過物理模型實(shí)驗(yàn)對岸礁破碎帶附近沿礁分布的波浪和水流進(jìn)行了詳細(xì)的測量, 實(shí)驗(yàn)在典型卷破波條件下測試了礁冠存在與不存在時的兩種礁形。分析表明, 在珊瑚岸礁上由于存在多重的波浪反射作用形成了不完全駐波, 但沿礁平均水位以及增水的產(chǎn)生并未受到駐波的影響; 破碎帶附近由于淺化作用產(chǎn)生高次諧波, 波浪破碎所耗散的能量主要由主頻波提供, 造成礁坪上透射波的組成中二次諧波與主頻波的能量相當(dāng); 礁冠的存在使破碎帶寬度減小、礁坪上增水變大以及礁坪上各次諧波變小; 在典型的卷破波條件下, 礁冠不存在時岸礁上波生流的沿礁分布特征與已有文獻(xiàn)中報(bào)道的平直海岸的情況相似, 而礁冠存在的作用類似于明渠中寬頂堰, 一定程度上阻礙了靠岸側(cè)水體向離岸側(cè)的回流, 進(jìn)而影響礁坪和遠(yuǎn)海之間的水體交換。研究成果有助于進(jìn)一步研究珊瑚礁附近的水動力過程和珊瑚沙的輸運(yùn)機(jī)理, 并為島礁工程的建設(shè)和維護(hù)提供一定的理論參考, 也可用于校核與波浪破碎相關(guān)的數(shù)值模型。

    [1] Vetter O, Becker J M, Merrifield M A, et al. Wave setup over a Pacific Island fringing reef [J]. Journal of Geophysical Research, 2010, 115: C12066.

    [2] Jago O K, Kench P S, Brander R W. Field observations of wave-driven water-level gradients across a coral reef flat [J]. Journal of Geophysical Research, 2007, 112: C06027.

    [3] Brander R W, Kench P S, Hart D E. Spatial and temporal variations in wave characteristics across a reef platform, Warraber Island, Torres Strait, Australia [J]. Marine Geology, 2004, 207(1-4): 169-184.

    [4] Lowe R J, Falter J L, Bandet M D, et al. Spectral wave dissipation over a barrier reef [J]. Journal of Geophysical Research, 2005, 110: C04001.

    [5] Yao Y, Becker J M, Ford M R, et al. Modeling wave processes over fringing reefs with an excavation pit [J]. Coastal Engineering, 2016, 109: 9-19.

    [6] 張慶河, 劉海青, 趙子丹. 波浪在臺階地形上的破碎[J]. 天津大學(xué)學(xué)報(bào), 1999, 32(2): 204-207. ZHANG Qinghe, LIU Haiqing, ZHAO Zidan. Wave breaking on a submerged step [J]. Journal of Tianjin University, 1999, 32(2): 204-207.

    [7] Yao Y, Huang Z, Monismith S G, et al. 1DH Boussinesq modeling of wave transformation over fringing reefs [J]. Ocean Engineering, 2012, 47: 30-42.

    [8] Yao Y, Huang Z, Monismith S G, et al. Characteristics of monochromatic waves breaking over fringing reefs [J]. Journal of Coastal Research, 2013, 29(1): 94-104.

    [9] 梅弢, 高峰. 波浪在珊瑚礁坪上傳播的水槽試驗(yàn)研究[J]. 水道港口, 2013, 34(1): 13-18. MEI Tao, GAO Feng. Flume experiment research on law of wave propagation in reef flat [J]. Journal of Waterway and Harbor, 2013, 34(1): 13-18.

    [10] 柳淑學(xué), 劉寧, 李金宣, 等. 波浪在珊瑚礁地形上破碎特性試驗(yàn)研究[J]. 海洋工程, 2015, 33(2): 42-49. LIU Shuxue, LIU Ning, LI Jinxuan, et al. Experimental researches on wave propagation characteristics on reefs terrain [J]. The Ocean Engineering, 2015, 33(2): 42-49.

    [11] 姚宇, 杜睿超, 袁萬成, 等. 珊瑚岸礁破碎帶附近波浪演化實(shí)驗(yàn)研究[J]. 海洋學(xué)報(bào), 2015, 37(12): 66-73. YAO Yu, DU Ruichao, YUAN Wancheng, et al. Experimental study of wave transformation around the surf zone over fringing reefs [J]. Haiyang Xuebao, 2015, 37(12): 66-73.

    [12] 姚宇, 袁萬成, 杜睿超, 等. 岸礁礁冠對波浪傳播變形及增水影響的實(shí)驗(yàn)研究[J]. 熱帶海洋學(xué)報(bào), 2015, 34(6): 19-25. YAO Yu, YUAN Wancheng, DU Ruichao, et al. Experimental study of reef crest’s effects on wave transformation and wave-induced setup over fringing reefs [J]. Journal of Tropical Oceanography, 2015, 34(6): 19-25.

    [13] Gourlay M R. Wave set-up on coral reefs. 1. Set-up and wave-generated flow on an idealised two-dimensional reef [J]. Coastal Engineering, 1996, 27: 161-193.

    [14] Hench J L, Leichter J J, Monismith S G. Episodic circulation and exchange in a wave-driven coral reef and lagoon system [J]. Limnology and Oceanography, 2008, 53: 2681-2694.

    [15] Brossard J, Perret G, Blonce L, et al. Higher harmonics induced by a submerged horizontal plate and a submerged rectangular step in a wave flume [J]. Coastal Engineering, 2009, 56: 11-22.

    [16] Longuet-Higgins M S, Stewart R W. Radiation stress and mass transport in gravity waves, with application to surf beats [J]. Journal of Fluid Mechanics, 1962, 13: 481-504.

    [17] Ting F C K., Kirby J T. Observation of undertow and turbulence in a laboratory surf zone [J]. Coastal Engineering, 1994, 24: 51-80.

    [18] Dr?nen N, Karunarathna H, Freds?e J, et al. An experimental study of rip channel flow [J]. Coastal Engineering, 2002, 45: 223-238.

    Experimental study of wave transformation and wave-driven current around the surf zone over coral reefs

    YAO Yu1, 2, TANG Zheng-jiang1, DU Rui-chao1, JIANG Chang-bo1, 2

    (1. School of Hydraulic Engineering, Changsha University of Science & Technology, Changsha 410004, China; 2. Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province, Changsha 410004, China)

    To investigate the wave transformation and wave-driven current around the surf zone over fringing reefs, a series of laboratory experiments are performed in a wave flume to measure the cross-reef water lever and flow in details. Experimental results are reported for two reef profiles with and without the reef crest under a typical plunging breaker condition. Data analysis shows that partial standing waves across the reef are formed due to multiple wave reflection, and higher harmonics around the surf zone are generated through wave shoaling. The energy dissipation in the surf zone results primarily from the dissipation of fundamental harmonics, and the resultant second harmonic waves are comparable to the fundamental harmonics on the reef flat. The presence of the reef crest reduces the surf zone width, increases the wave-induced setup, and decreases the harmonic wave amplitudes on the reef flat. The variation of the wave-driven current across the fringing reef without the reef crest is similar to that for the plane beaches. The existence of the reef crest causes the retention of the mean-flow on the leeside of the crest to some extent.

    higher harmonic waves; wave-driven current; surf zone; fringing reef

    TV139.2

    A

    1000-3096(2017)02-0012-08

    10.11759/hykx20160313001

    2016-03-13;

    2016-05-06

    國家自然科學(xué)基金資助項(xiàng)目(51309035,51239001);水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室開放研究基金(2014491011);高等學(xué)校博士學(xué)科點(diǎn)科研基金新教師類(20134316120004)

    姚宇 (1982-), 男, 湖南湘潭人, 副教授, 博士, 主要從事近海水動力學(xué)和環(huán)境流體力學(xué)研究, E-mail: yaoyu821101@163.com

    Mar. 13, 2016

    [National Natural Science Foundation of China, No.51309035, No.51239001; the Open Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, No.2014491011; the Specialized Research Fund for the Doctoral Program of Higher Education, No.20134316120004]

    (本文編輯: 李曉燕)

    猜你喜歡
    坪上波浪諧波
    在路上(外一首)
    詩歌月刊(2023年3期)2023-05-30 14:51:08
    少年毛澤東曬谷坪上的風(fēng)波
    波浪谷和波浪巖
    梅映青山 小鎮(zhèn)悠然
    中國公路(2022年4期)2022-04-13 05:17:34
    波浪谷隨想
    去看神奇波浪谷
    小饞嘴
    虛擬諧波阻抗的并網(wǎng)逆變器諧波抑制方法
    基于ELM的電力系統(tǒng)諧波阻抗估計(jì)
    基于ICA和MI的諧波源識別研究
    電測與儀表(2015年8期)2015-04-09 11:50:00
    最黄视频免费看| 亚洲情色 制服丝袜| 午夜91福利影院| 一级av片app| 国产高清三级在线| 午夜福利在线观看免费完整高清在| 日韩精品有码人妻一区| 校园人妻丝袜中文字幕| 亚洲国产成人一精品久久久| 久久精品国产亚洲av涩爱| 91久久精品电影网| 亚州av有码| 久久午夜福利片| 欧美丝袜亚洲另类| 精品人妻熟女毛片av久久网站| 亚洲精品视频女| 在线观看三级黄色| videossex国产| 国产欧美日韩一区二区三区在线 | 一个人看视频在线观看www免费| 精品国产露脸久久av麻豆| 国产永久视频网站| 五月伊人婷婷丁香| 在线观看三级黄色| 91精品伊人久久大香线蕉| 久久午夜福利片| 一本一本综合久久| 久久久久国产精品人妻一区二区| 日本黄色片子视频| 精品亚洲乱码少妇综合久久| 亚洲第一av免费看| 18+在线观看网站| 人妻一区二区av| 久久女婷五月综合色啪小说| 中文字幕人妻丝袜制服| 精品人妻偷拍中文字幕| 美女主播在线视频| 黑人高潮一二区| 精品一区在线观看国产| 久久97久久精品| 在线 av 中文字幕| 赤兔流量卡办理| 边亲边吃奶的免费视频| 男女国产视频网站| 国产精品久久久久久久久免| h日本视频在线播放| 国产一区二区在线观看日韩| 精品国产国语对白av| 丝瓜视频免费看黄片| 性高湖久久久久久久久免费观看| 搡女人真爽免费视频火全软件| 一区二区三区四区激情视频| 亚洲精品乱码久久久v下载方式| 久久久国产精品麻豆| 少妇丰满av| 久久97久久精品| 亚洲av不卡在线观看| 不卡视频在线观看欧美| 日本wwww免费看| 卡戴珊不雅视频在线播放| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 99热国产这里只有精品6| 国产高清有码在线观看视频| 国产精品伦人一区二区| 成人免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 日韩亚洲欧美综合| 国产又色又爽无遮挡免| 欧美精品一区二区免费开放| 91在线精品国自产拍蜜月| 国产一级毛片在线| 亚洲欧美日韩东京热| 3wmmmm亚洲av在线观看| av国产久精品久网站免费入址| 国产乱人偷精品视频| 免费av不卡在线播放| 99久久精品热视频| 美女内射精品一级片tv| 精品亚洲成国产av| 伊人亚洲综合成人网| 欧美少妇被猛烈插入视频| 最近的中文字幕免费完整| 午夜视频国产福利| 午夜日本视频在线| 亚洲欧美日韩卡通动漫| 高清在线视频一区二区三区| 精品国产一区二区三区久久久樱花| 国产成人91sexporn| 亚洲va在线va天堂va国产| 又黄又爽又刺激的免费视频.| 免费观看在线日韩| 搡老乐熟女国产| a级毛色黄片| 久久久久久久久久成人| 一区二区av电影网| 好男人视频免费观看在线| 亚洲国产精品国产精品| 又大又黄又爽视频免费| 国语对白做爰xxxⅹ性视频网站| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区| 亚洲精品国产成人久久av| 久久人人爽人人片av| 日韩成人伦理影院| 亚洲怡红院男人天堂| 午夜91福利影院| 国产精品麻豆人妻色哟哟久久| 大话2 男鬼变身卡| 一级二级三级毛片免费看| 欧美丝袜亚洲另类| 夜夜看夜夜爽夜夜摸| 精品午夜福利在线看| 中文乱码字字幕精品一区二区三区| 一级毛片aaaaaa免费看小| 汤姆久久久久久久影院中文字幕| 18禁在线无遮挡免费观看视频| 国产日韩欧美亚洲二区| 久久99一区二区三区| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| www.av在线官网国产| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 日韩精品免费视频一区二区三区 | 久热这里只有精品99| 婷婷色麻豆天堂久久| 色视频在线一区二区三区| 99热这里只有精品一区| 天天操日日干夜夜撸| 欧美少妇被猛烈插入视频| 色视频www国产| 大香蕉久久网| 一边亲一边摸免费视频| 美女视频免费永久观看网站| 精品熟女少妇av免费看| 午夜免费鲁丝| 国产日韩欧美在线精品| 街头女战士在线观看网站| 亚洲色图综合在线观看| 精品亚洲成国产av| 国产精品一区二区在线观看99| 街头女战士在线观看网站| 女人精品久久久久毛片| 精品一区在线观看国产| 中文字幕精品免费在线观看视频 | 国产男女内射视频| 亚洲高清免费不卡视频| 国产精品久久久久久精品古装| 久久婷婷青草| 欧美日韩一区二区视频在线观看视频在线| 免费观看a级毛片全部| 黑人猛操日本美女一级片| 2021少妇久久久久久久久久久| 熟妇人妻不卡中文字幕| 久久鲁丝午夜福利片| 国产精品成人在线| 亚洲精品国产成人久久av| 成人二区视频| 国产熟女欧美一区二区| 青青草视频在线视频观看| 国产在线一区二区三区精| 51午夜福利影视在线观看| 国产高清videossex| 免费在线观看影片大全网站| 欧美在线一区亚洲| 亚洲 欧美一区二区三区| 国产亚洲av片在线观看秒播厂| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频| 一个人免费看片子| 成年人黄色毛片网站| 国产成人精品久久二区二区免费| 日韩 欧美 亚洲 中文字幕| 久久性视频一级片| 不卡av一区二区三区| 久久天堂一区二区三区四区| 一本色道久久久久久精品综合| av免费在线观看网站| 国产欧美亚洲国产| 欧美中文综合在线视频| 黑人巨大精品欧美一区二区蜜桃| 最近最新免费中文字幕在线| 搡老熟女国产l中国老女人| 欧美黑人精品巨大| 精品久久久久久电影网| 妹子高潮喷水视频| 亚洲一码二码三码区别大吗| 亚洲中文日韩欧美视频| 国产91精品成人一区二区三区 | netflix在线观看网站| 久久精品aⅴ一区二区三区四区| 午夜精品久久久久久毛片777| 国产精品久久久久成人av| 精品乱码久久久久久99久播| 成年女人毛片免费观看观看9 | 天堂8中文在线网| 国产亚洲精品第一综合不卡| 日韩大码丰满熟妇| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区三 | 国产亚洲一区二区精品| 成年美女黄网站色视频大全免费| 黑人欧美特级aaaaaa片| 精品一品国产午夜福利视频| 免费在线观看黄色视频的| 国产精品久久久久久精品电影小说| 日本黄色日本黄色录像| 黑人猛操日本美女一级片| 好男人电影高清在线观看| 久久久久久久久久久久大奶| 久久精品国产综合久久久| av视频免费观看在线观看| 最黄视频免费看| 国产精品久久久久久精品电影小说| 在线天堂中文资源库| 成人三级做爰电影| 97在线人人人人妻| 一区在线观看完整版| 女性被躁到高潮视频| 国产在线视频一区二区| 黑人猛操日本美女一级片| 黑人巨大精品欧美一区二区蜜桃| 日本av免费视频播放| 国产在视频线精品| 99国产精品免费福利视频| 各种免费的搞黄视频| 国产一级毛片在线| a级毛片在线看网站| 亚洲avbb在线观看| 自线自在国产av| 中亚洲国语对白在线视频| 大片免费播放器 马上看| 1024香蕉在线观看| 亚洲精品国产一区二区精华液| 18禁观看日本| 亚洲色图 男人天堂 中文字幕| 热99国产精品久久久久久7| 日韩三级视频一区二区三区| 国产精品 欧美亚洲| av线在线观看网站| 国产福利在线免费观看视频| 男女高潮啪啪啪动态图| 人人妻,人人澡人人爽秒播| 免费观看av网站的网址| 一级片'在线观看视频| 999久久久精品免费观看国产| 人人妻人人澡人人爽人人夜夜| 久久ye,这里只有精品| av在线app专区| 大香蕉久久成人网| 免费在线观看视频国产中文字幕亚洲 | av超薄肉色丝袜交足视频| 99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| av一本久久久久| 国产精品影院久久| 欧美人与性动交α欧美精品济南到| av不卡在线播放| 国产一区二区三区综合在线观看| 正在播放国产对白刺激| 搡老岳熟女国产| 99久久国产精品久久久| 亚洲三区欧美一区| 亚洲欧美精品自产自拍| 无限看片的www在线观看| 十八禁网站网址无遮挡| 国产精品一区二区免费欧美 | 午夜免费鲁丝| 亚洲第一青青草原| 男女高潮啪啪啪动态图| 一区福利在线观看| 可以免费在线观看a视频的电影网站| 国产高清视频在线播放一区 | 欧美+亚洲+日韩+国产| 视频在线观看一区二区三区| 老司机在亚洲福利影院| 亚洲欧美一区二区三区黑人| 欧美日韩国产mv在线观看视频| 男女床上黄色一级片免费看| 成人三级做爰电影| 免费少妇av软件| 久久这里只有精品19| 一区二区三区激情视频| 又大又爽又粗| 侵犯人妻中文字幕一二三四区| 俄罗斯特黄特色一大片| 久久精品国产综合久久久| 国产在线视频一区二区| 黄色怎么调成土黄色| 国产主播在线观看一区二区| 91麻豆av在线| 免费少妇av软件| 极品人妻少妇av视频| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 老司机亚洲免费影院| 黑丝袜美女国产一区| 亚洲色图综合在线观看| 亚洲国产日韩一区二区| 久久青草综合色| av视频免费观看在线观看| 亚洲第一欧美日韩一区二区三区 | 青春草视频在线免费观看| 日本wwww免费看| 亚洲国产欧美日韩在线播放| 两性夫妻黄色片| 一级a爱视频在线免费观看| 肉色欧美久久久久久久蜜桃| 精品少妇内射三级| 大香蕉久久网| 啦啦啦视频在线资源免费观看| 国产91精品成人一区二区三区 | 亚洲专区国产一区二区| 男女边摸边吃奶| 超碰成人久久| 久久精品国产亚洲av高清一级| 丝袜美足系列| 日韩,欧美,国产一区二区三区| 久久ye,这里只有精品| 青春草视频在线免费观看| 成人国语在线视频| 久久久国产欧美日韩av| 免费不卡黄色视频| 嫩草影视91久久| 香蕉国产在线看| 97精品久久久久久久久久精品| 欧美xxⅹ黑人| 老司机亚洲免费影院| 久久人妻福利社区极品人妻图片| 搡老熟女国产l中国老女人| 性高湖久久久久久久久免费观看| 中文字幕高清在线视频| 亚洲精品美女久久久久99蜜臀| 国产日韩一区二区三区精品不卡| 久久天堂一区二区三区四区| 亚洲av电影在线观看一区二区三区| av欧美777| 中文精品一卡2卡3卡4更新| 新久久久久国产一级毛片| 精品国产一区二区久久| 日本五十路高清| 国产成人av教育| 999久久久国产精品视频| 午夜精品久久久久久毛片777| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 另类精品久久| 美女主播在线视频| 伊人久久大香线蕉亚洲五| 午夜精品久久久久久毛片777| 不卡一级毛片| 国产欧美日韩一区二区精品| 午夜视频精品福利| 亚洲人成电影免费在线| 国产深夜福利视频在线观看| 两人在一起打扑克的视频| 99国产极品粉嫩在线观看| 久久精品久久久久久噜噜老黄| 婷婷成人精品国产| 精品国产乱码久久久久久小说| 婷婷色av中文字幕| 成人黄色视频免费在线看| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 老汉色∧v一级毛片| 国产一区二区三区av在线| 亚洲欧美激情在线| 一区二区三区精品91| 侵犯人妻中文字幕一二三四区| 亚洲精华国产精华精| 午夜激情av网站| 日韩人妻精品一区2区三区| 欧美日韩精品网址| 亚洲精品自拍成人| 韩国高清视频一区二区三区| 国产熟女午夜一区二区三区| 亚洲精品国产av蜜桃| 免费在线观看黄色视频的| 午夜精品国产一区二区电影| 国产成人av教育| 99久久人妻综合| 日本撒尿小便嘘嘘汇集6| 亚洲少妇的诱惑av| 欧美一级毛片孕妇| av欧美777| 欧美另类一区| 国产伦理片在线播放av一区| 一区二区三区四区激情视频| 免费高清在线观看日韩| 超碰成人久久| 亚洲精华国产精华精| 亚洲精品第二区| 亚洲精品自拍成人| 欧美日韩黄片免| 午夜精品久久久久久毛片777| 999久久久国产精品视频| 亚洲国产av新网站| 在线观看人妻少妇| 亚洲国产中文字幕在线视频| 亚洲,欧美精品.| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 亚洲色图综合在线观看| 大陆偷拍与自拍| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品大桥未久av| 夜夜骑夜夜射夜夜干| 黄频高清免费视频| 一二三四在线观看免费中文在| 午夜福利在线观看吧| 秋霞在线观看毛片| 美女主播在线视频| av在线app专区| 精品久久蜜臀av无| av免费在线观看网站| 日韩大码丰满熟妇| 女人被躁到高潮嗷嗷叫费观| 国产精品秋霞免费鲁丝片| 91成人精品电影| 国产男人的电影天堂91| 亚洲视频免费观看视频| 黄色怎么调成土黄色| 国产三级黄色录像| 午夜福利影视在线免费观看| 久久久久久久久久久久大奶| 久久久久久免费高清国产稀缺| 69精品国产乱码久久久| 午夜免费成人在线视频| 色综合欧美亚洲国产小说| 国产精品久久久人人做人人爽| 搡老乐熟女国产| www.999成人在线观看| 91精品国产国语对白视频| 久久女婷五月综合色啪小说| 老司机亚洲免费影院| 亚洲激情五月婷婷啪啪| 悠悠久久av| 亚洲欧美成人综合另类久久久| 国产野战对白在线观看| 国产精品欧美亚洲77777| 欧美黄色片欧美黄色片| 美女福利国产在线| 男女免费视频国产| 久久国产精品影院| 后天国语完整版免费观看| 久久国产精品男人的天堂亚洲| 极品人妻少妇av视频| 欧美av亚洲av综合av国产av| 伦理电影免费视频| 国产在线视频一区二区| 一本综合久久免费| 1024视频免费在线观看| 欧美 日韩 精品 国产| 三级毛片av免费| 午夜影院在线不卡| 国产日韩一区二区三区精品不卡| 国产免费现黄频在线看| 欧美激情久久久久久爽电影 | 中国国产av一级| 欧美激情极品国产一区二区三区| 又黄又粗又硬又大视频| 午夜福利在线观看吧| 老汉色∧v一级毛片| 日韩一卡2卡3卡4卡2021年| 少妇精品久久久久久久| 精品一区在线观看国产| 国产在线免费精品| 国产成+人综合+亚洲专区| 亚洲黑人精品在线| 美女视频免费永久观看网站| 少妇人妻久久综合中文| 久久久精品免费免费高清| 精品国产乱子伦一区二区三区 | 九色亚洲精品在线播放| 亚洲熟女精品中文字幕| 国产亚洲精品第一综合不卡| av一本久久久久| 19禁男女啪啪无遮挡网站| 又黄又粗又硬又大视频| 嫁个100分男人电影在线观看| 女人被躁到高潮嗷嗷叫费观| 精品久久蜜臀av无| 午夜福利,免费看| 精品少妇内射三级| 久久久久久久久久久久大奶| 亚洲国产成人一精品久久久| 自拍欧美九色日韩亚洲蝌蚪91| 精品乱码久久久久久99久播| 午夜91福利影院| 国产精品 国内视频| 成人av一区二区三区在线看 | 久久天堂一区二区三区四区| 亚洲成人手机| 电影成人av| 少妇粗大呻吟视频| 91麻豆av在线| 一二三四社区在线视频社区8| 一区二区三区乱码不卡18| 国产精品久久久久久人妻精品电影 | 老熟妇仑乱视频hdxx| 亚洲免费av在线视频| 久久毛片免费看一区二区三区| 国产成人免费无遮挡视频| 亚洲av日韩在线播放| 日韩中文字幕欧美一区二区| 极品人妻少妇av视频| 久久精品国产亚洲av香蕉五月 | 中国国产av一级| 亚洲av电影在线进入| a级毛片黄视频| 亚洲天堂av无毛| 久久久国产欧美日韩av| 国产一区二区 视频在线| 99国产综合亚洲精品| 亚洲欧美精品综合一区二区三区| 欧美日韩一级在线毛片| 男人添女人高潮全过程视频| 色视频在线一区二区三区| 少妇裸体淫交视频免费看高清 | 纵有疾风起免费观看全集完整版| 亚洲精品中文字幕在线视频| 欧美一级毛片孕妇| 国产欧美日韩精品亚洲av| 超碰成人久久| 波多野结衣一区麻豆| 搡老岳熟女国产| 精品一区二区三区av网在线观看 | 免费高清在线观看视频在线观看| 亚洲少妇的诱惑av| 91大片在线观看| 国产精品av久久久久免费| 少妇 在线观看| 亚洲国产av影院在线观看| 欧美变态另类bdsm刘玥| 少妇猛男粗大的猛烈进出视频| 伊人亚洲综合成人网| 久久久久久久国产电影| 午夜成年电影在线免费观看| 色94色欧美一区二区| 侵犯人妻中文字幕一二三四区| 国产一区二区在线观看av| 欧美午夜高清在线| 国产亚洲精品一区二区www | a级毛片在线看网站| www.自偷自拍.com| 一级毛片电影观看| 黑丝袜美女国产一区| 一边摸一边做爽爽视频免费| 国产97色在线日韩免费| 国产精品二区激情视频| 不卡av一区二区三区| 久久香蕉激情| 午夜免费观看性视频| 一区二区三区四区激情视频| 美女大奶头黄色视频| 日韩制服丝袜自拍偷拍| 搡老乐熟女国产| 欧美日韩亚洲综合一区二区三区_| 性色av一级| 亚洲国产欧美在线一区| 国产欧美日韩精品亚洲av| 美女国产高潮福利片在线看| 99久久国产精品久久久| 蜜桃在线观看..| 国产免费一区二区三区四区乱码| 免费人妻精品一区二区三区视频| 精品亚洲成国产av| 精品国产乱子伦一区二区三区 | 国产一级毛片在线| 熟女少妇亚洲综合色aaa.| 国产免费现黄频在线看| 国产精品国产三级国产专区5o| 十分钟在线观看高清视频www| 亚洲成av片中文字幕在线观看| 色婷婷久久久亚洲欧美| 一进一出抽搐动态| 国产不卡av网站在线观看| 一本大道久久a久久精品| 免费观看a级毛片全部| 久久亚洲精品不卡| 热re99久久国产66热| 麻豆av在线久日| 精品人妻一区二区三区麻豆| 国产一区二区 视频在线| 狠狠狠狠99中文字幕| 久久久久久人人人人人| 国产片内射在线| 国产精品av久久久久免费| 蜜桃国产av成人99| 法律面前人人平等表现在哪些方面 | 久久国产精品人妻蜜桃| 少妇裸体淫交视频免费看高清 | 久久精品成人免费网站| 在线观看人妻少妇| 国产99久久九九免费精品| 男女边摸边吃奶| 波多野结衣av一区二区av| 窝窝影院91人妻| 又大又爽又粗| 天天影视国产精品| 亚洲成人手机| 两性夫妻黄色片| 午夜福利乱码中文字幕| 精品久久久精品久久久| 亚洲美女黄色视频免费看| 亚洲精品久久成人aⅴ小说| 亚洲av片天天在线观看| 国产一区二区三区av在线| 日韩欧美一区视频在线观看| 一本久久精品| 久9热在线精品视频| 国产免费av片在线观看野外av| 老司机深夜福利视频在线观看 | 美女午夜性视频免费|