JPS6322275B2 - - Google Patents

Info

Publication number
JPS6322275B2
JPS6322275B2 JP17115080A JP17115080A JPS6322275B2 JP S6322275 B2 JPS6322275 B2 JP S6322275B2 JP 17115080 A JP17115080 A JP 17115080A JP 17115080 A JP17115080 A JP 17115080A JP S6322275 B2 JPS6322275 B2 JP S6322275B2
Authority
JP
Japan
Prior art keywords
display
display device
attenuation
output
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
JP17115080A
Other languages
English (en)
Other versions
JPS5794671A (en
Inventor
Hidehiro Shiraishi
Hiroshi Iino
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Furuno Electric Co Ltd
Original Assignee
Furuno Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Furuno Electric Co Ltd filed Critical Furuno Electric Co Ltd
Priority to JP17115080A priority Critical patent/JPS5794671A/ja
Publication of JPS5794671A publication Critical patent/JPS5794671A/ja
Publication of JPS6322275B2 publication Critical patent/JPS6322275B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/96Sonar systems specially adapted for specific applications for locating fish

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】 この発明は魚群等の反射体から帰来した反射波
に基づいて、魚群の散乱強度を図示することを目
的とする。
従来この種装置としては、例えば特公昭48−
26299号公報に見られるように、反射波強度、船
の航走距離、振動子の指向角に基づき魚群の密度
並びに魚群総量をデイジタル表示する技術が存在
する。しかし、この装置は計算結果を数値で表示
するため統計資料として直観的に利用することは
困難である。
本発明は反射体の密度及びその深度分布等の解
析を迅速に行なえる装置を提供するものである。
具体的には微少深度範囲ごとの反射波の積分結
果を各深度に対応してグラフ表示しようとするも
のである。
以下本発明を魚群反射波の解析に実施した場合
について説明する。
本発明の原理 船から超音波を魚群に向けて発射し、その反射
波を分析すれば、何らかの形で魚量に関係する情
報が得られる。特に反射波のレベルは魚群密度を
代表する重要な情報を有する。しかし、反射波レ
ベルは音源から反射物体までの距離の相違によつ
て大きく異なる。このため第1に反射波の減衰補
正を行なう必要がある。従来この減衰補正の方法
としては、単に指数関数的に出力電圧が変化する
関数発生器により利得が制御される増巾器が用い
られているが、定量的な魚群密度測定を行なうに
は不充分である。
一般に、超音波の水中における減衰特性TLは、 TL=20log r+α/1000r ……(1) で与えられる。
ここで、 r:音源からの距離でr=Δr×in α:吸収減衰係数で、周波数依存性を持ち、次
の様な実験式が得られている。
α=43.52/7000+2+0.000332(dB/Km)……
(2) :周波数(kHz) (1)式で第1項は拡散による減衰を示し、第2項
は吸収による減衰を示す。通常第1項の拡散減衰
に対する補正は時定数を利用した充電又は放電曲
線関数発生器を利用して補正を行なつているが、
第2項の吸収減衰に対しては、厳密な補正を行な
つていない。
本発明の一目的は、(1)式に示される減衰特性に
従う厳密な補正を行なうことである。
魚群密度の測定 今第1図に示すように、反射係数がTSの魚が
単位体積当りn尾の密度で分布する散乱層Fを考
える。又この散乱体の散乱強度をSVとする。一
般に音源から距離rの点における音波の強度Iiは Ii=I0b(θ,φ)1/r2e-r ……(3) 微少単位体積ΔVの散乱体からの反射強度Irは Ir=IisvΔV=IitsnΔV ……(4) 音源における反射強度は ΔI=Irb(θ,φ)1/r2e-r ……(5) =I0b2(θ,φ)1/r4e-2rsvΔv ……(6) ここで I0:音源から単位距離の点の強度 b:指向性強度 β:α/1000r=10logeβr ……(6A) 微少体積ΔVを音源から見込む立体角をΔΩとし、
音速C、パルス巾τとすると ΔV=r2ΔΩcτ/2 ……(7) (6)式から ΔI=I0b2(θ,φ)1/r2e-2r svcτ/2ΔΩ…
…(8) 従つて、振動子の指向角内に含まれる散乱体から
の反射強度は I=∫I0b2(θ,φ)1/r2-2r svcτ/2dΩ…
…(9) =I01/r2e-2r svcτ/2 40∫b2(θ,φ
)dΩ ……(10) =I01/r2e-2r svcτ/2ψ ……(11) =I01/r2e-2r ntscτ/2ψ ……(12) ここでψは、ビーム巾で指向性定数である。
10logI=10logI0−(20logr+2αr/1000) +10log sv+10log(cτ/2ψ) ……(13) 一方音源の送波強度SLは次のように定義され
る。
SL=10logI0/Iref ……(14) ここで Iref:基準強度 受信エコーレベルをEL=10log1/Irefとすると EL=SL−(20logr+2α/1000r)+sv +10log(cτ/2ψ) ……(15) =SL−(20logr+2α/1000r)+TS +10logn+10log(cτ/2ψ)……(16) 魚群密度 振動子受波感度をER=20logErとするとき、振
動子端における出力V0を求める。
受波点における音圧Pは 振動子の出力電圧Vは入射音波の音圧Pに比例
し、 ここで P0=音源電圧 Pref=基準音圧(1μbar) (18)及び(19)式には、魚群量に関する情報と
してsv又はntsを含んでいる。又対数表現すると、 20logV=SL+ER−(20logr+2α/1000r) +10log(cτ/2ψ)+sv ……(20) と表わす事が出来る。
式(18)に示す様に、魚群エコーの受信出力電
圧は散乱強度svの√ に比例している。
又、散乱強度svと魚鮮密度nの関係は sv=nts=nσ/4π ……(21) であるから、一尾当りの反射係数ts又は散乱断面
積σがわかれば、出力電圧Vから、魚群密度nを
求める事が出来る。
式(18)又は(20)において、 P0/Pref(SL),Er(ER)及c,τ,ψはすべ
てシステム定数で既知である。従つて、音波の伝
搬減衰に関する項、1/re-rを補正すれば、出力 電圧Vを推定する事により、距離rに無関係にsv
を求めることが出来る。。
1/re-rの補正は一般にTVG(TIME VARIED GAIN)と呼ばれ、距離に対応し、re〓rの利得
特性を持つ可変利得増巾器により行われる。
そのTVG増巾器の特性をEg re〓rとすれば増巾
器の出力電圧V0は、 従つて SV=V0 2/(Eg・Er・Po/Pref)cτ/2ψ ……(23) =V0 2/10{(EG+ER+SL)/10}・cτ/2ψ ……(24) 又魚群密度n(尾/m3)は n=V0 2/10{(EG+ER+SL+TS)/10}・cτ/2ψ……(25) EG=受信アンプゲイン dB ER=受波感度 dB SL=ソースレベル dB TS=ターゲツトストレングス dB τ=パルス巾 sec c=音速 m/sec ψ=指向性定数 EG=20logEg 0dB=1V/1μBar 0dB=1μBar/1m 例えば、 EG:0dB ER:−75dB SL:120dB τ :3msec ψ :0.062 c :1500m/sec とし、ある魚群エコーとしてV=1Vrmsの出力
が得られたとする。
この魚群を構成する魚単体の反射係数T sが
−45dBであれば、この魚群の散乱強度及び密度
は式(23)と(25)より SV=−36.5dB n=7.2尾/m2 と求められる。
実際の魚量測定においては、発射パルス毎のデ
ータにはバラツキがあるので、何マイルか毎の平
均値を求める。積分の範囲を水平方向(航走区
間)l、垂直方向(深度区間)R1〜R2とすれば、
その区間の平均散乱強度は、 =〔∫L 0(∫r2 r1sv(r,x)dr)dx〕/L(r2
−r1
……(26) =10log ……(27) となる。
本発明は式(27)の計算結果をグラフ化するこ
とを目的とするものであつて、以下実施例に従い
これを説明する。
第2図で、1は超音波帯域のパルス波を周期的
に発生する送信器、2はパルス波の周期及びパル
ス巾τを設定する送信制御器、3は送受切換器、
4は超音波送受波器、5は減衰器、6は増巾器7
の利得を距離に対し直線的に増加させる制御電圧
を発生するTVG電圧発生器である。8は増巾器
9内の受波信号から海底反射波を選別する海底検
出器で、後述のゲート手段12を制御して海底波
到来後の受波信号データを処理しないようにここ
で遮断する。10は入力アナログ信号をそのレベ
ル情報も含めてデジタル化するA―D変換器であ
る。11はA―D変換器の出力データを記憶する
記憶器、13は入力が予め設定したレベルを越え
る飽和信号である時には出力を送出して表示ラン
プ14を点灯する飽和検出手段。15はレベル設
定器16で設定されたレベル以下の信号を除外す
るスライス手段である。
17は減衰補正手段であつて、第8頁後半で説
明した1/re-rの項を補正するためのものであ
る。18は2乗回路であつて、魚群密度を算出す
る(24)式の分子に相当する出力を作成する。
19は魚群探索毎に受波されるデータの加算手
段である。
次に、20は船速信号発生器で、例えば超音波
によるドツプラ効果を利用したもの等を使えばよ
い。21は送信器1から発射周期に応じた信号を
作成する周期信号発生器である。22は船速信号
と周期信号との掛算手段、23は船の探索距離範
囲Lの設定手段で単位マイルである。24は掛算
手段22の算出結果が設定手段23の設定値に達
したときに出力を送出する比較手段である。25
は比較手段24の出力発生時の加算手段19の出
力をラツチするラツチ手段である。
26は送信制御器2の出力からパルス巾信号を
作成するτ信号発生手段、27は(24)式におけ
るシステム定数EG,ER,SL、を設定する定数
信号設定手段、28は同じく指向性定数ψの設定
手段、29は音速信号Cの設定手段である。30
はこれら所定の定数から(24)式の分母K K=10{(EG+ER+SL)/10}・cτ/2ψ ……(28) を算出するK算出手段である。
なお、本実施例ではK算出手段30には減衰器
5から減衰係数Gが入力される。したがつて、結
局 K=10{(EG+ER+SL+G)/10}・cτ/2ψ……(29) を得る。
31はラツチ手段25の出力をKで割る割算手
段、32はさらにその出力を探索範囲Lで割つて
平均散乱強度を算出する割算手段、33は割算出
力の対数を取る対数演算手段である。
この対数演算手段33の出力は、右側に示す記
録ブロツクの比較手段34に供給されている。比
較手段34の他の入力端にはレベル信号発生器3
5からの出力電圧が供給されている。レベル信号
発生器35は、レベルカウンター36が出力する
毎に一段づつレベルが増加する直流信号を発生す
る。レベルカウンタ36は比較手段24の出力に
よつて0にリセツトされ、後述の記録器37のキ
ーパルスが線路Aから到来する毎に登算して所定
入力数毎に出力をレベル信号発生器に供給する。
38はアドレスカウンターであつて、記録紙3
9の有効紙巾に対する深度区間のアドレス指定を
行なう。このため、カウンタ38は記録ペン40
が記録紙39の始端に達したときに発生されるキ
ーイングパルスでリセツトされ、記録ペンの進行
に従つて発生されるクロツクパルス発生器41の
出力を計数する。
アドレスカウンタ38の出力は左上方の記憶手
段11のアドレスを指定する。
37は記録器であつて、矢印方向に低速で移動
する記録紙39と、記録紙を横断走行する無端ベ
ルト42とこのベルトに取付けられた記録ペン4
0を有している。43は記録ペンに信号を与える
通電レール、44,45はベルト走行させるプー
リー、46は記録ペン40が記録紙39の始端に
来た時にパルスを発生するようにベルトに設けら
れたマグネツト、47は信号発生用リードスイツ
チである。48はタイマーであつて一分毎に記録
信号を送出して記録接片49に信号を供給する。
41はクロツクパルス発生器であつて記録ペン4
0の走行と同期してパルスを発生し、アドレスカ
ウンタ38を駆動する。
減衰補正手段 減衰補正は前述したように拡散減衰(20logr)
の補正と、吸収減衰(α/1000r)の補正に大別出 来る。
従来のTVG電圧による補正は主として拡散減
衰を補正するものであるが、深度に換算して約
300m程度までが限度であつた。本発明の実施例
ではこの従来のTVG技術に加えて電子計算機に
よる演算を行ない、300m以上の深度における精
密補正を可能とした。これを次に説明する。
再び第2図において、50は上述TVG電圧発
生器6の補正深度係数inおよびその限界値iiを設
定する補正限界設定手段である。一方別に周波数
設定手段51が設けられ、その出力によりα算出
手段52で(2)式によりαが算出される。53はさ
らに(6A)式によりe〓rを算出する計算手段であ
る。54は測定深度がTVG有効範囲内である場
合(r≦Δrii)は計算手段53の出力をそのまま
減衰補正手段17に供給するが、もし測定深度が
TVG有効範囲を越えている場合(r>Δrii)は、
深深度用TVG電圧発生手段56を駆動してその
出力と計算手段53の出力の積を減衰補正手段1
7に供給する判定手段である。56は深深度用
TVG電圧発生手段であつて深度inに対し in/ii ……(30) なる補正係数信号を発生する。
ここで、 ii:補正深度限界値 減衰補正手段17では、スライス手段15の出
力である、距離に応じて変化する信号と判定手段
54又は深深度用TVG電圧発生手段56の出力
との積出力 VEge〓rr を送出する。
装置の動作 上述のような装置において、第3図に示すよう
に船SPが矢印方向に速度Uで区間Lの海域を探
索する場合を考える。船は海底方向に周期τで周
波数Fの超音波パルスを1,2,3……mと云う
ように発射し、その反射波を受波して信号処理し
て行く。この場合深度r方向をΔrづつに区切つ
て、多数の横断層を考え、この横断層毎に反射信
号を進行方向に積分して行くのが本発明の基本的
な信号処理である。
第2図で送信器1が高周波パルスを発生する
と、このパルスは送受切換器3を介して超音波送
受波器4を励振しパルス巾τ、指向角θの超音波
パルスを水中に発射する。各深度に散在する魚群
からの反射波は浅深度のものから順次帰来して超
音波送受波器4、送受切換器3を経て減衰器5に
至る。減衰器5では信号が後段の増巾器で飽和し
ないようにレベル調整される。さらに信号は増巾
器7で増巾されるが、この増巾器の増巾度は
TVG電圧発生器6から加わる電圧によつて時間
の経過と共に増加するように設定されている。し
たがつて(1)式の第1項に示される(20log r)の
減衰分は補正される。この補正深度限界はTVG
電圧の最大値で決定され通常深度にして約300m
程度でああ。信号は次に増巾器9に至るが、ここ
ではレベル差を利用した周知の海底選別回路8に
よつて線路Bに海底信号が検出される。次に信号
はA―D変換器10でデイジタル信号に変換され
記憶器11に上記Δr区間毎に記憶される。書込
アドレス信号はアドレスカウンタ38から入力さ
れている。
さて、以後のデータ処理は中央演算処理装置
(以下CPUと略す)が行なうが、その手順を説明
する。
第2図の記録器37において、記録ペン40が
記録紙39の始端に到来した時にリードスイツチ
47が発生するパルスによつてアドレスカウンタ
38が駆動される。このカウンタの読出し指令に
よつて記憶器11のデータが読み出される。これ
等データのうち海底反射波以後の反射波はゲート
手段12で除去される。さらに信号が所定レベル
を越えて飽和すると正しい情報処理ができないの
でこれを飽和検出手段13で検出して表示ランプ
14で報知する。低レベルの信号は処理の対象外
とするためにスライス手段15で除去される。こ
れ等信号は減衰補正手段17で前述の如く超音波
の伝播減衰成分が補正される。この時点で(22)
式を代表する出力が得られたことになる。すなわ
ち深度i、第j番目のパルス周期におけるΔr深
度区間の正しい魚群情報が得られたことになる。
さらにこのデータは2乗手段18で2乗され
る。その結果(23)式又は(24)式の分子が得ら
れたことになる。
以上の演算が各深度区間Δr毎に船の航走方向
にくり返し行なわれる。これを39に表示する。
云うまでもなく新データの記憶器への書込みと
この記憶器からの読出しとは同期して処理する必
要はなく直接関係ない。したがつて、上述の演算
動作は、船から超音波が発射される毎に行なう必
要はない。
さて、くり返し演算の航走方向の範囲の決定で
あるが、本実施例では、探索パルス回数の総和と
船速との積が所定値Lに達した場合に比較手段2
4によつてラツチ手段を働かせ、くり返し演算を
停止せしめるようになされる。これは、例えば船
速が1/2に変化しても、範囲L内の探索パルスの
数は2倍になり結局その積は一定となる。すなわ
ち船速の影響を受けずに範囲測定を可能とするた
めである。
一方生の魚群信号は記憶回路11から読み出さ
れてD―A変換された後記録ペン40に与えら
れ、記録紙にF1,F2、として示すように通常の
魚群探知器の記録が行なわれる。
比較手段24の出力によつて、レベルカウンタ
ー36がリセツト状態となる。そしてその後に到
来するリードスイツチ47によるパルスを計算す
る毎に一段づつ値が大きくなるレベル信号がレベ
ル信号発生器35から比較手段34に送出され
る。
他方演算手段ではK算出手段30から(24)式
の分母を得て、割算手段31でラツチ手段出力の
V0 2を割る。したがつて、割算手段31の出力は
結局(24)式の散乱強度svを示すこととなる。さ
らに割算手段32でsvをLで割つて(26)式の平
均散乱強度を得る。さらに対数演算手段33で
svの対数を取つて(27)式のを得る。この出
力は比較手段34に供給される。
結局比較手段34は、積分範囲Lを経過した直
後の記録位置(直線P)を基準線としてペン一走
査毎に図中左方向に増大する平均散乱強度基準レ
ベルとペン存在位置に該当する深度区間Δrの対
数演算結果とをペン走査と同期して比較すること
になる。そうして両入力が一致した場合に記録ペ
ン40に出力を送出する。したがつて、記録紙3
9には平均散乱強度svが各深度区間毎に直線Pを
基準線とした左方向への偏位グラフQとして表示
されることとなる。
なお、記録紙上方の空欄には分時マークM,M
……が表示され、記録紙送り速度を船速に応動す
るようにしておけば記録紙上で船速変化を識別す
ることができる。これは探索パルスの航走方向に
対する粗密度を知る上で役立つ。
以上説明したように本発明によれば、探知海域
の平均反射強度を深度方向の分布も含めて直観的
に観察することができ、海域の魚類資源量の把握
に便ならしめることができる。
なお、上記実施例では、表示器として記録紙を
採用し、現示手段として一本ペンを採用したが、
この現示手段としては、この種記録紙に汎用され
ているマルチペンを利用できる。この場合は、マ
ルチペンの本数をすなわち深度分割区間数とすれ
ば便利である。
【図面の簡単な説明】
第1図は本発明の原理を説明するための説明
図、第2図は本発明の実施例を示すブロツク図、
第3図は本発明の動作を説明するための説明図で
ある。

Claims (1)

  1. 【特許請求の範囲】 1 探索波を発射しその反射波を時間軸対距離軸
    (深度軸)を有する表示器に表示する装置におい
    て、上記反射波の強度を上記探索波の伝播による
    減衰関数に従つて補正する補正手段と、上記距離
    軸を複数の時間区間に分割する手段と、該時間区
    間別に上記反射波を所定回数分積分する演算手段
    と、上記表示器に設けられ且つ上記距離軸方向に
    所定速度で主走査される現示主段と、上記時間軸
    方向に上記表示器の表示結果を副走査させる副走
    査手段と、上記時間区間毎の演算出力値を上記所
    定回数経過後の所定の主走査部位を基準とした上
    記時間軸方向への偏位信号に変換する手段と、該
    偏位記号を上記現示手段に供給する手段とを備え
    たことを特徴とする反射情報の表示装置。 2 上記補正手段が、拡散減衰分を補正する
    TVG電圧発生手段と、吸収減衰成分を補正する
    演算手段と、両手段の和手段を有したことを特徴
    とする請求範囲1に記載した反射情報の表示装
    置。 3 上記演算手段は、超音波の伝播減衰補正手段
    の出力を2乗する2乗手段と、その値を 10{(EG+ER+SL)/10}・Cτ/2ψ で除する手段とを有したことを特徴とする請求範
    囲1に記載した反射情報の表示装置。 4 上記表示器が一定方向に副走査する記録紙を
    有してなり且つ上記現示手段が上記記録紙を横断
    するごとく並べられた針状記録ペンよりなること
    を特徴とする請求範囲1に記載した表示装置。 5 上記表示器がブラウン管表示器よりなり上記
    現示手段が上記ブラウン管面を掃引する電子ビー
    ムより成ることを特徴とする請求範囲1に記載し
    た表示装置。
JP17115080A 1980-12-03 1980-12-03 Indicator of reflective information Granted JPS5794671A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17115080A JPS5794671A (en) 1980-12-03 1980-12-03 Indicator of reflective information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17115080A JPS5794671A (en) 1980-12-03 1980-12-03 Indicator of reflective information

Publications (2)

Publication Number Publication Date
JPS5794671A JPS5794671A (en) 1982-06-12
JPS6322275B2 true JPS6322275B2 (ja) 1988-05-11

Family

ID=15917909

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17115080A Granted JPS5794671A (en) 1980-12-03 1980-12-03 Indicator of reflective information

Country Status (1)

Country Link
JP (1) JPS5794671A (ja)

Also Published As

Publication number Publication date
JPS5794671A (en) 1982-06-12

Similar Documents

Publication Publication Date Title
US4376301A (en) Seismic streamer locator
US4389893A (en) Precision ultrasound attenuation measurement
US4532617A (en) System for locating a towed marine object
GB2207757A (en) Ultrasonic rangefinder
JPS5812548B2 (ja) 液体内の含有気体泡測定装置
US3381264A (en) Submarine topography
US4399702A (en) Method of measuring stress distribution in a solid body
US4319348A (en) Method and apparatus of surveying nodular targets on the sea floor
US7363177B2 (en) Apparatus and method for performing the time delay estimation of signals propagating through an environment
JPH07140243A (ja) 海洋等の測深方法
RU2451300C1 (ru) Гидроакустическая навигационная система
RU2272303C1 (ru) Способ определения глубин акватории и устройство для его осуществления
US7768875B2 (en) Underwater sounding apparatus capable of calculating fish quantity information about fish school and method of such calculation
RU2178880C2 (ru) Способ обработки сигналов, характеризующих волны, которые отражены или пропущены объектом, для исследования и анализа структуры объекта и устройство для исследования и анализа структуры объекта
JPS6322275B2 (ja)
RU178905U1 (ru) Многолучевой научный эхолот для учёта водных биоресурсов
KR101928799B1 (ko) 재질정보 생성을 위한 측심 소나의 후방산란 음압 정보 처리 장치
RU2123191C1 (ru) Эхолот
JPH021273B2 (ja)
Sathishkumar et al. Echo sounder for seafloor object detection and classification
JPS6312544B2 (ja)
JP2572452B2 (ja) 水中探知装置
JP2828259B2 (ja) 魚群探知機
RU98254U1 (ru) Многочастотный корреляционный гидроакустический лаг
Sawada et al. Precision calibration of echo sounder by integration of standard sphere echoes