臧 濤 海 偉

（1.海軍駐武漢701所軍代室 武漢 430061）（2.海軍裝備研究院 北京 100161）

淺海海底聲學(xué)參數(shù)反演方法研究＊

臧 濤1海 偉2

（1.海軍駐武漢701所軍代室 武漢 430061）（2.海軍裝備研究院 北京 100161）

海底聲學(xué)參數(shù)反演問題近年來已經(jīng)成為水聲界的研究熱點(diǎn)之一。論文針對不同的海底模型和海底參數(shù)提出了兩種反演方法。淺海多途效應(yīng)導(dǎo)致聲場在距離－頻率平面形成干涉譜，利用距離域干涉譜提取出干涉準(zhǔn)周期作為反演匹配物理量，提出一種基于無限均勻沉積層海底模型的海底反演方法。針對傾斜海底，提出一種基于特征聲線傳播時間的海底反演方法。

淺海；海底參數(shù)反演；干涉譜；聲線傳播時間

Class NumberTN929

1 引言

淺海環(huán)境下，海底的聲學(xué)特性對聲傳播起著決定性的作用，海底的聲速、密度、衰減系數(shù)和分層特征等地聲屬性對聲場預(yù)報、匹配場定位［1］等問題都具有重要意義。但是海底的聲學(xué)參數(shù)很難通過直接測量而獲取，而且直接測量獲取海底參數(shù)的方法比較耗時費(fèi)力。利用反演技術(shù)，所期望的海底聲學(xué)參數(shù)就有可能快速、低成本地從海洋波導(dǎo)聲傳播的直接測量中提取出來。因此，海底聲學(xué)參數(shù)反演方法的探索和研究近年來已成為國內(nèi)外水聲學(xué)家研究的熱點(diǎn)問題，并取得了大量成果［2～7］。

2 海底模型

根據(jù)反演的目的，海底參數(shù)反演可大致分為兩類［8］：一是用盡可能簡單的海底等效模型，反演出一組等效的海底參數(shù)；二是盡可能準(zhǔn)確地反演出海底的地質(zhì)特性。兩類反演分別針對聲納應(yīng)用和水聲物理實(shí)驗(yàn)現(xiàn)象的解釋有不同的側(cè)重點(diǎn)。具體應(yīng)用中加以區(qū)分。

圖1給出幾種典型水平不變的海底模型。

從圖1可以看出，模型越復(fù)雜，待反演的參數(shù)越多，對反演的壓力也越大，甚至在參數(shù)過多的情況下，會導(dǎo)致計算失敗。

圖1 幾種常見的海底分層模型

3 利用距離域干涉譜的海底參數(shù)反演

3.1 距離域干涉譜

當(dāng)聲源和接收器之間的距離r滿足遠(yuǎn)場條件時，無指向性聲源的接收壓力場接收壓力場可以表示為有限階簡正波的累加。若聲源頻率為w，聲源和接收器的深度分別為zs和zr，那么，接收到的聲強(qiáng)場為［9］

其中

式（1）中ψn（z）為n號簡正波的本征函數(shù)；kn和βn分別為n號簡正波的水平波數(shù)和衰減系數(shù)。第一項(xiàng)為非相干項(xiàng)，隨距離和頻率緩慢變化；第二項(xiàng)為相干項(xiàng)，由于各號簡正波之間相互干涉，相干項(xiàng)會隨距離出現(xiàn)振蕩。當(dāng)聲源頻率較低時，激發(fā)的簡正波的號數(shù)較少，相干項(xiàng)會在距離－頻率平面上表現(xiàn)出明顯的干涉條紋。從式（1）進(jìn)行分析，對于某一固定頻率，聲場隨著距離的分布，有一定的周期性，即

式中，δr表示目標(biāo)聲源產(chǎn)生的聲場聲強(qiáng)值在距離坐標(biāo)上的周期值，當(dāng)然這一周期對應(yīng)第m號和第n號簡正波的干涉效果，不同的簡正波模態(tài)對應(yīng)不同的周期分布，把這種干涉周期分布稱為距離域干涉譜。對于某一頻率點(diǎn)對應(yīng)的距離域干涉譜，相鄰條紋距離差的平均值定義為距離域干涉準(zhǔn)周期。利用距離域干涉譜提取干涉準(zhǔn)周期作為反演匹配物理量。

3.2 敏感性分析

只有當(dāng)反演參數(shù)與采用的匹配物理量之間具有較好的敏感度才能進(jìn)行反演。在淺海等聲速梯度下，利用Kraken簡正波模型［10］仿真得到干涉準(zhǔn)周期隨沉積層密度、壓縮波聲速和衰減系數(shù)變化的曲線，如圖2所示。

圖2 干涉準(zhǔn)周期與密度、聲速以及衰減系數(shù)的關(guān)系圖

