基于最小二乘法和支持向量機(jī)的海洋內(nèi)孤立波傳播速度反演模型*

梁可達(dá) 劉滕飛 常哲 張猛 李志鑫 黃松松 王晶

(中國海洋大學(xué)物理與光電工程學(xué)院,青島 266100)

傳播速度是內(nèi)孤立波的重要參數(shù)之一,如何通過光學(xué)遙感圖像準(zhǔn)確快速地獲得內(nèi)孤立波速度是目前需要解決的重要科學(xué)問題.本文在實(shí)驗(yàn)室建立了模擬內(nèi)孤立波光學(xué)遙感成像系統(tǒng),開展了系列綜合實(shí)驗(yàn),獲取實(shí)驗(yàn)數(shù)據(jù)庫,并利用最小二乘法和支持向量機(jī)兩種方法分別建立基于單景光學(xué)遙感圖像的內(nèi)孤立波速度反演模型.利用中國南海的MODIS 圖像和GF-4 圖像數(shù)據(jù)對速度反演模型進(jìn)行精度檢驗(yàn).研究結(jié)果表明: 最小二乘法內(nèi)孤立波速度反演模型能夠給出回歸方程,物理意義更為直觀,且在300—399 m 水深范圍精度更高,而支持向量機(jī)內(nèi)孤立波速度反演模型在水深400—1200 m 和83—299 m 的范圍內(nèi)精度高.因此兩種內(nèi)孤立波速度反演模型各有優(yōu)勢,都可以應(yīng)用于真實(shí)海洋中內(nèi)孤立波速度的反演.

1 引言

海洋內(nèi)波是指在海水穩(wěn)定層化的海洋中產(chǎn)生的、最大振幅出現(xiàn)在海洋內(nèi)部的波動,是一種常見的海洋現(xiàn)象[1].海洋內(nèi)孤立波是周期短的非線性內(nèi)波,具有傳播的單向性和波形的不變性等特點(diǎn)[2].海洋內(nèi)孤立波會對海上石油、天然氣等資源的勘探以及開采設(shè)施的安全構(gòu)成破壞,也會改變水下聲信號傳播路線,對潛艇航行產(chǎn)生影響等等[3?6].同時海洋內(nèi)孤立波會影響海洋表面溫度的變化[7],對海洋的生態(tài)環(huán)境具有重要影響[8].因此,研究內(nèi)孤立波特性具有重要的科學(xué)意義和良好的應(yīng)用前景.

目前研究內(nèi)孤立波的方法有很多,其中利用衛(wèi)星遙感圖像是研究內(nèi)孤立波的重要方法之一[9,10].其中光學(xué)遙感具有時間、空間分辨率較高和圖像豐富等優(yōu)點(diǎn),常被學(xué)者用來研究內(nèi)孤立波的時空分布及傳播特性[11?13].海洋內(nèi)孤立波的特征參量主要包括半寬、振幅和速度等,這些參量的大小決定了內(nèi)孤立波在海洋混合中的作用量級,因此遙感圖像反演內(nèi)孤立波參數(shù)引起了研究人員的關(guān)注.而內(nèi)孤立波的速度作為重要的傳播特性,一直以來也被許多學(xué)者利用不同手段開展研究.目前測量內(nèi)孤立波傳播速度的常用方法主要是通過現(xiàn)場實(shí)測或利用遙感圖像計算內(nèi)孤立波速度.現(xiàn)場實(shí)測可以通過兩個或多個系泊裝置觀測內(nèi)孤立波計算傳播速度,該方法難度大、費(fèi)用昂貴.Alford[14]等利用2007 年的系泊數(shù)據(jù)研究了中國南海的非線性內(nèi)波的速度和演化.Huang 等[15]利用水下系泊陣列進(jìn)行了“南海內(nèi)波實(shí)驗(yàn)”.報道了2013 年12 月4 日的一次極端內(nèi)孤立波,其最大振幅為240 m,并對傳播速度進(jìn)行了分析.此外,呂海濱等[16]利用 2009 年 6 月24 日 15 時 40 分至 25 日 16 時 40 分“科學(xué)一號”考察船在東沙島東北部陸架上K106 站進(jìn)行的長達(dá) 25 h 的 X 波段雷達(dá)、溫度鏈、聲學(xué)多普勒流速剖面儀同步觀測,在該海域利用 Radon 變換技術(shù)獲取了本次觀測到的內(nèi)孤立波的傳播速度.雖然現(xiàn)場實(shí)測可以測量內(nèi)孤立波平均速度,但大規(guī)模大范圍的測量難于實(shí)現(xiàn).

