高　飛，潘長(zhǎng)明，李勝全，孫　磊，陳符森，2

（1.海軍海洋測(cè)繪研究所，天津，300061；2.解放軍理工大學(xué) 氣象海洋學(xué)院，江蘇 南京，211101）

基于HMG法的被動(dòng)聲吶淺海探測(cè)距離預(yù)報(bào)

高飛1，潘長(zhǎng)明1，李勝全1，孫磊1，陳符森1，2

（1.海軍海洋測(cè)繪研究所，天津，300061；2.解放軍理工大學(xué) 氣象海洋學(xué)院，江蘇 南京，211101）

基于HMG法并將實(shí)測(cè)數(shù)據(jù)與模型仿真相結(jié)合，以評(píng)估海洋“水聲環(huán)境效應(yīng)”對(duì)聲吶探測(cè)距離的影響。文章首先基于聲吶方程闡述 HMG法預(yù)報(bào)原理，給出其聲吶參數(shù)組合表達(dá)式；然后結(jié)合東海淺水區(qū)域水聲調(diào)查構(gòu)建海洋環(huán)境背景場(chǎng)，分析其水體及邊界聲學(xué)參數(shù)分布；最后，利用實(shí)測(cè)數(shù)據(jù)驗(yàn)證Kraken水聲傳播數(shù)值模型及風(fēng)生噪聲經(jīng)驗(yàn)?zāi)Ｐ?，并根?jù)HMG法模擬預(yù)報(bào)海面風(fēng)速、聲速剖面對(duì)被動(dòng)聲吶探測(cè)距離的影響。研究結(jié)果表明: HMG法可較好模擬“水聲環(huán)境效應(yīng)”對(duì)聲吶探測(cè)距離的影響；傳播損失和環(huán)境噪聲隨海面風(fēng)速增加而增大，從而降低了聲吶探測(cè)距離；聲速剖面通過影響聲能邊界損失作用于探測(cè)距離，垂向均勻的聲速剖面環(huán)境中探測(cè)距離最大。文中研究可為被動(dòng)聲吶水下探測(cè)及水下航行器的隱蔽航行研究提供參考。

水下航行器；HMG法；被動(dòng)聲吶探測(cè)；水聲環(huán)境效應(yīng)；水聲調(diào)查；聲速剖面

0　引言

聲波在水下環(huán)境中所表現(xiàn)出的普遍、異常變化特征和規(guī)律及其影響即水聲環(huán)境效應(yīng)［1］?！八暛h(huán)境效應(yīng)”對(duì)聲信道聲能傳輸特征、聲吶系統(tǒng)水下探測(cè)影響顯著，且已被大量研究所證實(shí)［2-7］。HMG（hydrology meteorology geology）法［8］用于系統(tǒng)評(píng)估“水聲環(huán)境效應(yīng)”對(duì)聲吶工作效能的影響，適用于模擬預(yù)報(bào)時(shí)空多變的海洋環(huán)境中被動(dòng)聲吶探測(cè)距離。

海洋環(huán)境影響著水下聲傳播特征規(guī)律進(jìn)而作用于聲吶探測(cè)效能。分析常見的海洋中尺度現(xiàn)象（中尺度渦［9］、內(nèi)波［10］、溫躍層［11-12］等）及海洋環(huán)境不確定性［13-14］多是從傳播損失角度研究其對(duì)聲吶探測(cè)效能的影響；對(duì)海面風(fēng)、浪［15-16］的研究則是分析海洋環(huán)境噪聲探索背景干擾級(jí)對(duì)聲吶檢測(cè)水下目標(biāo)信號(hào)的影響；對(duì)海底底質(zhì)分類［17］、地形［18］、粗超度［19］的研究則是從傳播損失、散射的角度評(píng)估其對(duì)聲吶工作效能的影響。

有關(guān)探測(cè)海洋環(huán)境對(duì)聲吶工作效能影響的研究也大量展開，Alec［20］以水聲信號(hào)處理為切入點(diǎn)，仿真分析海洋環(huán)境參數(shù)擾動(dòng)對(duì)聲吶水聽器能量接收機(jī)波動(dòng)峰值所造成的變化；Paul［21］、Wei［22］、Li［23］等通過構(gòu)建不確定海洋環(huán)境變化規(guī)律模型，與聲吶水下探測(cè)相結(jié)合，用以研究海洋環(huán)境不確定性對(duì)聲吶檢測(cè)精度確定性的影響；李玉偉［24］、岳雷［25］從信號(hào)處理模擬研究海洋環(huán)境噪聲對(duì)聲吶水下檢測(cè)概率、方位角偏離的影響；趙志允［26］、He［27］通過構(gòu)建聲吶檢測(cè)概率模型，綜合利用已有的水聲傳播、環(huán)境噪聲級(jí)、混響級(jí)預(yù)報(bào)模型綜合評(píng)估海洋環(huán)境背景場(chǎng)對(duì)聲吶探測(cè)距離和預(yù)報(bào)精度的影響；楊小小［28］、戴明強(qiáng)［29］則是從聲吶本身性能參數(shù)及人為因素，利用 ADC（application dependability capability）法、蟻群算法、層次分析法等對(duì)聲吶探測(cè)進(jìn)行評(píng)估。

