• <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天堂久久9| 亚洲情色 制服丝袜| 亚洲午夜理论影院| 国产一区二区三区视频了| 久久狼人影院| 国产精品久久久久久人妻精品电影| 两性夫妻黄色片| 男女高潮啪啪啪动态图| 人人妻人人添人人爽欧美一区卜| 老司机亚洲免费影院| 欧美性长视频在线观看| 日韩中文字幕欧美一区二区| 老汉色∧v一级毛片| 91字幕亚洲| 欧美在线黄色| 琪琪午夜伦伦电影理论片6080| 免费观看人在逋| 午夜福利在线观看吧| 在线观看日韩欧美| 1024视频免费在线观看| 欧美日韩视频精品一区| 亚洲欧美激情综合另类| 色尼玛亚洲综合影院| 看免费av毛片| 国产乱人伦免费视频| 黄色 视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 琪琪午夜伦伦电影理论片6080| av网站免费在线观看视频| 国产欧美日韩一区二区精品| 中出人妻视频一区二区| 亚洲精品久久午夜乱码| 老司机深夜福利视频在线观看| 国产99白浆流出| 亚洲精品一区av在线观看| 久久天堂一区二区三区四区| 99国产综合亚洲精品| 十八禁人妻一区二区| 97超级碰碰碰精品色视频在线观看| 国产成人啪精品午夜网站| 窝窝影院91人妻| 国产乱人伦免费视频| 91字幕亚洲| 久热这里只有精品99| 国产一区二区在线av高清观看| 51午夜福利影视在线观看| 久9热在线精品视频| 少妇的丰满在线观看| 久久久国产一区二区| 免费观看精品视频网站| www日本在线高清视频| 国产99白浆流出| 国产在线观看jvid| 日韩三级视频一区二区三区| 一级,二级,三级黄色视频| 视频区图区小说| 美女扒开内裤让男人捅视频| 亚洲自偷自拍图片 自拍| 激情视频va一区二区三区| 乱人伦中国视频| av视频免费观看在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成av片中文字幕在线观看| 久久人妻福利社区极品人妻图片| 精品熟女少妇八av免费久了| 精品久久久久久成人av| 免费一级毛片在线播放高清视频 | 日韩欧美在线二视频| 我的亚洲天堂| 在线观看66精品国产| 91成人精品电影| a级片在线免费高清观看视频| 欧美日韩亚洲综合一区二区三区_| 欧美成狂野欧美在线观看| 国产区一区二久久| 丰满迷人的少妇在线观看| 美国免费a级毛片| 亚洲成人免费电影在线观看| 欧美在线黄色| 少妇裸体淫交视频免费看高清 | 国产欧美日韩一区二区三| 亚洲精品一二三| 亚洲性夜色夜夜综合| 久久九九热精品免费| 91麻豆av在线| 天堂俺去俺来也www色官网| 精品乱码久久久久久99久播| 人人妻人人澡人人看| 看黄色毛片网站| 新久久久久国产一级毛片| 中文字幕另类日韩欧美亚洲嫩草| www.精华液| 一级黄色大片毛片| 日韩欧美三级三区| 两性夫妻黄色片| 黄色毛片三级朝国网站| 久久久久亚洲av毛片大全| 亚洲性夜色夜夜综合| 亚洲第一欧美日韩一区二区三区| 黑人欧美特级aaaaaa片| 成人特级黄色片久久久久久久| 99国产精品免费福利视频| 一级作爱视频免费观看| 91精品国产国语对白视频| 日韩欧美国产一区二区入口| 99国产综合亚洲精品| 91麻豆av在线| 色哟哟哟哟哟哟| 欧美激情极品国产一区二区三区| 欧美成人免费av一区二区三区| 成年版毛片免费区| 18禁国产床啪视频网站| 久久天躁狠狠躁夜夜2o2o| 久久热在线av| 日韩免费av在线播放| 久久久精品欧美日韩精品| 国产成人欧美在线观看| 亚洲午夜精品一区,二区,三区| 精品人妻在线不人妻| 日韩大码丰满熟妇| 黄色a级毛片大全视频| 午夜视频精品福利| 午夜福利在线免费观看网站| 夫妻午夜视频| 999久久久国产精品视频| 俄罗斯特黄特色一大片| 欧美日韩亚洲国产一区二区在线观看| 老司机靠b影院| 淫秽高清视频在线观看| 十分钟在线观看高清视频www| 免费在线观看视频国产中文字幕亚洲| 久久精品国产清高在天天线| 免费不卡黄色视频| 日韩精品中文字幕看吧| 免费在线观看完整版高清| 久久香蕉激情| 热99国产精品久久久久久7| 在线观看一区二区三区| 久久精品91无色码中文字幕| 黄色毛片三级朝国网站| 亚洲成人精品中文字幕电影 | 看免费av毛片| 久久狼人影院| 操出白浆在线播放| 韩国av一区二区三区四区| 美女午夜性视频免费| 18美女黄网站色大片免费观看| 久久人妻av系列| 久久人人97超碰香蕉20202| 久久人妻熟女aⅴ| 亚洲男人的天堂狠狠| 亚洲三区欧美一区| 女生性感内裤真人,穿戴方法视频| 亚洲人成伊人成综合网2020| 欧美不卡视频在线免费观看 | 亚洲欧美日韩高清在线视频| 国产av在哪里看| 国产av精品麻豆| 欧美丝袜亚洲另类 | 亚洲成人国产一区在线观看| 精品国产一区二区久久| 亚洲自偷自拍图片 自拍| 在线永久观看黄色视频| 波多野结衣一区麻豆| 黑人操中国人逼视频| 国产xxxxx性猛交| 桃红色精品国产亚洲av| 免费一级毛片在线播放高清视频 | 满18在线观看网站| 成人影院久久| 国产av在哪里看| 亚洲精华国产精华精| 久久青草综合色| 精品久久久久久,| 最近最新中文字幕大全免费视频| 黄色毛片三级朝国网站| 一级a爱片免费观看的视频| 日韩高清综合在线| 色播在线永久视频| 日韩精品免费视频一区二区三区| 亚洲人成电影免费在线| 欧美 亚洲 国产 日韩一| 丰满的人妻完整版| 一进一出抽搐动态| 欧美国产精品va在线观看不卡| 精品国产美女av久久久久小说| 中文字幕人妻熟女乱码| 超碰成人久久| 久久久水蜜桃国产精品网| 久久久国产成人精品二区 | 国产精品秋霞免费鲁丝片| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 精品无人区乱码1区二区| 色老头精品视频在线观看| 国产精品二区激情视频| 亚洲免费av在线视频| 国产精品1区2区在线观看.| 在线天堂中文资源库| 欧美黄色淫秽网站| 午夜两性在线视频| 国产精品自产拍在线观看55亚洲| 丰满人妻熟妇乱又伦精品不卡| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 久久影院123| 国产欧美日韩一区二区三区在线| 免费人成视频x8x8入口观看| 国产精品野战在线观看 | 黄色a级毛片大全视频| 亚洲成人精品中文字幕电影 | 国产精品98久久久久久宅男小说| 亚洲九九香蕉| 免费观看精品视频网站| 国产不卡一卡二| av福利片在线| 又紧又爽又黄一区二区| 黑丝袜美女国产一区| www.熟女人妻精品国产| 欧美+亚洲+日韩+国产| 亚洲精华国产精华精| 丝袜美足系列| 国产精品一区二区三区四区久久 | 99国产精品一区二区蜜桃av| 在线国产一区二区在线| 他把我摸到了高潮在线观看| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 嫁个100分男人电影在线观看| 久久精品国产综合久久久| 91麻豆精品激情在线观看国产 | 一区二区日韩欧美中文字幕| 欧美丝袜亚洲另类 | 琪琪午夜伦伦电影理论片6080| 法律面前人人平等表现在哪些方面| 91字幕亚洲| 大香蕉久久成人网| 久久精品亚洲精品国产色婷小说| 国产三级黄色录像| 国产激情久久老熟女| 欧美+亚洲+日韩+国产| 视频在线观看一区二区三区| 日韩精品中文字幕看吧| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 搡老岳熟女国产| tocl精华| 欧美乱码精品一区二区三区| 久久久水蜜桃国产精品网| 亚洲国产中文字幕在线视频| 午夜影院日韩av| 如日韩欧美国产精品一区二区三区| www日本在线高清视频| 国产一区二区激情短视频| 成人18禁在线播放| 亚洲欧美日韩无卡精品| 精品熟女少妇八av免费久了| 亚洲精品在线美女| 久久久久久久精品吃奶| 国产成人精品久久二区二区免费| 日本精品一区二区三区蜜桃| 国产aⅴ精品一区二区三区波| 久久国产乱子伦精品免费另类| 久久久国产成人免费| 欧美不卡视频在线免费观看 | 亚洲精品一卡2卡三卡4卡5卡| 久久久精品欧美日韩精品| 人人澡人人妻人| 精品国产一区二区三区四区第35| 久久久久久亚洲精品国产蜜桃av| 亚洲全国av大片| 日韩大码丰满熟妇| 亚洲国产精品合色在线| 亚洲狠狠婷婷综合久久图片| 国产精品野战在线观看 | 在线国产一区二区在线| 国产成年人精品一区二区 | 啦啦啦 在线观看视频| 男女高潮啪啪啪动态图| 丁香欧美五月| 很黄的视频免费| 老汉色∧v一级毛片| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 五月开心婷婷网| xxx96com| 国产欧美日韩一区二区三| 91麻豆av在线| 免费少妇av软件| 日韩一卡2卡3卡4卡2021年| 色哟哟哟哟哟哟| 在线观看午夜福利视频| 女性被躁到高潮视频| 一个人免费在线观看的高清视频| 黄色毛片三级朝国网站| 人人妻人人爽人人添夜夜欢视频| 天堂影院成人在线观看| 国产主播在线观看一区二区| 人妻久久中文字幕网| 乱人伦中国视频| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美激情在线| 亚洲五月天丁香| 日韩国内少妇激情av| 久久精品人人爽人人爽视色| 看黄色毛片网站| 久久久水蜜桃国产精品网| 丰满迷人的少妇在线观看| 亚洲精品一区av在线观看| 在线观看免费午夜福利视频| 亚洲国产欧美日韩在线播放| 真人一进一出gif抽搐免费| 欧美不卡视频在线免费观看 | 亚洲欧美激情在线| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 99精国产麻豆久久婷婷| 日日爽夜夜爽网站| 99国产精品免费福利视频| 中文字幕最新亚洲高清| 黄色女人牲交| 高清在线国产一区| 婷婷六月久久综合丁香| 成人三级做爰电影| 丁香欧美五月| 欧美日韩亚洲国产一区二区在线观看| 国产又色又爽无遮挡免费看| 亚洲精品国产一区二区精华液| 国产精品久久久av美女十八| 精品午夜福利视频在线观看一区| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 久久国产乱子伦精品免费另类| 日本vs欧美在线观看视频| www.自偷自拍.com| 人人妻人人澡人人看| 国产精品综合久久久久久久免费 | aaaaa片日本免费| 三上悠亚av全集在线观看| 如日韩欧美国产精品一区二区三区| 日韩中文字幕欧美一区二区| 香蕉国产在线看| 黑人巨大精品欧美一区二区mp4| cao死你这个sao货| 曰老女人黄片| 母亲3免费完整高清在线观看| 欧美国产精品va在线观看不卡| 欧美日韩乱码在线| 91九色精品人成在线观看| 久久久久久大精品| 91精品国产国语对白视频| 亚洲精品粉嫩美女一区| 国产99久久九九免费精品| 精品乱码久久久久久99久播| 亚洲欧美日韩高清在线视频| 99精品欧美一区二区三区四区| 两个人看的免费小视频| 午夜福利,免费看| 国产人伦9x9x在线观看| 另类亚洲欧美激情| 色尼玛亚洲综合影院| 久久精品国产亚洲av香蕉五月| 熟女少妇亚洲综合色aaa.| 国产又色又爽无遮挡免费看| 免费av中文字幕在线| 国产主播在线观看一区二区| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 久久久久精品国产欧美久久久| 久久亚洲精品不卡| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 日韩欧美在线二视频| 久热这里只有精品99| 久久伊人香网站| 国产蜜桃级精品一区二区三区| 精品久久久久久成人av| 80岁老熟妇乱子伦牲交| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三| 国产精品一区二区在线不卡| 精品乱码久久久久久99久播| 99精品在免费线老司机午夜| 五月开心婷婷网| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女 | 亚洲精华国产精华精| 一级a爱视频在线免费观看| 亚洲五月天丁香| 精品少妇一区二区三区视频日本电影| 91国产中文字幕| 亚洲精品在线美女| 麻豆国产av国片精品| 精品电影一区二区在线| 国产麻豆69| 91成人精品电影| 可以免费在线观看a视频的电影网站| 天天影视国产精品| 午夜成年电影在线免费观看| 午夜福利,免费看| 免费观看精品视频网站| 久久精品人人爽人人爽视色| 欧美 亚洲 国产 日韩一| 最近最新中文字幕大全免费视频| 久久国产精品影院| 在线观看免费日韩欧美大片| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 巨乳人妻的诱惑在线观看| 亚洲avbb在线观看| 日韩欧美一区二区三区在线观看| 人人妻人人添人人爽欧美一区卜| 90打野战视频偷拍视频| 久久精品亚洲精品国产色婷小说| 久久久久久久久久久久大奶| 热re99久久精品国产66热6| 波多野结衣一区麻豆| 又紧又爽又黄一区二区| 欧美成人性av电影在线观看| 麻豆一二三区av精品| aaaaa片日本免费| e午夜精品久久久久久久| 热99国产精品久久久久久7| 日韩欧美一区二区三区在线观看| 校园春色视频在线观看| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| 97碰自拍视频| 麻豆国产av国片精品| 亚洲精品成人av观看孕妇| 国产精品久久久久久人妻精品电影| 国产伦一二天堂av在线观看| 国产精品九九99| 50天的宝宝边吃奶边哭怎么回事| 日韩欧美一区二区三区在线观看| 欧美中文日本在线观看视频| 丝袜人妻中文字幕| 亚洲色图av天堂| 99riav亚洲国产免费| 9热在线视频观看99| 又大又爽又粗| 大型黄色视频在线免费观看| 久久草成人影院| 国产精品久久久人人做人人爽| 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 久久久久久久久中文| 一夜夜www| 亚洲欧美精品综合一区二区三区| 后天国语完整版免费观看| 99久久人妻综合| 日本撒尿小便嘘嘘汇集6| 国产在线观看jvid| 老司机福利观看| av视频免费观看在线观看| 黄色怎么调成土黄色| 午夜福利在线免费观看网站| 手机成人av网站| 99riav亚洲国产免费| 国产欧美日韩一区二区三区在线| 中文字幕色久视频| 精品人妻在线不人妻| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 久久久久国产精品人妻aⅴ院| 18禁裸乳无遮挡免费网站照片 | 婷婷丁香在线五月| 久久中文看片网| 真人做人爱边吃奶动态| 天堂中文最新版在线下载| 一级a爱视频在线免费观看| 一区在线观看完整版| 亚洲伊人色综图| 精品久久久久久成人av| 老汉色av国产亚洲站长工具| 91九色精品人成在线观看| 电影成人av| 999精品在线视频| 很黄的视频免费| 精品一品国产午夜福利视频| 精品国内亚洲2022精品成人| 丝袜美腿诱惑在线| 欧美成人午夜精品| 精品无人区乱码1区二区| 男女下面进入的视频免费午夜 | 麻豆国产av国片精品| 亚洲精品一区av在线观看| 午夜成年电影在线免费观看| 欧美激情极品国产一区二区三区| 老司机亚洲免费影院| 大香蕉久久成人网| 欧美日本亚洲视频在线播放| 国内久久婷婷六月综合欲色啪| 久久久水蜜桃国产精品网| 中文字幕精品免费在线观看视频| www.www免费av| 韩国精品一区二区三区| 国产精品 欧美亚洲| 欧美黑人欧美精品刺激| av电影中文网址| www.www免费av| 成人国产一区最新在线观看| 午夜成年电影在线免费观看| 国产三级黄色录像| 亚洲色图av天堂| av福利片在线| 在线观看一区二区三区| 亚洲精品av麻豆狂野| 满18在线观看网站| 成人av一区二区三区在线看| 99久久综合精品五月天人人| 免费高清视频大片| 亚洲国产精品合色在线| 中文字幕人妻熟女乱码| 在线观看一区二区三区| 丁香六月欧美| 老鸭窝网址在线观看| 三上悠亚av全集在线观看| 国产亚洲欧美精品永久| 国产一区二区激情短视频| a级毛片黄视频| 99香蕉大伊视频| 午夜福利欧美成人| 中文字幕av电影在线播放| 久久久久久久精品吃奶| 午夜精品国产一区二区电影| 亚洲国产中文字幕在线视频| 成人免费观看视频高清| 亚洲va日本ⅴa欧美va伊人久久| 一区二区三区精品91| 亚洲美女黄片视频| 亚洲在线自拍视频| 亚洲欧美一区二区三区黑人| 999久久久精品免费观看国产| 日本撒尿小便嘘嘘汇集6| 午夜免费观看网址| 一本综合久久免费| 一区二区三区精品91| 无遮挡黄片免费观看| 精品久久久久久久久久免费视频 | 亚洲av片天天在线观看| 1024视频免费在线观看| 亚洲少妇的诱惑av| 久久国产精品人妻蜜桃| 男女之事视频高清在线观看| 老司机靠b影院| 熟女少妇亚洲综合色aaa.| 国产亚洲欧美在线一区二区| 久热爱精品视频在线9| 美女 人体艺术 gogo| 国产精品久久久人人做人人爽| 色播在线永久视频| 嫁个100分男人电影在线观看| 久久久国产欧美日韩av| 亚洲情色 制服丝袜| 亚洲精品中文字幕一二三四区| 成人av一区二区三区在线看| 欧美中文综合在线视频| 国产一区二区三区在线臀色熟女 | 日韩一卡2卡3卡4卡2021年| 国产91精品成人一区二区三区| 日韩免费高清中文字幕av| 青草久久国产| 亚洲第一av免费看| 国产1区2区3区精品| 一边摸一边抽搐一进一小说| 国产野战对白在线观看| 怎么达到女性高潮| 亚洲国产欧美一区二区综合| 国产男靠女视频免费网站| 精品久久久久久久久久免费视频 | 精品久久久精品久久久| 久久午夜亚洲精品久久| 啦啦啦免费观看视频1| 91精品国产国语对白视频| 在线观看午夜福利视频| 在线视频色国产色| 日日爽夜夜爽网站| 久久久久国内视频| 亚洲精品一区av在线观看| 夜夜躁狠狠躁天天躁| 狠狠狠狠99中文字幕| 国产精品综合久久久久久久免费 | 一级毛片高清免费大全| 一级,二级,三级黄色视频| 久久人人精品亚洲av| 看免费av毛片| 久久精品91无色码中文字幕| 精品久久久久久,| 日日爽夜夜爽网站| www.www免费av| 一级a爱片免费观看的视频| 看片在线看免费视频| 一边摸一边抽搐一进一小说| 国产成人系列免费观看| 一级毛片高清免费大全| 在线观看免费日韩欧美大片| 一a级毛片在线观看| 午夜91福利影院| 一区在线观看完整版| 午夜福利一区二区在线看| 在线观看舔阴道视频| 波多野结衣一区麻豆| av视频免费观看在线观看| 又黄又粗又硬又大视频| av福利片在线| 宅男免费午夜| 热99国产精品久久久久久7| 国产成人av教育| 国产精品二区激情视频| av免费在线观看网站| 久久精品国产99精品国产亚洲性色 | 一a级毛片在线观看|