遙感觀測手段測量內(nèi)孤立波速度主要是利用多時間圖像法,即采用同一天不同時刻的包含同一內(nèi)孤立波的遙感圖像,利用傳播路程與傳播時間的比值估算內(nèi)孤立波的傳播速度.Hong 等[17]利用衛(wèi)星數(shù)據(jù)估算南海淺水區(qū)內(nèi)孤立波相速度,基于環(huán)境衛(wèi)星高級合成孔徑雷達(dá)(ASAR)和中分辨率成像光譜儀(MODIS)在很短的時間間隔獲得不同的位置圖像,進(jìn)行內(nèi)孤立波速度計算.Meetei 等[18]利用2002 年12 月 至2016 年5 月 的ERS2 和Envi Sat 衛(wèi)星的不同圖像對和MODIS 的光學(xué)圖像,以及可見紅外成像輻射計來調(diào)查和估計安達(dá)曼海內(nèi)孤立波的傳播速度.Sun 等[19]利用MODIS Terra/Aqua 衛(wèi)星2014 年1 月至2018 年12 月采集的195對影像對估算安達(dá)曼海內(nèi)孤立波傳播速度.多時間圖像法計算內(nèi)孤立波速度相比現(xiàn)場實(shí)測手段更加便捷,但很難獲取某一海域任意時刻的兩景包含同一內(nèi)孤立波的光學(xué)遙感圖像,且該方法計算的速度為內(nèi)孤立波平均速度，這些都是該方法的不足之處.然而,發(fā)展一種只用單景圖像反演內(nèi)孤立波速度的方法顯然是非常重要的.

本文通過實(shí)驗(yàn)?zāi)M光學(xué)遙感內(nèi)孤立波成像,獲得了大量實(shí)驗(yàn)數(shù)據(jù),采用最小二乘法和支持向量回歸兩種方法建立利用單景光學(xué)遙感圖像反演內(nèi)孤立波速度的模型,為內(nèi)孤立波速度的研究和獲取提供了新方法.無疑,這對于海洋內(nèi)孤立波傳播特性的研究起到了推動作用.

2 實(shí)驗(yàn)設(shè)計及數(shù)據(jù)處理

2.1 實(shí)驗(yàn)設(shè)計

實(shí)驗(yàn)在長寬高分別為5 m×0.8 m×0.35 m 三維長直水槽中進(jìn)行.為模擬海洋中的層結(jié)現(xiàn)象,根據(jù)流體力學(xué)相似性原理,在水槽中鋪設(shè)上下兩層水體,下層水體模擬高鹽度海水,上層水體模擬低鹽度海水.長直水槽為玻璃材質(zhì),可從水槽側(cè)面觀測水體內(nèi)部的內(nèi)孤立波演化過程,從而獲取內(nèi)孤立波的半波寬度、振幅和速度等參數(shù).為了仿真光學(xué)遙感成像,將水槽中的水體表面視為海表面,在水槽一端上方放置LED 面光源,調(diào)整好角度,發(fā)出白光模擬太陽光;在水槽另一端上方放置CCD1,模擬衛(wèi)星的可見光傳感器,LED 光源發(fā)出的光照射到水體表面,經(jīng)反射后在CCD1 形成仿真光學(xué)遙感內(nèi)孤立波圖像;在水槽側(cè)面放置CCD2,與上下層水體的分層位置等高,用來采集記錄內(nèi)孤立波圖像,獲取內(nèi)孤立波特征參量.兩臺CCD 的視場邊界相同,并連接一臺電腦同步采集圖像,獲得了內(nèi)孤立波一一對應(yīng)的仿真光學(xué)遙感圖像和波形圖像,光學(xué)遙感圖像參數(shù)與內(nèi)孤立波特征參數(shù)完全匹配[20].實(shí)驗(yàn)平臺結(jié)構(gòu)如圖1 所示.