HMG法側(cè)重評(píng)估“水聲環(huán)境效應(yīng)”對(duì)聲吶探測(cè)效能的影響，文章從聲吶方程出發(fā)闡述HMG法原理，并以東海某一淺水調(diào)查海域?yàn)檠芯繀^(qū)域，采用實(shí)測(cè)數(shù)據(jù)構(gòu)建海洋環(huán)境背景場(chǎng)，進(jìn)而利用實(shí)測(cè)傳播損失數(shù)據(jù)、環(huán)境噪聲數(shù)據(jù)驗(yàn)證水聲模型，并模擬預(yù)報(bào)單因子環(huán)境參數(shù)（風(fēng)速、聲速剖面等）變化對(duì)被動(dòng)聲吶探測(cè)距離的影響。

文中的研究貼近真實(shí)海洋環(huán)境，結(jié)果可靠性強(qiáng)，可為被動(dòng)聲吶水下探測(cè)及水下航行器的隱蔽航行研究提供參考。

1　HMG法預(yù)報(bào)原理

聲吶方程通?？煞譃橹?被動(dòng)形式［30-31］

式中:TL，BL，SL，TS，AG，DT和NL分別表示傳播損失、背景干擾級(jí)、聲源級(jí)、目標(biāo)強(qiáng)度、聲吶陣列的增益、檢測(cè)閾和環(huán)境噪聲級(jí)。其中，BL取值與NL和RL（混響干擾級(jí)）相關(guān)

式中: SL，TS，AG和DT為聲吶及水下目標(biāo)的固有屬性，不受海洋環(huán)境影響；TL和BL則反映水聲信道聲波傳輸能力，受海洋環(huán)境影響巨大，是“水聲環(huán)境效應(yīng)”的集中體現(xiàn)。當(dāng)聲吶系統(tǒng)性能已知，能否探測(cè)到水下目標(biāo)物取決于“水聲環(huán)境效應(yīng)”的影響，即TL和BL的大小。

水聲信道傳輸介質(zhì)通?？煞譃樗w、海底及海面，三者時(shí)空特征差異造成TL和BL的不同分布特征，進(jìn)而影響聲吶有效探測(cè)距離。故可將“水聲環(huán)境效應(yīng)”概括為水文環(huán)境效應(yīng)（hydrology）、氣象環(huán)境效應(yīng)（meteorology）及海底環(huán)境效應(yīng)（geology）之和，分別利用聲學(xué)模型定量研究3種具體環(huán)境中TL和BL的分布特征，即 HMG法整體構(gòu)成。用組合聲吶參數(shù)表示如下

式中: 對(duì)于被動(dòng)聲吶，l=1，m=1，n=0；對(duì)于噪聲干擾背景的主動(dòng)聲吶，l=2，m=1，n=0；對(duì)于混響干擾背景的主動(dòng)聲吶，l=2，m=0，n=1。

結(jié)合式（2）、式（4）可得被動(dòng)聲吶檢測(cè)域

下文分析中，假設(shè) AG=0。當(dāng)聲源級(jí)已知時(shí)，利用水聲傳播數(shù)值模型、噪聲模型分別模擬具體海洋環(huán)境中的傳播損失和環(huán)境噪聲譜級(jí)，則基于式（4）可得被動(dòng)聲吶HMG值，進(jìn)而利用式（5）預(yù)報(bào)其有效探測(cè)距離。

2　海洋環(huán)境背景場(chǎng)

文章選取東海淺水海域水聲調(diào)查中某一實(shí)測(cè)調(diào)查航線為海洋環(huán)境背景場(chǎng)［12，32］（見圖1）。

圖1　水聲調(diào)查航線位置及邊界要素分布Fig.1　Location of underwater acoustic survey route and distribution of boundary parameters

2.1水聲調(diào)查概況

水聲調(diào)查以走航觀測(cè)與定點(diǎn)觀測(cè)相結(jié)合的方式進(jìn)行。調(diào)查項(xiàng)目含聲傳播損失、海洋環(huán)境噪聲、GPS數(shù)據(jù)、溫度、鹽度、密度、水深、海流、海面氣象和海表風(fēng)浪等。

