Sonar For Practising Engineers PDF
Sonar For Practising Engineers PDF
Sonar For Practising Engineers PDF
Third Edition
This page intentionally left blank
Sonar for Practising
bngineers
Third Edition
A. D. Waite
8
JOHN VVILEY & SONS, LTD
Copyright 02002 by Ashley Waite Published by John Wiley & Sons Ltd,
Baflhs Lane, Chichester,
West Sussex PO19 lUD, England
National 0 1243 779777
International (+44) 1243 779777
e-mail (for orders and customer service enquiries): [email protected]
Visit our Home Page on http://www.wiley.co.uk
This work is based on Sonar for Practising Engineers, Second Edition, published and distributed by Thales
Underwater Systems Limited (formerly named Thomson Marconi Sonar Limited), Ocean House,
Templecombe, Somerset, BA8 ODH (www.tms-sonar.com), 1998
All Rights Reserved No part of this publication may be reproduced, stored in a retrieval system, or
transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, scanning or
otherwise, except under the terms of the Copyright, Designs and Patents Act 1988 or under the terms of a
licence issued by the Copyright Licensing Agency Ltd, 90 Tottenham Court Road, London, UK W 1P OLE
without the permission in writing of the Publisher with the exception of any material supplied specifically
for the purpose of being entered and executed on a computer system for exclusive use by the purchaser of
the publication.
Neither the authors nor John Wiley & Sons, Ltd accept any responsibility or liability for loss or damage
occasioned to any person or property through using the material, instructions, methods or ideas contained
herein, or acting or refraining from acting as a result of such use. The authors and publisher expressly
disclaim all implied warranties, including merchantability or fitness for any particular purpose. There will
be no duty on the authors or publisher to correct any errors or defects in the software.
Designations used by companies to distinguish their products are often claimed as trademarks. In all
instances where John Wiley & Sons, Ltd is aware of a claim, the product names appear in capital or all
capital letters. Readers, however, should contact the appropriate companies for more complete information
regarding trademarks and registration
Ashley Waite has asserted his right under the Copyright, Designs and Patents Act 1988 to be identified as
the author of this work.
A catalogue record for this book is available from the Brisish Library
Preface xiii
About the Author xv
Introduction xvii
1 Sound 1
1.1 Wave motion 1
1.2 Sound pressure 1
1.3 Reference intensity 3
1.4 Source level 4
1.5 Radiated power 5
1.6 Limitations to sonar power 5
1.7 Cavitation 6
1.8 Interaction 8
1.9 Changes to arrays 8
1.10 Projector sensitivity 9
I . 11 Hydrophone sensitivity 9
I . I2 Spectrum level 10
1.13 Sound in air and in sea water 11
1.14 Problems 12
2 Arrays 13
2.1 Need for projector arrays 13
2.2 Need for hydrophone arrays 13
2.3 Beam patterns 14
2.4 Directivity of a dipole 15
2.5 The general line array 18
2.6 Line array: beam pattern vs. steer angle 20
2.7 Broadside array: length and spacing 20
2.8 Beam pattern for a continuous line 21
2.9 Shading 24
2.10 Shaded arrays: transmit source levels 27
2.1 1 Directivity index 28
... Contents
Vlll
4 Target Strength 67
4.1 Definition 67
4.2 Formulae 67
4.3 Measurement 68
4.4 Dependence on pulse type and duration 69
4.5 TS of a sphere 69
4.6 TS of some simple shapes 70
4.7 TS of small targets 72
4.8 Mine target strength 72
4.9 Torpedo target strength 73
4.10 Submarine echoes 74
4.1 1 Beam aspect target strength 74
4.12 Bow aspect target strength 75
4.13 Submarine target strengths 75
Contents iX
6 Reverberation 103
6.1 Sources of reverberation 103
6.2 Scattering and reflection 103
6.3 Boundary roughness 104
6.4 Classes of reverberation 105
6.5 Backscattering strength 106
6.6 Reverberation target strength 106
6.7 Volume reverberation 107
6.8 Boundary reverberation 110
6.9 Scattering layers 111
6.10 Volume scattering strength 112
6.1 1 Sea surface scattering strength 113
6.12 Bottom scattering strength 114
6.13 Variation with fi-equency 116
6.14 Reverberation under ice 117
6.15 Problems 11:
9 Activesonar 161
9.1 Range, pings and doppler shift 161
9.2 Pulse types 163
9.3 CW processing 163
9.4 FM processing 164
9.5 Active sonar equations 165
9.6 Reverberation index 168
9.7 FM pulses 169
9.8 CW pulses 171
9.9 Reverberation rejection by CW pulses 174
9.10 Reverberation and Target Echoes in the main lobe 175
9.1 1 Reverberation and Target Echoes in the sidelobes 178
Contents xi
Conclusion 277
Solutions to Problems 281
Index 293
Preface
Most books on sonar - the use of underwater sound for the detection, classifica-
tion and location of underwater targets and for communications and telemetry -
have been written by physicists and mathematicians. They are not always easily
understood, nor are they immediately useful for solving the problems met by
engineers and technicians.
The aims of this book - written by a practising engineer for practising
engineers - are to provide an understanding of the basic principles of sonar and to
develop formulae and rules of thumb for sonar design and performance analysis.
No prior knowledge of sonar is assumed, and the physical principles and mathe-
matics will be readily understood by engineers and technicians.
The earlier editions were produced to be supplied as back-up material to a short
sonar course given by the author. This edition has been extensively rewritten to
facilitate its use by an individual reader. Several new topics have been included:
Many chapters contain worked examples and most chapters conclude with a few
problems for the reader to solve; solutions are given at the end of the book. I hope
these problems will be particularly useful to lecturers and students.
The book can be divided into three parts:
Practical sonar systems: this part develops passive and active sonar systems; it
looks at passive broadband, narrowband, intercept and communications sonars;
and it considers active sonars for the detection of submarines, mines and
torpedoes.
About the Author
Ashley Waite retired from the UK Ministry of Defence in 1990 after 40 years of
Underwater Warfare Research at Portland.
Primarily engaged in Sonar Research and Development, his experience spans
Surface Ship Sonars, Submarine Sonars, Surface Ship Torpedo Defence and a
lesser involvement with Helicopter and Minehunting Sonars. He has also worked
on Submarine Command Systems and the Underwater Aspects of Surface Ship
Command Systems.
xvi About the Author
Many methods of detecting the presence of underwater targets in the sea have
been investigated. Here are some non-acoustic methods which have had varying
degrees of success:
Magnetic
Optical signatures
Electric field signatures
Thermal detection (infrared)
Hydrodynamic changes (pressure)
Basic passive systems, however, give no information about the range of a target;
a signal may belong to a close, quiet target or a noisy, distant target. More complex
passive systems estimate range by the following methods:
All of these methods are fundamentally dependent on the accuracy of the bearing
measurements and therefore demand large arrays and large separations to achieve
useful range estimates.
Active sonar uses a projector (an underwater loudspeaker) to generate a pulse
of sound which travels through the water to a target and is returned as an echo to a
hydrophone, often the same device as the projector and in this context more
commonly known as a transducer. The echo now has to be detected against a
background of noise and reverberation (unwanted echoes from the sea surface
and sea bed and from scatterers within the volume of the sea). Because the time
between transmission of a pulse and reception of an echo can be measured and the
speed of sound in the sea is known, the range of the echoing target is simply
calculated. Active sonars are sometimes known as echo ranging systems.
To survive in sonar, newcomers must become familiar with the decibel. They
will hear old hands discussing the design and performance of sonar systems in
exchanges where almost every other word seems to be deebee (dB or decibel).
What is this decibel?
First of all, the be1 is inconveniently large so it has been divided by 10 to
become the decibel. It simply compares the power or intensity of the sound at one
point in a system with that at another. The decibel defined:
many parameters all expressed in decibel form. And to avoid serious errors, the
decibel must be carefully defined and correctly applied to all of the terms in the
equations.
Power ratios in sonar systems are frequently very large numbers. Calculations
are greatly simplified when very large numbers are expressed in logarithmic form,
so that values can be added instead of multiplied. Interestingly, although radar
engineers may use decibels to describe the terms in a radar system (e.g. antenna
gain, receiver noise factor), the radar equation is more commonly written in linear
form rather than logarithmic form:
The individual terms may be given their linear values and multiplied, or expressed
in decibels and added.
Compare with an active sonar equation:
All the examples, here and in the rest of the book, will use credible values for
these parameters, and the reader should thereby learn to recognize and question
suspiciously large or small values. (None of the terms in the sonar equations is
likely to have a value outside the range - 100 dB to +250 dB, at the very most.)
The speed of sound in water is about 1500 m/s (3355 mph), much faster than its
speed in air of 340 m/s (760 mph) at sea level, but very much slower than the
speed of light (electromagnetic radiation). Velocity, frequency and wavelength are
related by
xx Introduction
where
The radar will detect the submarines periscope and, assuming that the electro-
magnetic energy spreads spherically, the inverse square law will apply and the EM
propagation loss (two-way) is 40 log 10 000 = 160 dB.
The sonar will detect the submarines hull and, again assuming spherical
spreading, the sound propagation loss is also 160 dB.
Therefore, given echoing surfaces above and below water, radar and sonar
systems have comparable performances. But if the submarine is completely
submerged, the EM losses over the final, say, 100 m are completely prohibitive at
any radar fiequency (Figure 0.1).
EM losses in the highly conductive sea water are given by 1400f1/* dB/km,
where f is in kilohertz and, for just a 100 m path in the sea, the losses at
2000 MHz are 200000 dB! Even at 30 kHz the losses equal 770 dB; the sea
effectively presents a short circuit to the EM energy.
A modern nuclear-powered submarine can remain completely submerged for
indefinite periods and during a mission may never expose any reflecting surface
above water, but will rely completely upon its sophisticated sonars to navigate, to
build up a complete acoustic scenario of its surroundings, to deploy counter-
measures, and finally to prosecute an attack. Radar systems, therefore, will never
be given opportunities to detect the submarine.
xxii Introduction
i_ I ) ) ) ;
Radiated noise at source = 120 dE3 c&
Radiated noise at hydrophone = 50 db
Background noise = 40 db
PL = 120 - 50 = 70 dB
20 log R = 70 dB
therefore R = 3000 m.
Practical passive sonars will have directional arrays and limited operating
bandwidths to reduce background noise, and they will integrate the signal (the
radiated noise from a target) over a period of time. These measures will all
improve the allowable PL significantly and greater ranges are possible. (But
120 dB represents quite a noisy submarine, so perhaps 3000 m is not too
unrealistic a range.)
Introduction xxiii
Transmittedpulse = 200 d
The simplest active sonar will have an omnidirectional projector and hydro-
phone (Figure 0.3). Its detection performance is given by the allowable two-way
propagation loss (2PL), i.e. from projector to target and back to the hydrophone. If
now the echo at the projector is 10 dB greater than the background noise, which
for a surface ship sonar might be 60 dB, say, then
therefore R = 3000 m.
Practical active sonars will have directional transmit and receive arrays, and a
knowledge of the pulse will allow the processing to be matched to the returned
echo. These measures will all improve the allowable PL significantly, and much
greater ranges will be achievable.
This page intentionally left blank
1
Sound
The speed of sound is not to be confused with the particle velocity, u, which refers
to the movement of the molecules in the medium.
2 Sound
-
Particle motion
The sound wave carries mechanical energy with it in the form of the kinetic
energy of the particles and the potential energy of the stresses in the medium.
Because the wave is propagating, a certain amount of energy per second will flow
across unit area normal to the direction of propagation.
This energy per second (power) crossing unit area is known as the intensity of
the wave (power per unit area). For a plane wave, the intensity is related to the
pressure by
1 I, = 0.67 X
1
W/m2 J 1
Intensities are often loosely stated as re 1 p a . This is clearly incorrect since the
micropascal is a unit of pressure not intensity (power per unit area). Strictly,
intensities should be stated as re the intensity due to a pressure of 1 @a.
4 Sound
SL = lOlog );(
where
propagation loss - may only increase the range by a small percentage. Attempts to
increase the radiated power are eventually limited by two effects: cavitation and
interaction.
6 Sound
1.7 Cavitation
When the power applied to a projector or an array is increased, bubbles form on the
surface and there is a resultant loss in power by absorption and scattering within the
bubbles, degradation of the beam pattern and a reduction in the acoustic impedance
into which the projector generates (resulting in a mismatch with the transmitters
supplying the projector array). Cavitation is a function of depth (pressure) and can
be avoided by not exceeding the cavitation threshold (Figure 1.2).
Figure 1.2 gives a conservative estimate of the cavitation threshold for fiequen-
cies up to about 10 kHz and for pulses of at least 5 ms duration. It predicts the
start of cavitation at some point before it spreads to the total radiating surface. In
practice, higher power intensities are possible, perhaps 3 or 4 times higher, but
these should be confirmed by experiment before finalizing a design.
Because cavitation takes a finite time to build up, albeit very short, the threshold
also changes with pulse length and fiequency:
Pulse length: the threshold is increased for very short pulses, up to about 3
times for a 0.5 ms pulse compared with its 10 ms value, but is sensibly constant
for longer durations.
Example 1.1
What is the maximum radiated power for a 100 ms pulse of frequency 10 kHz to avoid
cavitation, for a cylindrical array of 1 m diameter and 1 m height, containing circular
elements closely packed and operating at a depth of 5 m?
The radiated power intensity may be increased 3 times for the short pulse and 5 times
for the higher frequency. Therefore
Maximum radiated power = 0.1 X 2 X 3 X 5 = 3 kW
The radiated power intensity may be increased 3 times for the short pulse and 20 times
for the higher frequency. Therefore
Note, however, that the power intensity is 3000 kW/m2 and for even moderate losses,
say 20 per cent, the array must dissipate some 600 kW of heat per square metre or, to
lapse into imperial, 6 kW over a 4 in square! Clearly the limiting factor at higher
frequencies andor greater depths is not cavitation but the power-handling capability
of the array.
8 Sound
1.8 Interaction
When a number of projectors are assembled together in an array and driven
electrically, the velocity of motion of the individual elements is not constant, but
varies from element to element in a complex manner due to the acoustic inter-
actions between them.
Unless this is compensated in the design, the interactive effects will reduce the
total power output; affect the transmitted beam pattern; and produce a mismatch
with the transmitters, which could damage transmitters, projectors or both.
Interactive effects may be reduced in three ways:
Separating the elements of the away produces a larger than optimum array and,
particularly if the spacing is much greater than A/2, deterioration in the transmit
beam pattern. This may be acceptable if the resultant transmitter design problem
is eased and total costs reduced.
Making the individual elements large so that their self-radiation impedances are
much greater than the mutual radiation impedances between elements.
Using individual amplifiers to drive each element at the correct amplitude and
phase to yield a uniform velocity of motion across the array.
sv = lOlog ;( );
= SL - 20logv
sw = S L - 10logP
Example 1.4
If SL = 200 dB re 1 pPa (strictly, re the intensity of a plane wave of pressure 1 @a)
and P = 30 W, then SW= 200 - 15 = 185 dB/W.
Example 1.5
If 20 log p = 80 dB re 1 Pa and v = 1 pV, then S h = - 120 - 80 = -200 dB/V. To
find the output voltage from such a hydrophone, we use
BL = SpL+ 10logAf
Example 1.6
If the spectrum level is 40 dE3 and Af = 3000 I-kq then BL = 40 + 35 = 75 dE3.
If the spectrum level is not flat, the band level can be obtained by integrating the
intensity over the complete band. For ambient sea noise, the intensity may fall by 6 dB
per octave. The total intensity is then given by
Example 1.7
What is the sound level (total intensity), BL, of sea noise in a band from 100 Hz to
10 kHz, given that the spectrum level at 100 Hz is 100 dB?
and in dB form, BL = 100 + 20 = 120 dB. This is 20 dB less than if the spectrum
were flat.
1.14 Problems
1.1 The pressure, p , of an underwater sound is 100 pPa. What is the intensity, I, of
the sound and its level expressed in decibels compared with the reference intensity, I,?
1.2 What is the source level, SL, of a projector radiating 40 kW of acoustic power and
with a directivity index, DIt, of 15 dB?
1.3 What is the safe maximum radiated power to avoid cavitation, given a 100 ms
pulse of frequency 20 kHz and a planar array of 2 m length by 1 m height, containing
circular elements closely packed and operating at a depth of 10 m?
1.4 The total intensity of the sound in an octave from 2000 to 4000 Hz is 80 dB
relative to the intensity resulting from a sound pressure of 1 pPa. Assuming the band
is flat, what is the spectrum level, SpL, of the sound?
2
Arrays
directivity of a hydrophone array is due to all signals impinging upon the array
from the same direction being in phase, and therefore reinforcing one another.
S2H = dsin0
2nd
cp = -sin0
I.
The sound pressure at the hydrophone is the vector sum of the pressures p1 and
p2, of the same amplitude and with a phase difference cp (Figure 2.2):
For a single source S1, the pressure is constant whatever the direction. When a
second source Sz is present, the pressure varies between 0 and 2 p as a function of 0:
cp
= 2pocos- = 2pocos
2
16 Arrays
PI +P2
/
/
/
/
/
/
/
/
P2
2po is the pressure, pa, on the axis of symmetry of the dipole (where 8 = 0) and
we can write
- = cos
Pa
G sin e)
The 3 dB beamwidth (the width in degrees at the half-power points of the dipole
response) for different spacing of the dipole elements is obtained by putting
[
2010g cos -sin8
(7 11 = 3 dB
Cd ) 7
cos -sin8 =-
nd . n
--slnO =-
;1 4
Directivity ofa Dipole 17
Table 2.1 gives the 3 dB beamwidth at four element spacings and Figure 2.3 shows
the beam patterns. When the spacing is exactly L/2, nulls appear at 90" and 270".
As the spacing is increased, secondaries (sidelobes) appear which are soon
unacceptable, and when d = L they totally change the beam pattern so that the
major responses now occur at the previous null angles.
U-
270" 800
b(8) = 20 log
{sin(nnd/A)sin 0)
{ n sin(xd/A)sin e)1
The main lobes of an array can be steered by introducing phase or time delays in
series with the elements. When this is done the beamwidth and sidelobe structure
of the array is changed. The beam pattern of the line array is then modified to
b(0) = 20log
sin nn(d sin O / L - d sin OJL)
n sin x(d sin 8/A - d sin &/A) 1
where 6, is the angle of steer at broadside.
This function is used to plot Figure 2.4, the 3 dB beamwidth of a line array at
various steer angles. In endfire (OS = 90") the beam pattern is searchlight shaped
and there is no lewright ambiguity. In a transition region close to 90" the 3 dB
beamwidth is ambiguous. Although the beamwidth varies as the beam is steered,
the directivity index of the array stays approximately constant. An approximate
formula for the 3 dB beamwidth of a line array, useful for steers up to about 60", is
where
L = length (m)
f = frequency (kHz)
0, = steer angle (degree)
I I:
2e3 = -
The General Line Array 19
3 dB beamwidth (degrees)
5 6 7 8 910 20 30 40 50
Array length (wavelengths)
Figure 2.4 Beamwidth of a line array of A/2 spaced elements at various steer angles
20 Arrays
OdB
- 20
OdB OdB
Figure 2.5 Line array, 10 elements spaced A/2: beam patterns at various steers
Figure 2.6 Effects of length and element spacing on a broadside array: (a) element spacing
constant at A / 2 , array length increasing; (b) array length constant at 3L, element spacing
increasing
pe
_.
-
-
sin{ (nL/il)sin S}
pa (nL/il)sin 8
b(0) = 20 log
sin{ (nL/iZ)sin 0)
(nL,/l)sin 8 1
22 Arrays
0
dB
- 10
-20
-30
(nLA) sin 0
There are nulls in the beam pattern where x = f n , f2n, f3n, . . ., and peaks in
the beam pattern where x = 0, fsn, f;n, . . ..
The sidelobe levels are given by
where rn is the sidelobe number. The results apply quite well to an array of spaced
elements provided their spacing, d, is not more than about A/2 and the array is at
least 3A in length.
Table 2.2
X x-' sin x 20 log(x-' sin x)
0 1 0
1.39 0.7 1 -3
2.27 0.32 - 10
2.85 0.1 -20
The numerical values in Table 2.2 can be used to derive some useful approxima-
tions for the shape of the main lobe. Beginning with the 3 dB beamwidth (Hz):
Beam Pattern-for a Continuous Line 23
nL
1.39 = -sni 03
A
Converting to radians, for small angles we obtain
1.39 = n-e3
L )'(
1
263 = 100/n
Example 2.1
What are the 3 dB, 10 dB and 20 dB beamwidths for a line array of length 6;1? What
are the bearings of the first-, second- and third-order sidelobes relative to the direction
of maximum intensity?
2.9 Shading
Except for the simple dipole with element spacings less than or equal to 1/2, all
the above beam patterns have significant sidelobe levels. The first-, second- and
third-order sidelobes are respectively only 13.5, 18 and 21 dB below the peak level
of the main lobe. Strong signals will be detected through the sidelobes of adjacent
beams as well as, correctly, in the main lobe of the beam at the bearing of the
signal. The resultant bearing ambiguity and additional, false, signals complicate
all further processes. The aim, therefore, should be to produce the narrowest
possible main lobe consistent with some reasonable level of sidelobes.
By using amplitude shading or weighting, i.e., varying the amplitudes of the
signals applied across the array (in transmit) or coming fiom the array (in receive),
it is possible to reduce the sidelobes, but always at the expense of widening the
main lobe and reducing the DI somewhat. In transmit, reducing the DI means that
the source level (SL) is reduced as a consequence. In an active system it is the
combined, transmit and receive, beam pattern which is important and it is hardly
ever desirable to hrther reduce the sidelobes at the expense of SL - far better to
adequately shade the receiver elements only.
All shading functions widen the main lobe by some factor; the factor is 1.36 for
the cosine shading hnction in Figure 2.8. Note also that the skirts of the main lobe
be (dB)
21 dB
.-
It 2n 3n
(nLA) sin 8
extend into the first sidelobes of the unshaded pattern, and therefore simply
quoting the level of the first sidelobe is an incomplete description and may be
misleading (Section 12.6).
Mathematical procedures, based on work by Dolph and Chebyshev, determine
the amplitude shading coefficients which yield the narrowest main lobe for a
specified level ofsidelobes (Figure 2.9). For example, the coefficients for a six-
element line array to yield -30 dB sidelobes are 0.30, 0.69, 1, 1, 0.69, 0.30. And
the main lobe is broadened, at the 3 dB points by a factor of 1.3.
-10 -
-20 -
-30 -
( J c L A )sin 8
Example 2.2
If the array of the previous example is (a) Dolph-Chebyshev 30 dB shaded and (b)
Hamming shaded, what are the 3 dB and 20 dB beamwidths?
From the previous example, the unshaded beamwidths are 8.3" and 17.3". Therefore
the shaded beamwidths are as follows:
Example 2.3
A Dolph-Chebyshev array designed to achieve -30 dB sidelobes will have its SL
reduced due to the reduction in transmit power, 3 dB, and due to the broadening factor,
1 dB, Le., a total reduction of 4 dB.
Check
Given a line array of six elements, spaced A/2 and each capable of a maximum
output of 1000 W omnidirectionally, we obtain
We will first derive the directivity index for a simple dipole and, avoiding complex
mathematical treatment, extend this result to cover the directivity indices for more
complex arrays.
DI of a Simple Dipole 29
If I is the average intensity at the surface of the sphere, the total sound flux
passing through its complete surface is
and therefore
and substituting
and
We would expect both of these results intuitively: two elements are twice (3 dB) as
good as one; and when the two become one (d = 0), the dipole is omnidirec-
tional and DI = 0 dB.
DI qf a Line Army 31
IDI= 10lognl
Clearly the elements in the line are only spaced A/2 at one frequency. Therefore,
for a wideband system the DI will change with frequency. As the frequency is
reduced, the line approximates to a continuous line of length L = n&/2, where A0
is the wavelength at fo, the design frequency for the line array and
Note that, from its definition, DI is always positive and cannot fall below 0 dB.
As the frequency is increased, the line can no longer be considered as
approximating to a continuous line, and excessive sidelobes appear as the spacing
approaches A. Therefore this equation should not be used for frequencies above,
say, 1.5f0.
Example 2.4
A line consists of 64 elements spaced A/2 at 1 kHz. What is its DI at 1 kHz, 100 Hz
and 1.5 kHz?
Note again that DI is always positive and cannot fall below 0 dB. As the fiequency
is increased, the planar array can no longer be considered as approximating to a
continuous surface, and excessive sidelobes appear as the spacing approaches A.
Therefore this equation, in common with the equation for a line array, should not
be used for frequencies above, say, 1.5fo.
Example 2.5
A flank array has 8 rows of 32 elements spaced A/2 at 1 kHz. What are its DI at
800 Hz and its 3 dB beamwidths at 1 kHz and at steer angles of (a) broadside, (b) 30,
(c) 60?
The length of the array is 161 and its height is 4L The horizontal and vertical
beamwidths are given by the same formula:
DI of a Planar Array 33
+
8 h = (76/24)( 1 900/4000) = 3.9"
(b) 0, = 12.7" (steer is only horizontal)
Element
Stave
Ring
where
h = height (m)
d = diameter (m)
fo = design frequency (kHz)
This is for a 'full' beam, i.e., a beam formed by using all the usefully
contributing elements on the surface of the cylinder (an aperture of about 120" of
the cylinder). The design frequency, fo, is the fkequency where the elements are
spaced A/2.
Example 2.6
A cylindrical array has height 1 m, diameter 2 m and design frequency 5 H z , so
DI = 1010g(5 X 1 X 2 X 25) = 24 dB.
If 'half-beams' are formed for cross-correlation or phase comparisons (see later), the
DI of each half-beam is 3 dB less than the DI of the corresponding full beam. Here are
some practical formulae for the beamwidths of a cylinder.
Horizontal beamwidth
& = - 88 = 8.8"
Vertical beamwidth
This is the same as the formula for line arrays, and we obtain
76
0 -- = 15.2"
"-1X5
36 Arrays
DI = 101og5hdf20
= 18.5 dB
Helicopter dipping sonars often use volumetric arrays. The array is mounted on
ribs which fold (umbrella fashion) for stowage and lowering &om the helicopter
into the sea. When the array reaches the required depth, the rib fastenings are
released and the array configures as a number of concentric cylinders. A separate
transmitting array of omnidirectional projectors (omnidirectional in the horizontal
plane) may complete the assembly.
38 Arrays
I.75 m
t
I
Vectors add in quadrature
I
+----
(0)
-+-
A Y
\
A ------ 4
Vectors cancel
/ Vectors add
2.20 Beamformers
Analogue
Analogue beamformers are unlikely to be used except for very simple arrays in
modern systems. They are assembled using lumped constant (LC) networks which
must have a linear phase relationship with frequency. They are difficult to modify
and cannot, at least not simply, be made to respond to sound speed changes.
(When the speed of sound changes, because its frequency cannot change, its
wavelength must change. Therefore the parameters of the beam will also change
unless the speed of sound is used in the beamforming process.)
Digital
Digital beamformers use shift registers or random access memory store (RAM) to
hold sampled data. The beams are then formed by addressing the samples. Shading
functions and sound speed variations are easily accommodated.
Narrowband systems
If the beamformer is for use with a narrowband system then phase delays, which
are only exactly correct at one frequency, may be used to form the required beams
instead of the more usual time delays which are constant with frequency and hence
suitable for broadband systems.
40 Arrays
Cylindrical array
\ \ \
Plane wave incident on dome
The elongated dome meets the hydrodynamic requirements. The array, however,
can only be close to the dome over about half its circumference. Baffle 1 reduces
the stern noise. Performance astern (known in sonar as the non-operational
bearings) will inevitably be impaired and this determines the position of baffle 1.
It cannot be too close to the array without affecting too many non-astern beams.
Problems 41
The resultant gap between the array and baffle 1 introduces unwanted reflections.
Therefore baffle 2 is introduced to eliminate the more significant reflections.
Reflection A, in a region where the array is close to the dome, is harmless but
reflection B (without the baMle) would contribute energy to several incorrect
beams.
2.22 Problems
2.1 What are the beamwidths to the half-power points of a line array of length 10 m
and frequency 4 kHz (i) normal (broadside) to the array axis and (ii) at a steer of 60?
2.2 A line array has a length of 6L. What is the level of the fourth-order sidelobes (in
dB relative to the main lobe) and where do their peaks occur in the beam pattern?
2.3 What is the DI of a baffled planar array of length 5 m and height 2 m at 8 kHz?
The elements are spaced A / 2 at 10 kHz. First calculate the DI using the numbers of
elements and then compare with the result using the dimensions and wavelength.
/ P L = lolog($) (dE3)l
3.2 Losses
A first approach to quantifling the likely PL is to consider it as the sum of a
spreading loss and a loss due to absorption. Other losses - not functions of range
- include scattering and refraction, and these will be considered later on.
44 Propagation of Sound in the Sea
Refer to Figure 3.1(a). When the source is located in an unbounded and lossless
medium, the power is radiated equally in all directions and the total power, P,
crossing spheres of increasing radii surrounding the source does not change with
range. Therefore, since power = intensity X area,
PL = lolog(;) = 1010gr2
or in logarithmic form,
where Y is in metres.
Refer to Figure 3,l(b). When the source is bounded by parallel planes separated
by h metres, the power - again unchanging - crossing cylindrical surfaces of
increasing radii surrounding the source, is given by
PL = lolog(;) = 10logr
or in logarithmic form,
I P L = 10logrI
where Y is in metres.
46 Propagation of Sound in the Sea
Viscosity: losses due to viscosity are present in fresh water and salt water. This
contribution is proportional to the square of the frequency and accounts for the
straight line for fresh water when a is plotted on a log-log graph like Figure 3.2.
Molecular relaxation: losses due to molecular relaxation are only present in salt
water. The mechanism is a reduction of molecules to ions induced by the
pressure of the sound. At very high frequencies (greater than about 500 kHz)
pressure changes are too rapid for the relaxation to take place and therefore no
energy is absorbed. Magnesium sulphate relaxation is dominant over the
frequency range 2 to 500 kHz.Below 2 kHz boric acid relaxation contributes to
the losses.
Extensive measurements of these losses have been made and several empirical
formulae exist which relate them to frequency, depth (pressure) and salinity. The
total loss is given as an attenuation coefficient, a, in dB/ km.
The attenuation coe$cient, a, increases rapidly with frequency and changes
with temperature. It also varies with depth and salinity, but less strongly. Figure
3.2 plots the variation of a with frequency at a temperature of 10C and a salinity
of 35 parts per thousand (ppt), based on the formulae of Francois and Harrison.
An approximation for a, usehl between 0.5 kHz and 100 HZ,in 'standard' sea
water is given by
Table 3.1 is nut based on the approximation; it offers a selection of values for a
that will be useful in performance comparisons and in assessing the effects of
frequency changes on a design. Linear interpolation will be sufficiently accurate
to determine cr. at other frequencies.
0.5 MIZ 1 kHz 2 kHz 5 kHz 10 Mz 20 kHz 50kHz lOOkHz 200 kHz 500 kHz
5 0.02 0.06 0.14 0.33 1.00 3.80 15 30 55 120
10 0.02 0.06 0.14 0.29 0.82 3.30 16 35 60 125
15 0.02 0.06 0.14 0.26 0.68 2.80 17 40 65 130
48 Propagation of Sound in the Sea
lPL=20logr+arx (dB)]
PL = 20 log + ar x 10-3
a fWz)
a 0.06 1
b 0.2 3
c 1.0 10
d 7 30
e 33 100
f 64 300
t = temperature ("C)
s = salinity @pt)
h = depth (m)
The speed of sound at 1O"C, at zero depth, and for a salinity of 35 ppt is 149Om/s.
Here are some approximate coefficients for sound speed valid for use with this
'standard' sound speed:
The result of these opposing tendencies is to produce SSPs which vary widely
within the first few hundred metres of depth, and these SSPs are further
complicated by diurnal changes as well as mixing of the surface layer by wind and
waves.
Figure 3.4 shows a typical deep sea SSP divided into four principal layers:
Surface layer (duct): a layer of isothermal water mixed by the action of wind on
the surface of the sea. Sound tends to be trapped in this layer by surface
reflections and upward refractions.
Main thermocline: little affected by the seasons. This is where the main increase
in temperature over the cold depths of the sea occurs. Although the pressure
-
52 Propagation of Sound in the Sea
Speed of sound ( d s )
Surface layer (duct)
100 m
Seasonal thermocline
-------
Main thermocline
1OOOm
2000 m
Deep isothermal layer
5
e
P
1
Figure 3.4 Typical deep sea SSP
increases with depth, the net effect of temperature and pressure changes is to
reduce the sound speed through this layer.
20 40 60 km
Depth
1500 m - 10" vertical beamwidth
Note that, particularly at high latitudes, the DSC may be close to the RAP
critical depth. A capability to steer the sonar beam vertically by, say, 5" (assuming
a 10" vertical beamwidth) as well as a capability to deploy the sonar at variable
depths will assist in reducing the shadow zones that would otherwise exist.
The great depths at which it is necessary to deploy both source and receiver - at
54 Propagation of Sound in the Sea
0 0 0 0 0 0
\ I
1000 1000 lo00
Figure 3.7 Critical depths for RAP: (a) tropical, latitude 20"; (b) temperate, latitude 50"; (c)
Mediterranean
least 500 m for DSC mode and 1000 m for RAP mode - impose severe engineer-
ing and mobility constraints on their use. The design of elements capable of
operating at depths in excess of, say, 1000 m is difficult if not impossible and
losses in the long cables will limit the achievable source levels.
below the layer, are therefore difficult to detect. As with all propagation in the sea,
no mode is perfectly described by a simplified ray trace. The shadow zone is an
area where the sound intensity is greatly reduced and the transition from the
surface duct is not abrupt.
Increasing the depth of the source so that it is below the layer (Figure 3.9) has
the effect of increasing the range of the start o f the shadow zone, but it may then
extend into the duct. As with DSC and RAP modes, the ability to manoeuvre in
depth - a variable depth sonar (VDS) deployed from a moving surface ship or a
hovering helicopter - offers significant operational advantages.
Rays projected at steeper angles of depression, 5" or more, are bent downwards at
first - producing a shadow zone. At greater depths the pressure bends these rays
upwards to form annuli of high intensity (Figure 3.10); each annulus is known as a
convergence zone (CZ). The water must be deep enough for upwards refraction to
prevent the rays from hitting the bottom. Typically the water depth must be in
excess of 3000 metres. Depending on the bottom depth, the first convergence zone
will occur at around 30-50 km and it will be 3-5 km wide.
The propagation is sphericd spreading (there are no boundaries to confine the
sound) and the absorption losses are those appropriate to a temperature of about
4C (the temperature over most of the path) and the focusing effect gives a
convergence gain of typically 3-6 dB.
Successive convergence zones exist at multiples of the original range. The
propagation losses (two-way) will prohibit the use of any but the first convergence
zone for active systems, but passive systems can make detections at the range of
the second or even third convergence zone.
Bottom Bounce Propagation 57
An active bottom bounce (BB) system will be large, very demanding in ships
power supplies and transmitters and requiring a large site for the installation of its
arrays. Successive range annuli exist at multiples of the first annulus. The
propagation losses (two-way) will prohibit the use of any but the first for active
systems, but passive systems can make detections at the ranges of the second
annulus.The range annuli for BB mode are significantly wider than for the focused
CZ mode. There must always be a shadow zone out to some minimum range, at
least 10 km, dependent on depth and the allowable angle of depression of the
sonar beam (allowable from the viewpoint of bottom loss).
58 Propagation of Sound in the Sea
Ray theory: this postulates wavefronts and the existence of rays which indicate
where the sound from the source is going. It does not provide a good solution
when the radius of curvature of the wavefront or the pressure changes
significantly over the distance of a wavelength. It is therefore restricted to short
wavelengths.
The two approaches are compared in Table 3.2. Wave theory will not be consid-
ered further. The Sonar Modelling Handbook (published by DERA, UK Ministry
of Defence) is a useful source of information on the many wave theory and ray
theory models and their practical limitations.
Snell j . law
Snells law is fundamental to ray theory models and describes the refraction of
sound rays in a medium where the speed of sound is changing. When the speed of
sound varies continuously with depth, the medium can be considered as a number
of thin layers of constant but different sound speeds within them. Snells law is
applied to the boundaries of the layers, and the sound ray is seen to be curved:
60 Propagation of Sound in the Sea
If speed continues increasing with depth, the ray will become horizontal and
c o s e = 1. The speed of sound where the ray is horizontal, cv, is known as the
vertex speed, and
This equation makes it possible to generate the path of a ray through the layers
into which the velocity profile has been divided.
Hodgson Example 61
February
The ray trace shows shadow zones between about 8 and 25 km and this is evident
in the PL curves for the two receiver depths of 50 and 200 m. Note that the simple
law of spherical spreading plus absorption follows the model predictions quite
accurately, particularly if the shadow zones are excluded. The reduced losses
centred on about 30 km are due to a CZ gain.
August
Severe downward refraction occurs near the surface, resulting in a shadow zone
from about 2 km until the CZ returns, centred on about 40 km.Note the first BB
returns and the second CZ returns (at greater than 110 dB each way, the losses will
be too great for an active system but passive will sometimes be possible).
Figures 3.13, 3.14 and 3.15 have been produced using the Wader-Hodgson
system. Intellectual property rights and copyright to the Hodgson Acoustic
Propagation Loss model are owned by the United Kingdom Secretary of State for
Defence. Intellectual property rights and copyright to the Wader Global Ocean
Information System are owned by Ocean Acoustic Developments.
62 Propagation of Sound in the Sea
PL (dB)
60
1 00
140
20
Range (km)
PL (dB)
60
100
140
1 I
10 20
(b) Range (km)
PL (dB)
60
100
140
10 20
(C) Range (kin)
Figure 3.15 Propagation loss graphs for the North Arabian Sea: (a) standard SSP, (b) modified
SSP, (c) modified SSP
66 Propagation of Sound in the Sea
3.19 Problem
3.1 What are the detection ranges achieved for a one-way propagation loss of 80 dB
assuming (i) spherical spreading plus absorption, (ii) cylindrical spreading plus
absorption? Give solutions for 5 kHz (a= 0.3 dB/km) and 20 kHz (a = 3 dB/km).
4
Target Strength
4.1 Definition
Target strength (TS) refers to the echo returned by an underwater target:
submarines, surface ships, torpedoes, mines, fish. It is defined as the log of the
ratio, in dB, of the reflected intensity referred to 1 m fiom the acoustic centre of
the target, to the incident intensity:
ITS = lolog($) 1
The use of an arbitrary reference distance of 1 m gives many underwater targets
positive values of TS. This does not imply that more sound is reflected from the
target than is incident on it. The acoustic energy appears to come from a
hypothetical point source within the scattering surface of the target, and for a large
target such as a submarine, this point may be some metres from the surface and
within the target.
4.2 Formulae
The correct value of TS to be used in the sonar equations should be chosen
carefully. In practice, TS is calculated using either the peak pressures of the
incident and reflected pulses or their total integrated energies. Here are the
resultant parameters:
68 Target Strength
PeakTS
Integrated TS
where pi and pr are the peak pressures of the incident and reflected pulses; pi(t)
and pr(t) are the time hnctions of the pulses; Tp is the duration of the incident
pulse and T, is the time extent of the target.
4.3 Measurement
It is simplest to measure the peak pressures of the incident and reflected pulses.
Therefore peak TS is usually determined regardless of the method of measurement
and is the parameter normally used in the active sonar equations.
TS is defined at 1 m from the acoustic centre of a target. It is clearly impractical
and often impossible to make measurements at this distance. Therefore measure-
ments are made at greater distances and reduced to the definition range.
A method well suited to the measurement of the TS of small objects such as
mines is to compare the echo levels from the target with the level fiom a reference
target such as a sphere. For larger targets such as submarines and torpedoes, which
because of their size must be measured at longer range, say 1000 m, the nec-
essarily large size of the reference target (a sphere of TS = 0 dB has a diameter of
4m), the reference target may be replaced by a calibrated transponder or an
alternative method may be used.
Most TS determinations, particularly for these larger targets, have been made
by measuring the peak pressure of the reflected pulse at long range and then
reducing it to its 1 m value. The appropriate sonar equation is
EL = SL - 2PL + TS
The echo level (EL) and the source level (SL) are measured, the propagation loss
(PL) is either calculated, assuming spherical spreading plus adsorption, or meas-
ured using a calibrated hydrophone at the target, and the equation is solved for TS.
TS of'a Sphere 69
Long CW: the small bandwidth and thus low resolution of this pulse type
ensures the PTS is approximately equal to the ITS.
Long FM: the large bandwidth and thus high resolution of this pulse type means
that the effective pulse length is short, of the order of a few milliseconds for an
antisubmarine pulse, and therefore the correct TS value will be given by the
PTS as measured by an equivalent short pulse.
4.5 TS of a Sphere
The meaning of target strength can be demonstrated by calculating the TS of a
sphere. Let a large sphere - large compared with a wavelength - be insonified by
a plane wave of intensity I,. If the sphere has radius a metres, the energy
intercepted by it from the incident sound is xa21j; where za2 is the cross-section
of the sphere, the scattering cross-section.
Assuming the sphere reflects this energy isotropically, the intensity of the
reflected wave at a distance Y metres from the centre of the sphere will be
I, = nu2II/4xr2.At the reference distance of 1 m this reduces to Ir/Ij = a2/4 and
the target strength of the sphere is given by
ITS = 1 0 1 0 g ( ~ )
Note that the scattering cross-section is known in radar as the radar cmss-
section (RCS or a). Expressions for the RCS of simple shapes are prevalent in the
radar literature. If helpful, these expressions can be translated into sonar TS
by dividing them by 4n.The target strengths of certain shapes at an angle 8 to
the normal may be determined with the aid of Figure 4.1, which plots
20 log(x-' sin x) dB against (L/A)sin 8.
Example 4.1
A cylinder has length L = 5 m, radius a = 1 m and A = 0.2 m. What is the TS at an
angle of 2" to the normal?
+
TS = 10 10g(aL~/2A)~ 20 log(x-' sin x ) + 20 log(cos 6 )
(L/A)sin 8 = 0.87
(Llh) sin 0
+
Therefore TS = 18 - 16 0 = 2 dB. At the first null, (L//l)sin 8 = 0.5, therefore
8 = 1.15" and TS is theoretically -moo.
In practice, perfect nulls do not exist; the cylinder will not be perfectly uniform and
the incident pulse will have a finite bandwidth, i.e., the wavelength will vary during
the pulse. Broadband pulses (FM or short CW) will have significant percentage
bandwidths, perhaps 5 per cent and often greater. Therefore, particularly for broad-
band pulses, the nulls will be blurred and the peaks broadened. Narrowband pulses
(long CW) will have less effect. The dashed curve is a realistic plot avoiding the
extreme nulls and peaks. If it is used for this example, we obtain the following results:
F o r $ = 1.15" T S = 18- 1 8 + 0 = 0 d B
For 8 = 2" TS=18-24+0=-6dB
Example 4.2
A rectangle has sides a = 5 m, b = 3 m and I. = 1 m. What is the TS at an angle of 7"
to the normal?
Now the largest dimension is a, replacing L, so we use (a/ll)sin 8 = 0.61 and, fiom
the dashed plot, 20 log(x-' sin x ) = - 19 dB. Therefore
TS = 18 - 19 + 0 = -3 dB
It would be unsafe to attach too much importance to the accuracy of these theoretical
TS values; real targets are seldom regular, and their aspects are not precisely known.
f i e important point is how rapidly TS falls for very small departures from normal
incidence.
where
a = radius
L = length
2xL
x = -sin 8
/i
As we have seen, the TS falls rapidly away from its value on the beam, but because
of the hemispherical ends it does not fall below -22.5 dB.
Again, the TS falls rapidly away from normal aspect (on the beam), but because of
the nose, it will not fall below - 18 dB, except perhaps astern.
Some torpedo designs have flat noses to reduce flow noise and therefore
improve homing sonar performance. TS on the (flat, circular) nose at normal
incidence is
Once again, this high TS value will fall off rapidly away from normal and a f i t
nose torpedo will be virtually undetectable by an active sonar at any aspects other
than close to beam or bow.
TS = 10 log
(3
- (at normal incidence)
at 5 kHz we have
)::;r
TS = lolog =42dB
glint only occurs at this intensity over a fraction of a degree. But since the hull
tends to taper and is seldom exactly cylindncal, this large value is reduced and
exists over greater angles (compare with the narrow beamwidths achieved by long
towed arrays of large apertures). The fact that the TS varies rapidly with aspect
helps to explain the extreme variability of practical measurements.
Thickness (U4)
A prototype thin tuned coating was used in the Alberich design employed by the
Germans during World War 11. It consists of an inner layer of rubber perforated
with a pattern of holes and covered by an outer layer of similar thickness (2 mm).
Resonances occur and the attenuation is high over limited bands. At search sonar
frequencies the resonant approach is essential, and in order to provide adequate
transmission loss, the sound absorber may be combined with a decoupling base
layer. The resonances depend on depth and temperature (Figure 4.4).
140 m, 13C
70 m. 21C
70 m, 13C
1 I I I I
20 40 60 80 kHz
Submarines
Others
4.17 Problems
4.1 What is the target strength of a sphere of radius 0.5 m? Above what frequency is
the result likely to be reliable?
4.2 A cylindrical mine has a length 2 m and radius 0.5 m. The ends are hemispheres.
What are the target strengths normal to the mine and end on? Calculate TS for
frequencies of 10 kHz and 100 kHz.
4.3 What is the target strength of an irregular plate of area 4 m2 and largest dimension
2 m at an angle of 5" to the normal and frequency 4 kHz? Use the practical dotted
curve of Figure 4.1 and take L to be 2 m, the largest dimension.
This page intentionally left blank
5
Noise in Sonar Systems
Thermal noise
0 Noise from the sea (ambient noise)
0 Vessel noise
84 Noise in Sonar Systems
where
This thermal noise is an absolute minimum and would only be observed in the
absence of any other noise source, i.e.,
If the water had no other agitation except thermal agitation and if it were
completely isolated from any source of sound (a dead sea)
If the receivers were perfect and added no noise (the noise factor of the receiver
is NF, = 0 dB).
Note the thermal noise plotted in Figure 5.1 where it forms a lower bound to the
family of sea state noise curves. Sea state, wave height and wind speed are given
in Table 5.1.
Table 5.1 Sea state, wave and wind speed
Sea state 0 1 2 3 4 5 6
Wind speed (knots) Sl 5 13 16 19 22 28
Wave height (m),crest to 0 0.05 0.4 0.7 1.3 2 3
trough
Noise.from the Sea 87
80
- - - NoRainrain
-
4
i
I
Y
Rainfall rate (d)
k! 60
m
2
L_ 30
'u
'u
_I
E
2
4 4
0
3 40
\
.
0.1 I
'
1
DSSl
1 10
Frequency (kHz)
Rain can significantly increase the ambient noise. The spectra are quite flat with
frequency, and Figure 5.2 plots theoretical values which agree quite well with
measurements.
Shipping noise is important, particularly at frequencies below about 500 Hz.
In harbours and harbour approaches, the dominant noise is from shipping -
particularly small craft whose noise extends to frequencies of several kilohertz -
and from industrial activity along the shore.
Biological noise is produced by a variety of marine life. Besides the individual
sounds of sea mammals, snapping shrimps are a notable noise source. They
produce quite a flat spectrum of noise between 500 Hz and 20 kHz, and it can be
as high as 70 dB re 1 p a .
Figure 5.3 plots spectrum levels for these noise sources (except biological).
Some broad practical conclusions can be drawn:
88 Noise in Sonar Systems
0.1 1
Heavy rain can dominate at all frequencies, but rainfall less than about 3 mm
per hour is not significant.
Noise resulting from the agitation of the sea due to wind and waves is dominant
above 500 Hz.It is difficult to relate the actual noise level to the prevailing wind
speed or sea state; the wind speed may be measured but the sea state is often
based on a subjective judgement of wave height, and although wind speed and
sea state are clearly related, rarely do their maximum or minimum values
coincide in time.
Shipping noise is unimportant above 1 kHz (radiated noise itself falls rapidly
with frequency and also falls with range due to absorption losses). The noise
below 1 kHz is likely to be dominated by discrete tones, but if there are many
ships then the noise may appear to have a continuous spectrum.
89
10 20 30
Speed (knots)
The designer will ensure that the internal noise of the sonar is negligible, nothing
can be done about the ambient noise of the sea, and therefore the design of the
total system - vessel, dome and array - must attempt to minimize the, usually
dominant, own vessel noise. A practical aim is to attempt to match the ambient
noise of the sea corresponding to, say, SS2 for surface ships up to operating speeds
of 15 knots and somewhat better for submarines.
5.6 Self-noise
The self-noise of a sonar is perceived as a noise voltage appearing at the output of
the receiver. It is more convenient, however, for both design and performance
analysis, to define it at the entry to the system, i.e., at the array/water interface.
The noise intensity is given in dB relative to the omnidirectional intensity of a
plane wave having a pressure equal to 1 @a in a 1 Hz band, the isotropic spectrum
level. If we assume that the vessel is in a calm sea and there are no parasitic noise
sources, the self-noise will be solely due to the combination of vessel and sonar.
We now discuss its components and indicate the principal methods for reducing
them.
0 Isolate machines, shafts and piping from the hull using anti-vibration mounts.
Isolate the array and dome from the hull; provide damping between the array/
dome and the hull; dampen the hull itself.
Separate arrays and machines; site the arrays well forward of the main
propulsion machinery; avoid placing any auxiliary machinery close to arrays.
Agouti: a device which emits air bubbles in the vicinity of the blades to replace
the water vapour bubbles created by cavitation. An improvement of the order of
10 dB is possible at high ship speeds (above 20 knots).
Bafles: these are mounted within the dome to shield the array from propeller
noise. They inevitably reduce the performance of the sonar over the stem arc
protected by the baffle; a useful rule of thumb is to assume the allowable
propagation loss is reduced by 10 dB over the stern arc.
5.13 Directivity
Self-noise is markedly directive. There is clearly a maximum over the stem arc,
even with an intervening baffle. Measurements are frequently made in octave
bands and averaged over operational bearings (the forward 270), thus avoiding
any bias which would result from including measurements astern.
Self-noise: the total noise measured by a sonar system which forms the back-
ground to detection of a wanted signal.
0 Radiated noise: the noise radiated by a vessel and measured at some distance
from the vessel, typically between 100 and 1000 m. It is the source of signals
for passive sonars (Chapter 8).
Self-noise and radiated noise have many common sources and are frequently
simply two different aspects of the same phenomenon. For example, propeller
noise may dominate radiated noise at all speeds but only be an important
contributor to self-noise at higher speeds; flow noise may dominate the self-noise
when its source is close to the array but be an insignificant contributor to the
radiated noise.
where PT is the total noise power and PI and P2 are the individual noise sources,
all in dB re 1 pPa.
When P1 = Pz then PT = P2 + 3 dB
When P1 = P2 + 6 dB then PT = P1 + 1 dB
94 Noise in Sonar Systems
Therefore if one noise level exceeds the other by at least 6 dB, the other makes a
negligible contribution to PT. The equation is plotted in Figure 5.5.
where Noutis the noise power at the output of the receiver and N m i n is the un-
avoidable output noise, i.e., the noise at the input of the receiver due to the thermal
noise of the generator connected to the receiver input multiplied by the receiver
gain.
Noise Factor of a Sonar 95
Receiver
+
where R = R m I$ is the total hydrophone resistance.
The noise factor of the receiver alone is NFr. The noise factor of the receiver
plus hydrophone, NFrh, is greater. This is because, in defining NF,, the thermal
noise due to R is unavoidable, whereas for the ensemble (receiver plus hydro-
phone) the only unavoidable noise is due to R , (which is smaller than R and
would apply to a perfect hydrophone where Rp = 0).
The noise power, Nmin, is therefore reduced by r j ~= R,/R, hence
Now consider the case where the sea state is not negligible. The noise factor
becomes NF, (at a given frequency), and
NF, =
isotropic spectral noise intensity - N,
-
+ Nmin
thermal noise (min) Nmin
96 Noise in Sonar Systems
and the resultant noise factor for the sonar, NF,,, for a given sea state, will be
and since
we obtain
NFres =
Nmin + Nrtr + Ns
Nmin
Using the expressions for NF, and NFrh obtained above, we can deduce the
following expressions for the overall noise factor for the sonar in a given sea state,
NFres
-4cceptable Receiver Noise Level 97
which is the minimum value of NF,,, for a given sea state. We lose, therefore, due
to the noise introduced by the sonar itself, in the ratio
Example 5.1
Consider a sonar operating at 5 kHz on a fi-igate whose self-noise is unlikely to fall
below the ambient noise, which is 50 dB re 1 @a for SS2.
We have
10 log NF, = 50 dB
Suppose that we only lose 1 dB due to the internal noise of the sonar, then
If the efficiency of the hydrophone is 40 per cent (-4 dB) then the noise factor of the
receiver itself will need to be less than 39 dB, which is easy to achieve.
The hydrophone sensitivity is given by Sh = 20logv - 201og p, and if S h =
-200 dB/V then
i.e., the spectrum level of the receiver noise, referred to the input, must be less than
10 nV.
Example 5.2
Now consider a sonar, again operating at 5 kHz, installed on a submarine whose self-
noise, in its quietest state, can be less than the ambient noise of the sea, and take the
ambient noise in sea state 0 as 35 dB re 1 p a .
We have
lOlogNF, = 35 dB
Again suppose that we only lose 1 dB due to the internal noise of the sonar, then
NFr
5 = 1.26 and -= 0.26 X 3162 = 29 dB
N F S 7
and for the same hydrophone, the noise factor of the receiver has to be less than
25 dB:
i.e., the spectrum level of the receiver noise, referred to the input, must be less than
2 nV.
Alternative Calculation 99
For 1 dB loss due to receiver noise, and again assuming 4 dB loss in the hydrophone,
the spectrum level of the receiver noise must be
50 - 6 - 4 = 40 dB
From Figure 5.5, 1 dB loss requires the receiver noise to be 6 dB down. Now we
substitute 40 dB into the sensitivity equation:
i.e., the spectrum level of the receiver noise, referred to the input, must be less than
10 nV:
For 1 dB loss due to receiver noise, and again assuming 4 dB loss in the hydrophone,
the spectrum level of the receiver noise must be
35 - 6 - 4 = 25 dB
100 Noise in Sonar Systems
Therefore
i.e., the spectrum level of the receiver noise, referred to the input, must be less than
2 nV
Table 5.2 Practical noise values for design and performance comparison
Isotropic spectrum level (dB)re the intensity of sound due to a
pressure of 1 Pa
The thermal noise of the sea is increasingly the dominant source of noise and
reverses the steady falls in both ambient and self-noise with frequency.
Prt,hEerns 101
5.2i Problems
5.1 A hydrophone with a sensitivity of - 170 dB/V receives over a 1000 Hz band
centred on 100 kHz.What would its output be in SS2?
5.2 The self-noise of a sonar is 55 dB. If the ambient noise of the sea is also 55 dB,
what would be the level in a third-octave band centred on 6 kHz measured at the
output of a receiver beam of DI = 20 dB?
This page intentionally left blank
6
Reverberation
Here Zr is the specific acoustic impedance of air and 2, is the specific acoustic
impedance of water. Almost all of the sound is reflected and only a tiny fiaction,
0.0005, escapes into the atmosphere. Note the phase change.
When a reflection occurs at the sea bed, the situation is much more complicated.
Because Z,, the specific acoustic impedance of the sea bed, is now much greater
than Zw, the specific acoustic impedance of water, the value of Cr is positive (but
still almost unity) and there is no phase change at the reflection. Once again, much
of the sound is reflected and only a fi-action penetrates the sea bed.
The sea bed, however, has highly variable acoustic properties due to its
composition and inhomogeneous nature, with the result that Z,, the specific
acoustic impedance of the sea bed, changes markedly with both depth of penetra-
tion and laterally. A sufficiently accurate and detailed knowledge of the bottom to
ensure useful predictions for scattering, and for propagation within the bottom
itself, is never realistically available.
Rp = khsin8
where h is the wave height, measured from trough to crest, 8 is the grazing angle
and k = 2 n / A is the wave number. In practice when 4 << 1 the surface is taken to
be smooth and when l$ >> 1 the surface is rough.
Figure 6.1 illustrates the reflection and scattering of sound from (a) perfectly
smooth boundaries and (b) rough boundaries. The boundary shown is the sea
surface; inversion of the figure represents the similar situation when the boundary
is the sea bed.
C1assc.s of Reverberation
i 0.0005
Smooth
surface
/' i \
1.0 / \, 0.9995
/
1 \
+
Rough surface
f,
L
Specular direction
Figure 6.1 Scattering of sound from a boundary: (a) perfectly smooth, (b) rough
For analytical purposes, sea surface and sea bottom reverberation may be lumped
together as boundary (area) reverberation.
1 S,," = 10 log -
I
;
;J
where Iscatis the intensity of the sound scattered (back to the source) by a unit area
or unit volume, referred to a distance of 1 m from its acoustic centre, and Ii is the
intensity of the incident plane wave.
where A and V are the total reverberating area or volume, respectively, as defined
by the 'two-way' beamwidth of the sonar - the effective beamwidths of the
projector plus hydrophone combination in both azimuth and elevation. For
example, the horizontal two-way beamwidth of the combination of an omnidirec-
tional, in azimuth, projector and a 10" , in azimuth, receive array is 10" (the
smaller of the two).
Volume Reverberation 107
The length of the cylinder, dV (its extent in range), is such that the reflections
from all scatterers within it return to the source simultaneously. (The scattering of
the front end of the pulse by the rear scatterers in dV returns to the source at the
same time as the scattering of the rear end of the pulse by the front scatterers
in dV)
The extent in range is therefore cT/2, where T is the pulse length and c is the
speed of sound. (Note the division by 2, a constant source of confusion between
active, two-way, and passive, one-way, sonars.) The volume can now be written
V=-R2
I, B,Btd6
108 Reverberation
and if R is large compared with the cross-section of the volume (say R > 4R@,
the volume is given by
cT nR2&Ov
V=--
2 4
where & and OV are the two-way horizontal and vertical beamwidths expressed in
radians. Figure 6.3 is a graphical representation of this simplified expression for
the reverberating volume.
Example 6.1
A sonar has two-way horizontal and vertical beamwidths of 12" (0.2 rad) and transmits
a pulse of length 1 ms. What is the TSR at 1000 m for S, = -80 dB?
We calculate
and
When the beam patterns intercept a boundary, i.e., the surface or the bottom of the sea
- which, except in deep water, occurs at quite close ranges (perhaps 1000 m) - the
volume will be smaller than given by the formula. It is then possible, but not usually
necessary, to work out the cross-sectional area from a knowledge of the depth of water
and the depth of the source and to use this to determine the volume.
with @, in radians.
Example 6.2
For the same sonar, what is the TSR at 1000 m for S, = -30 dB?
We calculate
Scattering Layers 111
and
In deep water the dominant reverberation will be from the sea surface.
In shallow water (< 200 m) the dominant reverberation will be from the sea bed
unless the wind speed is high (see later).
I S[ = s, + 10 log H 1
and this value must be used in the expression for boundary reverberation.
112 Reverberation
\
\
'\
/ Night
/
/
\
/
\+L
\
\ 100 m
/
/
Frequency 5 1<Hz
Volume scatterers are not uniformly distributed in the sea and tend to be
concentrated in the deep scattering layer (DSL). This source of reverberation is
overwhelmingly biological in nature and is a complex mix of different organisms
resulting in scattering strength versus depth profiles that change with frequency,
location and time. The representative S, profiles in Figure 6.5 show vertical
migration of the DSL with time of day. Estimates of thickness are shown and
Grazing an@e
ss (dB)
-40
-50
- 60
-30
-40
-50
Sea state 6
Wind speed 30 knots
2 4 6 8 10
(b) Frequency (kHz)
The roughness of the sea surface, and the possibility of bubbles trapped just
beneath it, results in significant scattering of sound. The scattering strength varies
with the grazing angle of the incident wave and with frequency. Empirical
formulae by Chapman and Harris fit the many measurements quite well, particu-
larly at frequencies between 1 and 10 kHz.Figure 6.6 plots these formulae for two
representative sea states. Note that S, increases markedly with grazing angle and
also with frequency (around 3 dB per octave at the lower grazing angles which are
of particular interest to the sonar designer). The dashed lines show recommended
values for use in sonar design, representing high sea states (-40 dB) and low sea
states (-50 dB).
114 Reverberation
Figure 6.7 shows two directivity patterns for the sound incident upon the bottom
at point 0. The radius of the patterns in any direction is proportional to the
intensity of the sound reflected in that direction.
Bottom scattering strengths are generally greater than those from the surface
and therefore, when the geometry is such that the sonar vertical beam patterns
intercept the bottom (which is always the case for shallow water at practical
ranges), then bottom reverberation is often the dominant background to detection
and classification.
Figure 6.8 illustrates the wide range of S b which may be encountered in shallow
waters. We shall see later that s b is an extremely important factor in determining
the Performance and design of an active sonar and if, during its design, an
incorrect value is chosen for the maximum Sb likely to be encountered, the sonar's
performance will not be optimum.
Over a uniform bottom the basic detection performance will be determined by
the average value of Sb, which is little affected by comparatively rare, submarine-
116 Reverberation
like false alarms. The bottom may be mud, clay, sand, shingle, pebbles or a mix of
these and the calculated & (derived From measurements of reverberation levels)
falls within the range -35 to -50 dB at the low grazing angles encountered in
shallow water.
Over a non-uniform bottom (rocky outcrops, shorelines, undulations) the
average value of Sb is greatly influenced by such local features and it is suggested
that the actual values of &,, used to determine the initial detection ranges of a
sonar, should be some 5-10 dB less than the average measured values. Thus,
values for &, which are greater than -35 dB, i.e., falling within the range -35 to
-25 dB, and typical of non-uniform sea beds, should be reduced by 5 dB but
never to less than -35 dB. Table 6.1 should make this clear.
Table 6.1 Reduce St, values by 5 dB but never to less than -35 dE3
&, measured (dB) -25 -27 -29 -3 1 -33 -35
Use (dB) -30 -32 -34 -35 -35 -35
This argument is supported by initial detection trials where, over such difficult
bottoms, final detection and classification are delayed by a multitude of false
alarms but post-trials analysis Frequently reveals missed detection opportunities at
greater ranges (consistent with postulating lower values for Sb).
4- ~
5-10 kHz Summer
_cc____ 1-3 kHz polar ice 1
- 40 /
/
0
rr--
------ "
/ -
Ice-free
Wind speed = 30 knots
/
5 10 15 20
Grazing angle (degrees)
The underside of the ice cover in polar regions has a very high backscattering
strength, S,, when it is rough and non-uniform (e.g., spring pack ice). When the
underside is smoother (e.g., summer polar ice), Si is significantly less. Most of the
Si measurements show an approximately 3 dB per octave frequency dependency
and are generally greater than S,measurements in ice-free water (Figure 6.9).
6.15 Problem
6.1 A sonar system has a horizontal beamwidth of 10" and transmits a pulse of
duration 100 ms. If the backscattering strength is -40 dB, what is the reverberation
target strength, TSR,at 4000 m? What is the TSR at twice the range?
This page intentionally left blank
7
The Sonar Equations
control. This process is often severely constrained by, for example, platform size
or equipment cost and is considerably aided by the intuition of an experienced
sonar designer.
SE = S - N - DT
Signal power is in the analysis bandwidth and noise power is in a 1 Hz band. This
threshold is defined such that a signal with S - N (at the output of the
beamformer) equal to DT, has a specified (often 50 per cent) probability of
detection ( p d ) for a required probability of false alarms (Pfa). The value of DT
depends on these probabilities and the sonar signal processing. Put another way, a
signal excess of zero corresponds to a pd of 50 per cent and a positive SE indicates
pd > 50 per cent.
Display/auto
Preset Pd for given Y,
system
T
MDS DT d
The detection threshoZd (DT) is the ratio of mean signal power to mean noise
power (expressed in dB form) measured after beamfoming, which results in the
preset P d and Pfa values. The signal power is in the analysis bandwidth and the
noise power is in a 1 Hz band. Detection index, d, is defined as
d=(
+
mean (S N) - mean (N)
std dev(N)
Figure 7.2 shows curves of probability density plotted for N alone and for
S + N.
+
The detection index, d, is equivalent to the ratio (S N ) / N of the envelope of
the signal processing output where the threshold, T, is established. The area under
+
the curve of S N above T is the probability that an amplitude in excess of T is
due to signal plus noise; this area equals P d . The area under the curve of N above
T is the probability that an amplitude in excess of T is due to noise only; this area
equals Pfa.These two probabilities therefore vary with the threshold setting and
depend on the value of d.
Pd 1 .o
0.5
0.2
levels of signal variability are shown, defined in terms of the fluctuation index (FI)
given by
FI =
o(S + N)
uN
+
FI = 1 corresponds to a stable signal where the variability in S N is entirely due
to N.Increasing FI corresponds to increasing the signal variability.
124 The Sonar Equations
7.8 ROCCurves
There is no one ROC curve applicable to all sonars. Figure 7.3, which assumes
Gaussian statistics, is applicable to sonars with large BT products (broadband
sonars). Table 7.1 gives values of 5logd taken from the Sonar ModeZZing
Handbook ROC curves together with their recommended usage. But here are some
useful rule-of-thumb values:
Broadband: 6 dB
Intercept and communications: 10 dB
Narrowband and active: 10 dB
7.9 Problem
7.1 The Pd of a sonar is to be 0.5. What will be the change in Pf, which results from
reducing 5 log d from 6 dB to 5 dB? Use Figure 7.3.
8
Passive Sonar
Figure 8.1 Radiated noise: continuum and tonals for (a) low speed and (b) high speed
For surface ships and submarines this frequency is between 100 and 500 Hz,
depending on speed and, for submarines, depth. For torpedoes the crossover
fiequency is higher, perhaps between 500 and 1500 Hz, and the tonals extend to
higher frequencies, say 3000 Hz, because of the higher speeds of machinery and
propellers.
Pructical Values 127
I60
Frequency (Hz)
(therefore they need ever longer towed arrays) to detect the remaining lower
frequency lines, which are difficult to eliminate.
An older submarine, for example, may have a tonal at 200 Hz of intensity
125 dB, but tonals near this frequency for a modern submarine will have been
reduced to a value below the broadband continuum which tends to level off in
intensity at, say, 300 Hz (and perhaps 1000 Hz for torpedoes).
Broadband and Narrowband 129
the signal bandwidth, B,, is greater than B, the signal is over resolved and there will
be a reduction in gain of 10 log B / B,. Electronic integration sums the energy from
each analysis cell over a number of time steps. The number of steps is known as the
integrationfactor (IF) and the gain is 5 log IF.
8.6 Normalization
There are significant differences in the mean noise levels with fiequency, bearing
(beam) and time, which must be reduced before display to the operator. Normal-
ization schemes all depend on using the information surrounding the cell to be
displayed to reduce these differences to the dynamic range of the display. Figure 8.4
shows an algorithm for normalizing a line of data in the frequency dimension before
display. The same process is repeated for all analysis cells.
Power
Cell i
Surface ship hull-mounted arrays will generally take the form of cylindrical arrays
mounted in bow or keel domes. Operating frequencies may extend down to about
100 Hz for narrowband tonals and up to about 10 kHz for broadband. The DI will
clearly be small at frequencies below about 1 kHz and the self-noise high.
Performance against submarines will be ineffective but torpedoes (apart from
quiet, electrically propelled, marks) will be detected at adequate ranges for the
deployment of countermeasures.
Line arrays
Line arrays towed from submarines or surface ships - by virtue of their length,
typically 32 or 64 wavelengths, and low operating frequencies, down to about
10 Hz - are very effective against both submarines and torpedoes. Separation
from the tow vessel assists in reducing the self-noise seen by the towed array,
particularly at bearings other than those of the ahead beams which receive the
stern radiated noise of the tow vessel in their main lobes.
133
As a frequency analyser, the human ear can resolve two sounds which differ by as
little as 50- 100 Hz over much of the audible band. Because of these character-
istics, the ear is still important in passive sonar, and perhaps unrivalled as the final
arbiter in the classification process.
The ear is not a simple broadband listening channel, but rather may be likened
to a contiguous comb of narrowband filters covering the entire audio spectrum.
The bandwidth of each of these hypothetical filters is known as the critical
bandwidth (Af,) of the ear and is between 50 and 100 Hz for frequencies between
300 and 2000 Hz and increasing somewhat outside this band.
The quantity AfC is the bandwidth such that an increase in spectrum level
(within reasonable limits), outside this band has no masking effect on a signal
Myithin the critical band (Figure 8.5). Therefore, the masking noise is given by
to C(f>and, from the concept of critical bandwidth, it follows that as the level of
a signal (which is in fact the wideband radiated noise emitted by the threat) is
increased, as soon as this level reaches the noise level at any frequency within the
receiver bandwidth, it is perceived (at P d = 50 per cent) as a signal at that
frequency, known as the thresholdfrequency, f t .
Reducing the listening bandwidth risks eliminating the threshold frequency (or
moving it to another value where the noise is higher) and therefore the audio
bandwidth should cover the entire frequency spectrum of the expected signal.
The ear will have detected the signal at the level represented by Figure 8.5(b),
but a simple broadband receiver will not detect the signal until the level is
increased to about that represented by Figure 8.5(c), a positive signal-to-noise
ratio.
Passive Displays 135
.' '
; '..-.
; *:
**
. 3
L
- \
-...
*
*.
*
*. *.
:
I
:
I
I
#
#
- - - - .
2
Time
*.
.
.
8.
9-
*
I
I
9
. } Beam 1
+
mean(S N) - mean(N) signal power
( std dev( N)
Broadhand Square Law Detector 137
and, substituting for DT, the sonar equation for this passive receiver, also called an
energy detector, becomes
The expression for DT assumes a decision based on a single sample (or line for a
practical multibeam receiver), displayed to the operator or passed to an automatic
system. In practice, operators and automatic systems make decisions based on a
number of successive time samples. The equation must therefore be modified by
including a visual integration term, 5 log n, where n is the number of successive
time samples used by the operator or automatic system in making the decision.
Putting SE = 0 gives
Note that the DI is now that of the half-beams, or 3 dB less than for the full
aperture. However, the factor of 2 in the 5 log 2 BT, term implies a net loss of
1.5 dB (3 - 1.5) compared with the energy detector, a small price to pay for the
improved bearing discrimination.
Narvo whand Processor 139
Here signal and noise apply to identical bandwidths, but for narrowband the signal
is the intrinsic bandwidth of the line, whereas the noise is referred to a 1 Hz band.
Both levels must be referred to a common bandwidth, the analysis bandwidth, B ,
of each frequency cell. The input noise to the narrowband processor is therefore
narrowband, and for the narrowband case
An integration factor, IF, which is the number of independent samples from the
signal processor summed before display.
DT=510gd-510g
(3
-
If the signal bandwidth, B,, is greater than the analysis bandwidth, B, the signal
will be over resolved - it will appear in more than one cell - and there will be a
reduction in gain of 10 log( B/ 73,). Note that any mismatch cannot increase
processing gain and therefore any positive value for 10 log( B / B,) must be reduced
to zero.
I40 Passive Sonar
I - lOlogB+ lOlog(B/B,)+5logn
141
Example 8.1
The submarine passive sonar receives using flank arrays of length 8 m and height 4 m.
If the design frequency of the array is 2000 Hz then /2/2 = 375 mm and the array will
have 10 rows of 20 elements.
Assume the self-noise of the submarine to be equivalent to SS2. The background noise
will then be ambient sea noise or SS2, whichever is the greater. The DI will be given
by
At 3000 HZ +
DI = 3 10 log200 + 20 log 1.5 = 29 dB
At 400 HZ +
DI = 3 23 - 2010g5 = 12dB
At 200 HZ +
DI = 3 23 - 2O10g 10 = 6 dB
At 8 0 H z DI=3+23-2010g25 = 3 dB
At 40 HZ +
DI = 3 23 - 2010g50 = 3 dB
142 Passive Sonar
Because the flank arrays are baffled by the hull, DI cannot be less than 3 dB. The
performance of the system is given by the following tables. The ranges have been
calculated assuming spherical spreading plus absorption.
Narrowband
Broadband
Note that all three threats are detected by their broadband radiated noise at long
range but, except for the noisy torpedo (threat A), cZassiJication by narrowband
tonals is only possible at very close range. This is because, at the lower
frequencies of the tonals, the background noise is much higher and the DI of the
array is low (12 dB at best).
Ever quieter submarines and torpedoes are the stimulus behind the development
of towed arrays, which can achieve a good DI down to the very low fiequencies of
the residual tonals of modern threats. The second example will show how, even
when towed by a surface ship, towed arrays significantly improve classification
ranges.
Example 8.2
The towed array will be 321 at all frequencies by selecting elements. The background
noise at the array will be the ambient sea noise or the radiated noise of the tow ship at
the array, whichever is the greater. Note that a good towed array can be expected to
achieve self-noise levels equivalent to SS2. Take the spectrum level of the tow ship
radiated noise to be
The radiated noise at the array will be reduced by the separation from the ship (cable
length) and by the main lobe to sidelobe ratio of the array. Take the cable length as
600 rn and the ratio as -20 dB. Then, from 40 to 400 Hz, radiated noise at the array is
144 Passive Sonar
Operational bearings for hull-mounted arrays are conventionally fiom R 135" around
to G135" (red for port, green for starboard), i.e., those bearings not directly affected
by stem noise. By the same token, operational bearings for towed arrays are from
R45"around to G45".
If the ambient sea noise is taken as SS2, the levels are as follows:
~~~ ~
32 X1500
= 1200 m (using c = fA)
40
The performance of the system at operational bearings is given by the following
tables. The ranges have been calculated assuming spherical spreading plus
absorption. (Figure 3.3).
Narrowband
We have
Broadband
We have
Note that all three threats are now detected by their broadband radiated noise at
long range and also clussijied by narrowband tonals at useful ranges. This is
because, even at the lower frequencies of the tonals, the DI of the array is
maintained at 18 dB.
Over operational bearings, performance is limited by the ambient noise of the
sea and therefore the performance would be the same if the tow vessel were a
submarine.
8.18 Triangulation
Given two well-separated arrays, the range can be estimated using simple
trigonometry. This technique, called triangulation, is not limited to sonar and is
best illustrated by a practical example. A submarine is in contact with a target at
range R from a flank array and from a towed array whose centre is S = 1000 m
behind the flank array. The maximum response angles (MRAs) are shown in
Figure 8.8.
From simple trigonometry, given one side and two angles, we have
The accuracy of the range estimate is proportional to the angle subtended at the
target by the two arrays. Increasing the separation, S, of the arrays andor reducing
the range of the target increases this angle and improves the accuracy.
The arrays must be large enough for the beamwidths to be small at the
frequencies of interest. Furthermore, increasing beamwidths as steer angles move
towards endfire will limit the usefulness of the technique to angles not too far
removed from broadside.
Practical limits to the technique might be
Max R / S = 5
Max beamwidth = 5"
Max steer = f60" from broadside
The separation (or baseline) for measurements may be greatly increased by using
own vessel movement. This will greatly increase the baseline for the bearing
measurements. But the measurements are now also separated in time, and in
practice the target will also be moving. Initial assumptions are made for target
course and speed, e.g., constant values. Many mathematical procedures are avail-
able to refine these initial assumptions and the resulting range estimate will
improve with observation time. The technique is known as bearings only analysis
(BOA) and is clearly an implicit form of triangulation.
Vertical Direct Passive Ranging 149
- - - - __ - --
R B
Platform and target depths, p and t, are much less than the water depth, H
* The sea bed is flat at the bounce site
The range of the target in terms of arrival angle is given by the following
approximate expressions:
I R = (223 - p - t)/tan61
150 Passive Sonar
SB (not shown)
Phase comparisons between vertical half-beams (formed using top and bottom
halves of the vertical staves)
Example 8.3
Suppose the signal is contained within an octave from 1000 to 2000 Hz and the height
of the array is at least 5 wavelengths over this band, then from Figure 2.4, the vertical
beamwidth for steers up to 530"will be at most about 12". If we assume amplitude or
phase comparisons will improve on this by a factor of 4, then measurements of the
vertical arrival angle will be up to f 3 " in error.
For a simple bottom bounce only, H = 1000 m, p = 100 m, t = 0 (surface ship) and
8 = 30", the range is R = 3290 m. The limits (for 8 = 27" and 33") will be 3728 m
and 2925 m, or approximately f12%.
This might be seen as an acceptable result for a passive ranging system. But it does
depend fundamentally on the accuracy of 8, which in turn implies a narrow vertical
beamwidth. Unfortunately, detection - which must of course precede ranging - is
increasingly only likely at low frequencies where beamwidths will be large. Therefore
ranging using 8 only is likely to be very inaccurate.
The time diflerence of arrival, dT, between the two paths is given approximately by
Ertical Direct Passive Ranging 151
and for this example 6T = 67 ms. 6T is now used as an initial value for the off-
set between two 'windows' of time series data from the two paths shown.
Cross-correlation is performed between the two time series and the correlation peak is
used to improve on the initial value of 67 ms.
Suppose the improved value is 80 ms; using the above formula, 8 = 37" and the
range can now be determined from the earlier formulae:
There is an alternative method which can also yield an estimate of target depth.
Simple geometry in the vertical plane gives the following expressions for the sound
transit times for the four possible paths that include one bottom bounce:
1
Tss = ;[IP + (2H + p -- f)2] "2
and for this example, provided p and t are small compared with H , then 6T
= TBs - Ts is approximately
I I
If 6T = 80 ms then for this example R = 2667 m. This is very close indeed to the
average value obtained above (2654 m), by simply using the angle derived fi-om the
measured 6T in the coarse range formulae.
152 Passive Sonar
S S
Having obtained the range from one time difference, the depth of the target can be
estimated from the time difference between two other paths. If we take the difference
between the two paths in Figure 8.10,dT = TSBS- TB,then we obtain
Suppose 6T = 200 ms and the range, already obtained from different paths, is 2667 rn
then t = 150 rn.
Horizontal Direcf Passive Ranging 153
Target
,.+ f
/
+
Axis of arrays A
R2 = p2 + ( R - dr)
R2 = (dsin8)2 + ( R - (rr)2
To find the range, we therefore need to determine dr and 8. For a distant target,
angle 8 is almost the same for all three arrays. In practice, determine 8 by
measuring and averaging the maximum response angles from the three arrays; dr
is given by
dr = BE - tCF
Target
fl
Example 8.4
Three arrays, A , B and C, are mounted on the flank of a submarine (Figure 8.12).
Signals from a target at 30" to the axis of the arrays are detected in an octave centred
on 3000 Hz at all three arrays. The delays measured by cross-correlation are
BE = 0.0067 c = 10.05 m
CF = 0.0133 c = 20 m
= 10.05 - 10 = 0.05 m
For a target on the beam (6 = 90") then, for the same 6v, we have I< =
(20 sin 90")'/(2 X 0.05) = 4000 m. Also, for the same bearing measurement errors,
the actual range would be 3988 m, i.e., a very small range error. Therefore, not only is
6 r largest for the same range when the target is on the beam, but also range accuracy
is least sensitive to bearing errors. The submarine can optimize the range estimate by
manoeuvring to put and keep the target on a bearing close to broadside.
In practice the arrays will not be equispaced and neither will they be aligned in the
remaining two axes. The full equations for range and bearing of a target will therefore
be three-dimensional and more complex, although based on the same simple
geometric relationships. The simple case for collinear, equispaced arrays has clearly
demonstrated the principle of HDPR and how a range estimate may be optimized.
I56 Passive Sonar
Hose
Tail
VIM ~ o o o o o o o o o o o o o o o o o o o o o o o o o o o o o o ~
Element 1 Element n
A typical towed array will comprise several hundred hydrophones, together with
electronic circuits to preampliQ, sample and digitize their outputs. The hydro-
phones are enclosed in a plastic hose, liquid filled to achieve neutral buoyancy in
the sea. The array is towed by a cable of a length decided by the speed of the tow
vessel and the desired depth of operation; it will be several hundreds of metres. A
vibration isolation module (VIM) reduces the tow vessel vibrations transmitted
down the cable, which would otherwise augment the noise of the array. The
assembly is completed by a tail whose purpose is to provide some drag to the
array so that it is maintained reasonably linear.
In practice the array will have significant curvature, which would affect the
beam shapes unless corrected in the beamformer. Heading sensors are placed at
intervals within the array and their readings used to correct for array curvature;
there will be at least three sensors (front, rear and middle) and possibly more in
very long arrays.
Towed arrays can have diameters as small as, say, 50 mm (thin arrays) and as
large as, say, 150 mm. The length of a towed array is determined by its fiequency
of operation and the desired DI. When used as a receive array in an active system,
the length might be from 10 to 50 m, whereas for a passive array the length might
be from 100 to 1000 m.
Passive towed arrays operate over a frequency range of several octaves. The
spacing between the elements is maintained at about A / 2 at the centre frequency
of each octave by a suitable choice among the available elements.
Example 8.5
Design a towed array covering 200 to 1600 Hz with a DI of about 16 dB at all
frequencies within this range.
DI = 10 log IZ = 16. Therefore n = 40 and each octave must use 40 elements spaced
ilJ2 at its centre frequency (Table 8.7). The complete array is shown in Figure 8.14.
The higher octaves are nested (a term possibly borrowed from computer program-
ming) within the lower octaves, and elements are chosen to meet the requirements of
each octave. To gain 3 dB in DI, or to operate at the next lower octave, the array must
be doubled in length.
Course change
If the tow vessel, and therefore the towed array, changes its heading, it is possible
to resolve the ambiguity as shown in Figure 8.15. The array heading change is not
instantaneous and, particularly for a very long array, can take a considerable time.
Nor is the target stationary. Nevertheless, the true target bearing is often quickly
resolved, but may need confirming by making some assumptions about the targets
motion.
Array heading
Change to array heading
No
Parallel twin towed arrays use the time delays between the signals arriving at the
two arrays to resolve the lefthight ambiguity. Maintaining the spacing between
two flexible arrays is a practical problem, particularly for very long arrays and
during a change of course. Precision, however, is not necessary provided some
horizontal spacing survives and the arrays do not cross over. If both arrays are also
used to form beams, an improvement in DI of up to 3 dB will result over a limited
frequency range (particularly usefid when the twin array is used as a receive array
in an active system which will have a comparatively small percentage bandwidth).
Triplets
8.23 Self-noise
Towed arrays are well separated from tow vessels and therefore the vessel radiated
noise is significantly reduced (given a cable length of 500 m, the value of 20 log R
is 54 dB) and, except for the ahead bearings of the towed array, this is further
reduced by the main lobe to sidelobe ratio of the beams (perhaps 20 dB). The
160 Passive Sonar
hydrodynamic noise of the towed array can be made negligible at normal tow
speeds of up to about 12 knots, and therefore the remaining and dominant noise is
the ambient noise of the sea.
8.24 Problems
8.1 The dimensions of the flank arrays of Example 8.1 are halved. How are the
narrowband and broadband performances affected?
8.2 Refer to Figure 8.9. The platform depth, p, is 300m, the target depth, t, is
unknown but small compared to the water depth, H, of 2000 m. If the time difference
between the paths shown is 100 ms, what is the range of the target?
8.3 Design a towed array to have a DI of 21 dB at the centre frequency of the octave
from 2000 to 4000 Hz. What would be the spectrum level of the radiated noise of a
target which is first detected at a range of 10 km against a background noise equivalent
to SS4? Assume broadband detection and 5 log d = 6 dB, T, = 10 s, n = 60. Propaga-
tion is spherical spreading plus absorption.
9
Active Sonar
The first and fourth components are normally removed by own doppler nuZZiJica-
tion (ODN). The total frequency shift due to the target doppler, S cos cp, is given
by
lAf 2s cos cp
=flC
where f is the operating frequency. The factor of 2 results from the two equal
doppler shifts due to target motion (the second and third components).
Practical formulae for doppler shift are, for c = 1500 m/ s,
Note that 1 m/s = 1.945 knots. And remember that doppler is a shift of relative
velocity not speed. For a passive sonar, only the third and fourth components are
relevant. The frequency shift is positive for an approaching target but may be
changed by heterodyning in the receiver.
C W Processing 163
Frequency modulation (FM): the frequency of the pulse changes during the T
seconds duration of the pulse. The bandwidth, B , is not now the inverse of the
pulse length.
9.3 CW Processing
CW processing is similar to narrowband passive processing where, after beam-
forming, the data is analysed using a fast Fourier transform (FFT) process and
then displayed to an operator andor input to an automatic detection and classifica-
tion system. Traditionally, active CW processing used a bank of analogue filters,
each matched to the bandwidth of the CW pulse. Sufficient filters were provided to
cover the expected frequency shifts due to doppler. In modern sonars the analogue
filters are replaced by an FFT processor.
164 Active Sonar
9.4 FM Processing
FM processing replaces the FFT by a process where the output of the beamformer
is correlated with a replica of the transmitted pulse. The signal processing for both
classes of pulse is essentially matched filter processing and the same expression
for signal processing gain therefore applies to both.
The signal-to-noise gain fiom replica correlation arises from a reduction in the
mean of the noise because, unlike the signal, it is uncorrelated with the replica
pulse shape. The DT equation assumes that the noise is completely uncorrelated
with the transmitted pulse and therefore the noise power is reduced by the number
of independent samples combined in the correlation process, which is BT Thus
IDT=510gd- 10logBTJ
Note that the gain is 10 log BT compared with 5 log BT for broadband passive
sonars. This is because replica correlation is a coherent process (complete know-
ledge of the signal) whereas broadband processes are incoherent (nothing known
of the signal).
+
where ( N - DI 10 log B ) is the in-beam noise over the h l l bandwidth and
(5 log d - 10 log BT - 5 log n) is the detection threshold. Here y1 is the number of
pings used by the operator or automatic system to make the decision; in a passive
system, n is the number of lines used.
Putting SE = 0 and combining terms:
CWpulse: If the analysis bandwidth is greater than the pulse bandwidth there will
be a loss (due to the increase in noise): If the analysis bandwidth is smaller than
the pulse bandwidth there will also be a loss (due to the loss of some of the signal).
The loss is given by a mismatch term:
either, - 10 log( B,/B,) when B,, > B, or, - 10 log( B,/B,) when B, < B,
LPFM Pulse: If the replica is extended (to allow for target Doppler) there will be
a constant loss (due to the increase in noise) for all dopplers of
166 Active Sonar
If the replica is not extended but there is a mismatch (an incomplete overlap) with
the signal return due to target doppler, there is a variable loss given by
where Bc is that part of the signal which overlaps the reference. Note that these
mismatch terms are always negative, so that - when added to the RHS of the
equation - they always reduce 2PL.
+
where (SL - 2PLR TSR) is the in-beam reverberation over the full bandwidth
and (5 log d - 10 log BT - 5 log n) is the detection threshold.
When the propagation paths for echo and reverberation are identical, 2PL
= 2PLR and these terms will cancel out.
TSR,the target strength for reverberation, is given by
The dominant reverberation source is usually at the sea surface, the sea bed or in a
horizontal layer within the water column (see Chapter 6) and therefore
The term 10 log( B / O h ) is known as the reverberation index (RI) and is a useful
measure for comparing performance against a background of reverberation.
For an FMpulse, B = 1/Tr, where T' is the resolved pulse length, or the time
between independent output samples from the correlator.
Example 9.1
CW pulse: T = 100 ms, 8, = 10"
RI = 1010g(10/10) = 0 dB
Example 9.2
FMpulse: B = 300 HZ 4 = 10"
The wider bandwidth of the FM pulse therefore results in a much improved perform-
ance against reverberation. Looking at it another way, the effective reverberating area
is determined by the bandwidth and not by the pulse length.
If the dominant factor is voEurne reverberation resulting fiom backscattering fiom that
part of the water column which is actually included in both transmit and receive beam
patterns, then the equation becomes
A wideband FM pulse which spreads the reverberation power over the band-
width, B, of the pulse. When the signal is correlated with a replica of the
transmitted pulse - effectively coherent matched filter processing - the pulse is
compressed into an equivalent time - the resolved pulse length, Tr, which is
equal to the inverse of B. The area of the surface contributing reverberation is
therefore reduced and performance against reverberation is improved as B is
increased, increasing B in 10 log( B/&). Furthermore, the performance of the
wideband pulse is independent of target doppler - provided there is still an
adequate match with a replica.
Time Time
The periodicity of the pulse changes linearly with time (i.e., its frequency changes
hyperbolically with time). Therefore the echoes, shifted in frequency by target
doppler, will always correlate with a suitably extended replica. The extended
replica will result in a constant loss, perhaps as much as 2 or 3 dB, for the entire
range of target doppler and it is often preferable to use a replica equal to the pulse
bandwidth and to accept a loss increasing with target doppler. No doppler
information can be extracted from this pulse.
Linear FM
The frequency of the pulse changes linearly with time. The echoes shifted by
target doppler will not now fully correlate with an extended replica and a set of
replicas (analogous to the comb of filters or FFT cells used to process CW echoes)
are now required to cover the expected target dopplers. The best replica match
indicates the target doppler.
170 Active Sonar
Target doppler moves the instant in time when h l l correlation of the echo and
replica takes place. There will only be an error when a single replica is used (the
LPFM case). To a first order, the time error is
At = (T/B)Af
cT
Range error = - -Af (note the active 2)
2 B
Example 9.3
If T = 500 ms and B = 250 Hz,what is the error for 10 knots of doppler at 10 kHz?
Af = 0.69 X 10 X 10 = 69 HZ
1500 0.5
Range error = - -X 69 = 104 m
2 250
The target doppler may be determined by range rate calculations from two or more
pings or from a simultaneous CW transmission.
cw Pulses 171
9.8 CW Pulses
A rectangular shaped pulse is defined in the time domain as
A sin2n: ft for f f j d T / 2
otherwise
and the frequency spectrum of the pulse, normalized to 1 s and 1 kHz, is shown in
Figure 9.3. The 3 dB bandwidth is given by
0 4 8 12 16
Frequency (Hz)
and again its frequency spectrum, normalized to 1 s and 1 kHz, is shown in Figure
9.4.
The 3 dB and 40 dB beamwidths are given by
-801 I I I I '
0 4 8 12 16
Frequency (Hz)
Note the significant reduction in sidelobe levels which is achieved, but at the
expense of a 50 per cent increase in the 3 dB bandwidth.
Figures 9.3 and 9.4 are obtained by analysis using a rectangular window. If
pulse and window have the same length, it is possible to exchange pulse and
window weighting functions without changing the final observed pulse spectrum.
In other words:
A shaped pulse analysed using a rectangular window has the same spectrum as
a rectangular pulse analysed using a shaped window.
CW Pzilses 173
Af3 =-
dB Of I
-70,
-80- I I I 1
(J
0 4 8 12 16
Frequency ( Hz)
Rj (a)
- 20
I s, 10 kH2 pulses
-40
-60
I I I I I I I
-1.2 -0.8 -0.4 0 0.4 0.8 1.2
Target doppler (knots)
For example, -40 dB of R j is reached for 0.4 knots of target doppler using a
1 s Hamming pulse and window, but this is never reached using a rectangular
pulse for realistic submarine dopplers.
Figure 9.6 portrays a stable and stationary platform in a stable environment. If
this is not the case, the background reverberation in the main lobe and sidelobes of
a beam will be modified as follows.
Reverberation and Target Echoes in the Main Lobe 175
It is only possible to fully nullie own doppler at one bearing within a beam. The
finite beamwidth has two effects:
A0 f O Zt0 0
$10.04 h0.6 -1-0.9 10
ct0.06 i0.9 il.4 15
-1-0.08 f1.2 -1-1.8 20
Ship motion
Roll and pitch have little effect on reverberation spread but contributions due to
yaw are significant. NWS 1000 (UK Naval Weapons Specifications Publication)
gives a figure for yaw in large frigates and destroyers as 1.75" s - ~peak. (Yaw is
unlikely to be important in submarines or for towed arrays, except perhaps for
very long pulses.) If the array is 60 m from the centre of gravity then the linear
acceleration at a bow dome is 2 m/s2 and, by inference from other NWS 1000
figures, the doppler spreads for 1 s and 250 ms pulses are as in Table 9.2. These
spreads occur both at transmit and receive.
176 Active Sonar
Environment
In shallow water, bottom reverberation can produce spreads of between 0.2 and
0.5 knots. In deep water, surface reverberation is important and depends on wind
speed. A useful empirical formula is
IDopplerspread = 0 . 0 7 ~ 0.31 +
with doppler and wind speed, v, in knots. And, substituting for the appropriate
wind speeds, we obtain Table 9.3.
Sea state 2 4 6
Doppler spread (knots) 1.2 1.6 2.3
Total reverberation
The total reverberation spread in the main lubes results from the pulse spectrum
plus the three factors we have looked at:
Beamwidth
Ship motion
Environment
v knots
t
Doppler Platform
Ahead 'y;
Beam 0 -+ -v
- 12 -- 8 -4 0 4 # 12
Target doppler (knots)
The performance of the two pulse types is then assessed for both full-beam and
half-beam processing and the following benefits are shown to result from half-
beam processing. The last two items are applicable to both FM and CW pulses.
9.14 Beamforming
Segments of circular or cylindrical arrays are effectively reduced to equivalent
linear or planar arrays by applying time delays to the elements, hence a short line
of 16 elements spaced A / 2 is a useful and simple array with which to demonstrate
the principles of full-beam and half-beam techniques.
Full beams are fonned by the simple addition of all 16 elements, after shading
and time delays, to form the required set of beams. (For a circular array the time
delays are constant and different sets of elements are used to step around the array
in azimuth.) The resultant full beamwidths will be of the order of 8" and spacing
between beams is, in practice, somewhat less than this, to minimize scalloping
losses. For simplicity, we will ignore scalloping loss and make the spacing equal
to the full beamwidth (8").
Pairs o f codirectional half-beams are formed using the 2 X 8 elements of the
line in a similar manner to the full-beam case. Both left and right half-beams now
have widths of 16" and the spacing between each pair is also 16". The beam shapes
are shown in Figure 9.9.
182 Active Sonar
half-beams
\ \
/
/ / \
120"
60"
-0"
- 60"
- 120"
Full beam
A \
f
Figure 9.9 Half-beam and full beam: amplitude and phase plots
FM Phase Binning Process I83
Phase bins
Select L/R full beam
Rig1
half.
16
Figure 9.10 LPFM signal processing for one pair of codirectional half-beams
The target appearing in one phase bin, or overlapping into two phase bins (and
here we will always assume two phase bins), is thus detected against a background
- noise or reverberation - reduced by the ratio of the half-beam width to the width
of two phase bins. The reverberation index is thus improved by 6 dB (8/2)
compared with the full-beam output that would otherwise be used.
Integration of the independent samples (every 1/B) from the phase binning
process will give a post-detection gain until the integration time exceeds the target
extent. If we assume the minimum target extent to be 15 m, we would sample and
reset the phase bins every 20 ms and
value, which is perhaps about 6 dB less, should be used. In view of the uncertainty
surrounding TS measurements, although a PD gain certainly occurs, it is perhaps
better to err on the side of caution and ignore it when making performance
calculations.
9.16 CW Processing
Each pair of half-beams is processed as shown in Figure 9.1 1. The half-beams are
frequency analysed to form N doppler channels. Left and right full beams are then
formed for each doppler channel by a similar bulk steer and add process. Phase
channel
Doppler
channel 2
Left
half-beam Doppler
channel
(filter) N
comparison is again used to select the correct full-beam outputs from each doppler
channel. These outputs are combined with the measured fine bearings (within a
beam) and passed to data processing and display at a rate depending on the CW
pulse length (e.g., every 250 ms for the 1 s pulse). There is no possibility of using
phase bins for post-detection integration because the length of the CW pulse is
already comparable with, or much greater than, the dimensions of the target. The
long CW pulse has effectively integrated the echo energy in the water and
therefore this time the integrated TS value should be used in the sonar equation.
However, there are still positive advantages in forming half-beams and deter-
mining target bearing by phase comparisons. The CW pulse will benefit from
improved bearing accuracy and the elimination of echoes from adjacent beams in
exactly the same manner as for the FM pulse, There is also an improvement in the
CW detection performance, i.e., it becomes noise limited, immediately the doppler
spectrum of the reverberation within the main lobe of a beam is exceeded, as
explained in Section 9.12
-33-
-36-
-39-
-42-
-45-
-48-
Figure 9.12 (a) Half-beam response and (b) phase difference plot
Detection Performance
The noise-limited performance is identical for both pulse types and is given by
2PL=210+10-50+12+0-10+3= 175dB
(The integrated TS value is used for both pulse types and, for the FM pulse, this is
justified by ignoring the PD gain due to the phase binning process.) If we assume
spherical spreading and absorption, the range is 9 km.
Using the parameter values from Table 9.4, the LPFM reverberation-limited
ranges are given in Table 9.5 for three representative values of s b .
Sb = -37 dB, 75 per cent of UK inshore waters have sb less than this value
Sb = -40 dB, a low bottom, high surface value for S b
188 Active Sonar
At useful detection ranges (greater than say 5 km), submarines never subtend
angles greater than I", hence any processing and display scheme based on full
beams - which are inevitably at least a few degrees wide - will have no fine detail
in the bearing dimension to help with detection and classification.
For noise and uniform reverberation, where bearing extents are similar to the
submarine's, detection is not improved by using half-beam processing. But the
greater bearing accuracy and elimination of sidelobe returns will simplifL and
improve all subsequent processes.
For discrete reverberation, where bearing and bearing extents can differ
markedly from those of a submarine, detection is improved by using half-beam
processing.
Full-beam processing
-30 -27 8
-37 -20 3
- 40 -17 2
Half-beam processsing
--30 -27 1.1
-37 -20 1 .o
-40 - 17 1 .o
dopplers which exceed the platform speed for targets on the beam, and twice
platform speed for targets ahead or astern.
Half-beam processing is superior to hll-beam processing of CW pulses.
Detections remain noise limited down to target dopplers of about 1 knot for a 1 s
pulse. There is little to gain - in performance against low dopplers - from using
longer pulses. This is because the factors other than pulse length which cause
reverberation spreading - finite beamwidth, platform motion and the environment
- then become dominant.
ttxi'(s.S+ a )
- 10 0 + 10 +20
Target strength (dB)
per cent. The comparison indicates the importance of attempting to ensure that the
sonar is always noise limited. Clearly, there will be constraints which will often
make this impossible, but it must be the underlying aim in the design of any active
system.
192 Active Sonar
The 3D representation shows the correlator output above a 2D surface (f, Sfo).
The intersection at the detection threshold P (half-power) defines an ambiguity
contour within which a target cannot be located unambiguously, since all (f, bfo)
combinations result in detections.
A plot of this contour versus Sfo and f is known as the ambiguity diagram for
the pulse waveform. The ambiguity diagram, therefore, indicates the accuracies
and resolutions in range and doppler which are achievable from a given pulse.
The treatment which follows has been stripped of the somewhat difficult
mathematics necessary for even approximate derivations of the ambiguity dia-
grams for CW, FM and PRN pulses and concentrates on results and practical
examples.
A rnbigui? Diagrums 193
CW pulses
Example 9.4
A sonar transmits a 4 kHz CW pulse of duration 0.5 s and is required to detect targets
with dopplers up to f 2 0 knots. What are the range and doppler resolutions of the
pulse and how many references (filters or FFT cells) are required?
We have
(sfo
2vfo
z --
0.88
z ___
C T
Therefore
0 . 8 8 ~ 0.88 X 3000
7Jz-=- = 0.7 knot
2Tfo 2 X 0.5 X 4000
194 Active Sonar
Note that here c must also be in knots, 3000. The doppler resolution (ambiguity) of
the pulse is 0.7 knot and therefore to cover f 2 0 knots requires 40/0.7 =
57 references.
We have
2R
f = 0.6T = -
C
Therefore
The range resolution is 226 m and we need to sample at least twice during this range
to ensure we capture the peaks of the signal, i.e., at least every 0.6 X 0.5 = 300 ms
Note that, for the simple CW pulse, the concept of ambiguity functions is not
necessary to determine the range and doppler resolutions. For a rectangular CW pulse
the 3 dB width is 0.91/T (compare with 0.88/T above), which will determine doppler
resolution, and we could expect to sample at least twice every pulse length, i.e., every
250 ms, or every 188 m (compare with 226 m above).
In linear period modulation (LPM) doppler changes the received signal in frequency
only, causing an overlap loss with the reference, and a slope mismatch does not
occur whatever the magnitude of BT. By equating the ambiguity fhction for the
LPM pulse to 0.5 (half-power) and putting first t = 0 and then dfo = 0, we obtain
The endpoints of the ambiguity diagram (Figure 9.16) are dfo = f 0 . 3 B and
t = f0.3T. And for any length of pulse, the area of the ellipse is approximately
0.6B(0.88/B)(n/4) = 0.4.
The doppler resolution (ambiguity) of the pulse is given by the maximum dfo
value of the diagram:
900 B
Sfo =-2vfo = 0.6B therefore v = - (knots)
C fo
Ambiguity Diagru m s
The range resolution of the pulse is given by the intersections of the contour with
the time axis:
0.88 - 2R 660
t=-- - - therefore R = __ (metres)
B c B
Example 9.5
A sonar transmits a 10 kHz LPM pulse of 500 Hz bandwidth and 1 s duration. What
are the range and doppler resolutions of the pulse? What is the maximum target
doppler for 3 dB correlation loss?
660 900 B
R=-= 1.32m u=-------
= 45 knots
B fo
Because the doppler resolution is derived from the 3 dB contour, then the maximum
target doppler for 3 dB loss is f 22.5 knots (half of the resolution). The small
correlation losses resulting from large target dopplers explain why the LPM pulse is
known as a doppler-invariant pulse.
Once again a simpler derivation is possible. The range resolution is simply the inverse
of the bandwidth, i.e., 2 ms or 1.5 m. There will be a 3 dB correlation loss when only
0.707 of the signal overlaps the reference, which is equivalent to a doppler shift of
f 1 4 7 Hz or f 147/6.9 = f 21 knots (from Af = f0.69 Hz per knot per kHz),
which is in good agreement with the previous result of f 22.5 knots.
196 Active Sonar
Note that, because of the assumption that slope mismatch is the limiting factor,
the equation for doppler resolution only gives valid results for large BT products.
For small BTproducts, a better result is obtained from the simpler derivation used
for the LPM pulse: calculating a doppler shift equivalent to 0.29 times the signal
bandwidth.
Example 9.6
A sonar transmits a 10 kHz LFM pulse of 500 Hz bandwidth and 1 s duration. What
are the range and doppler resolutions of the pulse? How many references are required
to limit correlation loss to 3 dB for target dopplers up to f 3 0 knots?
A mhiguity Diagrams 197
We have
R = 6 6 0 / B = " ~ 0 . 6m
6
ZI = 5200/BT = 5200/500 = 10.4 knots
The area of the contour is 2.4fo/TB2 = 0.1, which is smaller than the corresponding
area of 0.4 for the LPM pulse. The smaller the area, the smaller the ambiguity of the
pulse (the better the resolution).
As we have seen, the smaller area of the contour for the LFM pulse indicates an
improved doppler resolution - about 10 knots instead of 45 knots for the LPM pulse.
The smaller ambiguity of the LFM pulse explains why it is sometimes known as a
doppler-intolerant pulse. Note, however, that while it may be intolerant from the
viewpoint of limiting correlation loss, the pulse only provides a poor measure of the
target doppler. (Compare this with a CW pulse of the same duration and frequency,
which would have a doppler resolution of 0.14 knot.)
Example 9.7
A sonar transmits 10 kHz LPM and LFM pulses of bandwidth 100 Hz and pulse
duration 1 s. What are the doppler resolutions of the pulses?
This is clearly an invalid result; the LFM pulse is apparently more tolerant of doppler
than the LPM pulse of identical bandwidth and duration. A better result is now
obtained from the equation used for the LPM pulse. Both pulses, then, have the same
doppler resolution of 9 knots.
Example 9.7 highlights the problem of how to decide which equation to use. For
practical values of B and T, a valid result will be obtained using
For intermediate values of BT, or indeed for any BT value, a recommended, robust
approach is to use both equations and take the result of smaller magnitude.
198 Active Sonar
PRN p ulses
at least twice in a period equal to the inverse of the bandwidth (at least every 1 ms
for a bandwidth of 500 Hz), and possible complications or compromises necessary
in the display of the data.
'The processing power requirement is unlikely to be a problem with current
technology but, as always, availability of an adequate bank of PRN trials data is a
strong argument for retaining combined CW and FM sonars.
Example 9.8
A sonar transmits a 5 kHz PRN pulse of 400 Hz bandwidth and 1 s duration. What are
the range and doppler resolutions of the pulse? How often must the signal be sampled
and how many references are needed for a maximum target doppler of &20 knots'?
We have
- 0.88 X 3000
v=- 0 . 8 8 ~- = 0.26 knot
2Tfo 2 X 5000
To cover f 2 0 knots, we need 40/0.26 = 154 references, each spaced 0.26 knot and
0.26 knot wide.
R = 660/B = 1.65 m
Statistically independent signal samples are separated by the inverse of the bandwidth.
(Note that the approximate solution using the ambiguity fimction is 0.88/B, very close
to the inverse of the bandwidth.) Here the inverse of the bandwidth is 2.5 ms, and to
ensure that peaks of the signal are not missed, sampling should occur at least twice
during this period, say, every 1 ms.
This area is very much smaller than for similar duration CW or FM pulses and
confirms the low ambiguity (good resolution) in both range and doppler of this class
of pulse.
200 Active Sonar
Table 9.7 Range and doppler resolution for different pulse types
Pulse type Range resolution (m) Doppler resolution (knots)
cw 450T
LPM 750/ B
LFM (BT < 100) 750/B
LFM (BT > 500) 750/ B
PRN 750/ B
Table 9.7 collates expressions for the range and doppler resolutions for the pulse
types considered. For intermediate BT values of LFM pulses, and preferably for
any values, calculate the doppler resolution using the two alternative formulae, as
discussed earlier, and use the smaller of the two results.
Example 9.9
To improve noise-limited performance, the term 10 log T in the active, noise-limited
sonar equation must be increased. Suppose a 10 s pulse is transmitted, then
10 log T = 10 dB. But if we assume the correlation loss to be 6 dB, the effective gain
over a 1 s pulse is only 4 dB. If the long pulse is replaced by 5 pulses, each of 2 s and
separately processed, then 10 log T = 3 dB but this can be followed by post-detection
integration, giving a further 5 log 5 = 3.5 dB gain. The effective gain over a 1 s pulse
is now 6.5 dB, an improvement of 2.5 dB over transmitting the single long pulse.
Operational Degradation Factor 20 I
This typical sonar uses half-beam processing. The 1 ms correlator samples are
placed in phase bins and integrated to 3 ms for classification (3 ms is the inverse
of the bandwidth and the minimum time for independent samples). The 1 ms
samples are integrated to 20 ms for surveillance, i.e., roughly matching the
minimum target dimension. The FM surveillance display must therefore display
the integrated contents of the phase bins for each beam as a function of range
(time).
RANGE
5 12 pixels
BEARING
5 12 pixels
32 X 16 phase bins
Range scale 4 8 16
Range increments 267 533 1067
Pixels per increment 2 1 0.5
204 Active Sonar
For the 4 km range scale, display the same information in two consecutive lines.
For the 16 km range scale, use one line to display two consecutive range
increments. (The maximum ranges displayed will be slightly less than the nominal
range scales.)
C W surveillance display
The CW processing provides independent samples every 133 ms. This is equiva-
lent to a 100 m range increment, which is already large compared to the minimum
target dimension, hence it requires no integration before display.
5 12 pixels
1 I
BEARING, DOPPLER
5 12 pixels
Range scale 4 8 16
Range increments 40 80 160
Pixels per increment 12 6 3
Note that, because of the small bandwidths associated with the CW pulses used
for submarine detection, range resolution is poor and, for all range scales, the
same information is displayed in several lines if the display is to be the same size
as the FM display. Alternatively, these surplus lines could profitably be used to
display ping history. Whatever solution is adopted, equal height displays are
clearly desirable to allow the operator to readily associate FM and CW echoes
from the same target.
The maximum number of increments in the horizontal dimension (bearing +
doppler) is 32 X 60 = 1920. This number must be reduced by ORing beams and/
or doppler channels to match the available number of pixels, e.g., 16
+
beams X (30 doppler channels 2 spaces) = 512. The doppler channels may be
varied in width, increasing with the magnitude of the target doppler, to reflect both
the improved performance and the lower probability of high magnitude target
dopplers.
Active Displa-ys 205
The range, bearing and doppler samples, however, will be passed to the
automatic detection and classification processes at their full resolutions.
History displa-vs
Active history displays are analogous to the BT and LOFAR passive displays, where
successive time samples are replaced by successive ping samples. Figure 9.20
shows the simplest form of history display, used by the early electromechanical
range recorders. A slope equivalent to the platform relative velocity will indicate
reverberation or a zero doppler contact. Other slopes indicate non-zero doppler
contacts. Noise will seldom correlate from ping to ping. Derivatives of this display
are still used successfidly in modern sonars and are not lightly disregarded.
1 1 J \ Ping n
Pingn - 1
Ping n - 2
Pingn - 3
Range _______+
display which will be easier to detect and class@. To avoid display whiteout after
several several pings, the individual pixel amplitude (brightness) values of the
samples from earlier pings must be decreased with time.
9.25 Bandwidth
For detection we need to match the minimum target dimension (d). Suppose for a
submarine this is 10 m, then the resolvable pulse length is f = 13 ms from
-
d = c t / 2 ; as we have already seen, t can be much less than the actual pulse length.
Therefore bandwidth = l / t = 75 Hz. For classification we need to resolve the
structures of the target and any non-targets (false alarms). This will require a
higher resolution, to say 1 m,giving bandwidth = 750 Hz.
If the target is small, e.g., a mine, the required resolution will be much higher,
say 0.01 m, giving bandwidth = 75 kHz. We could use either a CW pulse of
duration 1/ B = 13p s or a longer, perhaps 10 ms, broadband (FM type) pulse.
Bandwidth 207
The point to emphasize here is that bandwidth is decided by the dimensions and
structure of the target and the false alarms. There is nothing to be gained by
increasing the bandwidth further. Note that if there are reflecting objects much
smaller than the target, they will be irrelevant because their target strengths will
be too small to compete with those of the target and Comparable false alarms.
208 Active Sonar
9.26 Beamwidth
Just as bandwidth determines the resolution in the range dimension, so does
beamwidth in the bearing dimension. Again, for detection we would like to match
the minimum target dimension, which is now the minimum angle subtended at the
array by the target. This is directly proportional to range, and for a submarine at a
range of at least 5 km this will always be less than about 1". Given practical arrays,
the smallest beamwidth will be at best about 4" and matching will only occur at
around 1-2 km, depending on submarine aspect. (Smaller beamwidths are possible
using towed arrays, but only at bearings within, say, normal f60".)
Figure 9.21(a) shows, in B scan format, submarine targets at two different
ranges and a patch of discrete reverberation which might, for example, represent a
rocky outcrop in shallow water. Only one beam of the display is shown: the
beamwidth is 4" and the range is from 0 to 5 km. Note that the width of the display
is proportional to range: at 5 km the width is about 270 m and therefore the
submarine at this range has a bearing spread of about lo, but at 1.7 km the width is
90 m and its bearing spread is much greater, about 3".
Figure 9.21 Comparison of bearings and bearing spreads for submarines and false alarms
Figure 9.21(b) shows one beam of a sonar display in the same B scan format
and produced using full-beam processing. Because fU-beam processing cannot
resolve bearing within a beam, the targets and the false alarm will all be given the
same bearing and no values at all for bearing spread. The submarine targets and
the reverberation patch, therefore, all appear at the same bearing and with the
same (unknown) spread in bearing. Note, however, that the differences in range
spread are displayed.
Beamwidth 209
Figure 9.21(c) shows one beam of a sonar display again in the same format but
produced using half-beam processing. We can now resolve bearing within a beam,
and the submarine targets and the reverberation patch are all displayed with their
correct bearings and bearing spreads. We now have, within the sonar's limits of
resolution, a fair approximation to the overall shape and position of contacts and
false alarms,
At 5 km range, the submarine is in one phase bin whereas the reverberation
occupies four phase bins. Since, at this range, a submarine cannot spread into
more than two bins, the bearing spread is a useful discriminant between targets
and f'alse alarms. At closer ranges, this discriminant may not be so useful. The
target spread, however, when considered together with other parameters such as
target motion, will still be a strong classification feature. For classification we
again need to resolve the structures of the target and any non-targets (false alarnis)
in the bearing dimension. Discrete reverberation (clutter) can have much greater
bearing spreads than submarines, and in order to exploit this difference for
classification (false alarm reduction), we require a higher bearing resolution up to.
say, 1/10 of the beamwidth. A practical minimum value for the composite
detection and classification task might be 0.5".
Full beams, using the greatest possible aperture of the array, cannot achieve
such small beamwidths, therefore detection and classification processing using
signals from full beams is not optimal against a background of high discrete
reverberation (clutter).
However, against a background of noise or uniform reverberation and only if
the other important advantages of the alternative, half-beam processing, are
ignored, then full-beam processing, which uses the highest possible array DI, is
marginally better: the 3 dB greater DI provides an extra 1.5 dB of allowable
propagation loss to use in the sonar equations.
In summary, half-beam processing can be used to effectively achieve small
beamwidths and therefore detection and classification processing using signals
from half-beams is optimal against a background of high discrete reverberation
(clutter). The effective beamwidth is the size of a phase or bearing bin.
The phase bins, as recipients of the output signals from the process, resolve the
structures of both targets and false alarms in both range and bearing dimensions,
to provide the best possible data for the composite detection and classification
task.
210 Active Sonar
9.27 CADAC
The output from the signal detection process is a time series of samples from
every beam of the sonar, and for modem high resolution sonars there are many
samples every ping. Suppose pd = 0.5 and Pfa = (this is for the simple
signal detection process only). For the typical sonar which was used to illustrate
active displays, the FM display has 512 X 500 = 2.5 X 105 pixels (samples).
Therefore there are 2.5 false alarms every ping, or 15 false alarms per minute!
What appears to be a very low Pfa has become an intolerably high rate of false
alarms. The post signal detection processes of CADAC (computer-aided detection
and classification) and display will help the operator to reduce this to an
acceptable false alarm rate to pass to the command system.
Levels qf'CADAC 21 1
Level 1
Level 1 is the lowest level; it detects signals within single cells; the signals may
originate from vessels, noise or reverberation. A fixed threshold is normally used,
at around 10- 13 dB above the average of the background (the threshold may be
lowered to dig deeper into the background in areas of interest, or increased to
inhibit inputs from improbable areas).
Level 2
If the sonar resolution is finer than the typical size of a contact (including false alarms
such as discrete reverberation) then returns will occur in several adjacent range and
bearing cells. Clustering these returns together forms the second level of classifica-
tion, since some estimate of the dimensions and shape of the contact can be formed
from the distribution of returned energy in the contiguous cells. This estimate is
commonly known as the individual weight of evidence of a contact. The parameters
used to make this estimate include cell amplitudes, total energy, range, range spread
and bearing spread. At level 2, additional processes can be applied to the cluster of
individual returns making up a contact. There are, broadly speaking, two ways of
doing this: statistical analysis and amplitude profiles (Sections 9.30 and 9.3 1).
Level 3
Level 3 examines how well contacts associate from one ping to another by seeking
consistent tracks. Goodness of track and individual weights of evidence of the
212 Active Sonar
contacts forming the track help to build up a parameter known as the weight of
evidence (WOE)of a contact. Courses and speeds of all contacts are also calculated
at this level (the speed of a contact is clearly an important classification clue).
Target doppler, if available, should be input to the WOEat this level.
Level 4
At level 4, information gathered from all other possible sources is used to assist in
the final classification of all contacts. Here are some possible sources:
Sonobuoys
Radar
Many of these sources will not help with the classification of the submarine itself.
They will, however, materially assist in the reduction of false alarms. The
reduction of false alarms can itself be viewed as classification. If all false alarms
are eliminated, surviving contacts must be submarines.
Level 5
Level 5, the final level, includes all the experience, knowledge of tactics, strategy,
and intelligence data possessed by the operator and the command, and it is only at
this, human, level that the classification can be positively confirmed. For this
reason, any automated classification system, no matter how clever its algorithms,
should be thought of as a classification aid not an automatic classifier.
echoes fall into regions of the reverberation spectrum where the reverberation
energy is low. The important classification clue of target doppler is therefore
always available simultaneously with CW detections.
For the FM pulses, the range resolution cell is small compared with the typical
extent of a contact and therefore it is possible to make improved estimates of the
dimensions and echo profile of the contact. Furthermore, half-beam processing, by
improving bearing resolution, provides additional information (bearing spread)
towards this estimate, which can help to discriminate between targets and false
alarnis at CADAC levels 1 and 2.
For both pulse types, half-beam processing improves bearing accuracy and
eliminates echo returns in adjacent beams, from which ping-to-ping association
and tracking processes will benefit (CADAC level 3).
which, after normalizing and clean-up processes, are compared with memorized
profiles of known contacts to give an indication of the contacts type.
The method has the advantage that it can learn new profiles from just a single
exposure to a particular contact, and can dynamically adjust its repertoire of stored
profiles in accordance with the prevailing conditions. Furthermore, it has the
important attribute of being able to indicate to the operator why it has made its
hypothetical choice between alternative classifications in terms that the operator
can understand, rather than as obscure statistical technicalities.
There is a case for providing both statistical analysis and amplitude profile aids to
classification. The statistical analysis classifier is probably superior in classifying
those targets with which it has been trained. But it is not robust to changes in proc-
essing, it is sensitive to the environment (e.g., using different algorithms for deep and
shallow water), and from an operator viewpoint it makes its decisions in the dark.
The amplitude profile classifier is perhaps somewhat less than optimum when
operating against the specific targets with which the statistical analysis classifier has
been trained. On the other hand, it is robust to changes in signal processing, it is
insensitive to the environment, it can rapidly learn new profiles from just a single
exposure to a new target, and because the profiles are displayed to the operator, it
makes its classification in terms that the operator can see and understand.
0 Point a on the target submarine returns signals to the source via two paths, A
and AA, for this simplified case (Figure 9.22). The same point therefore puts
returns in two separated cells along the profile. Signals along paths AA and B
arrive simultaneously but from different points ( a and h) on the target, and are
detected against a background of reverberation returning along path C and also
arriving at the same time. Different point sources therefore put returns in the
same cell of the profile.
An echo from point b returned via the bottom will extend the length of the
profile. The magnitude of this extension will increase with range - as path
differences and the number of reflections from the boundaries increase.
Multipath propagation increases the total energy returning to the source from the
target without increasing the reverberation background, and therefore detection is
enhanced. On the other hand, classification or at least those features used by
classification which rely on the fidelity of the profile, will suffer to an extent that
depends on the energies of the misplaced returns.
The likely magnitudes of possible profile extensions and highlight displace-
ments will now be estimated for the most relevant propagation modes, i.e., surface
duct in deep water and in shallow water.
216 Active Sonar
hd = duct depth
c = speed of sound
g = velocity gradient
112
Skip distance = [8 X (E)
0.017
X 1001 = 8400 m
From simple geometry, the extra path length over the direct path for the limiting
ray for each surface reflection is about 2.5 m (Figure 9.23). The limiting ray gives
the maximum difference as 0.3 m per km. So if the range of the target submarine
is 20 km, for example, its extent could be increased by up to 20 X 0.3 = 6 m (and
highlights misplaced by up to similar amounts along the profile). Given a large
submarine (at least 5 0 m long, say) the extension is unlikely to seriously affect
detection and echo shape recognition, but classification could suffer if it relies on
profile details.
t
Duct
Suppose 8 = 10" and h = 200 m. The distance between bottom reflections will
be 2268 m. For the two bottom reflections in Figure 9.25, R = 3403 m and the
indirect path is 3455 m. The path difference is 52 m. This is already quite large
and will increase linearly with the range, R. Bottom losses at each reflection,
however, reduce the amplitudes of the multipath arrivals. (Surface reflection losses
are much smaller and are neglected.)
Table 9.10 shows differences in path lengths for different angles at representa-
tive ranges of around 11 km and 23 km. Note that the large differences are
accompanied by large total reflection losses and can therefore be neglected.
Furthermore, bottom reflection losses increase rapidly at grazing angles greater
than about lo", so even if vertical half-beamwidths are greater than lo", there
should be no need to consider rays at angles greater than this.
9.37 Problems
9.1 A shipborne sonar transmits a pulse at 5 kHz and receives an echo at 5.025 kHz
from a target bearing G45". If the ship is moving at 8 knots, what is the doppler, or
relative velocity, of the target?
9.2 An LPFM pulse has a bandwidth of 400Hz centred on 4000Hz. Given own
doppler nullification, what are the target dopplers for processing losses of 1 and 3 dB?
9.3 An active sonar uses a Hamming shaped CW pulse to improve the reverberation
rejection, Rj, due to target doppler. How much target doppler would be needed to
achieve X, = -40 dB for a pulse of 1 s duration at 3000 Hz, 1000 Hz and 300 Hz?
This page intentionally left blank
10
Echo Sounding and Side Scan
Sonars
is on the vertical axis, and time, or distance travelled over the sea bed, is on the
horizontal axis.
The echogram in Figure 10.1 displays some 30 ping returns and therefore, at a
ships speed of 10 knots (5m/s) the distance travelled to build up the echogram is
about 40 m. The school of fish can be seen to extend about 5 m in plan, and about
40 m in depth. The brightness or colour of the display will indicate the strength of
the echoes.
There are no fundamental differences between the design of echo sounding and
side scan sonars and the design of a submarine detection sonar. The differences
are only of scale - frequency, range and target size. Because we are here
concerned only with short-range or vertical transmissions, refraction of sound,
which severely limits the performance of long-range antisubmarine sonars, is
unimportant and the propagation losses are well described by spherical spreading
plus absorption.
Echo Sounder: Design Example 223
The propagation loss will be given by spherical spreading plus absorption (quite
accurately in this instance, due to the lack of boundaries, rather than the usual
working approximation).
A solid angle beamwidth of 5" will result from a square array of side L, or a
224 Echo Sounding and Side Scan Sonars
DI = 2010g(F) = 2010g(~X030;:8) = 26 dB
The pulse length would have to be very short in order to resolve individual fish.
Suppose the individual fish to be separated by 0.1 m, the pulse length would need
to be
This may be too short to meet the propagation loss requirement but let us stay with
it for the present. Let the background noise, N, be the ambient noise at SS2, or
35dBat40kHz.Let5logd= 10dBand5logn=3dB.
Then to solve the noise-limited active sonar equation for the required source
level, SL, all we require is TS. The sea bed is a very strong target at normal
incidence and will be quite unmistakable. We will arbitrarily assign TS = 10 dB
to the sea bed. Fish, individually or in schools, will have very much smaller target
strengths, say -50 dB, but the depths will be less, say 100 m where PL = 4 1 dB.
The noise-limited active sonar equation is
SL = 118 - 10 + 35 - 26 + 39 + 10 - 3 = 163 dB
SL = 82 + 50 + 35 - 26 + 39 + 10 - 3 = 187 dB
SL,,, is therefore 187 dB. Is this practical? We have
Therefore
lOlOgP=187-171 - 2 6 z - l O d B
Figure 10.2 'Vertical' beams: combined echo sounder and side scan sonar
As the transmissions from the side scan transducer glance across the sea bed, each
ping builds up a single line on the display (Figure 10.3). As the ship moves,
therefore, successive pings build up a two-dimensional map of the sea bed. If the
bottom is smooth, the display will simply show a characterless, noise-like picture.
If, however, the sea bed has features, such as peaks and valleys, the picture will be
quite different; the peaks will backscatter strongly and the valleys, shielded by the
peaks, will display as shadows. Objects on the sea bed - wrecks, mines and
bottomed submarines - are frequently detected by their shadows. This mechan-
ism, the so-called shadowgraph eBect, is very important for the detection and
classification of bottom mines. Because of their very low target strengths, echoes
from mines are very difficult to detect against a background of bottom reverbera-
tion; the shadowgraph can, however, indicate the presence of a mine even when its
echo is hidden by the reverberation.
250 tK Ship's
track
I I
I I
I I
I I
I 1
I I
I 50m
4----'1
I I
I I
I I
2
Range (&)
A true geographic representation of the sea bed or a target requires that the
scales of both display axes shall be the same. This is achieved by choosing a
platform speed, V, which will match the track dimension to the maximum range
dimension of the selected range scale:
where
But because the horizontal beamwidth is finite, albeit small, the beam patterns
from successive pings could overlap significantly (depending on PRI and ship's
speed) and therefore the data on successive lines may not be independent, except
at close range.
Example 10.1
Suppose the display has 250 lines, each one displaying a ping's worth of data, then
The pings or lines are therefore spaced by 1 m in the track dimension. Since the linear
beamwidth increases with range - from 0 at zero range to 9 m at 250 m range for the
2" wide beam in Figure 10.3 - the beam patterns from successive pings overlap
significantly; at 250 m, for example, the data from 8, i.e., 9 - 1, successive pings fully
overlaps, and completely independent data is only available every 9th ping.
A compromise, avoiding the storage and display of redundant data, might be to use the
data from every 3rd ping, thus ensuring independent data at about one-third of
maximum range and reasonable overlaps at greater ranges. The data from every 3rd
ping would then be displayed on three consecutive lines.
Side Scan Sonar: Design Example 229
Transmit and receive on separate arrays of identical size and tilted downwards by
45". To achieve the required combined beamwidths, each array must have double
the combined beamwidths and the dimensions are given by L = 76/f&. There-
fore, the width is 76/(100 X 2) = 0.38 m and the height is 76/(100 X 90)
230 Echo Sounding and Side Scan Sonars
= 0.0084 m. Hence each array will have two rows of 50 elements, all spaced by
112.
DI, in both transmit and receive, is given by
DI = 10 log(4Lh/A2)
Baffling will add 3 dB to this, making DI = 21 dB. Once again, the pulse length
would have to be very short in order to resolve individual fish. Suppose the
individual fish to be separated by at least 0.1 m; the pulse length would need to be
2 x m - 2XO.l
T=- - -= 133
C 1500
This is probably shorter than necessary to map the sea bed or to display
submarines and wrecks but it may well meet the propagation loss requirement, so
let us stay with it for the present.
Let the background noise, N, be the ambient noise at SS4, equal to 30 dB at
100 kHz, and let 5 log d = 10 dB and 5 log n = 3 dB. Then to solve the noise-
limited active sonar equation for the required source level, SL, all we require is
TS. Backscattering from the sea bed at high grazing angles will be strong and its
target strength will be given by
where
cT
A = -R(3h ((3h in radians)
2
= (750 X 0.000 133)(200 X 0.017)
= 0.34 m2
At shorter ranges, the reverberating area will be smaller but the backscattering
strength will be greater (higher grazing angle), therefore the TSR may not change
Problem 23 1
much. This is very much smaller than the TS of 10 dB used in the echo sounder
design, but remember that was for a normal incidence reflection rather than
backscattering. Again fish, individually or in schools, will have small target
strengths, say -50 dB.
The noise-limited active sonar equation is
SL = 106 + 40 + 30 - 21 + 39 + 10 - 3 = 201 dB
and to detect fish out to 200 m range,
SL = 106 + 50 + 30 - 21 + 39 + 10 - 3 = 211 dB
SL,,, is therefore 21 1 dB. Is this practical? We have
I S L = lOlogP+ 170.8+DIt(
Therefore
10logP = 21 1 - 171 - 21 = 19 dB
giving P = 79 W.
The acoustic power intensity is 79/(0.38 X 0.084) = 2.5 kW/m2 and again
cavitation will not be a problem. The array design is therefore practical and the
remaining parameters may be inserted in Table 10.2.
10.6 Problem
10.1 An echo sounder is to be used to survey the deep ocean floor. The depth
requirement is 12 km to cope with the deepest ocean trenches. Given similar
parameters to Section 10.3 (but take SL = 2 10 dB, T = 10 ms and TS of the deepest
ocean floor as 40 dB), what is the highest possible operating frequency? What range is
achievable at 30 kHz?
This page intentionally left blank
11
Mine Hunting Sonars
1 1.1 Overview
Side scan sonars, operating at similar frequencies to the example in Chapter 10,
would appear to have the necessary performance to be used for mine hunting. The
short pulse lengths and narrow beams would provide sufficient discrimination in
the range and bearing dimensions to detect and classify mines. Operationally,
however, it is clearly preferable to detect a mine ahead of the ship - before it
becomes a danger - and the time taken to search an area using side scanning
techniques is a major liability.
A mine hunting sonar, then, is usually very similar to a hull-mounted anti-
submarine sonar. The array will be housed in a keel or bow dome and, in
surveillance, will have a wide arc of cover centred on the ahead bearing.
Classification of mine-like contacts is difficult, particularly bottom or close
tethered mines, where the background to detection and classification is frequently
high discrete reverberation (clutter) &om the sea bed. Target strengths are low and
therefore, to discriminate against this reverberation, narrow beamwidths and short
(actual or resolved) pulse lengths are required. The short pulse length may be
either a CW pulse, typically less than 1 ms, or a broadband FM pulse, in which
case the actual pulse length can be much longer. This will improve noise-limited
performance, but now the bandwidth can also be quite large; at least 10 kHz is
desirable, which gives a resolved pulse length of 0.1 ms.
Modern mines are very difficult acoustic targets; they are shaped and clad to
ensure very low target strengths. In extremis, detection and classification may only
be possible by closely approaching the mine using a remotely operated vehicle
(ROV) housing a very similar array.
234 Mine Hunting Sonars
Bottom or close tethered this class of mine must be detected against a background
of noise and unavoidable bottom reverberation. It is by far the most difficult of
the two and definite classification may sometimes only be possible using an ROY
A mine hunting sonar, therefore, will benefit from having more than one fan of
beams available. A versatile design would include several fans, all or some of
which should be displayed simultaneously. Alternatively, split-beam processing in
the vertical plane may be used to achieve the same result.
beam has a linear width of 1.7 m, acceptable for detection perhaps but too large to
resolve structure for classification. The use of a broadband FM pulse and phase
binning will improve the resolution; it is perhaps unrealistic to expect any better
than 0.25" from the phase binning, which would give about 0.5 m resolution at
100 m and 0.05 m at 10 m. Hence the suggestion that an ROY which can get close,
might be needed for the final classification of the most difficult targets.
PL =20logvfar X (dB)
Transmit array
To cover all depths and ranges, the transmit array must be capable of being steered
in the vertical plane. It must therefore have a number of separate elements. The
required number of elements will be given by 8 = 100/n. The transmit array is
therefore a single vertical stave of 20 elements at spacing L / 2 . To reduce the
maximum steer needed, both transmit and receive arrays will be mounted at an
angle of depression of perhaps 30". Wavelength L is only 7.5 mm and therefore the
radiating surface is about 75 X 4 = 300 X m2. To avoid cavitation the maxi-
mum power (from Figure 1.2) is 2 kW/m2.
Because of the high frequency and short pulses, this may safely be increased to,
say, 10 kW/m2 or P = 3 W. The array will be baffled and therefore DIt =
3 + lOlog20 = 16 dB and
Receive array
Again, to cover all depths, azimuths and ranges, the receive array must be capable
of being steered in both the vertical and horizontal planes. The beamformer will
produce fans of beams to meet the requirements of surveillance (the initial sweep)
followed by a more detailed look with narrower beams for classification. The array
must be large enough to form the narrowest beams; wider beams will be formed
by subsets of elements or by combining beam outputs. The required number of
elements will be given by 8 = 100/n, and since the narrowest beams are 1" in
both dimensions, n = 100. The receive array is therefore a square of side 0.75 m.
Again, the array will be baffled and therefore
Pulses
The major problem for a mine hunting sonar is cZassiJication, particularly for
mines on the sea bed, which have to be detected and classified against a back-
ground of bottom reverberation. Therefore the pulses must have sufficient
bandwidth to be able to discriminate against discrete reverberation (clutter). We
will propose two pulse types and compare their performances:
FMpuZse: bandwidth 80 kHz, duration 10 ms. The pulse needs to be quite short
to avoid excessive dead range; a 10 ms pulse will give a dead range of 7.5 m.
The target will not have any doppler, hence a linear FM pulse will suffice. (Note
that the dead range is not absolute; as soon as the system goes into receive
mode, processing of the pulse begins but the full gain is not achieved until a
complete pulse is processed, i.e., 10 ms after the start of the receive period.)
11.8 Performance
The noise-limited active sonar equation is
0 + +
FMpulse: 2PL = 192 - 30 - 30 43 - 20 - 10 3 = 148 dB therefore PL =
74 dB. Using this pulse, the requirement for PL = 65 dB is easily met.
0 + +
CWpulse: 2PL = 192 - 30 - 30 43 - 39 - 10 3 = 129 dB therefore PL =
65 dB. Using this pulse, the requirement for PL = 65 dB is now onlyjust met.
* F-M pulse: 1 0 l o g R = 4 9 + 3 0 - 1 1 - 3 0 - 1 0 + 3 = 3 1 d B
therefore R =
1300 ni. The reverberation-limited range requirement is easily met with this
pulse.
Note that a shorter CW pulse, about 50 ps, would meet the reverberation-limited
range requirement, but not the noise-limited range requirement. More importantly,
however, even a 50 ps duration CW pulse has poorer range resolution than the FM
pulse. We therefore conclude:
The best solution for the composite detection and classification task is to use the
FM pulse.
If this conclusion appears to conflict with some in-service mine hunting sonars, it
is because of technological limitations when they were designed.
240 Mine Hunting Sonars
11.9 Classification
Classification demands the best possible discrimination in both range and azimuth.
In the range dimension, the resolved pulse length for the FM pulse is 12.5 p,
giving a runge resolution of about 9 mm. This is, of course, very small compared
to the dimensions of a mine and there will be many returns in the range dimension
to help with classification.
In azimuth, the linear resolution is proportional to range - at 100 m the
resolution is 1.7 m for the 1" wide receive beams. Only at very close ranges is the
resolution in azimuth comparable with the range resolution. Phase binning can
improve this, but realistically not more than a factor of 4 (0.25"). So the best we
can expect at 100 m is an azimuth resohtion of about 0.5 m. This resolution,
although quite large compared with the dimensions of most mines, means that,
dependent upon aspect, several returns from different bearings are still possible
from larger mines every ping.
An important aid to the classification of mines is the acoustic shadow projected
onto the sea bed by bottom objects. The effect is illustrated in Figure 11.2. At the
beginning of the receive period, the background is just noise. When the lower edge
of the vertical beam intercepts the sea bed, reverberation adds to the noise as the
pulse sweeps along the sea bed. Once the pulse has reached the mine then echoes
- which may not exceed the background - will be added to the background. The
The hypothetical sonar used for comparisons between full- and half-beam proces-
sing will be used again here. Because we need the smallest possible beamwidths,
phase binning is used and a CW pulse will not be considered. The bandwidth
available from the omnidirectional projector is assumed to be half an octave, or
5 kHz. Table 11.2 shows the parameters, modified where appropriate (modified
parameters are in italics).
Detection performance
If we assume spherical spreading and absorption, the range (from Figure 3.3) is
1000 m.
The reverberation-limited performance is given by
1010gR=37+3O- 11 - 3 0 - 1 0 + 3 = 19dB
The range necessary for reliable mine avoidance is unlikely to be achieved at the
very high frequencies used for mine hunting. The frequency chosen must, how-
ever, be high enough to provide sufficient bandwidth for a good reverberation-
limited performance. We have already seen that an FM pulse gives superior
performance in all environments. (A CW pulse can only be better given target
doppler, which is not normally the case with mines.)
Let us begin with an outline specification (italic in Table 11.3), for a mine
avoidance sonar to combat the threat posed by a mine of TS = -30 dB, and then
build up a more complete specification. The arrays are hull mounted and the
maximum range requirement is 1000 m. The operating fiequency is 80 kHz and
the bandwidth 40 kHz (half an octave).
Mine Avoidance Sonars 245
horizontal 2"
PL 83 dB at 1000 m
R1 (= 10 log(40 000/2) 43 dB
Directivity index, Dl
transmit 10 dB
receive 45 dB
Source leve, SL 193 dB
Pulse length, T 100 ms
10 log T -10 dB
Ambient sea noise, N 35 dB
TS -30 dB
5 log d 10 dB
5 log n 3 dB
The propagation loss is calculated using spherical spreading plus absorption (not
very accurately because the sea surface and the sea bed will present boundaries for
much of the ping).
PL = 201og Y + aY X (dB)
Y = 1000 m, a = 23 dB/km therefore PL = 83 dB
Transmit array
The height of the transmit array is given by 0 = 20" = 1 OO/n. The transmit array
is therefore a single vertical stave of 5 elements at spacing d/2 and mounted at an
246 Mine Hunting Sonars
Receive array
To cover all depths and azimuths, the receive array is an incomplete cylinder made
up of a number of staves spaced at 1/2 (Figure 11.3). Each stave will have
100/2 = 50 elements. The diameter of the array is given by &, = 88/df. There-
fore diameter d = 88/(2 X 80) = 550 mm. Each full beam, or codirectional pair
of half-beams, is produced using 120" of the array. To provide 90" of azimuth
cover, the staves must occupy 240" around the periphery of the cylinder. There will
therefore be a total of (nX 550 X 2)/(3 X 9.375) = 120 staves. The ahead beam,
for example, will use the 60 staves in segments B and C. Adjacent beams are
formed by stepping around the array by one stave (the staves have an angular
spacing of 2"). Several fans of beams in the vertical dimension will be formed to
cover the required depths and range. The near-horizontal fans will avoid bottom
reverberation until quite long ranges, well beyond avoidance range.
The directivity index of the array is given by the formula for a cylindrical array
used previously:
3
Ahead beam uses the 60 staves in segments BandC
FM pulse
The FM pulse has bandwidth 40 kHz and duration 100 ms. The pulse needs to be
quite short to avoid excessive dead range; a 100 ms pulse will give a dead range
of 75 m.The target will not have any doppler, hence a linear FM pulse will do.
Detection performance
The noise level is SS4; if we had taken SS2, the noise level would have been
30 dB.
The pulse length could be increased (but watch the dead range); increasing T to
300 ms would increase PL by 2.5 dB.
Classification
Classification demands the best possible discrimination in both range and azimuth.
In the range dimension, the resolved pulse length for the FM pulse is 25 p,giving
a range resolution of about 18 mrn. This is very small compared to the dimensions
of a mine, and there will be many returns in the range dimension to help with
classification.
In azimuth the linear resolution is proportional to range. Classification can be at a
closer range than initial detection; if we assume 300 m, the resolution is 10 m for the
2" wide receive beams. The resolution in azimuth is never remotely comparable with
the range resolution. Phase binning can improve this, but realistically not more than
a factor of 8 (0.25"). So the best we can expect at 300 m is an azimuth resolution of
about 1.3 m. This resolution is comparable with the dimensions of most mines,
therefore the majority of returns from a mine will be at the same bearing and a
complete picture of a mine will not be available. However, even this resolution is
useful for classification because it will discriminate against large discrete patches of
bottom reverberation, which can extend over much greater dimensions.
11.12 Problems
1 1.1 How would the mine hunting sonar of Section 1 1.7 perform if TS were -25 dB?
11.2 The mine avoidance sonar design of Section 11.1 1.2 took the PL at the centre
frequency of the pulse (80 kHz). At the highest frequency of the pulse, what SL is
necessary to achieve a range of 800m? And assuming SL does not change with
frequency over the bandwidth of the pulse, what equalization is needed to maintain a
flat response at the receiver?
12
Intercept and Communications
Sonars
Carrier frequency
Bandwidth
Pulse type (FM, CW, PRN, etc.)
Pulse duration
Pulse shape
Pulse interval
Multipath structure
With the exception of multipath structure, all serve to help identify the transmitter
and hence the likely platform.
The carrier frequency may be between, say, 1-100 Wz. A typical receiver to
cover this spectrum might do its surveillance processing in octave bands
(1 -2 kHz, 2-4 kHz, etc.,), but should also have the capability to vary the exact
bands given some a priori knowledge of likely threats.
Bandwidths may be between, say, 1 Hz (a long CW pulse) and 1000 Hz (a short
CW pulse or broadband FM or PRN pulses). The intercept sonar should at least
Function of an Intercept Sonar 2s 1
measure start and finish frequencies and ideally it would completely analyse the
amplitude and frequency structure of the pulse.
Pulse interval, and any change in the interval, are important clues to the likely
range of the intercepted sonar. Suppose pulse interval = 15 s (corresponding to a
maximum range of about 10 km) the transmitter will be within this range,
assuming it is in contact. If the pulse interval changes to say 8 s, this is a strong
indication that the transmitting platform is closing range and is within 5 lun of the
intercept sonar.
Time delays between direct and bottom bounce arrivals in deep water can also
indicate range, but again only very broadly given that the simple broadband
intercept receiver is unlikely to be able to make a good measure of the vertical
angle of arrival. The large planar or conformal arrays fitted to submarines,
however, can measure vertical angles of arrival to sufficient accuracy to give a
good estimate of range, either from time delays between direct and indirect
arrivals of the pulse, or from the vertical angle of maximum response (compare
with VDPR). Given a knowledge of the original pulse length of a transmission, it
may be possible to infer the range of its source from any stretching of the pulse
resulting from multipath. The range estimate may be no more than that if the pulse
is stretched then the source is distant, and if the pulse is not stretched then the
source is close.
252 Intercept and Communications Sonars
SE = (SL - PL) - N + DI - DT
This is the basic passive sonar equation: SL is no longer the radiated noise of a
target but is the source level of the active transmission.
Because nothing can be assumed about the signal to be detected, the detection
process for the simplest intercept sonar is based on the received signal-to-noise
ratio without any processing gain, therefore
Both signal and noise are referred to the same bandwidth, hence the term 10 log Br
where B, is the bandwidth of the receiver.
Combining these equations and putting SE = 0, we obtain the intercept sonar
equation:
Worked Examples 253
PL = 2010g Y + 0.9~
X = 138 dB
Intercept range = 50 km
How does this compare with the detection range from the ship?
2PL=220+10-60+20+0- 10+3=183dB
PL = 20 log Y + 0 . 9 ~
X = 92 dB
Detection range = 12 km
Therefore the presence of the ship can be known to the submarine long before it
becomes an actual threat.
Example 12.2
An intercept receiver mounted on a submarine processing an octave band from 30 to
60 W z detects a signal of bandwidth 1000 Hz centred on 50 kHz with a ping interval
of 0.5 s. What are the probable source and range of the signal?
Assume DI = 3 dB and N = 30 dB. The frequency and ping interval indicate that the
source is probably a torpedo in its active homing phase. The SL may therefore be
about 200 dB.
254 Intercept and Communications Sonars
PL = 20logr + 16r X = 83 dB
But because the ping interval is only 0.5 s (corresponding to a range scale of about
400 m), the submarine may well infer that the torpedo is well within this maximum
range.
e Reducing SL
Increasing frequency to limit range by absorption losses
e Increasing the bandwidth
Using sector transmissions to confine the energy within azimuths of interest
If increased bandwidth is achieved by reducing the pulse length, as with CW, then
detection performance will also suffer. If it is achieved without reducing the pulse
length, detection performance will not suffer, but since the intercept receiver
bandwidth is likely to be larger still, the performance of the intercept sonar will
not suffer either.
Clearly, with the possible exception of sector transmissions, all these measures
severely prejudice the detection performance of an active sonar. But a simple
example will suffice to show that LPI may not be usefully improved even by the
use of sector transmissions.
Example 12.3
An intercept receiver mounted on a submarine processing an octave band from 4 to
8 kHz detects a signal centred on 5 kHz from a sonar with a known sector SL of
230 dB. Given spherical spreading and absorption, what is the maximum intercept
range of the sector transmission and the reduced intercept range for bearings outside
the sector, for transmission sidelobes of -30 dB?
So the intercept range outside the sector (although only about half the range within the
sector) is still very long and much greater than the probable maximum range of the
intercepted sonar, perhaps 20 km.A useful goal would be to reduce the out-of-sector
intercept range to this value, i.e., 20 km. What level of sidelobes would achieve this?
Therefore the sidelobes would need to be 150 - 92 = 58 dB down. This sidelobe level
is quite impractical and, even theoretically, it could only be achieved by severely
compromising the sector SL and beamwidth.
Since all measures to reduce LPI compromise the performance of an active sonar, it is
reasonable to conclude:
Virtually any sonar transmission will be intercepted at a range far in excess of its
own detection range.
The emphasis must therefore fall upon LPE. Note, however, that mine hunting sonars,
which operate against close targets using high frequency transmissions, have a very
low probability of intercept. This is because absorption severely reduces the range
where the sound intensity becomes negligible.
the incidental advantage that it can often use an existing search sonar array,
transmitters and receivers.
Secure encoded message channels must be of short duration to reduce intercept
opportunities. Therefore they need more bandwidth than a simple voice channel.
A typical bandwidth is 10-50 W .The actual operating fiequencies are corre-
spondingly high and the range much reduced.
SE = (SL - PL) - N + DI - DT
This is the basic passive sonar equation: SL is no longer the radiated noise of a
target but is the known source level of the active transmission.
The detection process for the simplest communications sonar is based on the
received signal-to-noise ratio, in the bandwidth of the communication channel,
and therefore
D T = 5 1 0 g d + 10logBc
Both signal and noise are referred to the same bandwidth, hence the term
10 log B, where B, is the bandwidth of the communication channel.
Combining these equations and putting SE = 0, we obtain the communications
sonar equation:
I PL = SL + DI - N - 5 10gd - lolog Bc I
Examples of Communications Sonars 259
I PL = SL + DI - N - 5 log d - 10 log B, I
Let DI = 20 dB, 5 log d = 10 dB, N = 45 dB (SS4 at 12 kHz, the vessels will
be stationary or slow moving and the background will be ambient sea noise),
B, = 3000 Hz. Assume spherical spreading and absorption:
PL = 2010g4000+ (1.6 X 4) = 78 dB
SL = 78 - 20 + 45 + 10 + 35 = 148 dB
This very modest SL results from propagation being one-way only, and the transmitter
power will need to be significantly reduced for this application of the array in order to
minimize the probability of intercept.
The dedicated system on the submarine will also have modest requirements. It may
only be necessary to provide a suitable projector and use an existing passive array
(possibly an intercept array). The SL of the submarine array must be adjustable to
minimize the probability of intercept, and in practice both vessels would reduce their
transmitted powers to levels consistent with intelligible communication.
Example S2.5
An encoded message channel uses a band fi-om 50 to 1OOkHz. What SL is now
required for a range of 4 km, the same range as in Example 12.4?
At50kHz
At 1OOkHz
The large difference in SL between the ends of the band is range dependent, even if a
source level of 274 dB were possible, it would be very difficult to equalize success-
fully. Therefore a channel of this bandwidth is only capable of much shorter ranges. A
more realistic design would aim for a range of 1 km,say, then:
At 50 kHz
At 1OOkHz
Now the difference in SL between the ends of the band is only 19 dB and easy to
equalize.
12.12 Problems
12.1 An intercept receiver mounted on a submarine processing an octave band fiom
10 to 20 IcHz detects a CW pulse of duration 1 s at a frequency of 15 kHz. Given
spherical spreading and absorption, what are the range limits of the source of the
signal, assuming it is fiom a sonar whose SL can be varied between 180 and 220 dB?
Assume that DI = 10 dB at 15 kHz and background noise N = 40 dB.
12.2 If the submarine has a TS = 15 dB and the DI of the intercepted sonar is 20 dB,
what are the intercepted sonars likely noise-limited ranges at these limits of SL?
13
Active Sonar Design
13.1 Introduction
This chapter will use the concepts of the previous chapters to design representative
active sonar systems for the detection of submarines and torpedoes. The inter-
actions between the initial requirements and the limitations imposed by platform
size and the environment will be demonstrated by considering practical systems.
Mounted in a keel dome or, more commonly for modern sonars, in a bow dome;
the bow site offers the lowest self-noise.
Operating frequencies between 3 and 15 kHz; the frequency is chosen (it almost
chooses itself) to give an acceptable balance between performance against noise
and reverberation backgrounds for the array size.
Detection ranges from 5 km( 15 kHz) to 20 km(3 kHz) but highly dependent on
the environment.
Representative Hull-Mounted Design 263
Transmit DI = 10log(2h/l) = 9 dB
The source level is SL = 10 log P + 171 + DIt
Ifwe specify SL = 224 dB then the total acoustic power is P = 25 kW
The array elements should be spaced no more than l / 2 = 136 mm. Make the
actual spacing 120 mm. A stave can therefore have 1 150/120 = 10 elements.
(Note all dimensions are approximate and would be modified by practical
considerations.)
At the crossover points of the beams there will be a 'scalloping loss' of 3 dB. To
limit this loss to about 1 dB, the staves must be spaced less than a beamwidth, say
6", so there will be 360/6 = 60 staves. There will therefore be a total of 600
elements in the array and, in transmit, each element must be capable of transmit-
ting 25 000/600 = 40 W of acoustic power. Suppose the radiating surface of each
element to be ~ ( 0 . 0 5 then
) ~ , the total radiating surface is 600r~(0.05)~
= 4.7 m2.
The radiated power intensity is therefore 25/4.7 = 5.3 kW/ m2 and Figure 1.2
indicates that cavitation could be a problem unless the array is at a depth of at
least 10 m. Practical arrays, however, exist and operate at about this intensity at
lesser depths, so it is safe to continue with this design.
264 Active Sonar Design
Based on the above analysis and the procedures and recommendations from
previous chapters, we are now in a position to attempt to complete Table 13.1, a
table of parameters for the sonar.
Pulse length
The pulse length, T, has been chosen to be 1 s for both FM and CW pulses.
Increasing T beyond this might increase the noise-limited range slightly but it
does have operational implications - the sonar would be blind out to 1500T/2
metres, or even longer for the first pulse if both pulses are transmitted sequentially.
The other argument for increasing T - improved doppler resolution and hence
possibly improved performance against very low target doppler - is particularly
important for a stationary platform (e.g., a helicopter dipping sonar). Chapter 9
discusses the CW target doppler problem in detail but note here that a 1 s shaped
CW pulse will have a bandwidth of about 2 Hz, which corresponds to a doppler
width of about 0.5 knots at an operating frequency of 5.5 kHz. This is already
comparable with spreads in the reverberation background, so increasing the pulse
length will not be useful.
Representcitive Hull-Mounted Design 265
FM bandwidth
Reverberation indices
Detection index
The detection index is 5 log d = 12 dB. This value, for an active sonar, has
P d = 0.9 and Pf, = A decision after n pings, made by the operator or an
automatic detector, results in an incoherent gain of 5 log n(dB). For n = 5 we have
3 dB.
Background noise
The background noise, 55 dB, is the spectrum level of sea state 4 (SS4) ambient
noise at the operating frequency of 5.5 kHz. It is likely to be the dominant source
of self-noise for a typical fkigate at speeds of up to about 15 knots.
Target strength
The target strength (TS) is an average value for a typical submarine. Should the
same value be used for both pulses? It is usual for the TS to have been measured
using a CW pulse of a duration which roughly matches the target extent, say
100 ms for a submarine, and it is the peak value which is quoted. For long CW
pulses, the peak TS equals the integrated TS and the suggested value of 10 dB is
therefore appropriate.
The large bandwidth of the FM pulse means that its effective pulse length is
short, 1.3 ms, and the peak TS does not equal the integrated TS. The correct value
to use for TS is now somewhat less (the peak TS as measured by an equivalent
short pulse, which is perhaps about 5 dB less). However, we can use the same,
266 Active Sonar Design
10 dB, value if we integrate the post-detection samples, but do not include a term
in the sonar equation for this. (In other words, the signal processing does the
integration and not the target.)
Performance
1010gR=20+33-41+10- 12+3=13dB
Pulse length
The pulse length, T, has been chosen to be 2 s for both FM and CW pulses. A 2 s
shaped CW pulse will have a bandwidth of about 1 Hz, which corresponds to a
doppler width of about 1.5 knots at an operating frequency of 1 kHz. This is larger
than possible spreads in the reverberation background, so increasing the pulse
length of the CW pulse may be useful. Note well that as frequency is reduced,
Representative Design 269
FM bandwidth
Reverberation indices
Target strength
Detection index
The detection index is 5 log d = 12 dB. This value, for an active sonar, has
P d = 0.9 and Pfa= A decision after n pings, made by the operator or an
automatic detector, results in an incoherent gain of 5 log n(dB). For n = 5 we have
3 dB.
Background noise
The background noise, 70 dB, is the spectrum level of sea state 4 (SS4) ambient
noise at the operating frequency of 1 kHz. The self-noise of a towed array is that
of the ambient sea noise at towing speeds up to perhaps 10 knots, except at ahead
bearings where the radiated noise of the tow ship may be dominant (see Example
8.2 1.
270 Active Sonar Design
Perjormance
10 log R = 20 + 33 - 41 + 10 - 12 + 3 = 13 dB
and the reverberation-limited range, using the FM pulse, is 20 km. It is the same
as the hull-mounted sonar because the RI is unchanged. (The smaller bandwidth is
compensated by the narrower horizontal beamwidth.)
The small, 1 Hz,bandwidth of the CW pulse results in a very low RI. and its
reverberation-limited range is negligible. As always the CW pulse is only useful
given sufficient target doppler to ensure that detections are noise limited. The
remarks on reverberation-limited performance made above for the hull-mounted
sonar apply equally here.
Representative Design 27 1
Here are some ways we could use to improve the noise-limited performance of
this sonar:
0 Increase SL = +3 dB
This is a total increase in 2PL of 9 dB, and PL becomes 93 dB. The noise-limited
range is now 35 km. Note that a similar increase in 2PL for the hull-mounted
sonar (which is probably only possible by resorting to sector transmissions) would
have increased its range to about 28 km.
A useful insight into the performance of any active system is obtained by
plotting the levels of echo, noise and reverberation against range:
Using the parameters fiom Table 13.2 and assuming spherical spreading and
negligible absorption, here is what we obtain:
EL = 234 - 40 logR
R L ~ ~ = 2 2 4 - 3 3 +lOlog(O.13R)-40l0gR= 182-3OlogR
272 Active Sonar Design
The echo, noise and reverberation levels are plotted against range in Figure 13.2.
The performance results using the plots are for a single ping and are therefore not
directly comparable with the results fiom the fbll sonar equations derived above.
Note in particular how excluding the 3 dB gain from multiple pings halves the
reverberation-limited range. The effects of modifying parameters (TS, Sb,
SL, B, N ) can be easily and quickly seen by drawing new lines parallel to the
appropriate existing lines. For example, increasing EL by 5 dB - by changing SL
and/or TS - results in equal noise- and reverberation- limited ranges of 30 km.
120
12 dB echo excess at
reverb-limited range
100
12 dB echo excess at
noise-limited range
80
60 -
Figure 13.2 Echo, noise and reverberation levels plotted against range
Low Frequency .4ctive: Beware! 273
The frequency shifts due to target doppler (target motion relative to the sonar
platform) are directly proportional to frequency. Target doppler is therefore
inversely proportional to frequency. This is an unfortunate and important result
for active sonars because it means that as the frequency is reduced to improve
noise-limited detection range, there is an increase in the magnitude of target
doppler needed to achieve noise-limited performance using a CW pulse.
(Remember that the performance of a long CW pulse against a background of
even low reverberation is always negligible.) Suppose a target echo needs a
frequency shift of 2 Hz to avoid the main lobe of the reverberation spectrum
(making any detection noise limited), this will be achieved by a target doppler
of about 1 knot at a frequency of 3000 Hz but not until about 10 knots at
300 Hz, a very unlikely target doppler indeed!
274 Active Sonar Design
Pulse length
The pulse length has been reduced to 100 ms. A 1 s pulse would have a dead range
of 750 m, which is obviously too long for this application, and the range accuracy
Torpedo Detection 275
would also suffer. The 100 ms CW pulse will have an acceptable dead range of
75 m and the range resolution of the pulse will be somewhat better than this, say
20 m.
Clcy bandwidth
Target strength
The target strength is the suggested value for a random aspect torpedo, - 15 dB.
Performance
0 +
For the FM pulse 10 log R = 20 33 - 41 - 15 - 12 + 3 = - 12 dB and the
reverberation-limited range is only 60 m
Neither of these results should be surprising. For the FM transmissions the short
range when compared with the submarine target is entirely due to the much lower
target strength of the torpedo (25 dB less than the submarine value), and for the
CW transmissions the negligible reverberation-limited ranges for either target are
entirely due to the much lower reverberation indices of the CW pulses.
The optimum transmission for torpedo detection and classification, then, is a
fairly long CW pulse (as long as possible consistent with dead range and
resolution requirements). But there must be some target doppler. Because of its
wide bandwidth, an FM pulse cannot use target doppler to avoid reverberation, so
unless the reverberation background is very low, an FM pulse will not detect a
torpedo at any target doppler.
Conclusion
Neither task is easy. The first would be helped by exploiting favourable propaga-
tion modes - for example, VDS and Bottom Bounce (not only in an attempt to
increase maximum range but also to fill in any intermediate shadow zones). The
second would be helped by a unified approach to detection and classification; the
careful choice of pulses and processing to best achieve this unified approach; and
thorough integration of data on contacts from all available sources. Submarine
detection needs to be optimised against both noise and reverberation, but without
compromising classiJication(which needs the best possible resolutions in range,
bearing and doppler).
It is the author5 firm belief that these requirements are best met by a
combination of an FM pulse (for best range resolution) and a CW pulse
(for best doppler resolution), and half beam processing (for best bearing
resolution).
Although much effort in recent years has gone into the design of alternative
pulses in attempts to meet the above requirements with one pulse, it has not
completely succeeded because, inevitably, any alternative pulse will compromise
the best performance obtained from combined FM and CW pulses. Furthermore,
careful study into the performances of alternative pulses against all possible
dopplers will frequently reveal significant ambiguities.
Despite the clear advantages of half beam processing, there are two arguments
frequently used against it:
(1) The half beam DI is 3 dB less than the full beam DI. Note, however, that the
bulk steer and add process after detection restores most or all of this reduction
(compare with passive broadband energy and cross correlator detection,
Sections 8.12 and 8.13). In any event the gains in the unified detection and
classification task from having additional information on the bearing spreads
of targets and reverberation will more than compensate for this loss.
(2) The success of the phase comparison process relies upon a good signal-to-
noise ratio. But this is inevitably the case for an active sonar system where
marginal (low SNR) detections are not a consideration: typically there will be
no discernible echo on one ping and a strong echo on the next, due mainly to
Conclusion 279
Finally, dear reader, I hope that you have enjoyed reading my book as much as I
have enjoyed writing it. I remain dedicated to the advancement of sonar - which,
in these cost conscious days, may simply be to arrest any backward move! - and
would welcome dialogue with interested readers, particularly those who point out
the inevitable errors for correction in future impressions.
This page intentionally left blank
Solutions to Problems
Problem 1.1
I = p 2 / p c = (10-4)2/(1.5 x lo6) = 0.67 x W/m2
1, = 0.67 X lo-'* W/m2
Problem 1.2
SL = 101og P + 170.8 + DIt
= IOl0g40000 + 171 + 15
= 232 dB
Problem 1.3
From Figure 1.2 the safe power density is 5 kW/m2. This may be doubled because of
the operating frequency. The radiating surface is 2 X 1 X (n/4) = 1.6 m2. Therefore
the total safe radiated power is 2 X 5 X 1.6 = 16 kW. The factor 3t/4 is needed to
calculate the actual radiating surface resulting from the use of circular elements.
Problem 1.4
BL = SpL + 10logAf
Af = 2000 Hz and BL = 80 dB
Therefore SpL = 80 - 10 log 2000 = 47 dB
282 Solutions to Pmblems
Problem 2.1
The 3 dB beamwidth of a line array is given by
Problem 2.2
The sidelobe levels are given by 2010g[2/n(2m + l)]. For m = 4 the sidelobe levels
are 23 dl3 below the main lobe. Their angular positions are given by
(nL/il)sin8 = f 9 n / 2 , and since L = 6il we have 8 = 3149". Note that the sidelobe
levels do not depend on the dimensions of the array. The higher-order sidelobes only
exist, however, if the array is large enough (in terms of wavelengths) to generate them
Problem 2.3
At lOkHz, il/2 = 75 mm. The array therefore has 2000/75 = 27 rows of
5000/75 = 67 elements. At 8 kHz we have
Problem 2.4
The DI of a baffled cylinder is given by
Halving the height of the array reduces the transmit DI by 3 dB. The number of
elements in the array is halved and therefore the transmitted power is also reduced by
3 dB. The source level is therefore reduced by a total of 6 dB.
283
Problem 3.1
(i) The equation for spherical spreading plus absorption is
At 5 kHz 80 = 20 log r + ( 0 . 3 X
~ therefore r = 7 km
A value of PL = 80 dB is realistic for many sonar systems, and so are the ranges
obtained by assuming spherical spreading plus absorption. The ranges obtained by
assuming cylindrical spreading plus absorption, however, are usually far too large in
practice. Such ranges would only be achieved using the deep sound channel (DSC)
mode where the sound is confined to cylindrical spreading by refraction alone. For all
other modes there will be losses at the boundaries, and these losses will significantly
reduce the range.
Problem 4.1
The target strength of a sphere is given by TS = 1010g(a2/4). Therefore TS =
10 log(OS2/4) = -12 dB. The minimum dimension of the sphere is 1 m, therefore the
maximum wavelength for a reliable value of TS = 0.2 m. The fi-equency must
therefore be at least c/A = 1500/0.2 = 7500 Hz.
Problem 4.2
The target strength of a cylinder is given by TS = 10 log(aL2/2il):
End on, the target strength is that of a sphere, TS = 1010g(a2/4). The TS does not
change with frequency (provided the diameter is at least 31) and therefore, for both
frequencies, TS = 10 log(OS2/4) = - 12 dB.
284 Solutions to Problems
Problem 4.3
The target strength of any plate at an angle 8 to normal is given by
TS = 1010g(5/0.375)~- 16 + 0 = 22 - 16 + 0 = 6 dB
Problem 5.1
At 100 kHz there will be two noise sources: thermal noise and ambient sea noise. The
thermal noise is given by - 15 + 20 logf (f in W ) ,therefore = 25 dB.
Ambient noise, from Figure 5.1, is also 25 dB. So the total noise is 25 + 3 = 28 dB.
The output voltage from the hydrophone is given by
To find the output in a 1000 Hz band, add 10 log 1000 = 30 dB, so 20 log v = 8 dB
and v = 2.5 pV.
Problem 5.2
The isotropic spectrum level of the noise at the array/water interface is 55 +
3 = 58 dB because self-noise = ambient noise = 55 dB.One third-octave is 2000 Hz,
therefore the isotropic band level is 58 + 10 log 2000 = 91 dB.The noise measured at
the output of the receiver beam is 91 - 20 = 71 dB.Because self-noise measurements
are traditionally made at the beam outputs of a sonar receiver, they are sometimes
misleadingly quoted without adding the DI of the beam. So beware of any suspiciously
low values for self-noise.
Problem 6.1
Start with
The TS of a vessel is constant with range. Suppose the echo needs to be 10 dB greater
than the background reverberation to be detected, then for a TS of 20 dB, the system
would be reverberation limited at 8000m. At greater ranges, because the TSR
increases, reverberation would mask the vessel and it would not be detected.
Problem 7.1
Using P d = 0.5, we can read off Pf,:
Therefore Pfa is increased 50 times. Note how a small reduction in the detection
threshold (of which 5 log d is one of the terms) results in a large increase in false
alarms. This will have important implications for an automatic detection system.
Problem 8.1
The only parameter to change is DI:
At 80 HZ DI = 3 + 17 2010g25 = 3 dB
-
At 40 HZ DI = 3 + 23 - 2010gSO = 3 dB
Because the flank arrays are baffled by the hull, DI cannot be less than 3 dB. The PL
and R values are then as follows.
Broadband
Narrowband
PL(dB) 89 86 59 56 64 61
R (km) 30 20 0.9 0.6 1.6 1.1
Broadband detection performance is still very good, even with the much smaller
flank arrays. Narrowband performance is still good for the noisy torpedo, poor but
unchanged for the submarine, and still further reduced for the quiet torpedo. Had the
size of the array been doubled, which would hardly be practical, performance would
still be unsatisfactory against the quiet targets, thus demonstrating the need for long
towed arrays.
Problem 8.2
The time difference between the two paths is given by
This solves to give R = 15490 m. Note that because neither path includes a
component due to a surface reflection from the target, t does not appear in the
equation. To find t one of the paths must include a surface reflection fiom the target.
Problem 8.3
DI = 10 log n = 21 dE3, therefore the array must have 126 elements. They will be
spaced i2/2 at 3000 Hz. A/2 = c/2f = 1500/6000 = 0.25 m, and the array length
is 126 X 0.25 = 31.5 m. At 3000 Hz, PL = 201og 10000 + (0.2 X 10000 X
= 82 dB. The sonar equation to use is
SL = 82 + 6 0 - 21 + 6 - 22 - 9 = 96 dB
This is the spectrum level of the radiated noise of a fairly quiet submarine, implying
that a quite short towed array may be adequate to detect submarines. The problem,
however, is classification and the need to discriminate against other, probably noisier,
targets such as surface ships, for which a capability at much lower frequencies is
essential.
287
Problem 9.1
The frequency shift due to ship (platform) motion is Af = 2Scos6 X (f/c). Be
carefid with the units: S must be in m/s if c is in m/s.
The total shift is 25 Hz, therefore the shift due to the target is 6.2 Hz:
Therefore S cos cp = 0.93 m/s or 1.8 knots. This is the relative velocity of the target
(its doppler); we cannot separate the target speed from its bearing without more
information.
Problem 9.2
There will be a 3 dB processing loss when the replica only matches the echo for 0.7 of
its bandwidth. This will result from a 120 Hz frequency shift:
Af = 2 X ~ ~()f / ~ )
(SCOSX
120 = 2 X (S cos cp) X (4000/1500)
S cos cp = 22.5 m/s or 45 knots
There will be a 1 dB processing loss when the replica matches the echo for 0.9 of its
bandwidth. This will result from a 40 Hz frequency shift:
Even 15 knots is a lot of doppler for a submarine (remember that doppler is relative
velocity, not target speed). Therefore there is seldom any point in using an extended
reference.
Problem 9.3
For a Hamming shaped pulse, Af40 = 3.5/T. Therefore the doppler shift must be
i l . 9 Hz to achieve Rj = -40 dB.
At 3000 Hz
At 1000 Hz
At 300 Hz
Problem 10.1
Make a guess and then iterate. Suppose f = 20 kHz. Propagation loss (dB) will be
given by
PL = 201og r + a r x 1 0 - ~
=20log12000+12a=82+12a
+
Again, equating to (i), 115 = 82 12a therefore u = 2.8 and the frequency is now
about 18 kHz. Note the significance of absorption at these higher frequencies: it is
not practical to achieve, say, 10 km range at frequencies in excess of about 15 M z ,
even given a very strong target like an ocean floor at normal incidence. At 30 kHz
we have
+
therefore r = 5700 m. The TS of the bottom is given by TSB = SB 10 log A. At
normal incidence SB will be at least -20 dB, and at a depth of 12 000 m and for a
beam of solid angle 5" then 10 log A = 60 dB. Therefore TSB = -20 +
60 = 40 dB.
But this target strength does assume the bottom to be flat (within 7.5 m for the
effective pulse duration of 10 ms) throughout the area A.
Problem 11.1
Noise limited
PL = 77 = 2010g r + ur X ( a = 50 dB/km)
giving r = 470 m.
therefore PL = 67 dB:
PL=67=201ogr+arX (a=50dB/km)
giving Y = 330 m.
290 Solutions to Problems
Reverberation limited
1 0 1 0 g R = 4 9 + 3 0 - 11 - 2 5 - 10+3=36dB
lOlogR=39+30- 11 - 2 5 - 10+3=26dB
Problem 11.2
Use the approximate formula for absorption attenuation, a and
= 0.05f'.4, assume
spherical spreading.
2PL = 168 = SL - 30 - 35 + 45 - 10 - 10 + 3
29 1
and SL = 205 dB. To flatten the response at the receiver, the attenuation at 60 kHz
must be 2 X (84 - 70) = 28 dB.
This equalization is only correct at 800 m. In practice it would probably be better
to equalize at, say, half maximum range (400m) and accept the smaller changes in
response at other ranges. SL will change with frequency and the equalization should
also take this into account.
Problem 12.1
The intercept sonar equation is
PL = SL + DI - N - 5 10gd - lolog Br
For SL = 180 dB we have
At 15 kHz, a = 1.6and
PL = 2 0 1 0 g ~ +1 . 6 X
~ lo-' = 100 dB
Y = 12000 m
PL = 2010g Y + 1 . 6 X
~ lop3 = 140 dB
Y = 31 000 m
Note that the pulse type and duration is irrelevant to intercept (provided its bandwidth
is less than B r .
Problem 12.2
The active sonar noise-limited equation is
2PL= 1 8 0 + 1 5 - 4 0 + 2 0 + 0 - 1 0 + 3 = 1 6 8 d B
PL=201ogv+1.6rX 1 0 - 3 = 8 4 d B
Y = 5000 m
292 Solutions to Problems