圖1 內(nèi)孤立波光學(xué)遙感仿真平臺結(jié)構(gòu)示意圖Fig.1.Schematic diagram of simulation platform of the internal solitary wave optical remote sensing.

實(shí)驗(yàn)采用重力塌陷法制造內(nèi)孤立波,即靠近水槽一側(cè)用擋板隔開的區(qū)域?yàn)樵觳▍^(qū),可以通過控制在造波區(qū)加入或抽取上層水體的多少來控制塌陷高度,獲得內(nèi)孤立波不同大小的振幅.針對不同水文參數(shù)和層結(jié)的海洋,根據(jù)流體力學(xué)相似性原理設(shè)計了不同水深、不同層結(jié)、不同密度差和不同塌陷高度的系列實(shí)驗(yàn),如表1 所列.

表1 內(nèi)孤立波實(shí)驗(yàn)設(shè)計表Table 1.Design table of internal solitary wave experiment.

2.2 數(shù)據(jù)處理

在采集的圖像中選取不同位置做時間序列圖,CCD1 采集圖像的時間序列圖模擬了光學(xué)遙感圖像,可從中提取亮暗間距和灰度差等遙感圖像特征參量.CCD2 采集圖像的時間序列圖可提取內(nèi)孤立波傳播速度、半波寬度和振幅等波要素參量.在提取同一組實(shí)驗(yàn)的參數(shù)時,根據(jù)標(biāo)尺選取表面與側(cè)面圖像相對應(yīng)的取樣線,利用MATLAB 程序做時間序列圖,如圖2 所示.圖2(a)為內(nèi)孤立波仿真光學(xué)遙感圖像,圖2(b)為內(nèi)孤立波波形圖,可提取內(nèi)孤立波振幅、半波寬度等參數(shù),圖2(c)為內(nèi)孤立波仿真光學(xué)遙感圖像灰度剖面圖,可提取灰度差和亮暗間距等遙感圖像參數(shù).

圖2 實(shí)驗(yàn)室內(nèi)孤立波參數(shù)示意圖 (a) 內(nèi)孤立波仿真遙感圖像;(b)內(nèi)孤立波波形圖;(c)內(nèi)孤立波灰度剖面圖Fig.2.Schematic diagram of the internal solitary wave parameters in the laboratory: (a) Simulated remote sensing images;(b) the internal solitary wave waveform diagram;(c) the gray scale profiles.

在實(shí)驗(yàn)中共提取1440 組實(shí)驗(yàn)數(shù)據(jù),每組實(shí)驗(yàn)數(shù)據(jù)包括上下層水深、上下層密度、振幅、速度、半波寬度、亮暗間距、灰度差等參數(shù).對獲得的1440 條實(shí)驗(yàn)數(shù)據(jù)進(jìn)行無量綱化處理,無量綱參量定義如下:

h1,h分別為水槽實(shí)驗(yàn)的上層水深和總水深,ρ1,ρ2分別為上層和下層水的密度,d,G1,G2分別為仿真遙感圖像內(nèi)孤立波條紋的亮暗間距、灰度最大值和灰度最小值,v為內(nèi)孤立波傳播速度,c0為線性相速度.建立了無量綱化的實(shí)驗(yàn)室內(nèi)孤立波數(shù)據(jù)樣本庫,利用樣本庫中的數(shù)據(jù)進(jìn)行內(nèi)孤立波速度反演模型的訓(xùn)練.