水聲探測(cè)浮標(biāo)系統(tǒng)通過海面浮標(biāo)下掛水聽器陣（含 8個(gè)不同深度的水聽器），下端以鉛魚固定。調(diào)查船拖曳近似無指向性人工聲源沿調(diào)查航線航行，并發(fā)射不同頻率的聲信號(hào)，其中580 Hz聲信號(hào)寬度為500 ms，聲源級(jí)約為200 dB。

水聲探測(cè)浮標(biāo)系統(tǒng)與調(diào)查船均配備 GPS定位系統(tǒng)，可實(shí)時(shí)觀測(cè)兩者位置及相對(duì)距離。實(shí)時(shí)控制拖曳聲源，使其保持在1/2倍水深的深度處。

2.2海洋環(huán)境背景場(chǎng)

以圖1所示的調(diào)查航線為例進(jìn)行分析。天氣晴，海面浪高2～3 m，東北風(fēng)5～6級(jí)。海底底質(zhì)以粉砂質(zhì)粘土為主，水深變化在53～65 m之間。

調(diào)查航線長(zhǎng)約 50 km，在其起點(diǎn)、終點(diǎn)處利用CTD （conductive temperature depth）進(jìn)行溫、鹽、深剖面觀測(cè)，走航調(diào)查過程中利用XBT（expendable bathythermograph）觀測(cè)溫、深剖面?；贏kima插值法［33］、GDEM（generalized digital environment model）聲速計(jì)算法［34］可得調(diào)查航線聲速斷面如圖2。

圖2　海洋環(huán)境背景聲速場(chǎng)Fig.2　The sound speed field of ocean environment background

聲速場(chǎng)整體呈現(xiàn)3層結(jié)構(gòu)，可分為淺層、中層和深層，淺層和深層聲速垂向變幅較?。恢袑訛槁曀佘S層，聲速大小向下迅速減小，分布于32～40 m深度間。

3　被動(dòng)聲吶淺海探測(cè)距離預(yù)報(bào)

首先，基于文中構(gòu)建的海洋環(huán)境背景場(chǎng)，采用實(shí)測(cè)數(shù)據(jù)驗(yàn)證水聲模型模擬效果，進(jìn)而模擬預(yù)報(bào)海洋環(huán)境背景場(chǎng)單因子變化時(shí)被動(dòng)聲吶的探測(cè)距離。

3.1模型配置及驗(yàn)證

分析文中的海洋環(huán)境背景場(chǎng)，調(diào)查航線水深較小，且海底地形變化緩慢，故可采用水平不變的 Kraken簡(jiǎn)正波理論模型［35］模擬其傳播損失分布特征。簡(jiǎn)正波理論采用快速有限差分法精確解算實(shí)數(shù)本征值，精度高、速度快，適于工程應(yīng)用。

海表風(fēng)、浪對(duì)傳播損失場(chǎng)的影響主要考慮海表粗糙度，忽略水體不均勻性的作用。文章參考文獻(xiàn)［36］海面粗糙度的計(jì)算公式

式中:κ為 vonkarman常數(shù)；z0為海面粗糙度，通常取0.35；C10為阻力系數(shù)，且

式中，U10為海面 10 m高處風(fēng)速，經(jīng)計(jì)算可得海底粗糙度計(jì)算參照文獻(xiàn)［37］。

水體聲吸收系數(shù)αw=0.000 1dB/λ。設(shè)粉砂質(zhì)粘土厚度10 m，ρs=1.650g/cm3，cs=1 700 m/s，往下為無限大基底層，

淺海水聲傳播損失垂向差異小，聲能匯聚特征不明顯，故將其垂向平均，得到的平均傳播損失隨水平距離的變化曲線與實(shí)測(cè)數(shù)據(jù)對(duì)比如圖3。分析圖3中實(shí)測(cè)數(shù)據(jù)可知，傳播損失隨距離變化總趨勢(shì)逐漸增加，局部存在波動(dòng)現(xiàn)象，同一距離不同水聽器觀測(cè)結(jié)果聲強(qiáng)差異在20 dB，Kraken模型輸出結(jié)果與實(shí)測(cè)數(shù)據(jù)基本吻合。

圖3　模型輸出與實(shí)測(cè)傳播損失對(duì)比Fig.3　Transmission loss comparison between acoustic model and experiment

