水中音響ソナー方程式探知距離シミュレーター

水中音響ソナー方程式（Urick 1983）からアクティブ／パッシブソナーの最大探知距離をリアルタイム計算するツールです。音源レベル・標的反射・環境ノイズ・指向性・検出閾値・周波数・水深・水温を変えると、音速・吸収係数・許容伝搬損失・探知距離・波長が同時に求まり、潜水艦探知や海洋調査の感度設計が直感的に把握できます。

パラメータ設定

ソナー方式

アクティブは往復 2·TL、パッシブは片道 TL

音源レベル SL

dB re 1μPa @ 1m

標的反射 TS

大型潜水艦 +25、小型機雷 −20 程度

環境ノイズ NL

dB re 1μPa/√Hz

指向性 DI

受波アレイの指向性ゲイン

検出閾値 DT

周波数 f

kHz

高周波ほど吸収 α が急増（f² 依存）

水深 D

水温 T

°C

計算結果

—

音速 c (m/s)

—

吸収係数 α (dB/km)

—

許容伝搬損失 TL_max (dB)

—

最大探知距離 (m)

—

探知距離 (km)

—

波長 λ (m)

—

海中・ソナー伝搬の概念図

中央のソナーから同心球状に音波（青）が広がり、潜水艦（標的）に当たって反射エコー（橙）が戻ります。下端は海底、上端は海面。色の濃さは現在の伝搬損失を表します。

伝搬損失 TL vs 距離 R

吸収係数 α vs 周波数 f

理論・主要公式

$$SL - 2\,TL + TS - (NL - DI) \geq DT,\qquad TL = 20\log_{10}R + \frac{\alpha\,R}{1000}$$

アクティブ・ソナー方程式（往復伝搬）。SL=音源、TL=伝搬損失（球面拡散+吸収）、TS=標的反射、NL=環境ノイズ、DI=指向性、DT=検出閾値。R は m、α は dB/km。

$$SL - TL - (NL - DI) \geq DT$$

パッシブ・ソナー方程式（片道伝搬、TS なし）。標的が放射する音を直接聴取するモード。

$$\alpha \approx 0.106\,\frac{f^{2}}{f^{2}+1} + 0.52\,\frac{f^{2}}{f^{2}+4100} + 4.9\!\times\!10^{-4}\,f^{2}\;[\text{dB/km}]$$

吸収係数 α（Francois-Garrison 簡易式、f は kHz）。第1項ホウ酸緩和、第2項硫酸マグネシウム緩和、第3項純水粘性。

$$c \approx 1449.2 + 4.6\,T - 0.055\,T^{2} + 0.016\,D\;[\text{m/s}]$$

音速の温度・深度依存（UNESCO 簡易式、塩分 35‰）。T は °C、D は m。

水中音響ソナー方程式 — 探知距離・伝搬損失

🙋

「ソナー」って、潜水艦が「ピーン…ピーン…」って音を出して相手を探すあれですよね？あの探知距離って、どう決まってるんですか？

🎓

そう、映画でよく聞くやつだ。あの 1 回の「ピーン」で何 km 先まで見えるかを決めるのが「ソナー方程式」というやつでね、Urick が 1983 年に整理した式が今でも教科書の定番なんだ。アクティブ式だと SL − 2·TL + TS − (NL − DI) ≥ DT。要するに「自分が出した音 SL から、行って戻ってくる伝搬損失 2·TL を引き、相手の反射 TS を足し、ノイズ NL から指向性 DI を引いたものが、検出閾値 DT を超えていれば見える」というだけのシンプルな引き算なんだよ。

🙋

引き算だけ！？意外と単純なんですね。じゃあ「伝搬損失 TL」が距離をどう食ってるのか、それさえ分かれば探知距離は出るってことですか？

🎓

そう、本質はそこ。TL は「球面拡散 20·log10(R)」と「吸収 α·R/1000」の和で、最初は距離の対数で増えるけど、遠くなるほど吸収の直線項がじわじわ効いてくる。右下の TL-R グラフを見てごらん。近距離（数百 m）は球面拡散が支配、数 km 越えると吸収が支配に変わる「肘」みたいなカーブが見える。許容 TL_max の水平線とこのカーブが交わる距離が、まさに最大探知距離 R_max なんだ。

🙋

