笪良龍，李韋華，韓 梅，程艷合

(海軍潛艇學院，山東 青島 266071)

隨機動態(tài)聲速重構(gòu)方法研究

笪良龍，李韋華，韓 梅，程艷合

(海軍潛艇學院，山東 青島 266071)

海洋環(huán)境對聲傳播具有重要的影響，海洋溫度、鹽度、壓力決定聲速的大小，研究隨機動態(tài)聲速剖面的最優(yōu)重構(gòu)可以為聲場模型的計算提供良好的基礎(chǔ)，為艦艇行動提供最優(yōu)的輔助決策。然而聲速剖面的重構(gòu)卻是一個復雜的過程，涉及到大量的時空分布不規(guī)則的數(shù)據(jù)的計算和處理，提取出合理的時空間系數(shù)能夠提高聲速重構(gòu)的精度。實驗證明，利用合理的時空分辨率，從海洋環(huán)境數(shù)據(jù)庫中選取數(shù)據(jù)可以提高海洋聲速剖面重構(gòu)的精度。

聲速；重構(gòu)；海洋環(huán)境；精度

我國周邊海域受到風浪、季節(jié)變化等的影響，海水的溫度變化較快，而溫度是影響海水聲速的重要因素，海水聲速的細微變化對聲納作用距離、匯聚區(qū)的大小影響重大，因此，對海水聲速的快速重構(gòu)直接影響著海軍作戰(zhàn)的決策效率。事實上，人們也在對海水聲速(或水溫、鹽度、密度等)進行大量的測試和分析，并試圖尋求對海水聲速的變化規(guī)律進行數(shù)學描述，建立足夠完備的海洋水文數(shù)據(jù)庫系統(tǒng)?，F(xiàn)實的海洋聲速場是時間和空間的函數(shù)，即使是間距很近的兩個海區(qū)或同一海區(qū)的不同時間內(nèi)，其聲速場(尤其是海面下的溫躍層)也可能存在較大的差異。因此，要建立這樣的數(shù)據(jù)庫，其容量會是十分可觀的，也是很難做到的，而有效地對聲速場變化特征進行數(shù)學描述，并減少數(shù)據(jù)庫的數(shù)據(jù)量，是十分必要的。根據(jù)我國現(xiàn)有的水聲數(shù)據(jù)庫，需要建立一個適合于我國周邊海域使用的聲速剖面重構(gòu)系統(tǒng)，通過對這些海域的水聲數(shù)據(jù)進行試驗，找出重構(gòu)不同位置點的最優(yōu)的時空分布是十分重要的。但以往海軍作戰(zhàn)使用的聲速數(shù)據(jù)都是基于月平均或季節(jié)平均的資料，掩蓋了某些實時的海洋環(huán)境信息，對指揮官的決策起到了誤導作用，隨著海軍信息化作戰(zhàn)的精確打擊的要求，建立以實測資料為基礎(chǔ)的海洋聲速重構(gòu)系統(tǒng)迫在眉睫。D.N.Fox[1]等將改進后的經(jīng)驗正交函數(shù)的方法融入到美國海軍MODAS數(shù)據(jù)同化系統(tǒng)中，實現(xiàn)剖面延拓、預報，具有良好的效果。彭臨慧[2]等也將該方法應用于深海，進行深海聲速剖面的預報。魏亮[3]等則采用該方法成功地對中國近海某點實現(xiàn)了預報，預報結(jié)果非常理想。在利用改進型經(jīng)驗正交函數(shù)算法的基礎(chǔ)上，結(jié)合1911—2008年存有的數(shù)據(jù)庫，實現(xiàn)了隨機動態(tài)的聲速剖面的預報，并對中國近海的兩個ARGO數(shù)據(jù)點進行了預報，與ARGO實測數(shù)據(jù)和月平均資料進行了對比，有效地提高了聲速剖面的預報精度。

1 隨機動態(tài)聲速重構(gòu)的基本原理

1.1 改進型經(jīng)驗正交函數(shù)算法