淺海環(huán)境噪聲源多樣，主要由海表風(fēng)浪、海流、船舶和海洋生物等混合形成，難以用一個(gè)“全能”數(shù)學(xué)模型來精確表述［38-39］。聲吶探測(cè)過程中可預(yù)先根據(jù)感興趣的水下目標(biāo)物輻射噪聲信號(hào)頻帶分布，選擇影響該頻帶區(qū)間的主要環(huán)境噪聲源，并基于對(duì)應(yīng)的噪聲源經(jīng)驗(yàn)公式進(jìn)行預(yù)報(bào)。

100～10 000 Hz頻帶內(nèi)的主要噪聲源為船舶、海上施工及刮風(fēng)下雨。以下預(yù)報(bào)仿真中主要考慮風(fēng)生噪聲，故文章引用 Hamson［40］提出的經(jīng)驗(yàn)公式，李整林等［41］對(duì)其進(jìn)行了修正，使其更適合于國(guó)內(nèi)近海風(fēng)生噪聲源級(jí)預(yù)報(bào)，則

式中: f為頻率；ν為風(fēng)速。

圖4為6級(jí)風(fēng)情況下的實(shí)測(cè)海洋環(huán)境噪聲譜級(jí)與經(jīng)驗(yàn)公式（8）、（9）計(jì)算結(jié)果對(duì)比，兩者基本相符，實(shí)測(cè)數(shù)據(jù)微弱偏大，這可能源于水聲調(diào)查海域有往來船舶影響。

圖4　模型輸出與實(shí)測(cè)環(huán)境噪聲譜級(jí)對(duì)比Fig.4　Spectral level of ambient noise comparison between acoustic model and experiment

3.2探測(cè)距離預(yù)報(bào)

設(shè)水下目標(biāo)物為155 dB的窄帶信號(hào)（400 Hz）水下航行器，置于航線起點(diǎn)30 m深度處（見圖1）。設(shè)檢測(cè)域 DT≥6 dB時(shí)，被動(dòng)聲吶可成功檢測(cè)到水下航行器。則基于第2節(jié)構(gòu)建的海洋環(huán)境背景場(chǎng)，下文分析單因子（風(fēng)速、躍層）變化時(shí)對(duì)聲吶探測(cè)距離預(yù)報(bào)的影響。

3.2.1不同風(fēng)速時(shí)被動(dòng)聲吶探測(cè)距離預(yù)報(bào)

圖5為風(fēng)速變化時(shí)被動(dòng)聲吶探測(cè)距離預(yù)報(bào)曲線，可知風(fēng)速對(duì)被動(dòng)聲吶影響巨大，風(fēng)速微弱時(shí)，探測(cè)距離可達(dá)近百公里；風(fēng)速為 15 m/s時(shí)（7級(jí)風(fēng)），理論上探測(cè)距離僅約為10 km。

圖5　被動(dòng)聲吶探測(cè)距離隨風(fēng)速變化曲線Fig.5　Curve of passive sonar detection range versus wind speed

事實(shí)上，風(fēng)速通過作用于傳播損失與海洋環(huán)境噪聲共同影響被動(dòng)聲吶的工作效能，圖6為海表不同風(fēng)速大小時(shí)，傳播損失隨距離的變化曲線，海表風(fēng)力級(jí)改變海面粗糙度，同時(shí)改變海洋淺層水體介質(zhì)的均勻性，進(jìn)而改變水聲信道的聲波傳輸能力。 風(fēng)速增長(zhǎng)，海洋環(huán)境噪聲級(jí)迅速增加（見圖7），掩蓋水下目標(biāo)物的聲信號(hào)，導(dǎo)致被動(dòng)聲吶有效探測(cè)距離減小。

圖6　不同風(fēng)速時(shí)傳播損失隨距離變化曲線Fig.6　Curves of transmission loss versus range at different wind speed

圖7　不同風(fēng)速時(shí)400 Hz環(huán)境噪聲級(jí)Fig.7　Ambient noise level at 400 Hz with different wind speed

3.2.2不同躍層時(shí)被動(dòng)聲吶探測(cè)距離預(yù)報(bào)

文中研究海域位于東海大陸架近海區(qū)域，該區(qū)域冬季出現(xiàn)逆溫躍層，春末夏初溫度垂向混合較均勻，夏秋季節(jié)為典型的淺海正溫躍層［42］。故設(shè)模擬躍層類型如圖8所示，分別為類型1（強(qiáng)負(fù)躍層型，即聲速向下迅速減小）、類型2（弱負(fù)躍層型）、類型3（無躍層型）、類型4（弱正躍層型）、類型5（強(qiáng)正躍層型）。

圖8　5種類型聲速剖面Fig.8　Five sound speed profiles