3 內(nèi)孤立波速度反演模型

從經(jīng)典的KdV 理論模型中可知內(nèi)孤立波傳播速度與水深、層結(jié)、密度差和振幅相關(guān),對于海洋內(nèi)孤立波的振幅可以通過現(xiàn)場實(shí)測和光學(xué)遙感圖像反演模型獲取,Wang 等[21]利用內(nèi)孤立波光學(xué)遙感圖像中條紋的灰度差和亮暗間距以及其他海洋水文參數(shù)反演了內(nèi)孤立波振幅.而本文是利用光學(xué)遙感圖像來進(jìn)行內(nèi)孤立波速度的反演,所以輸入?yún)⒘恐邪伺c內(nèi)孤立波振幅相關(guān)的光學(xué)遙感圖像特征參量亮暗間距和灰度差,以這兩項(xiàng)遙感圖像參數(shù)反映內(nèi)孤立波振幅的作用.因此建立內(nèi)孤立波速度反演模型時將水深、水體分層結(jié)構(gòu)、密度差和遙感圖像的亮暗間距、灰度差作為速度反演模型的自變量.下文分別利用最小二乘法和機(jī)器學(xué)習(xí)兩種技術(shù)建立內(nèi)孤立波速度反演模型.

3.1 最小二乘法的內(nèi)孤立波速度反演模型

最小二乘法是解決曲線擬合的常用方法,其原理是尋找最小化的誤差平方和來確定與數(shù)據(jù)對應(yīng)的最佳匹配函數(shù),進(jìn)而可以直觀表達(dá)因變量與自變量的函數(shù)關(guān)系,物理意義也就一目了然,所以最小二乘法在多個領(lǐng)域研究中被廣泛地應(yīng)用.本文利用所獲得的大量內(nèi)孤立波實(shí)驗(yàn)數(shù)據(jù),采用最小二乘法尋找最佳的回歸方程,建立基于最小二乘法的反演內(nèi)孤立波速度模型.將無量綱數(shù)據(jù)庫中的H,△ρ,D,△G作為模型自變量,將V作為因變量.考慮到4 個自變量是相互獨(dú)立的,所以在回歸方程中將其作為獨(dú)立項(xiàng)進(jìn)行擬合.由KdV 方程等理論可知速度與自變量具有非線性關(guān)系,回歸方程中應(yīng)包含有自變量的高次方項(xiàng).因此通過最小二乘法擬合的二次回歸方程和三次回歸方程如下:

用實(shí)驗(yàn)數(shù)據(jù)檢驗(yàn)回歸方程反演值的結(jié)果如圖3 所示,圖3(a)為2 次方的回歸方程的結(jié)果,方程反演值與實(shí)驗(yàn)值的相關(guān)系數(shù)為0.927,平均絕對誤差為1.64 cm/s,平均相對誤差為8.25%.圖3(b)為最高次為3 次方的回歸方程的結(jié)果,方程反演值與實(shí)驗(yàn)值的相關(guān)系數(shù)為0.928,平均絕對誤差為1.61 cm/s,平均相對誤差為8.13%.從圖中可知,2 個回歸方程都有效,它們的相關(guān)系數(shù)和誤差大小接近,差異不大.

圖3 最小二乘法回歸方程反演結(jié)果散點(diǎn)圖 (a)最高冪次為2 次方;(b)最高冪次為3 次方Fig.3.Scatter plot of the inversion results of the least squares regression equation: (a) With the highest power of 2;(b) with the highest power of 3.

3.2 支持向量回歸內(nèi)孤立波速度反演模型

支持向量機(jī)中的支持向量回歸(support vector regression,SVR)是一種應(yīng)用于函數(shù)回歸預(yù)測領(lǐng)域的機(jī)器學(xué)習(xí)方法,SVR 相較于其他機(jī)器學(xué)習(xí)方法而言更適合于小樣本的回歸預(yù)測,且能夠處理樣本數(shù)據(jù)線性不可分的情況,所以本文采用SVR 的方法對內(nèi)孤立波速度反演模型進(jìn)行訓(xùn)練.選擇局部性強(qiáng)且有將訓(xùn)練樣本的特征空間從低維映射到無窮維性能的高斯徑向基函數(shù)(radial basis function,RBF)作為模型核函數(shù).