なるほど〜！じゃあ周波数を下げれば吸収が減って、めちゃくちゃ遠くまで届きそうですね。なんで全部のソナーが低周波じゃないんですか？

🎓

いい着眼点だね。確かに 1 kHz と 100 kHz では吸収が 400 倍くらい違う（α の右下のグラフ参照）。だから対潜哨戒ソナーは数 kHz〜十数 kHz の低周波を使い、収束帯まで含めて 30 km 級を狙う。一方で漁業のフィッシュファインダや海底地形図の Multibeam（FURUNO や Kongsberg EM124）は分解能を稼ぐために 100〜400 kHz を使い、見える範囲は数百 m に絞る。「遠く＝低周波／細かく＝高周波」というトレードオフで、用途ごとに帯域を選ぶんだ。

🙋

なるほど…じゃあ NL（環境ノイズ）や DI（指向性）はどう効いてくるんですか？SL を上げるだけじゃダメなんですか？

🎓

SL を上げるのは一番手っ取り早いけど、上限がある。海生哺乳類への影響規制（NOAA や IMO MEPC）で、安易に音圧を上げられない海域も増えてきた。そこで効くのが「DI を稼ぐ＝大きなアレイで指向性を絞る」「NL を下げる＝静かな環境やフィルタで雑音を切る」「DT を下げる＝信号処理（相関、整合フィルタ、ビームフォーミング）で SN を引き出す」の 3 つ。実艦のソナーシステムは、SL を上げずに DI と DT で 20 dB 以上稼いでいるんだ。左のスライダーで NL を 5 dB 下げてみて。同じ SL でも探知距離がぐっと伸びるのが分かるはずだよ。

🙋

最後に1つ。映画だと「サーモクライン（水温躍層）に隠れる」って話が出てきますが、これは式のどこに入ってますか？

🎓

鋭いね。実はこの単純な式（球面拡散+吸収）には入っていない。音速は水温・深度・塩分で変わって、海中には音速の谷（SOFAR チャネル、深度 1000 m 付近）や急変層がある。そこに音線が曲げられて陰影帯（Shadow Zone）が生じたり、逆に収束帯（CZ、約 30 km 周期）が生じたりする。この「音線追跡」を厳密に解くのが Bellhop や RAM などのモデルで、本ツールはあくまで第一近似。ただ、第一近似でも SL や NL を 5〜10 dB 変えたときの感度はかなり正確に出るので、システム設計の最初の当たりづけにはちょうどいいんだ。

よくある質問

Urick (1983) によりまとめられた水中音響の基本式で、アクティブソナーでは SL − 2·TL + TS − (NL − DI) ≥ DT、パッシブソナーでは SL − TL − (NL − DI) ≥ DT で表されます。SL は音源レベル（dB re 1μPa @ 1m）、TL は伝搬損失（球面拡散 20·log10(R) と吸収 α·R/1000 の和）、TS は標的反射強度、NL は環境ノイズ、DI は受波指向性、DT は検出閾値です。このツールは左辺を最大化する R を反復解で求めます。

海水中の吸収はホウ酸 B(OH)3 と硫酸マグネシウム MgSO4 の分子緩和、および純水の粘性の3項からなり、いずれも f² に比例する項を含みます。低周波では分子緩和が、100 kHz 以上では粘性吸収が支配的になります。Francois-Garrison (1982) の式が標準で、1 kHz で約 0.07 dB/km、10 kHz で約 1 dB/km、100 kHz では 30 dB/km を超えます。本ツールはこの式の主要3項を実装しています。

アクティブは自分で音を出して反射を聞くため伝搬損失が往復で 2·TL となり、許容 TL が約半分になります。一方パッシブは標的が出す音を片道で聞くだけで TS は登場しません。同じ SL を仮定すると一般にパッシブの方が距離が伸びますが、対象が静粛化された潜水艦などでは SL が低く成立しません。本ツールでは Active/Passive を切り替えると分母（SE_factor）が 2/1 で切り替わり、TL_max が再計算されます。

本ツールは球面拡散と吸収だけの単純なモデルです。実海域では音速分布により SOFAR チャネル（深度 1000m 付近の音速極小）や収束帯（約 30km 周期で再収束）が生じ、伝搬損失が大きく変動します。長距離（>30 km）の評価には Bellhop や RAM などの音線・PE モデルが必要で、本ツールはあくまで第一近似（システム設計の初期検討、感度の傾向把握）として使ってください。