躍層類型通過影響傳播損失作用于被動(dòng)聲吶探測(cè)距離。水下聲傳播中，聲線總是彎向聲速減小的方向，在深?？蓪?dǎo)致匯聚區(qū)與聲影區(qū)的形成。淺海躍層越強(qiáng)，對(duì)聲線傳播方向作用越大，導(dǎo)致聲波與邊界作用次數(shù)越多，邊界損失越大。圖 9所示，躍層類型 3聲吶有效探測(cè)距離最遠(yuǎn)，約47 km，傳播損失最小（見圖10），該類型聲速剖面（均勻聲速剖面）中，聲能邊界損失最小；類型1、類型5聲速剖面環(huán)境中被動(dòng)聲吶探測(cè)距離最小，傳播損失最大（圖10）；類型2、類型4為弱躍層型，被動(dòng)聲吶探測(cè)距離與傳播損失居中。

圖9　被動(dòng)聲吶探測(cè)距離隨躍層類型的變化曲線Fig.9　Curve of passive sonar detection range versus different thermocline

圖10　不同躍層時(shí)傳播損失隨距離變化曲線Fig.10　Curves of transmission loss versus range at different thermocline

已有的研究工作認(rèn)為，負(fù)躍層聲吶探測(cè)距離最小，正躍層探測(cè)距離最大，與文中研究結(jié)論存在一定偏差。源于其將海表假設(shè)為彈性釋放界面，不考慮表層聲能衰減較小所致。文中模擬的是海表5～6級(jí)風(fēng)時(shí)不同躍層類型對(duì)聲吶探測(cè)距離的影響，海面風(fēng)速較大導(dǎo)致上邊界傳播損失增大，形成與底邊界相當(dāng)?shù)乃p作用，故類型3被動(dòng)聲吶探測(cè)距離最大。

3.2.3其他因素對(duì)被動(dòng)聲吶探測(cè)距離的預(yù)報(bào)

當(dāng)被動(dòng)聲吶工作海域內(nèi)存在往來船舶時(shí)，應(yīng)根據(jù)其噪聲源級(jí)經(jīng)驗(yàn)公式計(jì)算其噪聲源級(jí)，并結(jié)合風(fēng)生噪聲級(jí)解算背景噪聲干擾級(jí)

式中，SLW，SLS分別為風(fēng)生噪聲、船舶噪聲源級(jí)。

海底底質(zhì)類型、沉積層厚度對(duì)傳播損失影響較大，利用HMG方法預(yù)報(bào)被動(dòng)聲吶探測(cè)距離時(shí)，應(yīng)基于其具體分布情形，調(diào)整水聲傳播模型內(nèi)的邊界參數(shù)。

海洋流噪聲影響頻率通常不超過100 Hz，文中模擬400 Hz窄帶水下目標(biāo)物，故不考慮其影響。

4　結(jié)束語

被動(dòng)聲吶水下探測(cè)效能受“水聲環(huán)境效應(yīng)”影響顯著，其有效探測(cè)距離與背景噪聲干擾級(jí)、傳播損失緊密相關(guān)。文章利用東海水聲調(diào)查某一航線，構(gòu)建海洋環(huán)境背景場(chǎng)；在此基礎(chǔ)上，結(jié)合水聲數(shù)據(jù)處理及聲吶方程進(jìn)一步闡述HMG法預(yù)報(bào)原理，給出聲吶組合參數(shù)表達(dá)式；最后，基于水聲調(diào)查實(shí)測(cè)數(shù)據(jù)對(duì)聲學(xué)模型模擬效果進(jìn)行驗(yàn)證，并模擬預(yù)報(bào)不同海面風(fēng)速、水體聲速剖面時(shí)的聲吶有效探測(cè)距離，給出其相應(yīng)的物理形成機(jī)制。文章研究?jī)?nèi)容緊密集合實(shí)際海洋環(huán)境，理論性強(qiáng)、可靠性高，對(duì)被動(dòng)聲吶水下探測(cè)及水下航行器的隱蔽航行具有應(yīng)用價(jià)值。

由于海洋水聲環(huán)境的不確定性及時(shí)空多變性等特點(diǎn)，對(duì)模擬真實(shí)海洋環(huán)境中被動(dòng)聲吶探測(cè)效能提出重大考驗(yàn)，也對(duì)水聲模型的選配及其適用性提出很高要求，該領(lǐng)域的工作還有待于進(jìn)一步深化。

［1］ 笪良龍.海洋水聲環(huán)境效應(yīng)建模與應(yīng)用［M］.北京: 科學(xué)出報(bào)社，2012: 1-2.

［2］ Harry D F.Predicting Sonar Performance Using Observations of Mesoscale Eddies［J］.Journal of the Acoustical Society of America，2006，120（2）: 3060-3069.