與最小二乘法相同,將無量綱數(shù)據(jù)庫中的H,Δρ,D,ΔG作為模型輸入?yún)⒘?將V/C0作為輸出參量.首先將所有數(shù)據(jù)輸入SVR 程序中運(yùn)行,去除誤差較大的數(shù)據(jù)后,將剩余的1383 條樣本數(shù)據(jù)隨機(jī)排序,按照大致4∶1 的比例劃分為模型訓(xùn)練集和測試集,訓(xùn)練集為1100 條樣本數(shù)據(jù),測試集為283 條樣本數(shù)據(jù).之后再次將數(shù)據(jù)輸入SVR程序進(jìn)行模型訓(xùn)練,得到模型后,記錄模型的懲罰系數(shù)c和核函數(shù)參數(shù)g的值,c的值為1,g的值為16.圖4 為SVR 內(nèi)孤立波速度反演模型測試集實(shí)驗(yàn)數(shù)據(jù)與反演數(shù)據(jù)散點(diǎn)圖,反演值與實(shí)驗(yàn)值的相關(guān)系數(shù)為0.888,平均絕對誤差為2.10 cm/s,平均相對誤差為10.28%,表明模型是有效的.

圖4 SVR 內(nèi)孤立波速度反演模型測試集散點(diǎn)圖Fig.4.Scatter plot of test set of solitary wave speed inversion model in SVR.

4 結(jié)果與分析

4.1 光學(xué)遙感圖像下載與處理

選取2021 年中國南海海域的遙感圖像和數(shù)據(jù)來進(jìn)行模型精度驗(yàn)證,在NASA 官方遙感影像數(shù)據(jù)下載網(wǎng)站LADDS DAAC 下載所需的MODIS遙感圖像,在中國資源衛(wèi)星應(yīng)用中心下載高分4 號衛(wèi)星光學(xué)遙感圖像.對遙感圖像進(jìn)行篩選和處理,提取參數(shù),對模型精度進(jìn)行驗(yàn)證.

采用多時間圖像法對內(nèi)孤立波速度進(jìn)行計算,從獲得有內(nèi)孤立波的遙感圖像中選取同一天不同時刻的圖像,找到同一個內(nèi)孤立波,根據(jù)不同遙感圖像中的內(nèi)孤立波位置計算其傳播距離,兩景遙感圖像的間隔時間為內(nèi)孤立波的傳播時間,傳播距離與傳播時間的比值就是傳播速度.同時利用ENVI軟件提取內(nèi)孤立波的亮暗間距和灰度差等遙感圖像參數(shù).圖5 為多時間圖像法示意圖.一般在采用多時間圖像法計算內(nèi)孤立波傳播速度時兩景圖像的時間間隔約在1—3 h,時間間隔較短,所以將計算的平均速度近似作為兩景遙感圖像中內(nèi)孤立波的傳播速度.

圖5 多時間圖像法示意圖 (a) 2021年5月25日10:43南海海域GF4光學(xué)遙感圖像;(b) 2021年5月25日13:40南海海域MODIS 光學(xué)遙感圖像Fig.5.Schematic diagram of the multi-time image method: (a) An optical remote sensing image of GF-4 in the South China Sea at 10:43 on May 25,2021;(b) the MODIS optical remote sensing image of the South China Sea area is at 13:40 on May 25,2021.

從光學(xué)遙感圖像中提取內(nèi)孤立波條紋的亮暗間距和灰度差,對數(shù)據(jù)進(jìn)行無量綱化處理,將輸入?yún)⒘看敕囱菽Ｐ?利用上述模型計算內(nèi)孤立波速度.與多時間圖像法計算出的內(nèi)孤立波傳播速度比較,驗(yàn)證模型精度.