実世界での応用

対潜哨戒・水上艦ソナー：米海軍 Arleigh Burke 級駆逐艦の AN/SQS-53C や P-8 Poseidon 哨戒機の投下式ソノブイ、海上自衛隊 P-1 哨戒機などは数 kHz の低周波アクティブと広帯域パッシブを組み合わせ、収束帯を利用して数十 km 級の探知を実現しています。本ツールの SL 220 dB・f=3〜10 kHz・DT 5〜10 dB といった設定が概略レンジに対応します。

海底地形・海洋調査：Kongsberg EM124（12 kHz）や EM2040（200〜400 kHz）といったマルチビーム測深機、Edgetech のサイドスキャンソナーは、周波数で分解能と探査範囲を切り替えています。本ツールの周波数スライダーで f を 100〜200 kHz に上げると吸収が急増し、探知距離が数百 m スケールに収束する様子が確認できます。

漁業・養殖：FURUNO・JFE アドバンテック・Simrad EK80 などの魚群探知機・科学魚探は 38 kHz / 70 kHz / 120 kHz / 200 kHz を切り替え、対象魚種の体長と表層・中層・深層の使い分けに利用します。TS（魚体反射）は 30cm の鯖で約 −37 dB、1m のマグロで約 −25 dB が代表値。

水中ロボット・パイプライン点検：Saab Sabertooth や Boston Engineering Bluefin など AUV/ROV は前方障害物回避（FLS）に 300〜600 kHz、海底地形把握に 100 kHz 帯、長距離通信に低周波音響モデムを使い分けます。海底パイプラインの腐食・座屈点検では、サイドスキャン＋サブボトムプロファイラの併用が主流です。Acoustic Stealth（潜水艦の SL 低減）や Convergence Zone 利用は、戦術設計の核心テーマです。

よくある誤解と注意点

まず最大の落とし穴は、「球面拡散+吸収だけで実海域の TL が決まる」と思い込むことです。本ツールの TL = 20·log10(R) + α·R/1000 は深海・等音速・無境界の理想モデルで、実際の海では音速分布による屈折で SOFAR チャネル・収束帯・陰影帯が生じ、距離 10〜30 km で TL が ±20 dB も変動することがあります。浅海域（<100m）ではさらに海面・海底反射による干渉や境界吸収が支配的になり、円柱拡散 10·log10(R) に切り替わることもあります。本ツールは「自由音場の第一近似」として、感度設計の入口に使ってください。実運用評価には Bellhop（音線）や RAM（PE 法）など専用ソルバが必要です。

次に、「dB の基準値を取り違える」こと。水中音響の SL や NL は大気中（dB re 20μPa）と異なり dB re 1μPa @ 1m が基準で、その差は約 26 dB あります。文献で大気中の dB をそのまま水中の式に入れると、見かけの SN が大きくずれます。また NL の単位は「dB re 1μPa/√Hz」のスペクトル密度で、バンド幅 B [Hz] に対しては +10·log10(B) を足してバンド総ノイズに直す必要があります。本ツールは式の中で B=1000 Hz を仮定していますが、実機では帯域に応じて再計算してください。

最後に、「TS や DT は固定値である」という誤解。標的反射 TS は方位（横腹は正反射、艦首は−20 dB 級）、周波数、コーティング（吸音材）で大きく変動します。検出閾値 DT もパルス長 τ、相関時間、誤警報率 Pfa、検出確率 Pd によって変わる「設計値」で、典型的には 0〜15 dB の幅を持ちます。さらに海生哺乳類保護（NOAA・IMO MEPC 1/Circ.833）では特定周波数帯の SL に上限規制があり、設計時に必ず参照すべきです。本ツールの数値は教育・初期検討用で、実機設計には ONR や Navy の標準手順（Wagner & Mylander, RAYMODE 等）の併用を推奨します。

水中音響ソナー方程式探知距離シミュレーター

水中音響ソナー方程式 — 探知距離・伝搬損失

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点

水中音響ソナー方程式 探知距離シミュレーター

水中音響 ソナー方程式 — 探知距離・伝搬損失

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点

水中音響ソナー方程式探知距離シミュレーター

水中音響ソナー方程式 — 探知距離・伝搬損失