［3］ 秦繼興，張仁和，駱文于，等.大陸坡海域二維聲傳播研究［J］.聲學(xué)學(xué)報(bào)，2014，39（2）: 145-153.Qian Ji-xing，Zhang Ren-he，Luo wen-yu，et al.Two Dimensional Sound Propagation over a Continental Slope［J］.Acta Acustica，2014，39（2）: 145-153.

［4］ 高飛，張新睿，孫磊，等.白令海西部小區(qū)域聲傳播特征研究［J］.聲學(xué)技術(shù)，2015，34（4）: 306-311.

Gao Fei，Zhang Xin-rui，Sun Lei，et al.The Analysis of Acoustic Propagation Characteristics in the Small Area West of Bering Sea［J］.Technical Acoustics，2015，34（4）: 306-311.

［5］ 于源，鄢社鋒，侯朝煥.混響背景下主動(dòng)探測(cè)聲納性能預(yù)報(bào)［J］.魚雷技術(shù)，2011，19（3）: 192-194.

Yu Yuan，Yan She-feng，Hou Chao-huan.Performance Prediction of Active Sonar Under Reverberation［J］.Torpedo Technology，2011，19（3）: 192-194.

［6］ David R B，Michael J B.Depth Dependence of Winddriven，Broadband Ambient Noise in the Philippine Sea［J］.Journal of the Acoustical Society of America，2013，133（1）: 62-71.

［7］ Bernard C，James D，William L，et al.Coherent Bottom Reverberation: Modeling and Comparisons with At-sea Measurements［J］.Journal of the Acoustical Society of America，2004，116（4）: 1985-1994.

［8］ 李勝全，李宏武，高飛，等.主動(dòng)聲納探測(cè)效能評(píng)估的HMG法［J］.海洋測(cè)繪，2015，35（3）: 67-70.Li Sheng-quan，Li Hong-wu，Gao Fei，et al.HMG Method for Evaluation Active Sonar Detection Efficiency［J］.Hydrographic Surveying and Charting，2015，35（3）: 67-70.

［9］ Mellerg L E，Robinson A R，Botseas G.Modeled Time Variability of Acoustic Propagation Through a Gulf Stream Mender and Eddies［J］.Journal of the Acoustical Society of America，1990，87（3）: 1044-1054.

［10］ Alexey A S，James F L，Lin Y T.Three-dimensional Coupled Mode Analysis of Internal-wave Acoustic Ducts［J］.Journal of the Acoustical Society of America，2014，135（5）: 2497-2512.

［11］ Martin S，MichaelI B P，Paul H，et al.Effects of Ocean Thermocline Variability on Nonecoherent Underwater Acoustic Communications［J］.Journal of the Acoustical Society of America，2007，12（4）: 1895-1908.

［12］ 潘長(zhǎng)明，高飛，孫磊，等.淺海溫躍層對(duì)水聲傳播損失場(chǎng)的影響［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2014，35（4）: 312- 318.Pan Chang-ming，Gao Fei，Sun Lei，et al.The Effects of Shallow Water Thermocline on Water Acoustic Transmission Loss［J］.Journal of Harbin Engineering University，2014，35（4）: 312-318.

［13］ 程廣利，張明敏.淺海不確定聲場(chǎng)的隨機(jī)多項(xiàng)式展開法研究［J］.聲學(xué)學(xué)報(bào)，2013，38（3）: 294-299.Cheng Guang-li，Zhang Ming-min.On Polynomial Chaos Expansion Method for the Method for the Uncertain Acoustic Field in Shallow Water［J］.Acta Acustical，2013，38（3）: 294-299.

［14］ Laurie T F，Michael D C，John S P.Source Localization in Noisy and Uncertain Ocean Environment［J］.Journal of the Acoustical Society of America，1997，101（6）: 3593-3599.

［15］ Lin J H，Gao T F.Model of Wind-generated Ambient Noise in Stratified Shallow Water［J］.Chinese Journal of Oceanology Limnology，2005，23（2）: 144 -151.

［16］ MichaelI J B.Theory of the Directionality and Spatial Coherence of Wind-driven Ambient Noise in a Deep Ocean with Attenuation［J］.Journal of the Acoustical Society of America，2013，134（2）: 950 -958.

［17］ 過武宏，笪良龍，趙建昕.地聲參數(shù)及傳播損失不確定性估計(jì)與建模［J］.應(yīng)用聲學(xué)，2015，34（1）: 71-78.Guo Wu-hong，Da Liang-long，Zhao Jian-xin.Estimation and Modeling of Geoacoustic Parameters and Transmission Loss Uncertainty［J］.Journal of Applied Acosutics，2015，34（1）: 71-78.