4.2 模型精度驗(yàn)證

在利用最小二乘法內(nèi)孤立波速度反演模型對遙感數(shù)據(jù)進(jìn)行反演時,兩種模型的反演值均小于實(shí)測值,存在系統(tǒng)偏差,如圖6 中黑色散點(diǎn)所示.實(shí)驗(yàn)仿真與真實(shí)海洋的差異可能是導(dǎo)致驗(yàn)證時產(chǎn)生系統(tǒng)偏差的原因.為了消除系統(tǒng)偏差,對模型進(jìn)行校正,使得系統(tǒng)偏差和其他誤差最小,校正后的模型如下:

圖6 最小二乘法內(nèi)孤立波速度反演模型精度散點(diǎn)圖 (a)校正前后最高冪次為2 次方的最小二乘法回歸方程的反演結(jié)果散點(diǎn)圖,黑色散點(diǎn)為校正前模型結(jié)果,紅色散點(diǎn)為校正后模型結(jié)果;(b)校正前后最高冪次為3 次方的最小二乘法回歸方程的反演結(jié)果散點(diǎn)圖,黑色散點(diǎn)為校正前模型結(jié)果,紅色散點(diǎn)為校正后模型結(jié)果.Fig.6.Scatter plot of the solitary wave speed inversion model accuracy within least squares.(a) The scatter plot of inversion results of least squares regression equation with the highest power of 2 before and after correction.The black scatters are the model results before correction,and the red scatters are the model results after correction.(b) The scatter plot of the inversion results of the least squares regression equation with the highest power of 3 before and after correction.The black scatters are the model results before correction,and the red scatters are the model results after correction.

(8)式為校正后的2 次方最小二乘法反演模型,常數(shù)項(xiàng)由1.25 校正為1.84,校正前后的實(shí)測值與反演值分別如圖6(a)中黑色散點(diǎn)和紅色散點(diǎn)所示.(9)式為校正后的3 次方最小二乘法反演模型,常數(shù)項(xiàng)由0.93 校正為1.71,校正前后的實(shí)測值與反演值分別如圖6(b)黑色散點(diǎn)和紅色散點(diǎn)所示.由圖6 可知,校正后的模型消除了系統(tǒng)偏差,反演效果良好,所以將校正后的模型作為最終模型.

依據(jù)4.1 節(jié)的方法,共獲得南海光學(xué)遙感觀測的44 個內(nèi)孤立波速度數(shù)據(jù),將所得的遙感觀測數(shù)據(jù)按水深劃分為400—1200 m,300—399 m,200—299 m 和83—199 m 四個范圍,用3 種模型對光學(xué)遙感數(shù)據(jù)進(jìn)行反演.反演結(jié)果如表2 和圖7所示.

表2 內(nèi)孤立波速度反演模型精度驗(yàn)證表Table 2.Precision verification table of internal solitary wave velocity inversion model.

圖7 內(nèi)孤立波速度反演模型精度驗(yàn)證圖 (a)—(d) 3 種反演模型對400—1200 m、300—399 m、200—299 m、83—199 m 水深范圍的內(nèi)孤立波的速度反演結(jié)果散點(diǎn)圖;(e) 3 種反演模型對83—1200 m 水深范圍的內(nèi)孤立波的速度反演結(jié)果散點(diǎn)圖;(a)—(e) 黑色散點(diǎn)為最高冪次為2 次方的最小二乘法反演模型結(jié)果,紅色散點(diǎn)為最高冪次為3 次方的最小二乘法反演模型結(jié)果,藍(lán)色散點(diǎn)為SVR 反演模型結(jié)果;(f) 三種反演模型對400—1200 m、300—399 m、200—299 m、83—199 m 水深范圍的內(nèi)孤立波的速度反演結(jié)果的平均絕對誤差柱狀圖Fig.7.Precision validation diagram of the internal solitary wave velocity inversion model: (a)–(d) Scatter plots of inversion results of three inversion models for internal solitary waves speed in water depths of 400–1200 m,300–399 m,200–299 m and 83–199 m,respectively;(e) scatter plot of inversion results of three inversion models for internal solitary waves speed in water depths of 83–1200 m;(a)–(e) the black scattered points are the results of least squares inversion model with the highest power of 2,the red scattered points are the results of least squares inversion model with the highest power of 3,and the blue scattered points are the results of SVR inversion model;(f) the average absolute error histogram of the inversion results of the three inversion models for internal solitary waves speed in the water depths of 400–1200 m,300–399 m,200–299 m and 83–199 m.