若研究的海洋環(huán)境因子為聲速，假設(shè)在某個海域有N個已知的聲速剖面樣本，每個剖面有k+1層，則這些對應的聲速剖面可表示成矩陣形式為：

改進型經(jīng)驗正交函數(shù)則是指原始聲速矩陣C在深度m與深度n上的協(xié)方差矩陣，其組成元素ri,m,n：

式中：σTi,mσTi,n由（3）式確定：

式中：bi,j（m,n）是位置點 i和 j之間的相關(guān)性系數(shù)。依據(jù)特征值絕對值大小由大到小選出較大絕對值所對應的特征向量，用絕對值最大的三個特征值對應特征向量ek作為EOF分解的空間函數(shù)，得到合成剖面如式（5）：

式中：αj為特征向量ek的權(quán)值，即EOF分解的時間函數(shù)，采用最簡單的使用與預報點相關(guān)系數(shù)最大的完整聲速剖面擬合出時間函數(shù)作為αj的值。

1.2 數(shù)據(jù)選取的原則

海洋環(huán)境數(shù)據(jù)庫中的溫度、鹽度、密度、深度、聲速這些量的分布都是由各個剖面的測點編號唯一決定的，每個剖面的測點編號中包含了所給剖面的經(jīng)緯度、時間的信息，對數(shù)據(jù)進行選取時可以按照測點編號進行選取，對經(jīng)緯度和時間的范圍進行規(guī)定，采用該區(qū)域的時空分辨率對數(shù)據(jù)進行篩選，有效地提高聲速剖面重構(gòu)的精度。對于數(shù)據(jù)時空的自相關(guān)函數(shù)可以用式（6）表示：

式中：Δxij，Δyij，Δtij分別表示不同位置點間的東西方向、南北方向和時間的分離值。而Ak，Bk，Ck表示在不同位置層的東西方向、南北方向和時間的去相關(guān)尺度，這三個參數(shù)因海域不同而不同，它反映出海洋現(xiàn)象的變化范圍，相關(guān)尺度計算的越好，數(shù)據(jù)選取的范圍就會近似，計算結(jié)果就越精確[4]。根據(jù)我國近海各個海底地形地貌選取合適的時空間尺度，有效提高聲速精度。

2 重構(gòu)使用的數(shù)據(jù)資料

聲速重構(gòu)系統(tǒng)使用的資料來源于1911—2008年的包括CTD資料、南森數(shù)據(jù)、ARGO數(shù)據(jù)構(gòu)成的實測海洋環(huán)境數(shù)據(jù)庫，用于對比的數(shù)據(jù)來自于2009年的ARGO數(shù)據(jù)，以及30'×30'的月平均資料。

3 重構(gòu)系統(tǒng)模塊設(shè)計與實現(xiàn)

3.1 重構(gòu)系統(tǒng)的流程與結(jié)構(gòu)

系統(tǒng)的控制模塊主要分為參數(shù)的設(shè)定、數(shù)據(jù)庫的檢索、聲速剖面的篩選、算法的運行、結(jié)果可視化實現(xiàn)。

整個系統(tǒng)軟件的實現(xiàn)的流程圖可如圖所示，系統(tǒng)軟件是采用C++Builder設(shè)計的，并連接到存有海洋環(huán)境信息的數(shù)據(jù)庫中，主要包括溫度、鹽度、聲速、深度。在編程實現(xiàn)從數(shù)據(jù)庫中查詢的合適的數(shù)據(jù)時，對以下幾個條件進行了設(shè)定：

⑴時間、空間的范圍和向外搜索的時空分辨率由各個海區(qū)決定。

⑵對各個剖面的深度進行篩選，由于中國近海是大陸坡地形，水深的變化明顯，根據(jù)要求選取合適的深度的數(shù)據(jù)進行重構(gòu)非常必要。

⑶已選取各個剖面的最大測量深度不一致，必須對各個剖面的最小深度進行計算，便于數(shù)據(jù)的網(wǎng)格化。