［18］ Tapper F S，Spiesberger J L，Wolfson M A.Study of a Novel Range-Dependent Propagation Effect with Application to the Axial Injection of Signals from the Kaneohe Source［J］.Journal of the Acoustical Society of America，2002，111（2）: 757-762.

［19］ 袁駿，張明敏，肖卉.雙基地聲納海底界面?zhèn)认蛏⑸浞治雠c仿真［J］.火力與指揮控制，2010，35（7）: 133-138.Yuan Jun，Zhang Ming-min，Xiao Hui.Analysis and Simulation of Bistatic Sonar Seafloor Interface Side E-Scattering［J］.Fire Control & Command Control，2010，35（7）: 133-138.

［20］ Alec J D，Robert D M，Amos L M.Predicting the Environment Impact of Active Sonar［J］.AIP Conference Proceedings，2004，728: 280-287.

［21］ Paul B，Jeffery K，Shawn K.Adaptive Sonar Detection Performance Prediction in an Uncertain Ocean［J］.Journal of the Acoustical Society of America，2003，113（1）: 2253-2263.

［22］ Sha L W，Loren W N.Effects of Environmental Uncertainties on Sonar Detection Performance Prediction［J］.Journal of the Acoustical Society of America，2005，117（4）: 1942-1953.

［23］ Fan L，Guo S M，Chen Y M.Sonar Detection Range Index Estimation Approach in Uncertain Environments［J］.AIP Conference Proceeding，2010，1272: 375-382.

［24］ 李玉偉，姜可宇，黃建波.探雷聲納水雷發(fā)現(xiàn)概率仿真研究［J］.聲學(xué)技術(shù)，2015，31（5）: 505-509.Li Yu-wei，Jiang Ke-yu，Huang Jian-bo.Simulation Research on the Detection Probability of Mine-detection Sonar［J］.Technical Acoustics，2015，31（5）: 505-509.

［25］ 岳雷.線性調(diào)頻脈沖串信號(hào)及其檢測(cè)性能分析［J］.魚雷技術(shù)，2015，23（4）: 285-290.Yue Lei.Pulse Trains of Linear Frequency-modulated Signal and Its Detection Performance［J］.Torpedo Technology，2015，23（4）: 285-290.

［26］ 趙志允，孫明太，任東彥.吊放聲納主動(dòng)工作方式探測(cè)域建模與仿真［J］.指揮控制與仿真，2013，35（6）: 49-54.Zhao Zhi-yun，Sun Ming-tai，Ren Dong-yan.Modeling and Simulation of Dipping Sonar Active Detection Tegion［J］.Command Control & Simulation，2013，35（6）: 49-54.

［27］ He X Y，Cai Z M，Lin J Y，et al.The Simulation Research on Forecasting the Detection Range of Active Sonar［J］.Journal of System Simulation，2003，15（9）: 1304-1308.

［28］ 楊小小，綦輝，陳磊.基于改進(jìn)ADC法的潛艇反潛作戰(zhàn)方案效能評(píng)估［J］.海軍航空工程學(xué)院學(xué)報(bào)，2014，29（3）: 285-292.

［29］ 戴明強(qiáng)，李旺，項(xiàng)浩.改進(jìn)蟻群算法在潛艇規(guī)避機(jī)載聲納探測(cè)的路徑優(yōu)化應(yīng)用［J］.兵工自動(dòng)化，2010，29（9）: 9-14.Dai Ming-qiang，Li Wang，Xiang Hao.Application of Improved Ant Colony Algorithm in Path Optimization of Submarine Evading Search Dipping Sonar on Aircraft［J］.Ordnance Industry Automation，2010，29（9）: 9-14.

［30］ Urick R J.Principles of Underwater Sound（Third Edition）［M］.Los Altos，California，USA: Peninsula Pub，1983: 104-114.

［31］ Ainslie M A.Principles of Sonar Performance Modeling ［M］.Berlin: Springer & Chichester: Praxis，2010: 53-122.

［32］ 高飛，潘長(zhǎng)明，孫磊，等.溫躍層及其變化對(duì)被動(dòng)聲納檢測(cè)概率的影響［J］.應(yīng)用聲學(xué)，2014，33（2）: 138-144.Gao Fei，Pan Chang-ming，Sun Lei，et al.The Effects of Thermocline and Its Variability on the Detection Probability of Passive Sonar［J］.Journal of Applied Acoustics，2014，33（2）: 138-144.

［33］ Akima H．A New Method of Interpolation and Smooth Curve Fitting Based on Local Procedures［J］.Journal of Association Computer Mathematics，1970，17（4）: 589- 602．