由圖7(a),(c),(d),(f)可知,SVR 內(nèi)孤立波速度反演模型在水深400—1200 m 和83—299 m 的范圍精度要高于最小二乘法的兩種模型.由圖7(b),(f)可知在水深范圍300—399 m 時,最小二乘法內(nèi)孤立波速度反演模型精度較高,其中3 次方的反演模型誤差最小,精度最高.從圖7(f)還可看到,在各個水深二次方和三次方的最小二乘法內(nèi)孤立波速度反演模型精度相差較小,所以在滿足精度要求時內(nèi)孤立波速度反演方程就可以采用二次方程描述.

5 總結(jié)

本文在實(shí)驗(yàn)室中搭建內(nèi)孤立波光學(xué)遙感仿真平臺,進(jìn)行了系列仿真實(shí)驗(yàn).提取并處理實(shí)驗(yàn)數(shù)據(jù),建立樣本庫.利用最小二乘法和機(jī)器學(xué)習(xí)SVR 兩種技術(shù)建立了內(nèi)孤立波速度反演模型.共建立了三種反演模型,分別是最小二乘法3 次方內(nèi)孤立波速度反演模型、最小二乘法2 次方內(nèi)孤立波速度反演模型和SVR 內(nèi)孤立波速度反演模型.

收集南海內(nèi)孤立波光學(xué)遙感圖像,提取遙感參數(shù)和水文參數(shù),利用多時間圖像法計算內(nèi)孤立波速度,作為內(nèi)孤立波現(xiàn)場實(shí)測速度,對反演模型精度進(jìn)行驗(yàn)證.最小二乘法建立的兩個反演模型對實(shí)測數(shù)據(jù)反演時存在系統(tǒng)誤差,所以利用實(shí)測數(shù)據(jù)對反演模型進(jìn)行了校正.之后將三種模型按不同水深范圍分別進(jìn)行精度驗(yàn)證,結(jié)果為: SVR 內(nèi)孤立波速度反演模型在水深400—1200 m 和水深83—299 m的范圍精度要高于最小二乘法的兩種反演模型,表明了SVR 非線性擬合能力更強(qiáng).在水深范圍300—399 m 時,最小二乘法內(nèi)孤立波速度反演模型精度較高,其中3 次方的反演模型誤差最小,精度最高.三種模型的反演精度和誤差基本都在合理范圍,在進(jìn)行真實(shí)海洋條件實(shí)際應(yīng)用時,可以將兩種方法的模型結(jié)合應(yīng)用,當(dāng)SVR 內(nèi)孤立波速度反演模型更適用于水深范圍為400—1200 m 或83—299 m,水深范圍在300—399 m 時,優(yōu)選3 次方最小二乘法內(nèi)孤立波速度反演模型進(jìn)行速度反演.兩種方法結(jié)合的方式,可以從不同角度理解內(nèi)孤立波的速度.

通過南海海域遙感數(shù)據(jù)的驗(yàn)證可知,最小二乘法內(nèi)孤立波速度反演模型和SVR 速度反演模型均具有較高精度,最小二乘法內(nèi)孤立波速度反演模型給出回歸方程進(jìn)行反演,物理意義更為直觀;SVR 內(nèi)孤立波速度反演模型則在全水深都具有較高精度.兩種模型都可應(yīng)用于南海的內(nèi)孤立波速度反演.