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

徐宋昀,許惠平,耿明會,關(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.