從圖中可以看出，干涉準(zhǔn)周期整體隨著密度、壓縮波聲速的增大而減小，隨著衰減系數(shù)的增大而增大。干涉準(zhǔn)周期與各參數(shù)敏感性較好。

3.3 代價函數(shù)的建立及最優(yōu)解搜索方法

文獻(xiàn)［11］中提出利用最小均方差準(zhǔn)則建立代價函數(shù)。

式中，N表示頻率點(diǎn)數(shù)；δidata表示某一頻率對應(yīng)的實(shí)驗(yàn)數(shù)據(jù)值；δipred表示某一頻率正向聲場計算值。最優(yōu)組合即為A取得最小值時對應(yīng)的沉積層密度、壓縮波聲速和衰減系數(shù)。

目前搜索方法有很多，傳統(tǒng)的優(yōu)化算法在實(shí)際應(yīng)用中有諸多的限制，如問題的復(fù)雜性、算法的并行性、計算速度等方面的致命缺陷，使得在一些大型的、高強(qiáng)度的優(yōu)化計算中逐漸被一些發(fā)展相對較晚的啟發(fā)式算法所取代。這類算法有：遺傳算法、模擬退火算法、粒子群算法、蟻群算法等。尤其是模擬退火算法、遺傳算法以及二者相結(jié)合的SAGA算法，在實(shí)際應(yīng)用中表現(xiàn)出各種優(yōu)勢，使得它們在目前的優(yōu)化計算領(lǐng)域成為主流。其中，模擬退火算法（SA）是通過隨機(jī)地擾動模型參數(shù)，并接受使目標(biāo)函數(shù)變小的模型參數(shù)和有概率地接受跳出局部小的模型參數(shù)，最終能夠較好地收斂到搜索空間里的最小值。

本文在單無限均勻海底模型的基礎(chǔ)上，以最小均方差原則建立代價函數(shù)，采用模擬退火法作為代價函數(shù)最優(yōu)解的搜索策略。

4 基于特征聲線傳播時間的海底參數(shù)反演

圖3 傾斜海底參數(shù)反演射線模型

如果海底沉積層下邊界傾斜，假定聲源坐標(biāo)為T（0，0，H），水聽器坐標(biāo)為R（r，θ，Z），沉積層厚度為h，沉積層聲速為cb，沉積層下界面傾斜角度為α，聲線在沉積層上界面的入射點(diǎn)坐標(biāo)為E（ρ1，β1，0），出射點(diǎn)坐標(biāo)為F（rcosθ＋ρ2cosβ2，rsinθ＋ρ2sinβ2，0），在沉積層下界面的反射點(diǎn)坐標(biāo)為M′（k1ρ1cosβ1，k1ρ1sinβ1，－（k1ρ1sinβ1tanα＋h））。根據(jù)入射聲線與反射聲線在水平坐標(biāo)系中的關(guān)系可得：

根據(jù)Snell定律，聲線在水層和沉積層邊界上滿足下列關(guān)系：

又因?yàn)槁暰€入射角余弦和反射角余弦相等，同時有入射聲線、反射聲線和法線共面，即→a·（→b×→c）＝0，所以：

聲線傳播時間等于聲波傳播路程除以聲速，這里假設(shè)海水中聲速為c0，則

上述特征聲線所滿足的七個方程中有ρ1、f0、k1、k2、β1、β2、cb、h、α、θ共10個未知參數(shù)。其中cb、h、α、θ為沉積層參數(shù)，不隨水聽器位置而變化；為聲線參數(shù)，將隨著水聽器位置的變化而變化。對于該沉積層模型，一根特征聲線滿足上述七個方程，為了求出沉積層參數(shù)及特征聲線參數(shù)，此時至少需要根據(jù)四根不同的特征聲線所滿足的28個方程進(jìn)行求解。方程組中將含有24個未知的聲線參數(shù)及四個沉積層參數(shù)。采用優(yōu)化算法求解此方程組便得到了沉積層的聲速、厚度、傾斜角度以及傾斜方向，同時也求解得到了滿足條件的特征聲線。當(dāng)從數(shù)據(jù)中獲取兩根特征聲線的幅度后，根據(jù)信號的幅度比及特征聲線軌跡長度差計算得到沉積層的衰減系數(shù)。因此對于傾斜沉積層海底模型，海底參數(shù)的反演至少需要獲取四個位置的聲場信息。

采用牛頓迭代法求解上述非線性方程組，此時對于初值的選擇可根據(jù)不同的ε≈｜sinθs（kw－kb）｜／；λ／L設(shè)置不同的初值區(qū)間，將這些區(qū)間中不同的初值組合代入牛頓迭代法中進(jìn)行計算，根據(jù)方差最小準(zhǔn)則來判斷結(jié)果是否滿足要求，如果滿足要求就使程序停止運(yùn)行，將結(jié)果顯示出來。如果沒有滿足要求且迭代次數(shù)達(dá)到上限，則重新代入初值計算。

