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

    南海東沙海域內(nèi)孤立波形態(tài)研究

    2017-01-10 06:12:01徐宋昀許惠平耿明會關(guān)永賢
    海洋學研究 2016年4期
    關(guān)鍵詞:環(huán)礁東沙波群

    徐宋昀,許惠平,耿明會,關(guān)永賢

    (1.同濟大學 海洋與地球科學學院,海洋地質(zhì)國家重點實驗室,上海200092;2.中國科學院 深海科學與工程研究所,海南 三亞572000;3.國土資源部 海底礦產(chǎn)資源重點實驗室,廣州海洋地質(zhì)調(diào)查局,廣東 廣州510760)

    ?

    南海東沙海域內(nèi)孤立波形態(tài)研究

    徐宋昀1,許惠平2,耿明會1,關(guān)永賢3

    (1.同濟大學 海洋與地球科學學院,海洋地質(zhì)國家重點實驗室,上海200092;2.中國科學院 深??茖W與工程研究所,海南 三亞572000;3.國土資源部 海底礦產(chǎn)資源重點實驗室,廣州海洋地質(zhì)調(diào)查局,廣東 廣州510760)

    本文基于反射地震數(shù)據(jù)和MODIS遙感數(shù)據(jù),對南海東沙海域內(nèi)孤立波及孤立波群的形態(tài)有了系統(tǒng)的認識。內(nèi)孤立波存在上升型和下降型兩種極性波,又根據(jù)波形分成“鐘形”、“平底形”和“碗形”三種類型,孤立波在波谷處的形狀與孤立波振幅有關(guān)。單個內(nèi)孤立波在傳播一段時間后,受到各種因素的影響,會發(fā)育成波群。東沙海域的不同位置存在兩種波群:“有序型”波群和“復雜型”波群,結(jié)合它們的位置及波群傳播過程,認為這兩種波群可能是孤立波從深海向陸架的整個傳播過程中的兩個階段,“有序型”波群在被東沙島阻礙后,受到各種海底地形、東沙環(huán)礁、波-波相互作用的影響,轉(zhuǎn)變成“復雜型”波群。

    東沙環(huán)礁;內(nèi)孤立波群;波形;傳播過程

    0 引言

    內(nèi)波屬于海洋的中尺度環(huán)節(jié),它造成的強烈混合[1]對海洋的物質(zhì)能量交換、海洋環(huán)境及活動在海洋內(nèi)部的物體有重要的影響[2]。東沙海域的大振幅內(nèi)孤立波是一種非線性內(nèi)波[3],海洋潮汐與起伏的海底地貌的相互作用、內(nèi)潮之間的非線性相互作用[4]和內(nèi)潮的非線性效應[5]是產(chǎn)生它們的主要原因。目前學者們普遍認為這些孤立波發(fā)源于呂宋海峽[6-8],而具體的生成過程非常復雜,存在多種可能的生成機制:背風波機制[9]、局地內(nèi)潮凹陷變陡形成孤立波[10-11]、黑潮的不穩(wěn)定性[12]、潮汐模式之間的能量耦合[13]等,而且內(nèi)孤立波的生成過程中,正壓潮流強度、水體層結(jié)程度、背景流的流向、呂宋海峽中存在的多個海脊等因素對產(chǎn)生的孤立波的振幅和傳播方向等都有重要的影響[14]。東沙海域的內(nèi)孤立波多數(shù)以孤立波群的形式傳播,波群中包含數(shù)個內(nèi)孤立波,首波的波峰線可達200 km[15-16],振幅達170 m[17-19],波群中的波長從9.7 km至500 m依次減小,波群的傳播速度大約為2.9 m/s[20],擁有很大的能量,可以傳播到幾百公里之外[21]。

    前人對東沙海域內(nèi)孤立波的研究主要集中在孤立波的傳播演變過程,水動力參數(shù)和孤立波特征參數(shù)的計算等方面。起初由于缺乏實測條件,無法給出孤立波的整個傳播過程, KORTEWEG和VRIES[22]提出了理論描述,即KdV方程,在兩層模型[23]的條件下模擬了內(nèi)孤立波從深海到淺海的傳播。隨著對孤立波研究的深入,為使模型更好地符合真實情況,發(fā)展出了二階KdV[24]等演化方程,從理論上展示了孤立波在傳播過程中波形變化與振幅、寬度的關(guān)系[25],極性轉(zhuǎn)變[26]和破碎[27]等。與此同時,對內(nèi)孤立波的觀測手段逐漸增多,許多研究者利用錨系、遙感和地震方法觀測,研究南海北部內(nèi)孤立波的分布情況并計算相關(guān)參數(shù)。曾侃[28]和范植松 等[29]在兩層流體模型的條件下結(jié)合KdV方程反演了振幅和波速。甘錫林 等[30]提出用傅里葉分析方法和經(jīng)驗模態(tài)方法進一步反演孤立波的波長和波向。董崇志 等[31]研究分析了地震數(shù)據(jù)后認為,地震疊加剖面上同相軸呈現(xiàn)的起伏變化反映了內(nèi)波的總體形態(tài)。WANG et al[32]和孫麗娜 等[33]還通過遙感圖像宏觀分析了內(nèi)波的分布情況。雖然前人對KdV模型的運用、特征參數(shù)的計算進行了大量研究,也是對孤立波形態(tài)的一定描述,但對于實際觀測到的內(nèi)孤立波形狀和波群形態(tài)的研究很少。

    本文的研究區(qū)域東沙海域位于南海北部陸架陸坡,水深在50~2 500 m之間,海底坡度由陸架到陸坡逐漸變大,地形比較平坦,是潮、浪、內(nèi)孤立波相互作用強烈的區(qū)域,而位于20.7°N,116.8°E的東沙島及其環(huán)礁,對內(nèi)孤立波向淺海傳播過程有重要的影響。圖1是2009年6月—9月從MODIS圖像中觀測到的東沙海域內(nèi)孤立波分布情況,從孤立波的密集程度可以看出這幾個月東沙海域的孤立波出現(xiàn)頻率較高,每天同一時段出現(xiàn)的位置幾乎相同。內(nèi)孤立波向西北傳播到東沙島附近時波峰線出現(xiàn)彎折,碰到東沙環(huán)礁后被分成南北兩部分繼續(xù)傳播,在東沙島的西北方向海域,南北兩部分內(nèi)孤立波產(chǎn)生交集。在這整個傳播過程中,內(nèi)孤立波與東沙環(huán)礁,波與波之間發(fā)生了許多復雜的相互作用,同時單個孤立波形狀和孤立波群形態(tài)還受到水深變化及海底地形的影響,呈現(xiàn)出多種形態(tài),因此東沙海域是研究內(nèi)孤立波形態(tài)的理想?yún)^(qū)域。

    本文主要研究東沙海域單個內(nèi)孤立波的形狀和孤立波群的形態(tài)及演變原因?;谇叭说难芯浚梅瓷涞卣饠?shù)據(jù)和MODIS遙感圖像,研究每個孤立波的幾何形態(tài),得到直觀的特征參數(shù);分析波群的整體形態(tài),比較不同類型波群的分布差異;探究內(nèi)孤立波群形態(tài)差異的影響因素。

    圖1 MODIS影像上觀測到的2009年6月—9月東沙海域內(nèi)孤立波分布情況Fig.1 Distribution of internal solitary waves in Dongsha are observed on MODIS images in June-September, 2009

    1 數(shù)據(jù)和方法

    1.1 地震數(shù)據(jù)處理

    地震剖面由廣州海洋地質(zhì)調(diào)查局提供,是2009年6月—10月采集的數(shù)據(jù)。原始數(shù)據(jù)經(jīng)過解編、觀測系統(tǒng)定義、噪聲衰減、共中心點(CMP)選排、速度分析、動校正和疊加得到疊加剖面。在此剖面上拾取反射同相軸信息,測量計算得到每個孤立波的特征參數(shù):寬度、振幅和半高寬,同時能直接觀察到孤立波的形狀。

    1.2 遙感數(shù)據(jù)處理

    MODIS數(shù)據(jù)從NASA網(wǎng)站上獲取,該數(shù)據(jù)有36個離散光譜波段,波段分辨率有250、500和1 000 m三種。選取2009年6月—9月每天未被云層遮擋且觀測到內(nèi)孤立波的影像,經(jīng)過輻射校正和幾何校正處理后,利用波段的真彩色合成更清晰地展示內(nèi)孤立波信號。內(nèi)孤立波信號在MODIS影像上表現(xiàn)為明暗相間的一系列條紋。由于遙感影像獲得的海表面信息受到風場、海面波浪、云層等的影響,在目視解譯內(nèi)孤立波及波群的分布時,真彩色合成圖像并不能非常清楚地顯示內(nèi)孤立波信號,因此需要對圖像進行高斯拉伸,處理后的遙感圖像會突出內(nèi)孤立波信號,以暗色條帶的形式顯示,比其他背景信號色彩深很多。

    1.3 二階KdV方程的應用

    目前二階KdV方程在數(shù)值模擬應用中是一個基礎方程,適用于模擬陸架上的大振幅內(nèi)波,它能從理論上比較好地描述單個孤立波的形狀,尤其是形狀特殊的波形,方程中有特定的參數(shù)能分別對應控制內(nèi)孤立波的極性和波谷處的形狀,使得理論模型更加貼近真實海水中的內(nèi)孤立波形狀。

    二階KdV方程[24]為:

    (1)

    式中:c是線性速度,η是振幅,β是非靜力頻散系數(shù)。還有兩個重要的參數(shù):二次非線性系數(shù)α和三次非線性系數(shù)α1,對孤立波的形狀有決定作用。在兩層流體模型的條件下,兩個非線性系數(shù)是[34]:

    (2)

    式中:h1和h2分別是上、下兩層流體的厚度,ρ1和ρ2分別是上、下兩層流體的密度。從式(2)中可以明顯看出α1<0,α的正負與兩層流體的厚度有關(guān)。

    根據(jù)方程能得到對應于孤立波的振幅解[35]:

    (3)

    式中:1/Γ是孤立波的有效寬度,V是相速度,A 和B是表示V、c、α和α1之間相關(guān)性的系數(shù)。孤立波的極性由二次非線性系數(shù)α決定,特別在兩層流體模型中,如果密度躍層靠近海面,即h1h2,孤立波則是上升型波。更重要的是孤立波的振幅η受到α和α1的限制,其絕對值不能超過以下極值[35]:

    (4)

    這個極值還影響到內(nèi)孤立波在波谷處的形狀,當孤立波的振幅越接近于這個極值,波形就越趨向于矩形(圖2),同時振幅的極值接近于0.15H,H為水深。

    圖2 二次非線性系數(shù)為正數(shù)時孤立波的形狀Fig.2 Form of the soliton for positive quadratic nonlinear coefficient

    1.4 海水分層處理

    真實的海水密度、溫度和鹽度變化是連續(xù)的,不存在嚴格界線分明的層次,而數(shù)值模擬中應用的理論模型往往使用簡化的兩層流體模型。為了更好地運用理論模型研究真實情況的內(nèi)孤立波波形,需要對真實海水進行分層,使之符合兩層流體模型的條件。兩層流體模型的設定是:存在一條涇渭分明的分界線,將流體分為密度完全不同的上下兩層,忽略兩層流體密度之間的過渡變化。根據(jù)上述模型和曾侃[28]對海水密躍層深度和厚度的定義,本文把真實海水分為3層,把兩層流體之間的密度躍層分出作為中間一層,此密度躍層的上一層厚度設定為h1,密度躍層下一層的厚度設定為h2,即在應用兩層流體模型時忽略密躍層的厚度。圖3為某個內(nèi)孤立波所在位置的浮性頻率,示意了密躍層深度的定義:浮性頻率的最大值處對應的深度值。圖中h1對應于兩層流體模型中的上層流體厚度,h2為下層流體厚度,Hp為密躍層深度,dHp為密躍層的厚度,Nm為浮性頻率的最大值。

    圖3 浮性頻率中密躍層深度的定義示意圖Fig.3 Definition of pycnocline depth from the buoyancy frequency

    2 結(jié)果

    2.1 單個內(nèi)孤立波的形狀

    內(nèi)孤立波在形成之初是單個的孤立子,從研究區(qū)的反射地震同相軸圖中可以直接觀察到。內(nèi)孤立波的形狀多樣,存在兩種極性的波:下降型波和上升型波。從測量計算得到的眾多內(nèi)孤立波的特征參數(shù)可知,下降型波振幅和寬度相對較大,最小振幅為30 m,最大振幅可達80 m,寬度最大有1 600 m;上升型波的寬度和振幅很小,一般寬度只有560 m,振幅只有12 m。統(tǒng)計兩種極性波的數(shù)量可知,下降型波的數(shù)量遠遠大于上升型波,且上升型波幾乎只出現(xiàn)在水深較淺的陸架區(qū)域。下降型波和上升型波具有完全不同的特征,XU et al[36]曾研究過這兩種極性波在KdV模型適用性上的差異。也有前人研究認為這兩種極性是內(nèi)孤立波淺化過程中的兩種形態(tài)[21,37],下降型波向陸架傳播時,隨著水深變淺,密躍層上下兩層的厚度發(fā)生變化,波的形狀逐漸演變成上升型。將實際海水簡化分為3層后,上層厚度h1約為19~31 m,密躍層厚度dHp約為97~210 m,發(fā)現(xiàn)下降型波處的水深H最小的為320 m,則下層厚度(H-h1-dHp)約為79~104 m,從理論模型上看屬于h1h2的情況,不僅與KdV方程的理論相符合,也證明了波的極性變化受到水深的影響。而反射地震圖中能同時觀測到兩個極性波的存在,在一條同相軸上,水較深處是下降型波,水較淺的陸架處是上升型波,兩者之間波形比較復雜,不再有完整的孤立波,進一步驗證了下降型波在水深變淺時演變成上升型波。

    在觀測到的眾多下降型波中,波形差異也很大。根據(jù)所有孤立波的振幅、寬度和半高度的數(shù)值,把這3個參數(shù)分別分成大、中、小3個等級(表1)。并結(jié)合視覺形狀,把觀測到的下降型波分為3個類型(圖4):第一類是“鐘形”內(nèi)孤立波,這是孤立波的普遍形狀,數(shù)量最多,大都出現(xiàn)在地形較平坦區(qū)域,振幅和寬度一般為同一尺度,同為大或中尺度,視周期約為17 min;第二類是“平底形”內(nèi)孤立波,它的波谷處非常平坦,如同平頂建筑,數(shù)量極少,在所有數(shù)據(jù)中僅發(fā)現(xiàn)了2個,所在處水深分別為307和402 m,振幅和寬度也都為中尺度,視周期約為8 min;第三類是“碗形”內(nèi)孤立波,它的波谷處呈大圓弧形,振幅非常小,一般只有30 m,而寬度非常大,約為1 530 m,兩個參數(shù)形成巨大的反差,視周期約10 min,這類波的數(shù)量也不多,所處的水深一般大于750 m。這3個類型的波僅從波形特征上進行了分類,傳播方向基本一致,都向西或西北方向傳播,它們的視周期是在觀測船上測到的相對周期。

    3種類型的孤立波形狀均可用二階KdV方程進行描述,公式(2)和(3)給出了兩層流體條件下的理論模型。“平底形”內(nèi)孤立波用二階KdV方程能更好地研究其波形特征,由公式(4)知,此種孤立波的振幅和形狀由兩個非線性系數(shù)決定。以圖4c中的孤立波為例,該處的上層深度為22.8 m,密躍層厚度為156.3 m,水深402 m,從理論上推斷此孤立波的振幅不能超過|α/α1|=50.4 m;反射地震數(shù)據(jù)測得的實際最大振幅為47.5 m,兩個數(shù)據(jù)非常接近,η/ηlim=0.94,對比圖2中振幅與波形的變化,理論推測此波的波形已經(jīng)開始向“矩形”演變,這與地震數(shù)據(jù)實際觀測到的波形相符合。此前對“平底形”內(nèi)孤立波僅有理論上的描述研究,這是首次從地震數(shù)據(jù)的實際觀測中展示這種波形的存在。

    另一種特殊波形“碗形”內(nèi)孤立波在形狀上與內(nèi)孔(internal bore)相似,從生成機制和傳播上比較,兩者雖然都屬于內(nèi)重力波,但還是有明顯的差別?!巴胄巍辈ㄊ且环N內(nèi)孤立波,傳播過程中形狀穩(wěn)定;內(nèi)孔表現(xiàn)為溫躍層的逐步變化,經(jīng)常伴有類似脈沖的短周期振蕩[38],且尺度比內(nèi)孤立波大。內(nèi)孔的傳播可以用兩層Boussinesq模型描述,這個模型還能描述出特殊情況:內(nèi)孔破碎后可能會逐漸演變成內(nèi)孤立波[39]。

    表1 3個特征參數(shù)的分級

    圖4 3類下降型波的示意圖和地震剖面圖Fig.4 Schematic diagram and seismic sections of three kinds of depression waves

    2.2 內(nèi)孤立波群的形態(tài)

    單個內(nèi)孤立波在傳播過程中逐漸淺化,受到各種環(huán)境因素的影響,傳播進入東沙海域時大都演變?yōu)楣铝⒉ㄈ旱男螒B(tài),結(jié)合遙感圖像和地震數(shù)據(jù)的觀測結(jié)果,在東沙海域的不同位置發(fā)現(xiàn)兩種典型的波群:一種是波形有序的孤立波群,包含了多個波形完整的孤立波,波形相似、類型一致(圖5a);另一種是更復雜的孤立波群,包含了多個波長,波形無序、不規(guī)則變化(圖5b)。

    圖5a中的“有序型”孤立波群中包含了7個完整的內(nèi)孤立波,每個孤立波形狀完整清晰,均屬于下降型中的“鐘形”孤立波。①號波是此波群的首波,在波群中具有最大的振幅70 m,且前3個孤立波的振幅和波長依次減小,④~⑦號波的波長很小,振幅的大小變化無序。圖5b中的“復雜型”孤立波群包含的信息非常復雜,整個波群的形態(tài)是由各種波形疊加而成,其中完整波形的內(nèi)孤立波類型多樣,同時存在下降型和上升型內(nèi)孤立波。此處的波群已經(jīng)傳播到非常淺的陸架區(qū)域,孤立波的振幅和寬度均很小,在海水深度200 m處存在上升型波,振幅大約為30 m,在海水深度283 m對應的水層中存在下降型波,振幅約為25 m。波群中還存在處于淺化過程中的孤立波,波形并不完整,此外,下降型波的左側(cè)存在類似尾流的小波動。

    圖5 東沙海域兩種典型內(nèi)孤立波群的地震剖面圖Fig.5 Seismic sections of two typical internal solitary wave packets in Dongsha area

    從遙感影像中也能清晰觀測到兩種波群的存在,上述兩種波群的地震數(shù)據(jù)采集時間有相對應的MODIS影像,圖6a是UTC 2009年7月7日02:55的影像,L1是整條地震測線的位置,呈SE—NW向,紅色部分是圖5b所在的位置,受到天氣影響,地震數(shù)據(jù)同一天的MODIS影像受到云層遮擋,無法觀測到內(nèi)孤立波信號,所以離地震數(shù)據(jù)時間最近的MODIS影像比地震數(shù)據(jù)早3 d零1 h。由于內(nèi)孤立波傳播的規(guī)律性,認為每天在同一位置的孤立波波群應該相似,同時根據(jù)測線位置和成像時間,可以認為該處已有的波群Ⅰ與圖5b的“復雜型”波群類似,同時可以推斷此處觀測到的波群可能會受到地形、東沙環(huán)礁的影響而更加復雜。圖6b是UTC 2009年7月29日05:30的影像,L2是整條測線位置,其中的紅色部分是圖5a所在的位置,地震數(shù)據(jù)采集到的“有序型”波群的時間比影像時間早2 h。按照內(nèi)孤立波一般傳播速度估計,地震采集到的波群為Ⅱ處的波,測線L2為SW—NE向,與內(nèi)孤立波的傳播方向垂直,這可能用來解釋波群中④~⑦號波振幅的無序變化:圖5a的位置位于測線的前段,水深為300~400 m,水深變化較緩,逐漸變淺,根據(jù)測線和波群的波峰線之間的夾角,觀測到的孤立波不在同一條傳播方向線上,使得地震剖面中顯示的振幅不按大小有序排列。

    圖6 東沙海域的MODIS影像Fig.6 MODIS images of Dongsha area

    根據(jù)上述地震剖面數(shù)據(jù)可知,Ⅰ處的波群為“復雜型”,Ⅱ處的波群為“有序型”,再結(jié)合遙感影像上波群信號的特征,期望能直接從遙感影像上分辨兩種波群類型。對比兩種波群的遙感信號可知,Ⅰ處的“復雜型”波群在遙感圖像上條紋縱橫交叉,明暗不一,波長無規(guī)則變化;Ⅱ處的“有序型”波群在遙感圖像上條紋清晰,波峰線有序排列,波長依次減小。按照這種特征對6月—9月的遙感數(shù)據(jù)進行統(tǒng)計后,發(fā)現(xiàn)“有序型”波群多存在于東沙環(huán)礁的東部,最遠在119.5°E附近也發(fā)現(xiàn)了“有序型”波群,這些孤立波群應該是由東南部深海盆向東沙臺地傳播的內(nèi)孤立波演變而來的;“復雜型”波群多存在于東沙島的西北部,波群之間常?;ビ薪患?。

    3 討論

    東沙海域的內(nèi)孤立波往往發(fā)源于呂宋海峽,而且從遙感圖像上看出兩種類型的波群都向西北傳播,波群的形態(tài)卻完全不同,綜合分析從遙感數(shù)據(jù)上獲得的兩種波群的分布位置差異、傳播過程以及從地震數(shù)據(jù)顯示的海底地形,波群形態(tài)可能受到以下3個因素的影響。

    3.1 海底地形

    內(nèi)孤立波自產(chǎn)生向西或西北方向傳播時,途徑南海東北部的深海盆、陸坡區(qū)域、東沙臺地、陸架區(qū)域,海水深度由4 000 m減小為200 m(圖1)。深水海盆區(qū)域地形平緩,坡度差異不大,內(nèi)孤立波傳播時幾乎沒有受到任何阻礙,且孤立波離海底較遠,即使有局部的海山,對個別孤立波波形產(chǎn)生的影響很微小[40],所以在傳播過程中內(nèi)孤立波形態(tài)幾乎不變;內(nèi)孤立波傳播到東沙海域時,海水深度驟減,孤立波受到海底地形影響變大,逐漸演變成“有序型”波群;陸坡至陸架區(qū)域有較小的坡度變化,總體地形起伏不平,在相對較短的一段距離中水下凸起多,水深較淺,波群在經(jīng)過水下凸起時,波谷處會產(chǎn)生微小變形,恢復之后形變的影響以孤立波尾波的形式存在[41],而海底深度變淺到一定程度時,部分內(nèi)孤立波處于向上升型波過渡階段,因此波群中波形信息復雜,波群形態(tài)不斷變化。

    3.2 東沙環(huán)礁

    東沙環(huán)礁對內(nèi)孤立波群的傳播過程產(chǎn)生了阻礙,對波群形態(tài)有重大的影響。從遙感圖像上看原波群的波峰線呈現(xiàn)向西凸出的圓弧形,中部波速最快,兩端較慢[42-43],而遇到東沙環(huán)礁后波峰線中部速度驟減,兩端波速不變,波峰線被分成南北兩個部分,且都屬于“復雜型”,靠近東沙環(huán)礁位置的一小段波峰線有很大的彎曲(圖1),而且南北兩部分波峰線不對稱,一般北部波峰線更長,弧度與原波群差不多,南部的波峰線很短,弧度不大。隨著波群向西推進,兩部分波峰線弧度逐漸增大,靠近東沙環(huán)礁的彎曲部分的彎曲程度也逐漸變大。同時波群被東沙島折射、衍射[44]后,靠近東沙環(huán)礁部分的波形發(fā)生巨大變化,有的孤立波甚至會破碎[37],南北兩端的波形變化較小,因此總體波群形態(tài)變得更加復雜。

    3.3 波-波相互作用

    內(nèi)孤立波群經(jīng)過東沙環(huán)礁以后,被分成兩部分的波群已經(jīng)形成兩個獨立且較完整的“復雜型”波群,兩個波群會在東沙島以西海域發(fā)生干涉、融合[44]等相互作用(圖6a)。內(nèi)孤立波的非線性使得波具有粒子性[45],波-波相互作用并非簡單的波形疊加之類的干涉,而是產(chǎn)生周相移動[46]等現(xiàn)象,相互作用還會產(chǎn)生更多的孤立子波,這些孤立子波又會融合形成大振幅孤立波,如圖1中最西側(cè)的幾條波峰線,中部的凹陷位置是原本兩個波群波峰線的交點,在繼續(xù)傳播的過程中凹陷處會逐漸變平,最終兩個波群又融合形成一條圓弧形波峰線。同一條波峰線上的波由于所處位置不同,波形變化不同,波長和波向的變化也不同,使得波群變得極為復雜。

    4 結(jié)論

    本文利用反射地震剖面和MODIS遙感影像,研究了東沙海域內(nèi)孤立波及波群的形態(tài)。通過反射地震剖面觀測到,內(nèi)孤立波存在上升型和下降型兩種極性波,根據(jù)寬度、振幅和半高寬等參數(shù)可以分成3種形狀類型,即“鐘形”、“平底形”和“碗形”,運用二階KdV方程詳細研究發(fā)現(xiàn)孤立波形狀與振幅有關(guān)。單個內(nèi)孤立波在傳播到東沙海域時演化為兩種波群,分別分布在不同的位置:包含多個波形完整、波形相似、類型一致孤立波的“有序型”孤立波群和波長、波形無序不規(guī)則變化的“復雜型”孤立波群,前者多存在于東沙環(huán)礁的東側(cè)海域,后者多存在于東沙環(huán)礁的西側(cè)及西偏北的區(qū)域。研究波群形態(tài)的影響因素之后發(fā)現(xiàn),“有序型”波群一般存在于未遇到東沙環(huán)礁之前的傳播過程中,“復雜型”波群一般存在于遇到東沙環(huán)礁之后。這兩種波群可能是孤立波從深海向陸架的整個傳播中的兩個階段,“有序型”波群在被東沙島阻礙后,受到各種復雜作用,轉(zhuǎn)變成“復雜型”波群。地震數(shù)據(jù)能從垂直方向上顯示內(nèi)孤立波的形狀,刻畫波形的細節(jié)特征,遙感圖像能從宏觀角度顯示內(nèi)孤立波的分布和群體特征,通過兩種數(shù)據(jù)在對應時間的相互印證和綜合分析,能更全面地描述內(nèi)孤立波及波群的形態(tài)特征,驗證了單個內(nèi)孤立波在淺化時發(fā)育形成波群,又受到海底地形、東沙環(huán)礁和波-波相互作用的影響逐漸變得復雜的傳播過程。

    [1] TIAN Ji-wei,LEI Zhou,ZHANG Xiao-qian. Latitudinal distribution of mixing rate caused by theM2internal tide[J].Journal of Physical Oceanography,2006,36(1):35-42.

    [2] LIU A K,CHANG Y S,HSU M K,et al. Evolution of nonlinear internal waves in the East and South China Seas[J]. Journal of Geophysical Research: Oceans (1978-2012),1998,103(C4):7 995-8 008.

    [3] CAI Shu-qun,XIE Jie-shuo,HE Jian-ling. An overview of internal solitary waves in the South China Sea[J]. Surveys in Geophysics,2012,33(5):927-943.

    [4] XIE Xiao-hui,SHANG Xiao-dong,CHEN Gui-ying. Nonlinear interactions among internal tidalwaves in the northeastern South China Sea[J].Chinese Journal of Oceanology and Limnology,2010,28(5):996-1 001.

    [5] XIE Xiao-hui,SHANG Xiao-dong,HAREN H V,et al. Observations of enhanced nonlinear instability in the surface reflection of internal tides[J].Geophysical Research Letters,2013,40(8):1 580-1 586.

    [6] RAMP S R,TANG T Y,DUDA T F,et al. Internal solitons in the northeastern South China Sea. Part I: Sources and deep water propagation[J]. IEEE Journal of Oceanic Engineering,2004,29(4):1 157-1 181.

    [7] ZHAO Zhong-xiang. A study of nonlinear internal waves in the northeastern South China Sea[M]. USA:ProQuest Dissertations and Theses,2005:1 280-1 285.

    [8] JACKSON C R,APEL J. An atlas of internal solitary-like waves and their properties[M]. Alexandria:Global Ocean Associates,2004:371-374.

    [9] SHAW P T, KO D S, CHAO S Y. Internal solitary waves induced by flow over a ridge: With applications to the northern South China Sea[J]. Journal of Geophysical Research,2009,114:1-11.

    [10] ZHAO Zhong-xiang,KLEMAS V,ZHENG Quan-an,et al. Remote sensing evidence for baroclinic tide origin of internal solitary waves in the north eastern South China Sea[J]. Geophysical Research Letters,2004,31(6):1-4.

    [11] ZHAO Zhong-xiang,ALFORD M H. Source and propagation of internal solitary waves in the northeastern South China Sea[J]. Journal of Geophysical Research,2006,111(C11):11-17.

    [12] YUAN Ye-li,ZHENG Quan-an,DAI De-jun,et al. Mechanism of internal waves in the Luzon Strait[ J]. Journal of Geophysical Research,2006,111(C11):11-17.

    [13] APEL J R,BADIEY M,CHIU Ching-Sang,et al. An overview of the 1995 SWARM shallow-water internal wave acoustic scattering experiment[J]. IEEE Journal of Oceanic Engineering,1997,22(3):465-500.

    [14] CAI Shu-qun,GAN Zi-jun,LONG Xiao-min. Features and evolution of internal solitarywaves in the northern South China Sea[J]. Chinese Science Bulletin,2001,46(15):1 245-1 250. 蔡樹群,甘子鈞,龍小敏.南海北部孤立子內(nèi)波的一些特征和演變[J].科學通報,2001,46(15):1 245-1 250.

    [15] LIU A K,RAMP S R,ZHAO Yun-he,et al. A casestudy of internal solitary wave propagation during ASIAEX 2001[J]. IEEE Journal of Oceanic Engineering,2004,29:1 144-1 156.

    [16] LI Xiao-feng,ZHAO Zhong-xiang,PICHEL W G. Internal solitary waves in the northwestern South China Sea inferred from satellite images[J]. Geophysical Research Letters,2008,35(13):1-7.

    [17] RAMP S R,YANG Y J,BAHR F L . Characterizing the nonlinear internal wave climate in the northeastern South China Sea[J]. Nonlinear Process in Geophysics,2010,17:481-498.

    [18] KLYMAK J M,PINKEL R,LIU Cho-teng,et al.Prototypical solitons in the South China Sea[J]. Geophysical Research Letters,2006,33(11):607.

    [19] BAI Yang,SONG Hai-bin,GUAN Yong-xian,et al. Nonlinear internal solitary waves in the northeast South China Sea near Dongsha Atoll using seismic oceanography[J]. Chinese Science Bulletin,2015,60:944-951. 拜陽,宋海斌,關(guān)永賢,等.利用地震海洋學方法研究南海東北部東沙海域內(nèi)孤立波的結(jié)構(gòu)特征[J].科學通報,2015,60:944-951.

    [20] ALFORD M H,LIEN R C,SIMMONS H,et al. Speed and evolution of nonlinear internal waves transiting the South China Sea[J]. Journal of Physical Oceanography,2010,40(6):1 338-1 355.

    [21] LIU A K,HO Chung-ru,LIU Chu-teng. Satellite remote sensing of South China Sea[M]. Taipei:Tingmao Publish Company,2008.

    [22] KORTWEG D J, VRIES D G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves[J]. Philosophical Magazine,1895,39:422-443.

    [23] CHEN Chen-yuan,HSU J R,CHEN Cheng-wu,et al. Wave propagation at the interface of a two-layer fluid system in the laboratory[J]. Journal of Marine Science and Technology,2007,15(1):8-16.

    [24] MICHALLET H,BARTHELEMY E. Experimental study of interfacial solitary waves[J]. Journal of Fluid Mechanics,1998,366:159-177.

    [25] XU Zhen-hua,YIN Bao-shu,HOU Yi-jun. Highly nonlinear internal solitary waves over the continental shelf of the northwestern South China Sea[J]. Chinese Journal of Oceanology and Limnology,2010,28:1 049-1 054.

    [26] XU Zhen-hua,YIN Bao-shu,HOU Yi-jun,et al. A study of internal solitary waves observed on the continental shelf in the northwestern South China Sea[J]. Acta Oceanologica Sinica,2010,29(3):18-25.

    [27] AGHSAEE P,BOEGMAN L,LAMB K G. Breaking of shoaling internal solitary waves[J]. Journal of Fluid Mechanics,2010,659:289-317.

    [28] ZENG Kan. Three aspects of studying oceanic internal waves by space-borne Synthetic Aperture Radar images[D]. Qingdao:Ocean University of China,2002. 曾侃.從衛(wèi)星合成孔徑雷達海表圖像研究海洋內(nèi)波的三個問題[D].青島:青島海洋大學,2002.

    [29] FAN Zhi-song,GENG Jian,ZHANG Yuan-ling,et al. A preliminary study of SAR remote sensing inversion of internal solitary waves in the north of the South China Sea[J]. Ocean University of China:Natural Science Edition,2005,35(6):885-888. 范植松,耿建,張遠凌,等.南海北部內(nèi)孤立波SAR遙感反演的初步研究[J].中國海洋大學學報:自然科學版,2005,35(6):885-888.

    [30] GAN Xi-lin,HUANG Wei-gen,YANG Jing-song,et al.A new method to extract internal wave parameters from SAR imagery with Hilbert-Huang transfrom[J]. Journal of Remote Sensing,2007,11(1):39-47. 甘錫林,黃韋艮,楊勁松,等.基于希爾伯特-黃變換的合成孔徑雷達內(nèi)波參數(shù)提取新方法[J].遙感學報,2007,11(1):39-47.

    [31] DONG Chong-zhi,SONG Hai-bin,HAO Tian-yao,et al.Studying of oceanic internal wave spectra in the Northeast South China Sea from seismic reflections[J]. Journal of Geophysics,2009(8):2 050-2 055. 董崇志,宋海斌,郝天珧,等.南海東北部海洋內(nèi)波的反射地震研究[J].地球物理學報,2009(8):2 050-2 055.

    [32] WANG Juan,HUANG Wei-gen,YANG Jing-song, et al. Study of the propagation direction of the internal waves in the South China Sea using satellite images[J]. Acta Oceanologica Sinica,2013,32(5):42-50.

    [33] SUN Li-na,ZHANG Jie,MENG Jun-min. Analysis of characteristics of spatial and temporal distributions of large amplitude internal waves in the South China Sea using MODIS[C]//Information technology in risk analysis and crisis response—Chinese disaster prevention association of professional committee of the risk analysis,Proceedings of the sixth annual meeting,F(xiàn)rance: Atlantis Press,2014:801-805. 孫麗娜,張杰,孟俊敏.基于MODIS數(shù)據(jù)的南海大振幅內(nèi)波時空分布特征分析[C]//風險分析和危機反應中的信息技術(shù)--中國災害防御協(xié)會風險分析專業(yè)委員會第六屆年會論文集,法國:Atlantis Press,2014:801-805.

    [34] APEL J R,OSTROVSKY L A,STEPANYANTS Y A,et al. Internal solitons in the ocean and their effect on underwater sound[J]. The Journal of the Acoustical Society of America,2007,121(2):695-722.

    [35] PELINOVSKII E N, POLUKHINA O E,LAMB K. Nonlinear internal waves in the ocean stratified in density and current[J]. Oceanology C/C of Okeanologiia,2000,40(6):757-766.

    [36] XU Zhen-hua,YIN Bao-shu,YANG Hong-wei,et al. Depression and elevation internal solitary waves in a two-layer fluid and their forces on cylindrical piles[J]. Chinese Journal of Oceanology and Limnology,2012,30:703-712.

    [37] OSBORNE A R,BURCH T L. Internal solitons in the Andaman Sea[J]. Science,1980,208:451-460.

    [38] HENYEY F S,HOERING A. Energetics of borelike internal waves[J]. Journal of Geophysical Research,1997,102(C2):3 323-3 330.

    [39] WANG Tao,GAO Tian-fu. Internal tides, solitary waves and bores in shallow seas[J]. Chinese Journal of Oceanology and Limnology,2001,19(2):103-111.

    [40] LI Qun,XU Zhen-hua,YIN Bao-shu,et al. Modeling the interaction of an internal solitary wave with a sill[J]. Acta Oceanologica Sinica,2015,34(11):32-37.

    [41] HSIEH C M,HWANG R R,HSU J R C,et al. Numerical modeling of flow evolution for an internal solitary wave propagating over a submerged ridge[J]. Wave Motion,2015,55:48-72.

    [42] WANG Cai-xia,PAWLOWICZ R. Propagation speeds of strongly nonlinear near-surface internal waves in the Strait of Georgia[J]. Journal of Geophysical Research: Oceans,2011,116(C10):1-15.

    [43] ZHANG Xiao-dong,WANG Jian,SUN Li-na,et al. Study on group velocity of oceanic internal solitary wave[Z]. France:Atlantis,2015.

    [44] ZHAO Yun-he.SAR remote sensing of nonlinear internal waves in the South China Sea[R].Md: Caelum Research Corp Silver Spring,2006.

    [45] BAI Ye-fei. Numerical study of the evolution process of internal solitary waves in two-layer fluids[D].Qingdao: Institute of Oceanology, Chinese Academy of Sciences,2005. 白曄斐.兩層流體界面孤立子內(nèi)波演變過程數(shù)值研究[D].青島:中國科學院研究生院(海洋研究所),2005.

    [46] HSU Ming-kuang,LIU A K,LIU Cheng. A study of internal waves in the China Seas and Yellow Sea using SAR[J]. Continental Shelf Research,2000,20(4):389-410.

    Study on the form of internal solitary wave in Dongsha area of the South China Sea

    XU Song-yun1, XU Hui-ping2, GENG Ming-hui1, GUAN Yong-xian3

    (1.StateKeyLaboratoryofMarineGeology,SchoolofOceanandEarthScience,TongjiUniversity,Shanghai200092,China;2.InstituteofDeep-seaScienceandEngineering,ChineseAcademyofSciences,Sanya572000,China;3.KeyLaboratoryofMarineMineralResource,MinistryofLandandResources,GuangzhouMarineGeologicalSurvey,Guangzhou510760,China)

    Based on multi-channel seismic profiles and MODIS images, the form of internal solitary wave (ISW) packets in Dongsha area of the South China Sea was studied. There are two polar waves: elevation and depression wave, and internal solitary waves can be divided into “bell shape”, “flat shape” and “bowl shape” according to the waveforms. The shape of wave trough is associated with amplitude. Influenced by various factors, the solitions will develop into packets within the propagation process. Two typical wave packets are observed at different locations in Dongsha area: rank-ordered ISW packet and more complex ISW packet, and combined with their locations and propagation, they are considered to be two stages of the propagation process from deep sea to shelf. A rank-ordered wave packet is influenced by submarine topography, Dongsha Reef and wave-wave interaction after Dongsha Island, then turns into a more complex wave packet.

    Dongsha Reef;internal solitary wave packet;waveform;propagation process

    10.3969/j.issn.1001-909X.2016.04.001.

    2016-09-24

    2016-11-24

    國家自然科學基金項目資助(41576047,41576054);國家自然科學基金重大計劃重點項目資助(91128205);中國科學院創(chuàng)新交叉團隊項目資助

    徐宋昀(1991-),女,浙江紹興市人,主要從事海洋遙感方面的研究。E-mail:xusongyun@163.com

    P731.24

    A

    1001-909X(2016)04-0001-09

    10.3969/j.issn.1001-909X.2016.04.001

    徐宋昀,許惠平,耿明會,等. 南海東沙海域內(nèi)孤立波形態(tài)研究[J]. 海洋學研究,2016,34(4):1-9,

    XU Song-yun, XU Hui-ping, GENG Ming-hui,et al. Study on the form of internal solitary wave in Dongsha area of the South China Sea[J].Journal of Marine Sciences,2016,34(4):1-9, doi:10.3969/j.issn.1001-909X.2016.04.001.

    猜你喜歡
    環(huán)礁東沙波群
    自行車
    基于ResNet與BiLSTM的心電QRS波群檢測方法
    《思考心電圖之166》答案
    《思考心電圖之162》答案
    密克羅尼西亞
    華夏地理(2018年3期)2018-08-06 03:59:24
    大洋洲的朗伊羅阿環(huán)礁
    印焦慮馬爾代夫向沙特賣礁
    玩具
    圣約瑟夫環(huán)礁
    中外文摘(2016年11期)2016-08-31 19:52:50
    求人問題
    国产熟女午夜一区二区三区| 亚洲,欧美精品.| a级毛色黄片| 亚洲综合精品二区| 日韩三级伦理在线观看| 国产极品粉嫩免费观看在线| 亚洲精品日韩在线中文字幕| 亚洲成av片中文字幕在线观看 | 全区人妻精品视频| 国产免费一级a男人的天堂| 国产日韩欧美亚洲二区| 亚洲情色 制服丝袜| 黑丝袜美女国产一区| 成年女人在线观看亚洲视频| 午夜福利视频在线观看免费| 天堂俺去俺来也www色官网| 国产色婷婷99| 午夜视频国产福利| 久久精品国产亚洲av涩爱| 另类亚洲欧美激情| 国产日韩欧美在线精品| 丝袜脚勾引网站| av在线播放精品| 一级毛片电影观看| 一级毛片 在线播放| 欧美精品国产亚洲| 亚洲欧美中文字幕日韩二区| 久久精品熟女亚洲av麻豆精品| 亚洲欧美精品自产自拍| 少妇的逼好多水| 精品久久久精品久久久| 侵犯人妻中文字幕一二三四区| 在线免费观看不下载黄p国产| 亚洲国产精品一区三区| 精品午夜福利在线看| 国产男人的电影天堂91| 国产精品一区二区在线不卡| √禁漫天堂资源中文www| 日韩成人av中文字幕在线观看| 一级黄片播放器| 欧美变态另类bdsm刘玥| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 日韩成人av中文字幕在线观看| 久久精品久久精品一区二区三区| 爱豆传媒免费全集在线观看| 久久精品久久久久久噜噜老黄| 91午夜精品亚洲一区二区三区| 纵有疾风起免费观看全集完整版| 这个男人来自地球电影免费观看 | 男的添女的下面高潮视频| 日韩精品有码人妻一区| 日韩伦理黄色片| 大话2 男鬼变身卡| 人体艺术视频欧美日本| 亚洲av电影在线进入| 亚洲成色77777| 国产精品 国内视频| 久久久国产一区二区| 久久久久网色| 欧美精品高潮呻吟av久久| 捣出白浆h1v1| 国产免费福利视频在线观看| 色网站视频免费| 99国产精品免费福利视频| 18+在线观看网站| 亚洲精品一二三| 亚洲人成网站在线观看播放| 日韩在线高清观看一区二区三区| 久久99精品国语久久久| 亚洲欧美日韩卡通动漫| 少妇人妻精品综合一区二区| 如何舔出高潮| 激情视频va一区二区三区| 亚洲色图综合在线观看| 亚洲av电影在线观看一区二区三区| 中国国产av一级| xxxhd国产人妻xxx| 国产免费视频播放在线视频| 啦啦啦啦在线视频资源| 国产成人精品在线电影| 大片电影免费在线观看免费| 国产精品一二三区在线看| 久久亚洲国产成人精品v| 91午夜精品亚洲一区二区三区| 色婷婷av一区二区三区视频| 国产亚洲午夜精品一区二区久久| 国产精品欧美亚洲77777| 各种免费的搞黄视频| 777米奇影视久久| 日韩av免费高清视频| 精品第一国产精品| 免费黄色在线免费观看| 国产探花极品一区二区| 国产麻豆69| 午夜91福利影院| 91成人精品电影| 国产有黄有色有爽视频| av播播在线观看一区| 国产麻豆69| 亚洲美女视频黄频| 日韩大片免费观看网站| 美女国产视频在线观看| 青春草亚洲视频在线观看| 视频中文字幕在线观看| 夜夜骑夜夜射夜夜干| 一级爰片在线观看| 国产黄色免费在线视频| 久久99热这里只频精品6学生| 中文字幕亚洲精品专区| 最近最新中文字幕大全免费视频 | 欧美精品一区二区免费开放| 亚洲欧美一区二区三区黑人 | 久久狼人影院| 人妻一区二区av| 日韩av不卡免费在线播放| 在线观看国产h片| 香蕉精品网在线| 国产一区二区三区综合在线观看 | www日本在线高清视频| 夜夜骑夜夜射夜夜干| 99九九在线精品视频| 成人无遮挡网站| 亚洲人与动物交配视频| 亚洲精品456在线播放app| 日本91视频免费播放| 久久精品国产亚洲av涩爱| 永久免费av网站大全| 国产精品秋霞免费鲁丝片| 丰满少妇做爰视频| 中国美白少妇内射xxxbb| 国产精品偷伦视频观看了| 国产国拍精品亚洲av在线观看| 综合色丁香网| 久久影院123| 午夜免费鲁丝| 国产成人a∨麻豆精品| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩另类电影网站| 乱码一卡2卡4卡精品| 欧美日韩一区二区视频在线观看视频在线| 欧美少妇被猛烈插入视频| 免费看av在线观看网站| 看十八女毛片水多多多| 婷婷色综合大香蕉| av又黄又爽大尺度在线免费看| 国产精品三级大全| 欧美少妇被猛烈插入视频| 精品少妇久久久久久888优播| 18禁在线无遮挡免费观看视频| 2021少妇久久久久久久久久久| 亚洲少妇的诱惑av| 女人久久www免费人成看片| 亚洲精品久久久久久婷婷小说| 免费观看a级毛片全部| 国产黄频视频在线观看| 最新中文字幕久久久久| 一级a做视频免费观看| 丰满少妇做爰视频| 久久精品国产鲁丝片午夜精品| 男女下面插进去视频免费观看 | 午夜福利视频在线观看免费| 免费大片黄手机在线观看| 国产免费现黄频在线看| 国产成人精品在线电影| 满18在线观看网站| 26uuu在线亚洲综合色| av天堂久久9| av卡一久久| 欧美丝袜亚洲另类| 亚洲成色77777| 免费久久久久久久精品成人欧美视频 | 人妻系列 视频| av有码第一页| 国产欧美日韩综合在线一区二区| 亚洲精品久久午夜乱码| 欧美 亚洲 国产 日韩一| 国产精品嫩草影院av在线观看| 性色av一级| av在线app专区| 尾随美女入室| 男女啪啪激烈高潮av片| 免费大片18禁| 在线观看一区二区三区激情| 男女下面插进去视频免费观看 | 女人被躁到高潮嗷嗷叫费观| 一级爰片在线观看| 18禁动态无遮挡网站| 在现免费观看毛片| 99热这里只有是精品在线观看| 精品国产一区二区三区四区第35| 天堂中文最新版在线下载| 26uuu在线亚洲综合色| 两个人看的免费小视频| 国产精品99久久99久久久不卡 | 亚洲欧美色中文字幕在线| 日产精品乱码卡一卡2卡三| 9191精品国产免费久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 女人久久www免费人成看片| 国产黄色免费在线视频| 免费高清在线观看视频在线观看| 国产一区二区在线观看av| 国产片内射在线| 日本爱情动作片www.在线观看| 一边摸一边做爽爽视频免费| 亚洲精品美女久久久久99蜜臀 | 丰满少妇做爰视频| 99热这里只有是精品在线观看| 亚洲精品成人av观看孕妇| 青春草国产在线视频| h视频一区二区三区| 国产成人精品无人区| 中文欧美无线码| 国产男人的电影天堂91| 在线亚洲精品国产二区图片欧美| 日韩人妻精品一区2区三区| 在线观看免费高清a一片| 亚洲精品久久成人aⅴ小说| 免费播放大片免费观看视频在线观看| 51国产日韩欧美| 欧美最新免费一区二区三区| 一级毛片 在线播放| 国产精品久久久久成人av| 97在线人人人人妻| 熟女av电影| 插逼视频在线观看| 亚洲精品视频女| 亚洲成国产人片在线观看| 国产成人av激情在线播放| 国产亚洲一区二区精品| 成人毛片60女人毛片免费| 校园人妻丝袜中文字幕| 丝瓜视频免费看黄片| 国产高清不卡午夜福利| 人人妻人人澡人人看| 国产激情久久老熟女| 啦啦啦中文免费视频观看日本| 在现免费观看毛片| 久久久国产欧美日韩av| 大陆偷拍与自拍| 如日韩欧美国产精品一区二区三区| 色哟哟·www| 亚洲精品日本国产第一区| tube8黄色片| 满18在线观看网站| 性色avwww在线观看| 欧美人与善性xxx| 大片电影免费在线观看免费| 免费黄频网站在线观看国产| 精品一区二区三区四区五区乱码 | 最黄视频免费看| 日韩大片免费观看网站| 视频中文字幕在线观看| 欧美成人午夜精品| 天天影视国产精品| 夫妻性生交免费视频一级片| 2018国产大陆天天弄谢| 国产亚洲一区二区精品| 少妇被粗大的猛进出69影院 | 一区二区日韩欧美中文字幕 | 女人精品久久久久毛片| 亚洲,一卡二卡三卡| 美女视频免费永久观看网站| 日韩人妻精品一区2区三区| 丰满少妇做爰视频| 亚洲成色77777| 青青草视频在线视频观看| 成年女人在线观看亚洲视频| 中国三级夫妇交换| 侵犯人妻中文字幕一二三四区| 国产探花极品一区二区| 99香蕉大伊视频| 国产国拍精品亚洲av在线观看| 一级a做视频免费观看| 毛片一级片免费看久久久久| 成人二区视频| 久久av网站| 母亲3免费完整高清在线观看 | 国产精品女同一区二区软件| 一级黄片播放器| 考比视频在线观看| 九草在线视频观看| a级毛片在线看网站| 免费不卡的大黄色大毛片视频在线观看| 日韩免费高清中文字幕av| 丝袜喷水一区| 大陆偷拍与自拍| 亚洲av.av天堂| 久久久久精品性色| 永久免费av网站大全| 久久精品国产综合久久久 | 91精品国产国语对白视频| 国产精品人妻久久久久久| 一区二区av电影网| 中文乱码字字幕精品一区二区三区| 精品人妻在线不人妻| 永久免费av网站大全| 免费在线观看黄色视频的| 校园人妻丝袜中文字幕| 欧美日韩av久久| 午夜视频国产福利| 亚洲精品乱码久久久久久按摩| 激情五月婷婷亚洲| 久久狼人影院| 国产成人免费观看mmmm| 啦啦啦视频在线资源免费观看| 蜜桃国产av成人99| 人妻 亚洲 视频| 寂寞人妻少妇视频99o| 中文精品一卡2卡3卡4更新| 97精品久久久久久久久久精品| 日韩视频在线欧美| av.在线天堂| 国产精品一区二区在线不卡| 国产又色又爽无遮挡免| 欧美精品av麻豆av| 一级爰片在线观看| 91aial.com中文字幕在线观看| 亚洲高清免费不卡视频| 久久精品国产鲁丝片午夜精品| 蜜桃在线观看..| 高清黄色对白视频在线免费看| 十分钟在线观看高清视频www| 午夜91福利影院| 久久精品熟女亚洲av麻豆精品| 欧美日韩一区二区视频在线观看视频在线| 国产男人的电影天堂91| 欧美丝袜亚洲另类| 91国产中文字幕| 在线免费观看不下载黄p国产| 边亲边吃奶的免费视频| 国产片特级美女逼逼视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区| 少妇人妻久久综合中文| 一级毛片我不卡| 欧美亚洲日本最大视频资源| 一区二区日韩欧美中文字幕 | 国产免费一级a男人的天堂| 日韩三级伦理在线观看| 精品久久久精品久久久| 国产成人一区二区在线| 色视频在线一区二区三区| 美女中出高潮动态图| 午夜av观看不卡| 欧美人与性动交α欧美软件 | 国产黄色免费在线视频| av又黄又爽大尺度在线免费看| 欧美人与性动交α欧美软件 | 大码成人一级视频| 水蜜桃什么品种好| 久久精品国产亚洲av涩爱| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美日韩在线播放| 国产极品粉嫩免费观看在线| 黑人巨大精品欧美一区二区蜜桃 | 亚洲第一av免费看| 亚洲内射少妇av| 久久狼人影院| 草草在线视频免费看| 国产69精品久久久久777片| 日本-黄色视频高清免费观看| 黑人欧美特级aaaaaa片| 国产爽快片一区二区三区| 亚洲久久久国产精品| 国产熟女午夜一区二区三区| 亚洲,一卡二卡三卡| 新久久久久国产一级毛片| av网站免费在线观看视频| 亚洲,一卡二卡三卡| 美女中出高潮动态图| 国产极品粉嫩免费观看在线| 高清欧美精品videossex| 国产伦理片在线播放av一区| 我的女老师完整版在线观看| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜制服| 亚洲情色 制服丝袜| 纯流量卡能插随身wifi吗| 考比视频在线观看| 日韩av不卡免费在线播放| 美女xxoo啪啪120秒动态图| 韩国av在线不卡| 国产在线视频一区二区| 有码 亚洲区| 国产精品国产三级国产专区5o| 赤兔流量卡办理| 成人国产av品久久久| 日韩av免费高清视频| 久久精品熟女亚洲av麻豆精品| 免费久久久久久久精品成人欧美视频 | 久久国产精品男人的天堂亚洲 | 在线免费观看不下载黄p国产| 尾随美女入室| 成人手机av| 啦啦啦视频在线资源免费观看| 91成人精品电影| 日产精品乱码卡一卡2卡三| 国产成人a∨麻豆精品| 少妇的逼水好多| av福利片在线| 久久久久久人人人人人| 久久精品久久久久久噜噜老黄| 熟妇人妻不卡中文字幕| 国产精品一国产av| a 毛片基地| 免费av不卡在线播放| xxx大片免费视频| 在线免费观看不下载黄p国产| 欧美激情 高清一区二区三区| 国产国语露脸激情在线看| 高清欧美精品videossex| 十八禁网站网址无遮挡| 欧美国产精品一级二级三级| 99热这里只有是精品在线观看| 美女脱内裤让男人舔精品视频| 亚洲精品久久午夜乱码| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 老司机影院毛片| 我要看黄色一级片免费的| 永久免费av网站大全| 少妇人妻精品综合一区二区| 日本色播在线视频| 97在线人人人人妻| 自拍欧美九色日韩亚洲蝌蚪91| 伊人久久国产一区二区| 五月伊人婷婷丁香| 国产精品一区www在线观看| 伦理电影大哥的女人| 国产永久视频网站| 久久久久久久精品精品| 在线观看人妻少妇| 在线观看免费日韩欧美大片| videossex国产| 久久精品国产亚洲av涩爱| av一本久久久久| 丁香六月天网| 午夜激情久久久久久久| 欧美激情 高清一区二区三区| 在线 av 中文字幕| 免费观看性生交大片5| 国产淫语在线视频| 国产成人免费观看mmmm| 欧美日韩av久久| 国产不卡av网站在线观看| 亚洲av福利一区| 久久精品国产亚洲av天美| 丝袜人妻中文字幕| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站在线观看播放| 精品国产乱码久久久久久小说| 一边摸一边做爽爽视频免费| 最近最新中文字幕大全免费视频 | 午夜老司机福利剧场| 久久人妻熟女aⅴ| 国产精品99久久99久久久不卡 | 91午夜精品亚洲一区二区三区| 久久久精品区二区三区| 精品一区二区免费观看| 午夜福利乱码中文字幕| 免费大片18禁| 国产伦理片在线播放av一区| 欧美少妇被猛烈插入视频| 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 欧美xxxx性猛交bbbb| 久久久久精品人妻al黑| 中文精品一卡2卡3卡4更新| 热re99久久精品国产66热6| 人体艺术视频欧美日本| 国产精品女同一区二区软件| 在线观看免费视频网站a站| 亚洲一级一片aⅴ在线观看| 97在线人人人人妻| 国国产精品蜜臀av免费| 免费大片黄手机在线观看| 最近手机中文字幕大全| 午夜日本视频在线| 九九爱精品视频在线观看| 国产男人的电影天堂91| 亚洲成av片中文字幕在线观看 | 热99国产精品久久久久久7| 18在线观看网站| 免费黄网站久久成人精品| 777米奇影视久久| 午夜福利,免费看| 精品国产国语对白av| 交换朋友夫妻互换小说| 晚上一个人看的免费电影| 各种免费的搞黄视频| 97超碰精品成人国产| 91在线精品国自产拍蜜月| 精品人妻熟女毛片av久久网站| 视频中文字幕在线观看| 欧美变态另类bdsm刘玥| 久久久久久伊人网av| 免费看av在线观看网站| 成人毛片60女人毛片免费| 少妇猛男粗大的猛烈进出视频| 色5月婷婷丁香| 一级片'在线观看视频| 久久精品国产a三级三级三级| 精品国产乱码久久久久久小说| 建设人人有责人人尽责人人享有的| 全区人妻精品视频| 美女福利国产在线| 九草在线视频观看| 久久久久网色| 国产极品天堂在线| av在线app专区| 欧美精品av麻豆av| 亚洲成人av在线免费| 少妇 在线观看| 搡女人真爽免费视频火全软件| 成年动漫av网址| 久久午夜综合久久蜜桃| 日韩成人伦理影院| 久久影院123| 欧美激情国产日韩精品一区| 在线观看国产h片| 国产色爽女视频免费观看| 亚洲欧美成人综合另类久久久| 香蕉精品网在线| 香蕉丝袜av| 日韩免费高清中文字幕av| 亚洲av电影在线进入| 卡戴珊不雅视频在线播放| 2018国产大陆天天弄谢| 水蜜桃什么品种好| 男女下面插进去视频免费观看 | h视频一区二区三区| 男女边吃奶边做爰视频| 午夜福利乱码中文字幕| 91在线精品国自产拍蜜月| 久久99一区二区三区| 欧美 日韩 精品 国产| 久久人人爽人人爽人人片va| 国产又色又爽无遮挡免| 久久影院123| 久久久久久久久久成人| 成人影院久久| 久久亚洲国产成人精品v| 欧美国产精品一级二级三级| h视频一区二区三区| 成人国语在线视频| 午夜福利视频在线观看免费| 色婷婷av一区二区三区视频| 国产欧美日韩一区二区三区在线| 大香蕉久久成人网| 亚洲国产欧美日韩在线播放| 日本av免费视频播放| 在线观看www视频免费| 国产精品偷伦视频观看了| 国产亚洲一区二区精品| 国产国语露脸激情在线看| 日韩电影二区| 91精品伊人久久大香线蕉| 久久久亚洲精品成人影院| 高清毛片免费看| 伦理电影免费视频| 最近中文字幕2019免费版| 韩国av在线不卡| 亚洲天堂av无毛| 18禁观看日本| 久久久久久伊人网av| 亚洲国产毛片av蜜桃av| 久久久久视频综合| av在线app专区| av片东京热男人的天堂| 丝瓜视频免费看黄片| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 久久99蜜桃精品久久| 蜜桃国产av成人99| 最黄视频免费看| 黑人巨大精品欧美一区二区蜜桃 | 热re99久久精品国产66热6| 日韩中文字幕视频在线看片| 久久久久国产网址| 亚洲欧美清纯卡通| 交换朋友夫妻互换小说| 国产日韩欧美视频二区| 亚洲美女黄色视频免费看| 国产亚洲一区二区精品| 婷婷色综合大香蕉| 男人舔女人的私密视频| 国产一区二区在线观看日韩| 亚洲精品456在线播放app| 亚洲色图 男人天堂 中文字幕 | 久久毛片免费看一区二区三区| 久久午夜综合久久蜜桃| 99国产精品免费福利视频| xxxhd国产人妻xxx| 亚洲经典国产精华液单| 午夜老司机福利剧场| 精品人妻偷拍中文字幕| 中文字幕免费在线视频6| 欧美精品av麻豆av| 亚洲第一区二区三区不卡| 亚洲国产成人一精品久久久| 五月玫瑰六月丁香| 女人久久www免费人成看片| 韩国av在线不卡| av又黄又爽大尺度在线免费看| 亚洲中文av在线| 日韩一区二区三区影片| 在线观看三级黄色| 成人亚洲精品一区在线观看| 欧美国产精品va在线观看不卡| 18在线观看网站| 18禁观看日本| 爱豆传媒免费全集在线观看|