［34］ Teague W J，Carron M J，Hogan P J.A Comparison Between the Generalized Digital Environmental Model and Levitus Climatologies［J］.Journal of Geophysical Research，1990，95: 7167-7183.

［35］ Porter M B，Reiss E L.A Numerical Method for Oceanacoustic Normal Modes［J］.Journal of the Acoustical Society of America，1984，76（1）: 244-252.

［36］ Wu H.Wind Stress and Surface Roughness at Air-sea Interface［J］.Geographical Reasearch，1969，21（5）: 707- 714.

［37］ 周艷蓮，孫曉敏，朱治林，等.幾種典型地表粗糙度計(jì)算方法與比較［J］.地理研究，2007，26（5）: 887-896.Zhou Yan-lian，Sun Xiao-min，Zhu Zhi-lin，et al.Comparative Research on Four Typical Surface Roughness Length Calculation Methods［J］.Geographical Research，2007，26（5）: 887-896.

［38］ 劉伯勝，雷家煜.水聲學(xué)原理［M］.第2版.哈爾濱: 哈爾濱工程大學(xué)出版社，2010: 248-249.

［39］ Martin S，MichaelI J B.Thirty Years of Progress in Application and Modeling of Ocean Ambient Noise［C］//AIP Conference Proceedings.Beijing，China，2012: 261-272.

［40］ Hamson R M.The Modelling of Ambient Noise Due to Shipping and Wind Sources in Complex Environments［J］.Applied Acoustics，1997，51（3）: 251-287.

［41］ 李整林，彭朝暉，何利.海洋環(huán)境噪聲源級(jí)經(jīng)驗(yàn)公式修正［J］.聲學(xué)技術(shù)，2010，29（6）: 36-37.Li Zheng-lin，Peng Chao-hui，He Li.Modification for the Empirical Relationship of Ambient Noise Source Level in Shallow Water［J］.Technical Acoustics，2010，29（6）: 36-37.

［42］ 葛人峰，郭景松，于非，等.黃、東海陸架海域溫度垂直結(jié)構(gòu)類型劃分與溫躍層分析［J］.海洋科學(xué)進(jìn)展，2006，24（4）: 424-435.Ge Ren-feng，Guo Jing-song，Yu Fei，et al.Classification of Vertical Temperature Structure and Thermocline Analysis in the Yellow Sea and East China Sea Shelf Sea Areas［J］.Advances in Marine Science，2006，24（4）: 424-435.

（責(zé)任編輯: 楊力軍）

Detection Range Prediction of Passive Sonar in Shallow-water Area Based on the HMG Method

GAO Fei1，PAN Chang-ming1，LI Sheng-quan1，SUN Lei1，CHEN Fu-sen1，2
（1.Naval Institute of Hydrographic Surveying and Charting，Tianjin 300061，China；2.College of Meteorology and Oceanography，PLA University of Science and Technology，Nanjing 211101，China）

A research is conducted by combining experiment data and simulation model based on the hydrology meteorology geology（HMG） method to evaluate the influence of the underwater acoustic environment effect on sonar detection range.Firstly，the HMG prediction theory is discussed based on the sonar equation，and a sonar parameters combination expression is proposed.Then，the ocean environment background field is constructed combining with shallow-water acoustic survey in the East China Sea，and its acoustic parameters distribution of the water and boundary are analyzed.Finally，the Kraken's numerical acoustic transmission model and the empirical model of wind-generated noise are verified with measured data，and the passive sonar detection range is simulated for different wind speeds and sound speed profiles with the HMG method.The results show that: 1） the HMG method can simulate the effect of the underwater acoustic environment effect on sonar detection range；2） the transmission loss and ambient noise increase with wind speed at sea surface，which decreases the sonar detection range；and 3） the sound speed profile influences detection range by affecting acoustic boundary loss，and the detection range gets the biggest in vertically uniform sound speed environment.This research may provide a reference for the underwater passive sonar detection and the stealth of an underwater vehicle.

underwater vehicle；hydrology meteorology geology（HMG） method；passive sonar detection；underwater acoustic environment effect；underwater acoustic survey；sound speed profile

TJ630.34；P733.23

A

1673-1948（2016）03-0177-07

10.11993/j.issn.1673-1948.2016.03.004

2016-02-29；

2016-04-21.

國(guó)家自然科學(xué)基金項(xiàng)目（41276088）；國(guó)家海洋公益性行業(yè)科研項(xiàng)目（2012050007-7）.

高飛（1988-），男，碩士，主要研究方向?yàn)楹Ｑ笏暛h(huán)境效應(yīng).