理論上，水平距離的變化越小，到達(dá)不同水聽器的特征聲線傳播時間則越接近，傳播時間計算誤差對反演精度的影響越大，從而導(dǎo)致反演結(jié)果誤差變大。對于接收深度的取值變化也是同樣道理。當(dāng)沉積層厚度增大時，特征聲線傳播時間受沉積層厚度的影響也越大，這樣就越有利于反演沉積層的厚度和聲速。應(yīng)盡可能多地使接收到的信號包含更多沉積層的信息。由于本模型中目標(biāo)函數(shù)對角度變量不敏感，這使得距離存在較小誤差的情況下，角度參數(shù)的反演結(jié)果就會存在較大的誤差。因此，反演結(jié)果對于特征聲線傳播時間的精度要求比較高。

5 結(jié)語

本文首先針對無限均勻沉積層海底模型，利用淺海距離域干涉譜相關(guān)理論，提出通過干涉準(zhǔn)周期這一物理量進(jìn)行匹配來反演地聲參數(shù)。這種方法對于反演沉積層密度以及壓縮波聲速效果較好，但對于海底衰減系數(shù)的反演效果一般。該方法在聲場預(yù)報方面有一定的使用價值。

而對于傾斜海底，本文提出一種基于特征聲線傳播時間的海底參數(shù)反演方法。理論上，應(yīng)用射線聲學(xué)模型可以快速獲取本地海底沉積層的聲速、衰減系數(shù)、厚度以及下界面的傾斜角度等參數(shù)信息。

［1］楊坤德，馬遠(yuǎn)良.基于扇區(qū)特征向量約束的穩(wěn)健自適應(yīng)匹配場處理器［J］.聲學(xué)學(xué)報，2006，31（5）：399－409.

［2］Gerstoft，P.Inversion of seismoacoustic data using genetic algorithms and a posteriori probability distribution［J］.J.Acoust.Soc.Amer.，1994（95）：770－782.

［3］Siderius，Nielsen，Sellschopp，Snellen，and Simons.Experimental study of geoacoustic inversion uncertainty due to ocean sound－speed fluctuations［J］.J.Acount.Soc.Amer.，2001（110）：769－781.

［4］劉建軍，李風(fēng)華，郭良浩.淺?；祉懙拇怪毕嚓P(guān)和海底反射損失與散射強(qiáng)度的反演［J］.聲學(xué)學(xué)報，2004，29（1）：49－56.

［5］Li Zhenglin，Zhang Renhe.Geoacoustic inversion by MFP combined with vertical reflection coefficients and vertical correlation［J］.IEEE J.Oceanic Engeneering，2004，29（4）：973－979.

［6］李整林，鄢錦，李風(fēng)華，等.由簡正波群延時及幅度反演海底參數(shù)［J］.聲學(xué)學(xué)報，2002，27（6）：487－491.

［7］楊坤德，馬遠(yuǎn)良.利用海底反射信號進(jìn)行地聲參數(shù)反演的方法［J］.物理學(xué)報，2009，58（3）：1798－1804.

［8］王占軍.垂直陣的海底參數(shù)反演研究［D］.哈爾濱：哈爾濱工程大學(xué)，2008.

［9］Jensen F B，Kuperman W A，Porter M B，et al.Computatinal ocean acoustics［M］.New York：AIP Press，1994：271－341.

［10］張歆，張小薊，李斌.基于Kraken簡正波模型的淺海聲場分析［J］.西北工業(yè)大學(xué)學(xué)報，2000，18（3）：405－408.

［11］K.D.Heaney.Rapid geoacoustic characterization using a surface ship of opportunity［J］.IEEE.J.Oceanic Eng，2004，29：88－99.

Inversion Method of Sea－bottom Acoustics Parameters in Shallowwater

ZANG Tao1HAI Wei2
（1.Navy Representative Office in 701th Research Institute，Wuhan 430061）（2.Navy Armament Academy，Beijing 100161）

Underwater acoustic association pays many attentions to the inversion of sea－bottom acoustics parameters in recent years.This paper puts forward two kinds of inversion methods about different sea－bottom model and sea－bottom parameters.Multipath propagation in a shallow－water environment leads to coherent structure in the

field，both spatially and temporally.In this paper，using range domain interference spectrum extraction cycle as a matching physical quantities，a seafloor parameters inversion method is put forward.With inclined seafloor，another seafloor parameters inversion method based on ray travel time is also put forward.

shallow water，inversion of sea－bottom parameters，interference spectrum，ray travel time

TN929DOI：10.3969／j.issn.1672－9730.2015.11.041

2015年5月10日，

2015年6月27日

臧濤，男，碩士研究生，工程師，研究方向：艦船總體質(zhì)量監(jiān)督。海偉，男，博士研究生，工程師，研究方向：艦船總體論證。