軟件設(shè)計的具體流程和基本結(jié)構(gòu)可見圖1、圖2。

圖1 重構(gòu)系統(tǒng)的流程圖

圖2 重構(gòu)系統(tǒng)操作過程圖

3.2 系統(tǒng)仿真的結(jié)果

為檢驗系統(tǒng)計算結(jié)果的效能，選取中國臺灣以東海域具有ARGO實測數(shù)據(jù)點的位置進行了聲速剖面的重構(gòu)。選取的位置點共2個，第一個剖面的位置點為 φ121.626°，λ21.863°，第二個位置點為 φ125.65°，λ24.578°，該位置點的ARGO數(shù)據(jù)的時間分別是2009年10月16日和2009年10月28日，具體位置見圖3。

圖3 選取位置點的地圖

采用聲速重構(gòu)方法構(gòu)造的兩個位置點的聲速剖面及與AGRO實測數(shù)據(jù)的比較如圖4～圖7所示。從圖4、圖5中可以看到，使用該重構(gòu)系統(tǒng)構(gòu)造的聲速剖面對ARGO實測數(shù)據(jù)的擬合比月平均資料要好，但在躍層位置處對聲速的重構(gòu)不夠好，這可能主要是由于黑潮的影響，使得重構(gòu)出的剖面最大誤差約為6 m/s。

圖4 第一個位置點聲速剖面圖

圖5 第一個位置點兩種聲速與ARGO的差值

從圖6、圖7中可以看到，第二個位置點用重構(gòu)系統(tǒng)對聲速剖面的重構(gòu)比月平均資料的擬合好很多，在此位置處各種海洋現(xiàn)象影響較小，對剖面的重構(gòu)起到了較好的作用，最大誤差不超過4 m/s。

圖6 第二個位置點聲速剖面圖

圖7 第二個位置點兩種聲速與ARGO的差值

4 結(jié)論

通過對隨機動態(tài)聲速重構(gòu)系統(tǒng)的研究，可以得到以下兩點結(jié)論：

⑴使用該系統(tǒng)可以隨機地對中國近海任意一點的聲速進行重構(gòu)，得到的結(jié)果比較好；

⑵通過對時空范圍最優(yōu)的選取，可以較好地選取數(shù)據(jù)，提高聲速重構(gòu)的精度。

[1]Fox D N,Teague W J,Barron C N,et al.The Modular Ocean Data Assimilation System(MODAS)[J].Journal of Atmospheric and Oceanic Technology,2002,19（2）:240-252.

[2]Peng L,Wang L,Qiu X F.Modal wave number tomography for South China Sea front[J].China Ocean Engineering，2003,17(2):289-294.

[3]韓梅，魏亮，周艷霞.改進型經(jīng)驗正交函數(shù)海洋聲速剖面預報方法[J].海洋科學，2009（1）:235-239.

[4]Chu P C,Gui W,Chen Y C.Japan Sea thermohaline structure and circulation[J].American Meteorological Society,2002,32:3596-3615.

Research on Stochastic and Dynamic Sound Speed Reconstruction System

DA Liang-long,LI Wei-hua,HAN Mei,CHENG Yan-He
(Navy Submarine Academy,Qingdao Shandong 266071,China)

Ocean environment has an important effect on the sound propagation.Temperature,salinity,pressure determine the sound profile.The research of stochastic and dynamic sound speed reconstruction can offer foundation for the computation of sound field and an optimize decision for the action of warship.But the reconstruction of sound profile is a complex process,including substantial data computing and processing that distributes in time and space anomalistically.Getting logical coefficient of time and space,the precision of sound profile will increase.Experiment proves that logical resolving power of time and space will increase the precision of sound profile.

sound speed;reconstruction;ocean environment;precision

P733

A

1003-2029（2011）03-0056-04

2011-04-11

總裝預研基金資助項目（9140A22050810JB1502）

笪良龍(1967-)，男，安徽桐城人，教授，主要研究領(lǐng)域為水聲環(huán)境效應技術(shù)。