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

    夏季太湖北岸的波浪組成及演變特征分析

    2022-04-16 11:58:45趙巧華徐嘉王健健
    人民長江 2022年3期
    關(guān)鍵詞:波高風(fēng)浪風(fēng)場

    趙巧華 徐嘉 王健健

    摘要: 波浪是驅(qū)動大型淺水湖泊中水生生態(tài)系統(tǒng)演變的關(guān)鍵動力過程?;谔辈拷秴^(qū)的波浪方向譜、同期風(fēng)速、風(fēng)向的觀測資料,分析了波浪組成及其演變特征。結(jié)果表明:太湖北部近岸區(qū)的波浪頻率譜主要為單峰或雙峰型,屬于中高頻波;風(fēng)是波浪形成的主要驅(qū)動力,風(fēng)向(或周邊地形)對有效波高的調(diào)制作用可用三角函數(shù)定量表征;風(fēng)浪的有效波高、譜峰頻率及其平均方向與風(fēng)速、風(fēng)向的關(guān)系證實了太湖北部近岸區(qū)的波浪由涌浪和風(fēng)浪組成,且涌浪占比較高。岸界的反射作用提高了北部近岸區(qū)涌浪的比例;涌浪的有效波高隨耦合風(fēng)速(風(fēng)速、風(fēng)向的耦合)不僅呈e指數(shù)增長,而且有快、慢分化的現(xiàn)象。研究成果有助于提高大型淺水湖泊中近岸區(qū)波浪估算的精度,確保水生生態(tài)修復(fù)效應(yīng)。

    關(guān) 鍵 詞: 波浪組成; 風(fēng)浪; 涌浪; 有效波高; 太湖

    中圖法分類號: ?P642

    文獻(xiàn)標(biāo)志碼: ?A

    DOI: 10.16232/j.cnki.1001-4179.2022.03.029

    0 引 言

    波浪不僅影響著航運,而且也影響水體中的生物-化學(xué)過程。波浪通常包含不同子波浪系統(tǒng)? [1-2] :① 風(fēng)浪。它的能量主要通過水-氣共振或平行流的不穩(wěn)定性機制攝取于風(fēng)場,因而風(fēng)浪的生成、發(fā)展與風(fēng)場密切相關(guān)? [3-4] ;② 涌浪。當(dāng)風(fēng)浪離開風(fēng)的直接作用區(qū)域,或在風(fēng)的直接作用區(qū)域內(nèi)風(fēng)速突然減小或風(fēng)向驟然改變,岸界對風(fēng)浪的反射會使風(fēng)浪演變成涌浪,因而涌浪散布的區(qū)域大于風(fēng)的作用區(qū)域??梢姡L(fēng)場與風(fēng)浪關(guān)系密切,而與涌浪的關(guān)系相對疏遠(yuǎn)? [5] 。

    基于實驗室分析和野外觀測,逄勇等通過對波浪切應(yīng)力、臨界切應(yīng)力的估算,驗證了波浪與太湖懸浮物濃度、透明度的關(guān)系? [6] ,證實了波浪是大型淺水湖泊中沉積物再懸浮、營養(yǎng)鹽內(nèi)源釋放的主要驅(qū)動力? [3-4] ,進(jìn)而可能也是水華頻發(fā),久治不愈? [7] 的關(guān)鍵原因之一。為此,諸多學(xué)者針對太湖的波浪要素(波高、周期等)進(jìn)行了大量的觀測和模擬分析,揭示了在開闊湖區(qū)波浪有效波高隨風(fēng)速呈e指數(shù)增長的規(guī)律? [4,8] 。但在近岸區(qū),由于岸界的反射、湖底地形的摩擦及風(fēng)區(qū)變化等作用,風(fēng)速與有效波高的關(guān)系弱化? [8] ,使得近岸區(qū)波浪的特征和演變規(guī)律遠(yuǎn)比湖心區(qū)復(fù)雜。范成新等的研究成果顯示:在湖流的長期作用下,太湖的沉積物主要分布于近岸區(qū)域(尤其是西岸),湖心區(qū)的沉積物分布相對較薄? [9] ??梢姡沂窘秴^(qū)的波浪特征及其演變規(guī)律對水生生態(tài)修復(fù)有重要意義。迄今為止,近岸區(qū)波浪與風(fēng)速之間的定量關(guān)系并未得到足夠的關(guān)注,導(dǎo)致其生態(tài)效應(yīng)常被忽略,進(jìn)而可能弱化了水生生態(tài)修復(fù)效果。

    就大型淺水湖泊而言,依據(jù)波浪的產(chǎn)生機制和演化規(guī)律分離近岸區(qū)的波浪組成(風(fēng)浪和涌浪),揭示近岸區(qū)波浪的特征和演變規(guī)律,提高近岸區(qū)波浪參數(shù)的估算精度,有利于水生生態(tài)系統(tǒng)的精準(zhǔn)修復(fù)。由于太湖是典型的大型淺水湖泊,它對氣象場的響應(yīng)迅速,因而其中波浪通常屬于中高頻波? [10] ,與風(fēng)浪頻率較為接近,增加了分離風(fēng)浪、涌浪的難度。本文擬基于太湖近岸區(qū)的波浪方向譜,結(jié)合同期風(fēng)速、風(fēng)向數(shù)據(jù),探究波浪的組成和演變特征,為太湖水生生態(tài)修復(fù)奠定基礎(chǔ)。

    1 研究數(shù)據(jù)與方法

    1.1 研究數(shù)據(jù)

    太湖是典型的大型淺水湖泊,面積為2 338 km? 2 ,平均水深不超過2 m。波浪的觀測點位于中國科學(xué)院太湖生態(tài)系統(tǒng)試驗站的棧橋頭(見圖1)。

    本文所用資料的觀測時段為2018年8月21日17:00至9月4日05:00。

    利用MIDASDWR方向波浪儀觀測波高、波周期等參數(shù),測定頻率為8 kHz,每30 min觀測一次,觀測時間持續(xù)10 min。

    同期的風(fēng)速、風(fēng)向數(shù)據(jù)來源于中國科學(xué)院太湖生態(tài)系統(tǒng)試驗站的觀測資料。風(fēng)場的觀測頻率為逐小時,其中正北風(fēng)向為0(rad),按順勢針方向旋轉(zhuǎn),正南風(fēng)向為π(rad)。

    為了研究風(fēng)場對波浪的影響規(guī)律,選取整點的波浪方向譜數(shù)據(jù)和風(fēng)速、風(fēng)向資料。

    1.2 波浪分離方法

    波浪是局地風(fēng)浪和涌浪的綜合,其能量是不同方向、不同頻率的波浪分量的總和。本文采用Hanson 和Phillips? [2] 介紹的方法實現(xiàn)波浪分離。 H,F(xiàn),Θ 分別表示有效波高、譜峰頻率及其平均方向,并用下標(biāo) w,ww,sw 分別表征觀測波浪(綜合波浪)、風(fēng)浪和涌浪。

    1.3 非線性分位擬合

    由于有效波高是刻畫波浪能量的關(guān)鍵參數(shù),是風(fēng)速、風(fēng)向等因子的綜合反映,在近岸區(qū),風(fēng)向?qū)Σɡ说挠绊懹葹橥怀?,但對它的定量描述卻一直未能得到有效的解決。為此,本文依據(jù)有效波高對風(fēng)向的敏感特征,采用分位擬合的方式定量表征風(fēng)向?qū)τ行Рǜ叩挠绊憽?/p>

    2 結(jié)果與分析

    2.1 風(fēng)場及有效波高的特征

    圖2(a)描述了有效波高的時序變化。該時序變化包含了3個主要過程(各子圖中均已標(biāo)注):① 時段為8月23日00:00至8月24日12:00,有效波高呈現(xiàn)三峰狀,中峰最強,且右肩峰強于左肩峰(簡稱Ⅰ過程);② 時段為8月30日03:00至8月31日23:00,同樣也存在三峰,右峰最強且窄(簡稱Ⅱ過程);③ 時段為9月3日10:00至9月4日05:00,基本呈兩峰(簡稱Ⅲ過程)。整個觀測時段內(nèi),最大有效波高為 0.26 ?dBar左右,平均為0.056 dBar。

    從10 km高處風(fēng)速( U ??10 )的時序變化來看,平均風(fēng)速為3 m/s,最大風(fēng)速為5.8 m/s(見圖2(b))。在該時段,受臺風(fēng)(蘇力)外圍的影響,風(fēng)速偏大。隨著臺風(fēng)逐漸靠近蘇南及隨后遠(yuǎn)離、消亡過程中,風(fēng)速整體偏大。與有效波高的過程相對應(yīng),每個過程的波高峰值均對應(yīng)風(fēng)速的峰值,但風(fēng)速的峰值大小并不能完全對應(yīng)波高峰值的強弱,如Ⅲ過程對應(yīng)的風(fēng)速顯著小于前兩個過程,但波浪峰值卻最強??梢姡L(fēng)速是影響波高的關(guān)鍵因子之一。

    從觀測時段開始至8月26日18:00,風(fēng)向均偏北,包含了Ⅰ過程(見圖2(b))。Ⅰ過程中左峰對應(yīng)的風(fēng)速大于最強的中峰所對應(yīng)的風(fēng)速;與此類似,8月24日 20:00 至8月26日06:00的風(fēng)速略大于Ⅰ過程對應(yīng)的風(fēng)速,但有效波高基本呈最低值(0.02~0.03 dBar左右),約為Ⅰ過程有效波高的1/10。二者之間的細(xì)微差別說明北風(fēng)或東北風(fēng)對有效波高影響偏弱,而西北風(fēng)的影響明顯(見圖2(b));Ⅱ過程的西南風(fēng)、東南風(fēng)對有效波高影響差別較小。結(jié)合觀測點的周邊地形可以發(fā)現(xiàn),北風(fēng)或東北風(fēng)(含東風(fēng))的吹程最小,波浪從風(fēng)場中獲得的能量最小,進(jìn)而說明該風(fēng)向?qū)τ行Рǜ哂绊戇h(yuǎn)小于西北風(fēng)和偏南風(fēng)??梢?,風(fēng)向是影響有效波高的另一關(guān)鍵因素。

    2.2 波能頻譜特征、風(fēng)速與有效波高的關(guān)系

    波能的頻譜時間序列刻畫了不同頻率( f )的波能隨時間的變化規(guī)律,圖3中色階表示波能頻率譜 ( S(f) ) 的對數(shù)值。為了便于比較,在圖3中標(biāo)注了波浪時序變化的3個主要過程。最大波能基本分布于0.5 Hz的頻率附近,而且可以看到在Ⅰ過程至Ⅱ過程過渡期間還有2個相對弱的過程(見圖3)。結(jié)合波能分布和風(fēng)矢量變化可以發(fā)現(xiàn),東北風(fēng)階段所對應(yīng)的波能很弱(見圖2~3);除位于0.5 Hz的波能主峰外,某些時段(如Ⅲ過程)位于高頻(1.0 Hz附近)處也呈現(xiàn)了次峰:如I過程的中峰及右肩峰;Ⅱ過程中最強的右峰;Ⅲ過程中的兩峰;在Ⅰ、Ⅲ過程中的更低頻率處,隱約呈現(xiàn)了弱波能分布(見圖3)。結(jié)合風(fēng)速、風(fēng)向的分布特征(見圖2)可知,風(fēng)速、風(fēng)向是影響波能頻率譜呈雙峰分布的關(guān)鍵因素。

    2.3 風(fēng)向?qū)τ行Рǜ叩挠绊?/p>

    水體波浪的能量主要來源于風(fēng)場,其中風(fēng)區(qū)的長短是波能從風(fēng)場中獲取能量多少的關(guān)鍵。本文的觀測站點位于太湖北部近岸區(qū),使得周邊地形對風(fēng)區(qū)的影響遠(yuǎn)強于湖心區(qū),可見風(fēng)向是通過影響風(fēng)區(qū)的方式修正風(fēng)速對波浪的驅(qū)動效應(yīng)。有效波高是定量表征波能的關(guān)鍵參數(shù),也是波浪對多因子(風(fēng)速、風(fēng)浪狀態(tài)、水深等)的綜合響應(yīng)。因而可通過非線性高分位回歸擬合法? [11] 篩選出對風(fēng)向最為敏感的有效波高樣本,進(jìn)而定量描述風(fēng)向?qū)τ行Рǜ叩挠绊憽?/p>

    由于采樣點位于梅梁灣和貢湖灣連接處的北部近岸區(qū),靠近太湖東岸(見圖1),從正北向(0 rad)開始,隨順時針旋轉(zhuǎn),即風(fēng)向轉(zhuǎn)向東,由于周邊地形作用使得風(fēng)區(qū)逐漸變短,導(dǎo)致風(fēng)能轉(zhuǎn)化為波能減弱,隨著風(fēng)向的持續(xù)增大,其風(fēng)區(qū)又逐步延長,風(fēng)向的限制作用減弱,一直持續(xù)到西北向,風(fēng)向的限制作用才再次逐漸顯現(xiàn)出來??梢婏L(fēng)向?qū)︼L(fēng)區(qū)的影響與角度密切相關(guān),且呈周期性,為此采用三角函數(shù)的高分位法(95%)擬合有效波高隨風(fēng)向的變化,擬合函數(shù)形式為

    H w(W d)=a sin (bW d+c)+d (1)

    式中: W d為風(fēng)向,a,b,c,d 為常數(shù)。

    通過高分位的方法篩選出有效波高對風(fēng)向敏感的樣本,采用三角函數(shù)形式擬合有效波高對風(fēng)向的依賴,擬合效果良好(見圖4),說明三角函數(shù)可以定量表征風(fēng)向?qū)τ行Рǜ叩挠绊憽?/p>

    2.4 風(fēng)場與有效波高的關(guān)系

    依據(jù)波浪的組成及其形成機制可知:風(fēng)浪是由局地風(fēng)場通過摩擦引起的不穩(wěn)定或共振方式形成的,因而風(fēng)浪的有效波高必定與風(fēng)速呈現(xiàn)密切的相關(guān)性,但涌浪通常與局地風(fēng)場無顯著相關(guān)性。

    采用指數(shù)函數(shù)形式? [3] 對綜合波浪的有效波高( H? w? )與風(fēng)速進(jìn)行擬合, H? w? 與風(fēng)速之間的確定系數(shù)為0.04,說明二者之間無顯著的相關(guān)性,可見除風(fēng)速外,另有其他因子影響風(fēng)速與有效波高關(guān)系。此類因素主要有兩種:① 風(fēng)向影響風(fēng)區(qū)的長短,進(jìn)而調(diào)制風(fēng)能轉(zhuǎn)化波能,從而降低了有效波高對風(fēng)速的依賴;② 波浪中的風(fēng)浪、涌浪與風(fēng)速的關(guān)系差異顯著,因而波浪中涌浪的比例是影響有效波高與風(fēng)速關(guān)系的另一關(guān)鍵因子。

    為此,采用風(fēng)向?qū)Σǜ叩恼{(diào)制形式定量表征耦合風(fēng)向的風(fēng)速 W? sd? (簡稱耦合風(fēng)速):

    W? sd =U? 10 [a sin (bW d+c)+d]/(a+d)

    進(jìn)而利用e指數(shù)函數(shù)形式擬合 H? w 隨W? sd? 的變化,發(fā)現(xiàn)二者的相關(guān)性得到顯著的升高,其決定關(guān)系 R? 2? 為0.45(見圖5(a))。風(fēng)速、耦合風(fēng)速分別與 H? w? 相關(guān)性的差別說明風(fēng)向是影響波浪關(guān)鍵要素,風(fēng)場為波浪能量的主要來源; W? sd? 可以有效表征風(fēng)向、風(fēng)速的耦合。

    為了探究耦合風(fēng)速對北部近岸區(qū)波浪組成比例的影響,選擇 W? sd? 大于2 m/s的樣本與對應(yīng)的有效波高進(jìn)行擬合,發(fā)現(xiàn)二者之間的關(guān)系急劇下降(見圖5(b)):結(jié)合風(fēng)向的調(diào)制作用, W? sd 越大,風(fēng)場對波浪的驅(qū)動效應(yīng)增強,綜合波浪可能受到邊界反射等影響增大,涌浪的比例增加,進(jìn)而導(dǎo)致其與局地W? sd 相關(guān)性急劇減弱??梢姡ɡ酥邪擞坷?W? sd 與H? w 的相關(guān)性及其隨W? sd 的變化證實了隨著W? sd? 的增大,涌浪比例升高。

    2.5? 風(fēng)浪有效波高、平均頻率及平均方向的特征

    風(fēng)場是太湖波浪產(chǎn)生的主要驅(qū)動力,使得波浪中涌浪比例提高,因而通過比較耦合風(fēng)速 (W? sd )與風(fēng)浪的有效波高(H? ww )、風(fēng)浪的譜峰頻率(F? ww )及風(fēng)浪的平均波向(Θ? ww ) 的相關(guān)性,探究太湖北部近岸區(qū)的波浪組成特征,驗證波浪分離方法的有效性。

    北部近岸區(qū)波浪的 H? ww 隨W? sd? 的變化與湖心區(qū)基本一致,即二者之間呈現(xiàn)e的負(fù)指數(shù)函數(shù)關(guān)系,其決定系數(shù)( R?? 2 )為0.65(見圖6(a));當(dāng) W? sd 大于2 ?m/s ,W? sd? 與相應(yīng)的 H? ww? 的決定系數(shù)( R? 2? )略有降低,其值為0.62(見圖6(b))。

    通過比較 W? sd 與H? ww 、W? sd 與H? ww? 的關(guān)系變化(見圖5~6),可以發(fā)現(xiàn)太湖北部近岸區(qū)波浪對耦合風(fēng)速的響應(yīng)迅速,風(fēng)場是太湖波浪的關(guān)鍵驅(qū)動因子;岸線的反射等作用使得北部近岸波浪中包含一定比例的涌浪,且涌浪的比例隨耦合風(fēng)速的增大而提高,從而使得當(dāng) W? sd? 大于2 m/s時, W? sd 和H w 的相關(guān)關(guān)系急劇下降;同樣條件下,風(fēng)浪的有效波高 (H? ww )和W? sd? 的相關(guān)關(guān)系證明分離出的風(fēng)浪中基本未受涌浪污染。波浪分離方法切實可行,且適用于太湖北部近岸區(qū)。

    風(fēng)浪的平均頻率( F? ww? )為1.37 Hz,最大為2.23 Hz,最小為0.42 Hz,并隨 W? sd? 的增大呈現(xiàn)遞減的趨勢(見圖7)。通常波浪的能量主要集中于波峰; F? ww? 遠(yuǎn)大于波浪的譜峰頻率,是對快變風(fēng)場特征的響應(yīng)? [12] 。通常風(fēng)浪的能量直接獲取于風(fēng)場,獲取的能量一方面彌補于高頻處的能量耗散(湍流、波浪破碎等);另一方面通過非線性作用向低頻傳輸,促進(jìn)波浪的發(fā)展,使得其頻率降低。誠然,當(dāng)從風(fēng)場中獲取的能量不足,低頻處的能量同樣補充高頻處能量損耗,導(dǎo)致波浪的減弱甚至消亡。因而, F? ww 隨W? sd? 的變化也從側(cè)面證實了波浪分離的有效性。

    風(fēng)浪是由風(fēng)直接驅(qū)動產(chǎn)生的,對風(fēng)場的響應(yīng)迅速,因而其波長相對較短,平均方向也與風(fēng)場方向較為接近。觀測點離岸界較近,且岸界基本圍繞在觀測點的正北到東南方向(見圖1),因而在岸界對波浪的反射作用下,觀測波浪中涌浪比例提高。波浪方向與風(fēng)場方向接近反向,因而從正北到東南方向,波浪方向與風(fēng)場基本成負(fù)相關(guān);隨著風(fēng)向角度的繼續(xù)增大,岸界離觀測點的距離逐漸加大,岸界對風(fēng)浪的反射作用基本消失,涌浪比例偏低,因而觀測的綜合波浪方向與風(fēng)向接近,二者呈正相關(guān)(見圖8(a))。從風(fēng)浪平均方向與風(fēng)向關(guān)系來看,二者呈線性關(guān)系,其決定系數(shù)為0.57,絕大部分樣本均圍繞在擬合直線周圍(見圖8(b)),但也存在少量的樣本偏離該直線。這些少量樣本對應(yīng)的風(fēng)向在正北(0 rad)附近,波浪方向也接近2π rad(見圖8(b)),其中0 rad和2π rad的方向?qū)嶋H上是圍繞在正北向附近,說明上述樣本的風(fēng)向和風(fēng)浪的平均方向也呈線性關(guān)系,即分離出的風(fēng)浪平均方向與風(fēng)向基本一致,這也證實了波浪分離方法適用于北岸區(qū)??傊瑥娘L(fēng)場與風(fēng)浪的譜峰頻率、平均方向及有效波高的關(guān)系3個方面,證實了耦合風(fēng)速為波浪的關(guān)鍵驅(qū)動因素,岸界的反射作用造成北部近岸區(qū)波浪中含有較高比例的涌浪,同時也從風(fēng)浪方向、頻率等方面證實了波浪分離的方法適用于太湖。

    2.6 涌浪特征及其與耦合風(fēng)速的關(guān)系

    波浪平均方向和風(fēng)向的關(guān)系證明北部近岸區(qū)的波浪組成中包含了涌浪,為此有必要探究涌浪的特征及其對耦合風(fēng)速的響應(yīng)規(guī)律。

    圖9描述了涌浪的有效波高 H? sw 隨W? sd? 的變化,該圖呈現(xiàn)了兩典型特征:① ?H? sw 隨W? sd? 呈快速升高(簡稱快變趨勢);② ?H? sw 隨W? sd? 呈緩慢增加(簡稱緩變趨勢)。為此,采用e指數(shù)形式的分位擬合,分位數(shù)分別為85%和35%,擬合效果良好。快變趨勢可能是由于風(fēng)向偏西南,即觀測點波浪受岸界反射疊加作用最強,隨耦合風(fēng)速的增大,岸界反射造成的疊加效應(yīng)增強,使得 H? sw 隨W? sd? 快速增長。當(dāng)風(fēng)向偏東或北時,岸界的反射作用對觀測點影響微弱,因而 H? sw 隨W? sd? 增長相對偏弱。另外,由于太湖涌浪頻率偏高,其涌浪傳播過程中的衰減也較低頻波快,因而反射傳輸路程越遠(yuǎn),衰減越大,進(jìn)而凸顯快、慢變化趨勢的分化??梢?,上述原因是導(dǎo)致 H? sw 隨W? sd? 呈增長趨勢分化的關(guān)鍵。

    H? sw 和W? sd? 之間存在線性關(guān)系,兩者之間的 R?? 2 為 0.158 ,遠(yuǎn)低于 H? ww 、H w分別與W? sd? 的關(guān)系。盡管涌浪的產(chǎn)生機制導(dǎo)致涌浪與局地風(fēng)場難以存在顯著的相關(guān)性,但由于太湖的波浪主要由風(fēng)場驅(qū)動,風(fēng)場作用于水體表面的吹程較短,使得由于岸界反射等作用形成北部近岸區(qū)的涌浪頻率偏大,在一定程度上也攜帶了風(fēng)場的信息,所以 H? sw 和W? sd? 也存在弱相關(guān)性。

    鑒于近岸區(qū)域波浪的生態(tài)意義極為重要,定量分析涌浪的比例有利于提高計算波浪切應(yīng)力的精度。 從 H? sw /H w 隨風(fēng)向的變化來看(見圖10),絕大部分? H? sw /H w? 的樣本均大于0.6。從其隨風(fēng)向分布來看,偏東風(fēng)的風(fēng)區(qū)較短,風(fēng)直接驅(qū)動產(chǎn)生的波浪很弱,因而波浪中涌浪占比接近于1。為此,在評價近岸區(qū)波浪的動力過程中,應(yīng)在風(fēng)浪上疊加基于耦合風(fēng)速估算的涌浪成分,提高估算波浪切應(yīng)力的精度。

    3 討 論

    太湖屬于典型的大型淺水湖泊,平均水深僅有 1.8 ?m。其譜峰頻率(0.5 Hz)遠(yuǎn)大于海洋的譜峰頻率( 0.075 ?Hz);波峰最大能量(0.001 m? 2 /Hz)遠(yuǎn)小于海洋的對應(yīng)值(2.5 m? 2 /Hz)? [13] 。在太湖北部開闊區(qū)域的波高隨風(fēng)速遵循e指數(shù)函數(shù)的形式,與本文中風(fēng)向?qū)τ行Рǜ叩挠绊懱卣飨辔呛?,但其參?shù)受周邊地形等因素的影響顯著? [8] 。與海洋相比,太湖水淺,加之面積不大,使得湖泊水體從風(fēng)場中獲取的能量較少,因而其波能遠(yuǎn)小于海洋。也正是因為水淺,使得波浪引發(fā)的切應(yīng)力可作用于湖底,成為導(dǎo)致底泥再懸浮、影響水生生態(tài)系統(tǒng)的關(guān)鍵因素。

    太湖屬于典型的內(nèi)陸湖泊,風(fēng)速和風(fēng)向的耦合與綜合波浪有效波高的關(guān)系證實了風(fēng)速(尤其是耦合風(fēng)速)是太湖的主要驅(qū)動力。內(nèi)陸湖泊風(fēng)向變化顯著,風(fēng)區(qū)長度的空間變化明顯,使得風(fēng)速、風(fēng)向?qū)Σɡ说尿?qū)動作用空間差異明顯,其中風(fēng)向通過三角函數(shù)的形式調(diào)制風(fēng)速對波浪起驅(qū)動作用。加上沿岸對波浪的反射等特征,導(dǎo)致波浪組成比例時空差異突出,因而近岸區(qū)的波浪特征遠(yuǎn)比開闊區(qū)復(fù)雜,且呈現(xiàn)出顯著的空間差異特征。

    本文所用數(shù)據(jù)的觀測時段超過13 d,期間經(jīng)歷了完整的臺風(fēng)(蘇力)影響過程,其風(fēng)速變化范圍寬泛(見圖2(b));風(fēng)向的變化也涵蓋了0~2π的范圍(見圖4),可見觀測到的波浪數(shù)據(jù)既體現(xiàn)了風(fēng)場(風(fēng)速、風(fēng)向)的驅(qū)動效應(yīng),也包含了岸界的反射信息,甚至也反映了風(fēng)場與岸界的耦合效應(yīng)。因而,本文的數(shù)據(jù)涵蓋了風(fēng)場和波浪的多個過程及其演變規(guī)律(見圖2(a)及圖3),具有良好的代表性。

    近岸區(qū)對波浪的反射作用會改變波浪的組成,提高涌浪的比例。這種改變效應(yīng)也受風(fēng)向及岸界的共同影響,因而近岸區(qū)的波浪中含涌浪比例很高,其比例最低大于0.6,最高大于0.9。由于在大型淺水湖泊中波浪對耦合風(fēng)速響應(yīng)迅速,加上近岸岸界反射作用使得涌浪在一定程度上攜帶了耦合風(fēng)速的信息,因而太湖中波浪與耦合風(fēng)速的關(guān)系較海洋中更為密切;但又和海洋一樣,波浪中包含了較高比例的涌浪,因而本次研究的進(jìn)展將提高近岸區(qū)波浪參數(shù)化方案的精度。

    4 結(jié) 論

    (1) 太湖北部近岸區(qū)波浪呈現(xiàn)中高頻特征,三角函數(shù)可以表征風(fēng)向?qū)Σɡ说亩坑绊憽?/p>

    (2) 風(fēng)速(尤其是耦合風(fēng)速)是太湖波浪的主要驅(qū)動因素,該區(qū)波浪既包含風(fēng)浪,也包含涌浪,其中涌浪的占比較高,且隨風(fēng)向變化。

    (3) 太湖北部近岸區(qū)的涌浪隨耦合風(fēng)速呈現(xiàn)快速和緩慢的e指數(shù)增長現(xiàn)象。

    致 謝

    感謝中國科學(xué)院太湖生態(tài)系統(tǒng)實驗站為本次研究提供相關(guān)的氣象數(shù)據(jù)。

    參考文獻(xiàn):

    [1]? GARCIA ?G W.Wave bimodal spectrum based on swell and wind-sea components[J].IFAC-Papers OnLine,2015,48-16:223-228.

    [2] HANSON J L,PHILLIPS O M.Automated analysis of ocean surface directional wave spectra[J].Journal of Atmospheric and Oceanic Technology,2009,18:277-293.

    [3] DING W,WU T,QIN B,et al.Features and impacts of currents and waves on sediment resuspension in a large shallow in china[J].Environmental Science and Pollution Research,2018,25:36341-36354.

    [4] 羅斂蔥,秦伯強.太湖波浪與湖流對沉積物再懸浮不同影響的研究[J].水文,2003,23(6):1-4.

    [5] 文圣常,余宙文.海浪理論與計算原理[M].北京:科學(xué)出版社,1984.

    [6] 逄勇,莊巍,韓濤,等.風(fēng)浪擾動下的太湖懸浮物實驗與模擬[J].環(huán)境科學(xué),2008,29(10):2743-2748.

    [7] 朱廣偉,許海,朱夢圓,等.三十年來長江中下游湖泊富營養(yǎng)化狀況變遷及其影響因素[J].湖泊科學(xué),2019,31(6):1510-1524.

    [8] 胡維平,胡春華,張發(fā)兵,等.太湖北部風(fēng)浪波高計算模式觀測分析[J].湖泊科學(xué),2005,17(1):41-46.

    [9] 范成新,劉元波,陳荷生.太湖底泥蓄積量估算及分布特征探討[J].上海環(huán)境科學(xué),2000,19(2):72-75.

    [10]? 許遐禎,陶蓉茵,趙巧華,等.大型淺水湖泊太湖波浪特征及其對風(fēng)場的敏感性分析[J].湖泊科學(xué),2013,25(1):55-64.

    [11] BISSINGER J E,MONTAGNES D J S,SHARPLES J,et al.Predicting marine phytoplankton maximum growth rates from temperature:Improving on the Eppley curve using quantile regression[J].Limnol.Oceanogr.,2008,53(2):487-493.

    [12] PORTILLA J,MONBALIU J.Spectral partitioning and identification of wind sea and swell[J].Journal of Atmospheric and Oceanic Technology,2009,26:107-122.

    [13] KPOGO-NUWOKLO KA,AILLIOT P,OLAGNON M,et al.Improvement ocean wave spectra estimation using the temporal structure of wave systems[J].Coastal Engineering,2015,96:81-91.

    (編輯:胡旭東)

    Analysis on wave composition and evolution characteristics in northern coastal ?area of Taihu Lake in summer

    ZHAO Qiaohua,XU Jia,WANG Jianjian

    ( School of Hydrology and Water Resources,Nanjing University of Information Science and Technology,Nanjing 210044,China )

    Abstract:

    Wave is the key dynamic process that drives the evolution of aquatic ecosystem in large shallow lakes.In this paper,the wave composition and its evolution characteristics were analyzed based on the observation data of wave direction spectrum,wind speed and wind direction in the northern coastal area of Taihu Lake.The results showed that the wave frequency spectrum in the near-shore area of northern Taihu Lake was mainly unimodal or bimodal,belonging to medium and high frequency wave.Wind is the main driving force for wave formation,and the modulation of wind direction or surrounding terrain on the significant wave height can be quantitatively characterized by trigonometric function.The relationship between the significant wave height,spectral peak frequency and its average direction of wind waves,and wind speed or direction confirmed that the wave in the northern coastal area of Taihu Lake were composed of surge and wind waves,and the surge proportion was relatively high.The reflection of the shoreline increased the surge proportion in the northern coastal area.The significant wave height of surge not only increased exponentially with the coupling wind speed (coupling of wind speed and wind direction),but also had the phenomenon of fast and slow differentiation.The research results can help to improve the accuracy of wave composition estimation in the near-shore area of large shallow lakes,and ensure the effect of aquatic ecological restoration.

    Key words:

    wave composition;wind wave;surge;significant wave height;Taihu Lake

    猜你喜歡
    波高風(fēng)浪風(fēng)場
    基于FHDI-GNWM 數(shù)據(jù)的全球超越概率波高宏觀分布特征分析
    坐在同一條船上
    風(fēng)浪干擾條件下艦船航向保持非線性控制系統(tǒng)
    基于FLUENT的下?lián)舯┝魅S風(fēng)場建模
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計有效波高研究
    海洋通報(2021年3期)2021-08-14 02:20:46
    非平整港池的多向不規(guī)則波試驗研究
    不怕風(fēng)浪的小船
    East Side Story
    “最美風(fēng)場”的贏利法則
    能源(2017年8期)2017-10-18 00:47:39
    側(cè)向風(fēng)場中無人機的飛行研究
    欧美一区二区精品小视频在线| 欧美在线一区亚洲| 亚洲人与动物交配视频| 制服人妻中文乱码| 亚洲av五月六月丁香网| 国产三级在线视频| 国产精品日韩av在线免费观看| 日本一二三区视频观看| 蜜桃久久精品国产亚洲av| 麻豆一二三区av精品| 日本黄色视频三级网站网址| 日韩欧美国产一区二区入口| 午夜免费观看网址| 久久久国产精品麻豆| av在线天堂中文字幕| 黄色片一级片一级黄色片| av天堂在线播放| 国产乱人伦免费视频| 首页视频小说图片口味搜索| 亚洲真实伦在线观看| 国产欧美日韩一区二区三| 99国产极品粉嫩在线观看| 日本与韩国留学比较| 波多野结衣高清无吗| 欧美又色又爽又黄视频| 亚洲成人精品中文字幕电影| 国产成年人精品一区二区| 美女大奶头视频| 精品久久蜜臀av无| 97超级碰碰碰精品色视频在线观看| 麻豆国产av国片精品| 这个男人来自地球电影免费观看| 国产精品久久视频播放| 51午夜福利影视在线观看| 久久中文字幕人妻熟女| 视频区欧美日本亚洲| 嫩草影院精品99| av在线天堂中文字幕| 我要搜黄色片| 免费高清视频大片| 女生性感内裤真人,穿戴方法视频| 久久久久久国产a免费观看| 成年女人永久免费观看视频| 在线观看美女被高潮喷水网站 | 精品久久蜜臀av无| 亚洲一区二区三区不卡视频| 亚洲自拍偷在线| 亚洲九九香蕉| 欧美成人免费av一区二区三区| 搞女人的毛片| 亚洲片人在线观看| 母亲3免费完整高清在线观看| 校园春色视频在线观看| 国产精品国产高清国产av| 97人妻精品一区二区三区麻豆| 国产一区二区三区在线臀色熟女| 国产亚洲欧美98| 国产亚洲精品av在线| 伊人久久大香线蕉亚洲五| 国产欧美日韩一区二区三| 精品国产乱码久久久久久男人| 国产一区二区三区在线臀色熟女| 波多野结衣高清作品| 在线观看一区二区三区| 成人三级黄色视频| 丁香欧美五月| 久久婷婷人人爽人人干人人爱| 国产精品,欧美在线| 国产一区二区激情短视频| 观看美女的网站| 18禁裸乳无遮挡免费网站照片| x7x7x7水蜜桃| 欧美成人免费av一区二区三区| 国产综合懂色| 一级毛片精品| 国内精品一区二区在线观看| 午夜视频精品福利| 搡老岳熟女国产| 噜噜噜噜噜久久久久久91| 欧美又色又爽又黄视频| 久久久国产成人精品二区| 亚洲国产色片| 99精品久久久久人妻精品| 女警被强在线播放| 亚洲第一欧美日韩一区二区三区| 欧美另类亚洲清纯唯美| 人妻丰满熟妇av一区二区三区| 国产1区2区3区精品| 国产亚洲精品av在线| 色老头精品视频在线观看| 国产午夜福利久久久久久| 精品久久久久久,| 国产一区二区在线av高清观看| 国产精品久久久久久人妻精品电影| 法律面前人人平等表现在哪些方面| 精品福利观看| 黄片小视频在线播放| 国产欧美日韩一区二区精品| 国产免费男女视频| 色噜噜av男人的天堂激情| a级毛片在线看网站| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 国产午夜精品久久久久久| bbb黄色大片| 欧美高清成人免费视频www| 日韩成人在线观看一区二区三区| 午夜福利在线观看吧| 黑人欧美特级aaaaaa片| 日韩欧美 国产精品| 九九热线精品视视频播放| 首页视频小说图片口味搜索| 亚洲人成网站在线播放欧美日韩| 国产av不卡久久| 在线免费观看不下载黄p国产 | 99在线人妻在线中文字幕| 亚洲黑人精品在线| 在线观看66精品国产| 亚洲欧美精品综合久久99| 国产精品久久电影中文字幕| 国产成人av激情在线播放| 欧美日本视频| 欧美三级亚洲精品| 97人妻精品一区二区三区麻豆| 日本精品一区二区三区蜜桃| 亚洲九九香蕉| 欧美成狂野欧美在线观看| 免费看日本二区| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 亚洲 国产 在线| 五月玫瑰六月丁香| 日本 欧美在线| 蜜桃久久精品国产亚洲av| 国产精品亚洲一级av第二区| 欧美日韩中文字幕国产精品一区二区三区| 成年女人毛片免费观看观看9| 亚洲av五月六月丁香网| 日本精品一区二区三区蜜桃| 国产v大片淫在线免费观看| 亚洲色图av天堂| 国产精品免费一区二区三区在线| 99久久精品国产亚洲精品| 岛国在线观看网站| 丰满人妻一区二区三区视频av | 在线看三级毛片| 成人高潮视频无遮挡免费网站| 黄色女人牲交| 精品无人区乱码1区二区| 少妇熟女aⅴ在线视频| 国产欧美日韩一区二区三| 国产真实乱freesex| 搡老熟女国产l中国老女人| 一区福利在线观看| 国产亚洲精品久久久com| 99国产精品99久久久久| av国产免费在线观看| 丝袜人妻中文字幕| 精品久久久久久久人妻蜜臀av| 亚洲自拍偷在线| 中国美女看黄片| 窝窝影院91人妻| 成人三级黄色视频| 俺也久久电影网| 小蜜桃在线观看免费完整版高清| 色噜噜av男人的天堂激情| 午夜福利视频1000在线观看| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av | 校园春色视频在线观看| 色哟哟哟哟哟哟| 最近视频中文字幕2019在线8| 精品99又大又爽又粗少妇毛片 | 在线看三级毛片| 国产乱人伦免费视频| 十八禁网站免费在线| 美女 人体艺术 gogo| 国产欧美日韩精品一区二区| 免费在线观看视频国产中文字幕亚洲| 国产毛片a区久久久久| 在线观看免费视频日本深夜| 精品一区二区三区视频在线观看免费| 最新中文字幕久久久久 | 亚洲一区高清亚洲精品| 国产av麻豆久久久久久久| 日韩中文字幕欧美一区二区| 精品国产三级普通话版| 亚洲欧美一区二区三区黑人| 欧美丝袜亚洲另类 | 国产乱人视频| a级毛片a级免费在线| 99热6这里只有精品| 精品福利观看| 偷拍熟女少妇极品色| 久久久久精品国产欧美久久久| 老熟妇仑乱视频hdxx| 日本a在线网址| 伦理电影免费视频| 色老头精品视频在线观看| 国产精品久久久久久亚洲av鲁大| 美女高潮喷水抽搐中文字幕| 最近视频中文字幕2019在线8| 最近在线观看免费完整版| 久久久精品大字幕| 999久久久精品免费观看国产| 88av欧美| 亚洲五月婷婷丁香| 99热这里只有精品一区 | 男女午夜视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产真实乱freesex| 日韩欧美在线乱码| 日本与韩国留学比较| 制服丝袜大香蕉在线| 中出人妻视频一区二区| 国内精品一区二区在线观看| 男女做爰动态图高潮gif福利片| 这个男人来自地球电影免费观看| 高清在线国产一区| 亚洲成人久久性| 丰满人妻一区二区三区视频av | 久久中文看片网| 国产精品久久视频播放| 国产高清videossex| 日韩中文字幕欧美一区二区| 亚洲av成人av| 麻豆国产av国片精品| 真实男女啪啪啪动态图| 精品一区二区三区视频在线观看免费| 精品乱码久久久久久99久播| 美女cb高潮喷水在线观看 | 国产精品av视频在线免费观看| 90打野战视频偷拍视频| 午夜精品在线福利| 久久久久久久久免费视频了| 757午夜福利合集在线观看| 我的老师免费观看完整版| 美女免费视频网站| 亚洲av熟女| 亚洲av电影不卡..在线观看| 国产又黄又爽又无遮挡在线| 成年版毛片免费区| 中文字幕人成人乱码亚洲影| 性色av乱码一区二区三区2| 色视频www国产| 亚洲自偷自拍图片 自拍| 欧美在线黄色| 精品久久久久久久末码| 亚洲av成人精品一区久久| 久久久久久人人人人人| 一个人免费在线观看的高清视频| 亚洲欧美激情综合另类| 久久久久免费精品人妻一区二区| 国产一区二区三区视频了| 天天一区二区日本电影三级| 一a级毛片在线观看| 久久九九热精品免费| 两人在一起打扑克的视频| 男女那种视频在线观看| 身体一侧抽搐| 日韩欧美精品v在线| 久久天躁狠狠躁夜夜2o2o| 99久久成人亚洲精品观看| 三级国产精品欧美在线观看 | 欧美中文综合在线视频| 国内精品美女久久久久久| 亚洲欧美一区二区三区黑人| 久久久精品大字幕| 日本成人三级电影网站| 成年女人永久免费观看视频| 免费在线观看亚洲国产| 少妇的逼水好多| 国产野战对白在线观看| 精品国产乱码久久久久久男人| 成年女人毛片免费观看观看9| 国产精品国产高清国产av| 国产蜜桃级精品一区二区三区| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻| 日本五十路高清| 欧美在线一区亚洲| 黑人欧美特级aaaaaa片| 国产高清视频在线播放一区| or卡值多少钱| 女人被狂操c到高潮| 亚洲欧美精品综合一区二区三区| 国产91精品成人一区二区三区| 可以在线观看的亚洲视频| 久久久久免费精品人妻一区二区| avwww免费| 久久九九热精品免费| 成人高潮视频无遮挡免费网站| 亚洲av电影在线进入| 久久人妻av系列| 一区福利在线观看| 欧美中文日本在线观看视频| 91字幕亚洲| 毛片女人毛片| 亚洲电影在线观看av| 嫩草影院入口| 国产乱人伦免费视频| 黄频高清免费视频| 国产亚洲精品综合一区在线观看| 91麻豆精品激情在线观看国产| 欧美日韩一级在线毛片| 国产伦一二天堂av在线观看| 亚洲 国产 在线| 久久久精品欧美日韩精品| 黄色片一级片一级黄色片| 免费一级毛片在线播放高清视频| 日韩免费av在线播放| 俄罗斯特黄特色一大片| 国产久久久一区二区三区| 国产99白浆流出| 偷拍熟女少妇极品色| 成人18禁在线播放| 色哟哟哟哟哟哟| 免费大片18禁| 色尼玛亚洲综合影院| 国产爱豆传媒在线观看| 国产高清有码在线观看视频| 制服人妻中文乱码| 99国产极品粉嫩在线观看| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 在线视频色国产色| tocl精华| 亚洲欧美日韩东京热| 19禁男女啪啪无遮挡网站| 淫妇啪啪啪对白视频| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| www日本黄色视频网| 成人鲁丝片一二三区免费| 神马国产精品三级电影在线观看| 成年人黄色毛片网站| 熟妇人妻久久中文字幕3abv| 免费无遮挡裸体视频| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 久久精品国产99精品国产亚洲性色| 黑人巨大精品欧美一区二区mp4| 欧美一区二区国产精品久久精品| 色av中文字幕| av在线蜜桃| 国产精品精品国产色婷婷| 亚洲av成人精品一区久久| 1024香蕉在线观看| 国产aⅴ精品一区二区三区波| 国产蜜桃级精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 婷婷六月久久综合丁香| 日韩欧美在线二视频| av欧美777| 特大巨黑吊av在线直播| www日本黄色视频网| 免费看日本二区| 欧美av亚洲av综合av国产av| 成熟少妇高潮喷水视频| av欧美777| 国产激情久久老熟女| 在线观看一区二区三区| 国产精品久久电影中文字幕| 国产欧美日韩一区二区三| 久久久久精品国产欧美久久久| 九九久久精品国产亚洲av麻豆 | 全区人妻精品视频| 免费观看的影片在线观看| 精品不卡国产一区二区三区| 国产亚洲精品久久久com| 动漫黄色视频在线观看| 1024手机看黄色片| av福利片在线观看| 噜噜噜噜噜久久久久久91| a级毛片a级免费在线| 精品免费久久久久久久清纯| 成人特级av手机在线观看| 亚洲国产欧美网| 香蕉国产在线看| 嫩草影院入口| 18禁裸乳无遮挡免费网站照片| 成人无遮挡网站| 好男人在线观看高清免费视频| 欧美成狂野欧美在线观看| 一个人免费在线观看的高清视频| 欧美绝顶高潮抽搐喷水| 9191精品国产免费久久| 女同久久另类99精品国产91| 亚洲欧美日韩无卡精品| 最近最新中文字幕大全免费视频| av视频在线观看入口| 亚洲国产精品成人综合色| 好男人在线观看高清免费视频| 国产精品美女特级片免费视频播放器 | 日本 欧美在线| 在线观看日韩欧美| 女人被狂操c到高潮| 亚洲av成人精品一区久久| 中文资源天堂在线| 国产一区二区三区视频了| 中文字幕人成人乱码亚洲影| 琪琪午夜伦伦电影理论片6080| 久久精品91无色码中文字幕| 97超级碰碰碰精品色视频在线观看| 中文字幕久久专区| 老熟妇乱子伦视频在线观看| 欧美日本亚洲视频在线播放| 日本精品一区二区三区蜜桃| 男女视频在线观看网站免费| 老司机午夜十八禁免费视频| 中文字幕av在线有码专区| 国产成+人综合+亚洲专区| av天堂中文字幕网| 久久久久亚洲av毛片大全| 国产精品九九99| 成人鲁丝片一二三区免费| 久久久久精品国产欧美久久久| 亚洲欧洲精品一区二区精品久久久| 高潮久久久久久久久久久不卡| 又黄又爽又免费观看的视频| 日本免费一区二区三区高清不卡| 97人妻精品一区二区三区麻豆| 五月玫瑰六月丁香| 亚洲欧美一区二区三区黑人| 99在线视频只有这里精品首页| 99久久综合精品五月天人人| 99热这里只有是精品50| 一个人看视频在线观看www免费 | 精品99又大又爽又粗少妇毛片 | 国产伦人伦偷精品视频| 久久久色成人| 欧美一级毛片孕妇| 成人欧美大片| 亚洲 欧美 日韩 在线 免费| 伊人久久大香线蕉亚洲五| 国产成人av激情在线播放| 国产1区2区3区精品| 999久久久国产精品视频| 哪里可以看免费的av片| 亚洲成av人片免费观看| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区| 久久精品国产清高在天天线| 亚洲精华国产精华精| 性欧美人与动物交配| 国产又黄又爽又无遮挡在线| 亚洲自偷自拍图片 自拍| 国产三级在线视频| 91av网站免费观看| 午夜福利成人在线免费观看| 在线观看舔阴道视频| 此物有八面人人有两片| 淫秽高清视频在线观看| 成人av在线播放网站| 亚洲在线观看片| 一区二区三区高清视频在线| 国产一区二区激情短视频| 无人区码免费观看不卡| 成人永久免费在线观看视频| 男女下面进入的视频免费午夜| 男人的好看免费观看在线视频| 最近最新中文字幕大全电影3| 欧美日韩黄片免| 国产午夜精品久久久久久| 日本撒尿小便嘘嘘汇集6| 国产免费男女视频| www日本黄色视频网| 日韩欧美免费精品| 成熟少妇高潮喷水视频| 亚洲欧美日韩东京热| 国产野战对白在线观看| 亚洲九九香蕉| 国产成人一区二区三区免费视频网站| 波多野结衣巨乳人妻| 波多野结衣高清无吗| 久久久精品欧美日韩精品| 无人区码免费观看不卡| 国产亚洲精品av在线| 啦啦啦免费观看视频1| 国产伦人伦偷精品视频| 欧美中文综合在线视频| 欧美日韩综合久久久久久 | 又爽又黄无遮挡网站| av福利片在线观看| 国产一区二区在线av高清观看| 免费在线观看成人毛片| 国产精品98久久久久久宅男小说| 19禁男女啪啪无遮挡网站| 国产97色在线日韩免费| 动漫黄色视频在线观看| 日韩欧美免费精品| 在线播放国产精品三级| 国产乱人视频| 欧美激情在线99| 99久久精品国产亚洲精品| 亚洲熟女毛片儿| 成在线人永久免费视频| 国产人伦9x9x在线观看| 男女午夜视频在线观看| 国产亚洲av嫩草精品影院| 超碰成人久久| 欧美黑人欧美精品刺激| 美女高潮的动态| 日韩中文字幕欧美一区二区| 欧美性猛交黑人性爽| 伦理电影免费视频| 欧美又色又爽又黄视频| 成人无遮挡网站| 国产伦一二天堂av在线观看| 丰满人妻一区二区三区视频av | 黄片小视频在线播放| 桃色一区二区三区在线观看| 午夜两性在线视频| 免费观看人在逋| 精品国产乱子伦一区二区三区| 午夜影院日韩av| 丝袜人妻中文字幕| 久久久久免费精品人妻一区二区| 国内揄拍国产精品人妻在线| 精品无人区乱码1区二区| 久久久久久九九精品二区国产| 美女cb高潮喷水在线观看 | 视频区欧美日本亚洲| 午夜福利成人在线免费观看| 国产乱人视频| 青草久久国产| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 91麻豆av在线| 又粗又爽又猛毛片免费看| 给我免费播放毛片高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 五月玫瑰六月丁香| www.精华液| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 午夜两性在线视频| 免费搜索国产男女视频| 成人精品一区二区免费| 亚洲,欧美精品.| 精品99又大又爽又粗少妇毛片 | 久久欧美精品欧美久久欧美| 午夜福利视频1000在线观看| 亚洲av五月六月丁香网| 俄罗斯特黄特色一大片| 免费一级毛片在线播放高清视频| 美女黄网站色视频| 亚洲avbb在线观看| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| av福利片在线观看| 18禁美女被吸乳视频| 人人妻人人看人人澡| 床上黄色一级片| 99热这里只有是精品50| 丰满的人妻完整版| 精品久久蜜臀av无| 免费看a级黄色片| 99久久国产精品久久久| 黄色 视频免费看| 在线观看午夜福利视频| 成人国产综合亚洲| 观看美女的网站| 9191精品国产免费久久| 午夜福利免费观看在线| 国产成人av激情在线播放| 一进一出好大好爽视频| 欧美另类亚洲清纯唯美| 91av网一区二区| aaaaa片日本免费| 日韩人妻高清精品专区| 国产真实乱freesex| 国产极品精品免费视频能看的| a在线观看视频网站| 亚洲av五月六月丁香网| 日日干狠狠操夜夜爽| 好看av亚洲va欧美ⅴa在| 欧美性猛交黑人性爽| 99国产精品99久久久久| av中文乱码字幕在线| 国产亚洲欧美98| 在线观看一区二区三区| 成熟少妇高潮喷水视频| 国产一区在线观看成人免费| av片东京热男人的天堂| 国产精品久久久久久亚洲av鲁大| 在线观看日韩欧美| 99久久久亚洲精品蜜臀av| 久久性视频一级片| 亚洲 欧美 日韩 在线 免费| 国产成人影院久久av| 黄色日韩在线| 欧美成人性av电影在线观看| 国产探花在线观看一区二区| 午夜影院日韩av| 99riav亚洲国产免费| 国产亚洲精品一区二区www| 97超视频在线观看视频| 毛片女人毛片| 久久人妻av系列| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| 两个人看的免费小视频| 精品国产美女av久久久久小说| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 亚洲avbb在线观看| 男女下面进入的视频免费午夜| 婷婷丁香在线五月| www日本黄色视频网| 国产精品久久电影中文字幕| 国产精品久久久人人做人人爽| 精品久久久久久久久久久久久| 美女黄网站色视频| 欧美最黄视频在线播放免费|