Waite A D Sonar For Practising Engineers
Waite A D Sonar For Practising Engineers
Waite A D Sonar For Practising Engineers
Third Edition
A. D. Waite
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
ISBN 0 471 49750 9
Typeset in 10/1 Ipt Times by Keytec Typesetting Ltd., Bridport, Dorset
Printed and bound in Great Britain by Antony Rowe Ltd., Chippenham, Wilts
This book is printed on acid-free paper responsibly manufactured from sustainable forestry,
for which at least two trees are planted for each one used for paper production.
Contents
Preface
About the Author
Introduction
xiii
xv
xvii
Sound
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
1.10
I . 11
I . I2
1.13
1.14
Wave motion
Sound pressure
Reference intensity
Source level
Radiated power
Limitations to sonar power
Cavitation
Interaction
Changes to arrays
Projector sensitivity
Hydrophone sensitivity
Spectrum level
Sound in air and in sea water
Problems
1
1
3
4
5
5
6
8
8
9
9
10
11
12
Arrays
13
2.1
2.2
2.3
2.4
2.5
2.6
2.7
2.8
2.9
2.10
2.1 1
13
13
14
15
18
20
20
21
24
27
28
...
Contents
Vlll
29
31
32
35
36
36
36
37
39
40
41
2.12
2.13
2.14
2.15
2.16
2.17
2.18
2.19
2.20
2.2 1
2.22
DI of a simple dipole
DI of a line array
DI of a planar array
DI of a cylindrical array
DI formulae based for simple arrays
Conformal arrays
Spherical arrays
Volumetric arrays
Beamformers
Domes and arrays
Problems
3.1
3.2
3.3
3.4
3.5
3.6
3.7
3.8
3.9
3.10
3.1 1
3.12
3.13
3.14
3.15
3.16
3.17
3.18
3.19
Propagation loss
Losses
Spreading losses
Absorption losses
Spherical spreading and absorption
Propagation in the real ocean
The speed of sound
Sound speed profiles
Deep sound channel
Reliable acoustic path
Surface duct propagation
Convergence zone propagation
Bottom bounce propagation
Propagation loss models
Ray theory and the Hodgson model
Hodgson example
Performance prediction
Multipath propagation
Problems
43
43
Target Strength
67
4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
4.9
4.10
4.1 1
4.12
4.13
Definition
Formulae
Measurement
Dependence on pulse type and duration
TS of a sphere
TS of some simple shapes
TS of small targets
Mine target strength
Torpedo target strength
Submarine echoes
Beam aspect target strength
Bow aspect target strength
Submarine target strengths
43
44
46
48
49
50
51
52
53
54
56
57
58
59
61
64
66
66
67
67
68
69
69
70
72
72
73
74
74
75
75
Contents
iX
4.14
4.15
4.16
4.17
Towed arrays
Target strength reduction
Practical values
Problems
77
78
80
81
83
83
84
86
89
90
90
90
91
91
92
92
92
93
93
93
94
95
97
99
100
101
103
Reverberation
6.1
6.2
6.3
6.4
6.5
6.6
6.7
6.8
6.9
6.10
6.1 1
6.12
6.13
6.14
6.15
Sources of reverberation
Scattering and reflection
Boundary roughness
Classes of reverberation
Backscattering strength
Reverberation target strength
Volume reverberation
Boundary reverberation
Scattering layers
Volume scattering strength
Sea surface scattering strength
Bottom scattering strength
Variation with fi-equency
Reverberation under ice
Problems
103
103
104
105
106
106
107
119
'7.1
110
111
112
113
114
116
117
11:
119
Contents
7.2
7.3
7.4
7.5
7.6
7.7
7.8
7.9
Passive Sonar
8.1
8.2
8.3
8.4
8.5
8.6
8.7
8.8
8.9
8.10
8.1 1
8.12
8.13
8.14
8.15
8.16
8.17
8.18
8.19
8.20
8.2 1
8.22
8.23
8.24
Radiated noise
Radiated noise: source level
Nature of radiated noise
Practical values
Broadband and narrowband
Normalization
A Note on Swaths
Passive arrays
Passive aural
Passive displays
Formulae for detection threshold
Broadband square law detector
Broadband cross-correlator detector
Narrowband processor
Narrowband amplitude detector processor
Worked examples
Passive ranging
Triangulation
Vertical direct passive ranging
Horizontal direct passive ranging
Towed arrays
Bearing ambiguity
Self-noise
Problems
Activesonar
9.1
9.2
9.3
9.4
9.5
9.6
9.7
9.8
9.9
9.10
9.1 1
119
120
120
120
121
123
124
124
125
125
125
126
127
129
131
131
132
133
135
136
137
138
139
140
141
146
147
149
153
156
157
158
159
161
161
163
163
164
165
168
169
171
174
175
178
Contents
9.12
9.13
9.14
9.15
9.16
9.17
9.18
9.19
9.20
9.2 1
9.22
9.23
9.24
9.25
9.26
9.27
9.28
9.29
9.30
9.3 1
9.32
9.33
9.34
9.35
9.36
9.37
xi
179
180
181
183
184
185
187
190
192
200
20 1
202
206
206
208
210
21 1
212
213
213
214
215
216
217
218
219
221
10
10.1
10.2
10.3
10.4
10.5
10.6
Common features
Echo sounders
Echo sounder: design example
Side scan sonar
Side scan sonar: design example
Problem
22 1
22 I
223
226
229
23 1
11
233
Overview
Two broad classes of mine
Backgrounds to mine detection
Range and bearing resolution
Design of a mine hunting sonar
11.6 The threats
11.7 Design example
11.8 Performance
11.9 Classification
1 1.10 Mine avoidance
233
234
234
235
235
235
236
239
240
242
11.1
11.2
11.3
11.4
11.5
Contents
xii
11.11 Mine avoidance sonars
11.12 Problems
Intercept sonars
Communications sonars
Function of an intercept sonar
Intercept sonar equation
Worked examples
Reduction in the probability of intercept
Reduction in the probability of exploitation
Effectiveness of intercept sonars
Communications sonars
Communications sonar equation
Examples of communications sonars
Problems
13
13.1
13.2
13.3
13.4
13.5
13.6
13.7
13.8
13.9
Introduction
Submarine detection
Hull-mounted surface ship sonar
Representativehull-mounted design
Longer ranges
Towed transmitters and towed array receivers
Representative design
Low frequency active: beware!
Torpedo detection
243
248
249
249
250
250
252
253
255
256
257
257
258
259
260
261
26 1
26 1
26 1
263
267
267
268
273
274
Conclusion
277
Solutions to Problems
281
Index
293
Preface
Most books on sonar - the use of underwater sound for the detection, classification 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 mathematics 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:
Echo sounding and side scan sonars for civil applications
Communications sonars
Low frequency active sonars
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:
xiv
Preface
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
Introduction
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)
Magnetic methods include self-generated fields or perturbations of the earths
magnetic field, known as magnetic anomaly detection (MAD). Research continues
into these methods but underwater sound is still unsurpassed, in spite of formidable difficulties facing its propagation through a highly variable medium and from
the noise and reverberation backgrounds to detection. Sonar (sound navigation
and ranging) - by analogy with radar (radio detection and ranging) - uses
underwater sound for the detection, classification and location of underwater
targets.
Passive sonar listens to the sound radiated by a target using a hydrophone, an
underwater microphone, and detects signals against a background of the ambient
noise of the sea and the seZfnoise of the sonar platform (an omnibus term to
describe any vessel or site possessing a sonar system). Passive systems can be
made directional, therefore the azimuth (horizontal bearing) of a signal is known.
The nature of the signal - its frequency spectrum and how it varies with time will help to classifL the target.
Introduction
xviii
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:
(2) I
(dB)
Introduction
xix
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:
2PL = SL
+ TS
N + DI + 10 log T - 5 log d
Introduction
xx
where
c = velocity in metres per second (
ds)
xxi
Introduction
The radar will detect the submarines periscope and, assuming that the electromagnetic 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 countermeasures, and finally to prosecute an attack. Radar systems, therefore, will never
be given opportunities to detect the submarine.
xxii
Introduction
&
c
i_ I ) ) ) ;
Radiated noise at source = 120 dE3
Figure 0.2
= 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.)
xxiii
Introduction
Transmittedpulse = 200 d
target strength
The simplest active sonar will have an omnidirectional projector and hydrophone (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
10 = (200 - 2PL
+ 10) - 60
PL = 70 dB
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.
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.
Sound
Particle motion
IP = ( P C ) U ]
where
1.5 X
Reference Intensity
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
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.
Sound
SL = lOlog
);(
= 1OlogP
+ lOlog(1.1846 X 10)
= 10logP+170.8dB
DIt = 10 log
(k)
lomni
where
+ 171.5 dB.)
+ 170.8 + 10 = 181 dB
and
SL = lOlog40000 + 170.8
+ 20 = 237 dB
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 fiequencies 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.
Frequency: the threshold is little changed up to about 10 kHz, but between
1OkHz and, say, 300 kHz the threshold increases roughly linearly with fiequency.
Cavitation
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?
Radiating surface area = (nX 1 X 1) X (n/4) = 2.5 m2
Maximum radiated power = 2.5 X 2 = 5 kW
And assuming DIt = 10 dB,the source level is
SL = lolog 5000
+ 171 + 10 = 218 dB
Example 1.2
What is the maximum radiated power at 50 kHz to avoid cavitation at a depth of 5 m
for a planar array of radiating surface 0.1 m2 and pulse length 0.5 ms?
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
And assuming DI, = 20 dB, the source level is
SL = lOl0g3000
+ 171 + 20 = 226 dB
Example 1.3
What is the maximum radiated power at 200 kHz to avoid cavitation at a depth of
50 m for a square array of radiating surface 0.01 m2 and pulse length 0.5 ms? At this
frequency, the array has a small size of 100 mm X 100 mm, which is 131 X 131.
The radiated power intensity may be increased 3 times for the short pulse and 20 times
for the higher frequency. Therefore
Maximum radiated power = 0.01 X 50 X 3 X 20 = 30 kW
And assuming DIt = 20 dB, the source level is
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.
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 interactions 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.
Hydrophone Sensitivity
sv = lOlog
;( );
= SL - 20logv
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.
Sh
=:
20 log(v/p) = 20 log LI
20 log p
(dB/V)
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
+ 100 = - 100 dB
v:= 10pv
10
Sound
Example 1.6
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.
It is seldom necessary to calculate BL for a non-flat spectrum. A practical sonar
receiver will either equalize @re-whiten) its input, or process (simultaneously) bands
narrow enough for the spectra to be assumed flat.
intensity
Pressure
I (W/m2)
p <N/m2>
Particle
velocity
24 (m/s)
Particle
displacement
u / o (m at
800 Hz)
Threshold of hearing
0 dB re 20 pPa
Conversation
60 dB re 20 pPa
Threshold of pain
120 dB re 20 pPa
10-12
2 x 10-5
5 x 10-8
lo-"
10-6
2 x 10-2
5 x 10-5
10-8
20
5 x 10-2
Sound level
DSS2 at 10 kHz
40 dB re 1 pPa
Typical reverb at 2 km
range 120 dl3 re 1 pPa
Array: source level
220 dB re 1 Clpa
intensity
I (W/m2>
Pressure
p (N/m2)
0.67 x 10-14
10-4
0.67 x
0.67 x 10-4
~-
Particle velocity
Particle
u (m/s)
displacement
U/UJ (m at
10 Hz)
6.7 x
10-11
1.1 x
10-15
6.7 x 10-7
1.1 x 10-11
10-5
6.7 X lo-*
1.1 x
(1 atm)
Note: it is interesting that 1.1 X 1Owl5 is only 3 millionths of the diameter of a hydrogen molecule (0.28 nm)
12
Sound
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
14
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.
15
Directivity of u Dipole
Figure 2.1
Simple dipole
S2H = dsin0
and the phase difference is
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
Arrays
16
PI +P2
/
/
/
/
/
/
/
/
P2
Figure 2.2
pa, on
- = cos
sin e)
Pa
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
(7 11
Cd ) 7
cos -sin8
nd .
--slnO
;1
= 3 dB
=-
=-
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.
Table 2.1 3 dB beamwidths
Element spacing, d
(wavelength units)
6 (degree)
3 dB beamwidth
(degree)
90
30
19.5
180
60
39
14.5
U-
270"
800
18
Arrays
b(8) = 20 log
{sin(nnd/A)sin 0)
{ n sin(xd/A)sin e)
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
where
L = length (m)
f = frequency (kHz)
0, = steer angle (degree)
At broadside, OS = O", this formula simplifies to
I I:
2e3 = -
19
3 dB beamwidth (degrees)
910
20
30
40
50
Figure 2.4
20
Arrays
OdB
- 20
0" steer
(broadside)
30" steer
OdB
60" steer
Figure 2.5
OdB
90" steer
(endfire)
21
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
pa
_
.
sin{ (nL/il)sin S}
(nL/il)sin 8
b(0) = 20 log
sin{ (nL/iZ)sin 0)
(nL,/l)sin 8
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
0
1.39
2.27
2.85
x-' sin x
20 log(x-' sin x)
1
0.7 1
0.32
0.1
0
-3
- 10
-20
The numerical values in Table 2.2 can be used to derive some useful approximations for the shape of the main lobe. Beginning with the 3 dB beamwidth (Hz):
23
nL
03
1.39 = -sni
1.39 = n-e3
L
)'(
1
if the elements are spaced by 1/2,then
263 = 100/n
Similar expressions can be produced for 10 dB and 20 dB beamwidths and the
results are given in Table 2.3. Note that the 20 dB beamwidth is about double the
3 dB beamwidth.
Table 2.3 Beamwidth expressions
50A/ L
1OO/ n
82A/ L
164/ n
104A/ L
208/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?
24
Arrays
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
Shading
25
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 sixelement 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
26
Arrays
Table 2.4 Commonly used weighting functions
Weighting function
Reduction in
transmit power
Broadening
factor, 3 dB
10 log BF3
(a)
Broadening
factor, 20 dB
(a)
None
Dolph-Chebyshev 30 dE3
Cosine
Gaussian
Hamming
Squared cosine
Binomial
0
3
1
1.3
3
3
3.5
3.5
1.4
1.4
1.5
1.6
2.0
1
1.5
1.5
1.6
1.6
1.5
2
2
3
1.7
I .9
2.2
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:
1.3 X 8.3" = 10.8"
(a) 1.5 X 17.3"= 26"
27
+ 10 log BF3
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
SL = lOl0g6000 + 171
+ 1010g6 = 217 dB
28
Am y s
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
30
Am y s
and substituting
and
31
DI qf a Line Army
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?
Length of line L = 32A = 32
1.5 = 48 m
At I kHz
At 100 Hz
At 1.5 kHz
DI = 18 + lolog 1.5 = 19 dB
= 8 dB
32
Arrays
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
Horizontal beamwidth
(a) Vertical beamwidth
33
= 76/(24 X 1) = 3.2"
= 76/(6 X 1) = 12.7"
4 = (76/24)(1 + 3600/4000)
(')
= 6.0"
0, = 12.7" (steer is only horizontal)
Arrays
34
good practical shape (Figure 2.11). The number of staves is simply related to the
beam spacing and is also, conveniently, a binary number. This makes the
beamforming simple because all beams are identical and provide 360" cover in
azimuth. The horizontal beamwidth is constant and the vertical beamwidth is a
function of the stave height.
Element
Stave
Ring
DI of a Cylindrical Array
35
where
h = height (m)
d = diameter (m)
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
Table 2.5
Formulae for DI
Line
lolog n
Planar
10 log rnn
Square
20 log n
Circular
20 log n
Cylindrical, baffled (see notes) 3 + 10 log rnn
10 log(2 L / A )
10 10g(4fi/1~)
20 log(2Ll.A)
20 log( 1.7742)
10 log(5hdfi) or 10log( 1 1 h d / A 2 )
Notes: m is the number of elements in one-third of a ring, n is the number of elements in a stave; fo is in
37
Volumetric Arruys
0.75
0.75 X 25)
= 18.5 dB
Arrays
38
I
I
.75 m
(0)
-+-
+----
Vectors cancel
Figure 2.12
A ------ 4
Vectors add
Beamformers
39
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
Figure 2.13
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.
2.4 What is the DI of a baffled cylindrical array of height 2 m and diameter 3 m at
5 kHz? How much would the source level of the array be reduced by halving its
height?
3
Propagation of Sound in the Sea
3.1
Propagation Loss
/ 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
PL = lolog(;)
= 1010gr2
or in logarithmic form,
where Y is in metres.
45
Spreading Losses
PL = lolog(;)
= 10logr
or in logarithmic form,
I P L = 10logrI
where Y is in metres.
46
3.4
Absorption Losses
When a sound wave travels through sea water, absorption losses occur through two
principal mechanisms:
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.
Absorption Losses
47
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.02
0.02
0.02
0.06
0.06
0.06
0.14
0.14
0.14
0.33
0.29
0.26
3.80
3.30
2.80
15
16
17
30
35
40
55
60
65
120
125
130
48
3.5
The sum of spherical spreading and absorption losses is a useful working rule for
initial design and performance comparisons:
lPL=20logr+arx
(dB)]
PL = 20 log
+ ar x
fWz)
0.06
b 0.2
c
1.0
d 7
e 33
f 64
a
1
3
10
30
100
300
10-3
49
at 30 kHz, a = 5 dB/km
EM. wave
50
3.7
The speed of sound in the sea depends on the temperature, pressure (depth) and
salinity. A variety of empirical formulae exist for its calculation; here is one due to
Leroy:
c = 1492.9
-4 X
10-2(t - 18)2
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:
Temperature
Salinity
51
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.
Seasonal thermocline: temperature decreases with depth. During sumrner and
autumn the thermocline is strong and identifiable. During winter and spring it is
weak and merges with the surface layer.
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
Speed of sound ( d s )
100 m
-------
Main thermocline
1OOOm
2000 m
Deep isothermal layer
e
P
1
Figure 3.4
increases with depth, the net effect of temperature and pressure changes is to
reduce the sound speed through this layer.
53
Range ___*
Speed
20" vertical beamwidth
Figure 3.5
20
Depth
1500 m
40
60 km
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
0
\ I
1000
1000
lo00
2000
2000
2000
3000
3000
3000
4000
4000
4Ooo-
Depth (m)
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 engineering 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.
55
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.
56
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.
57
Figure 3.11
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
Table 3.2
Wave theory
Ray theory
59
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
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
Figure 3.13
Hodgson Example
63
64
65
PeTformance Prediction
PL (dB)
60
1 00
140
20
Range (km)
PL (dB)
60
100
140
1
10
20
(b)
Range (km)
PL (dB)
60
100
140
10
(C)
20
Range (kin)
Figure 3.15 Propagation loss graphs for the North Arabian Sea: (a) standard SSP, (b) modified
SSP, (c) modified SSP
66
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($)
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
Peak TS = 20 log( pr/ pi)
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 measurements 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 necessarily 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 measured using a calibrated hydrophone at the target, and the equation is solved for TS.
69
TS of'a Sphere
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 ( ~ )
Therefore a sphere of radius 2 m has TS = 0 dB.
Target Strength
70
TS, dl3
Sphere
Convex
surface
10 log(A/A)'
Plate of
any shape
Rectangular I0 log(ab/A)2
10 log(ab/A)2 20 log(x-' sin x) + 20 log(cos 6)
plate
where x = (2na/A)sin 6
10 log(na2/A)2
Circular
plate
10 log(aL2/22)2
Cylinder
+ 20 log(x-' sin x) + 20 log(cos 6)
10 10g(aL~/22)~
where x = (2nL/A)sin 6
Incidence
Notes
h Y
Normal to
surface
Normal
a is radius
a1 and a2 are
principal radii
A is area
aandbare
Normal
6 to normal sides a 3 b
Normal
a is radius
u is radius L
Normal
6 to normal is length
Note that the scattering cross-section is known in radar as the radar cmsssection (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
angle of 2" to the normal?
= 0.2
m. What is the TS at an
(L/A)sin 8 = 0.87
and from the plot we obtain
201og(x-' sinx) = -16 dB
71
(Llh) sin 0
Figure 4.1
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?
TS = 10 log(ab/A)2
+ 20 log(x-'
sin x )
+ 20 log(cos 6 )
72
Target Strength
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.
71
13
where
a = radius
= length
2xL
/i
x = -sin 8
On the beam (normal to the axis) this reduces to
TS on the ends
4)/(2
0.15)) = 3 dB
1010g(0.15~/4)= -22.5 dB
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.
0.262/0.15) = 1.5 dB
Taeet Strength
74
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
at 5 kHz we have
(3
-
)::;r
TS = lolog
=42dB
75
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.
Table 4.2
Name
KILO
209
RUBIS (nuclear)
Mini-sub
Mini-sub
3GST9
R2 Mala
Mass (tomes)
3000
1300
2700
150
90
30
1.4
Dimensions (m)
Country
LARGE SUBMARINES
74 X 10 X 6.6
56 X 6.2 X 5.5
74 X 7.6 X 6.4
Russia
Germany
France
SMALL SUBMARINES
27
20 X 2 X 1.6
10
Italy, Maritalia
North Korea
Italy, Maritalia
Former Yugoslavia
76
Target Strength
Figure 4.2
Towed Arrays
77
Target Strength
78
Thickness (U4)
Figure 4.3
79
140 m, 13C
70 m. 21C
70 m, 13C
20
40
Figure 4.4
I
60
I
80 kHz
Target Strength
Submarines
Aspect
Beam
Intermediate
Bow/stern
Small
Large: clad
Large
5
3
0
10
25
15
8
5
10
Others
Target
Aspect
Surface ship
Beam
Off-beam
Beam
Off-beam
Random
Beam
Dorsal
Dorsal
Mines
Torpedoes
Towed array
Whale, 30 m
Shark, 10 m
Iceberg
h Y
TS (a)
+25
+15
0
-10 to -25
- 15
0 (max)
+5
-4
10 (min)
Problems
81
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.
5
Noise in Sonar Systems
5.1
Sources of Noise
Noise is the background against which sonars, active and passive, must detect
signals fiom targets. For an active sonar, noise is augmented by reverberations
from unwanted sources (Chapter 6) and the signal is an echo from the target. For a
passive sonar, the signal is also noise - the radiated noise of the target.
This chapter will first discuss the sources of noise in sonar systems and then
consider the total noise as a background to detection - the self-noise of a sonar
system. There are three sources of noise to be considered:
Thermal noise
0
Vessel noise
84
5.2
Thermal Noise
In common with any electrical receiving system, a sonar receiver adds its own
noise to the signals it receives. The designer must ensure that the noise introduced
from this source is negligible compared with the noise coming from the sea itself.
All the results established for the classical case of radio are valid for sonar when
aerial or antenna is replaced by hydrophone.
Any resistance, R, is the source of a thermal noise EMF (electromotive force)
resulting from the thermal agitation of its electrons. The value of this EMF, en
volts, is given by
where
J/K
Thermal Noise
85
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).
The thermal noise is given by
86
5.3
The thermal noise of the sea can only be the dominant background to a sonar at
high frequencies - at least 30 kHz and for all practical purposes at least 100 kHz,
where it equals the ambient noise expected from sea state 2 (SS2).
When the sea is not dead (Le., perfectly isolated from all sources of sound and
only subject to thermal agitation), even though it may appear perfectly calm, it is
subject to an agitation much greater than the thermal noise, particularly at the
lower frequencies (less than 30 kHz). Figure 5.1 plots the mean isotropic spectrum
levels as a function of frequency for various sea states. Departures from the
average of up to A10 dB are quite common, particularly in shallow water and
close to the surface. The ambient noise falls at about 5 dB per octave for
frequencies above 500 Hz.
6.2
0.5
10
20
50
100
Frequency (kHz)
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
Wind speed (knots)
Wave height (m),crest to
trough
0
Sl
0
2
13
0.4
3
16
0.7
19
1.3
5
22
2
6
28
3
87
80
- - - NoRainrain
I
Y
Rainfall rate (d)
k! 60
m
2
30
L_
'u
'u
_I
4 4
40
'
DSSl
0.1
1 10
Frequency (kHz)
Figure 5.2
Rain noise
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
0.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
* The propellers
Flow noise
Figure 5.4 shows the relative contributions of these sources of self-noise at
operational bearings for a typical surface ship at 10 kHz, particularly one where
little has been done to reduce its radiated noise. Modern frigates can improve
significantly on these figures, perhaps achieving SS2 at speeds up to 10 knots. The
self-noise of a sonar is the sum of all noise sources, ambient and vessel, that are
present at the array/water interface.
20
10
30
Speed (knots)
Figure 5.4
90
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.
Flow Noise
91
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.
92
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.
93
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.
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
When P1 = P2
+ 3 dB
+ 6 dB then PT = P1 + 1 dB
94
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.
Figure 5.5
where Noutis the noise power at the output of the receiver and N m i n is the unavoidable 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.
95
Receiver
Now consider the case where the sea state is not negligible. The noise factor
becomes NF, (at a given frequency), and
NF, =
+ Nmin
Nmin
96
and the resultant noise factor for the sonar, NF,,, for a given sea state, will be
and since
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
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
NFres
NFS
1.26 and
dB
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
20 log v = -200
=
-200
+ 20 log p
+ 40 = -160
v = 10nV
(dB)
dB/V
98
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 selfnoise, 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
since Ns >>
Nmin.
Therefore
Again suppose that we only lose 1 dB due to the internal noise of the sonar, then
5 = 1.26
and
NFr
-= 0.26 X 3162 = 29
N F S
dB
and for the same hydrophone, the noise factor of the receiver has to be less than
25 dB:
20 log v = -200
= -200
+ 20 log p
+ 25 = -175
(dB)
dB/V
v=2nV
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:
20 log v = -200
= -200
V =
+ 20 log p
+ 40 = -260
(dB)
dB/V
10nV
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
Therefore
20 log v = -200
+ 25 =
175 dB/V
= 2nV
i.e., the spectrum level of the receiver noise, referred to the input, must be less than
2 nV
75
65
84
74
70
60
78
68
65
55
72
62
60
50
66
60
55
45
60
50
50
40
54
44
45
35
48
38
40
30
42
32
Prt,hEerns
5.2i
101
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?
6
Reverberation
Reverberation
104
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 penetration 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
0.0005
Smooth
surface
/'
1.0
0.9995
\,
Rough surface
f,
Specular direction
Figure 6.1
106
Reverberation
For analytical purposes, sea surface and sea bottom reverberation may be lumped
together as boundary (area) reverberation.
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 omnidirectional, in azimuth, projector and a 10" , in azimuth, receive array is 10" (the
smaller of the two).
107
Volume Reverberation
where
Br
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
Reverberation
108
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
V=
1500 X
2
n X lo6 X 0.22
4
=2
x 104
109
W u m e Reverberation
and
TSR = -80
dB
When the beam patterns intercept a boundary, i.e., the surface or the bottom of the sea
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.
110
Reverberation
with @, in radians.
Figure 6.4
Boundary reverberation
Example 6.2
For the same sonar, what is the TSR at 1000 m for S, = -30 dB?
We calculate
111
Scattering Layers
and
TSR = -30
dB
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).
0
I S[ = s, + 10 log H 1
and this value must be used in the expression for boundary reverberation.
Reverberation
112
'\
/
/
/
\
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
St (day) = -76
Se (night) = -74
+ 10 log 80 = -57 dB
+ lolog 100 = -54 dB
I13
ss (dB)
-40
-50
- 60
(a)
6
Frequency (kHz)
10
Grazing angle
S, (dB)
-30
-40
-50
Sea state 6
Wind speed 30 knots
2
(b)
Figure 6.6
6
8
Frequency (kHz)
10
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, particularly 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).
Reverberation
114
Figure 6.7
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.
Pattern A is characteristic of a smooth bottom, where the specular reflection in
direction OR is large and the scattering in other directions including the
backscattering back to the source in the direction 01, at a grazing angle, 8,is
small.
Pattern B is characteristic of a rough bottom, where the specular reflection in
direction 01 is much greater. (The backscattering strength, Sb, is much greater.)
Sound is also absorbed where it penetrates the bottom in direction OK
115
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
-25
-30
-27
-32
-29
-34
-3 1
-33
-35
-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).
117
Problem
4~
- 40
/
rr--
------
"
10
15
/ -
Ice-free
Wind speed = 30 knots
20
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?
7
The Sonar Equations
120
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.
ISE
= (SL - PL) -
= DT(
121
Display/auto
system
T
MDS
DT
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=(
122
Figure 7.2 shows curves of probability density plotted for N alone and for
S + N.
Figure 7.2
Pd,
123
7.7
Pd
1 .o
0.5
0.2
+N
levels of signal variability are shown, defined in terms of the fluctuation index (FI)
given by
FI =
o(S
+ N)
uN
124
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
pfa= 10-5
Pd
= 0.5
6
6
12
10
9
8
8
10
Pd
= 0.9
8
>>8
20
15
13
9
12
12
USE FOR
pfa=
7
7
13
11
10
9
9
11
Broadband
Narrowband power
= 1
=4
=1
=2
=4
=8
Narrowband amplitude
Active, intercept and communications
Broadband
>>8
21
16
14
11
13
13
FI
FI
IF
IF
IF
IF
Narrowband power
FI
FI
IF
IF
IF
IF
=1
=4
=1
=2
=4
=8
Narrowband amplitude
Active, intercept and communications
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
8.1 Radiated Noise
Radiated noise is the noise emitted by a vessel and received by a hydrophone, or
an array of hydrophones, at some distance from the vessel. Radiated noise is the
source of signals for passive sonars which are designed to detect radiated noise
against a background of ambient and self-noise.
Analysis bandwidth
(Hz)
0-100
100- 1000
1000- 3000
0.3
1
4
126
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.
127
Pructical Values
I60
Frequency (Hz)
Figure 8.2
128
Passive Sonar
(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).
129
130
Passive Sonar
iz
4 Note on Swaths
131
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. Normalization 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
Figure 8.4
132
Passive Sonar
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 characteristics, 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
(b)
f;
lncreasing signal level
Figure 8.5
(c)
.fi
134
Passive Sonar
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.
135
Passive Displays
136
Passive Sonar
.'
;
-...
*.
*:
**
*
*
*. *.
- - - - .
:
:
I
9-
*.
Time
'..-.
Own vessel
Contact
'
.
8.
#
#
} Beam 1
mean(S N) - mean(N)
std dev( N)
signal power
137
BT(DT)2 or
and, substituting for DT, the sonar equation for this passive receiver, also called an
energy detector, becomes
SE = (SL - PL)
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
IPL = SL - N
Note that Te - the time that a signal may be expected to remain in a beam - is the
time between successive display samples (lines) and is greater than l/B.
138
Passive Sonar
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.
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.
So Te = IF/B; substituting for Te we obtain
DT=510gd-510g
(3
-
+ 10 log B
5 log IF
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.
Passive Sonar
I40
+ 10 log B - 10 log(B/B,)
5 log IF
+ DI
5 log d - 10 log B
DT = lolog
[?+
l.05(&)1/2]
+ 10logB
lOlog(B/B,)
PL=SL-N+DI-lOlog
p';id+ ):(
- 1.05
lOlogB+ lOlog(B/B,)+5logn
1'2]
141
140
145
110
115
100
105
155
165
125
135
.-
130
140
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
At 400 HZ
At 200 HZ
At 8 0 H z
At 40 HZ
+
+
+
+
= 29 dB
= 12dB
= 6 dB
= 3 dB
= 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
The narrowband equation is
(PL = SL - N
+ DI - Slogd
lolog B
Threat A
Tonal
SL
N
DI
5 log d
B
Bs
lolog B
10 log( B l B s )
IF (= BT,)
5 log IF
n
5 log n
PL
R (km)
200
165
73
6
10
2
4
3
-3
4
3
30
7
92
40
400
155
68
12
10
4
8
6
-3
8
5
30
7
92
40
Threat B
Threat C
400
40
200
135
125
140
68
85
73
6
12
3
10
10
10
0.1
2
4
0.2
4
8
-10
3
6
-3
-3
-3
1
8
4
5
0
3
30
60
30
7
9
7
64
62
62
1.6
1.3
1.3
80
130
80
3
10
0.2
0.1
-7
0
2
2
60
9
61
1.1
Broadband
The broadband sonar equation is
IPL = SL - N
Worked Exumples
143
Table 8.4 Broadband system performance for
Example 8.1
SL
N
DI
5 log d
B
re
5 log BT,
n
5 log n
PL ( a = 0.2)
(km)
Threat A
Threat B
140
54
29
6
2000
2
18
30
7
134
> 100
110
54
29
6
2000
2
18
30
7
104
50
Threat C
100
54
29
6
2000
10
21
60
9
99
40
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
140 dB from 40 to 400 Hz
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
Passive Sonar
144
140 - 20 log 600 - 20 = 64 dB at operational array bearings
0
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:
~~~
Hertz
Decibels
80
80
40
85
200
73
1500
58
400
68
32
1500
= 1200 m
40
(using c = fA)
Narrowband
We have
[PL = SL - N
145
Worked Examples
Table 8.5 Narrowband system performance for Example 8.2
Threat A
Tonal
SI,
N
D1
5 log d
H
23,
10 log B
10 log(B/ B s )
IF (= BT,)
5 log IF
ri
5 log 12
PL
R (km)
200
165
73
18
10
2
4
3
-3
4
3
30
7
104
160
400
155
68
18
10
4
tl
6
-3
8
5
30
7
98
80
Threat B
200
135
73
18
10
2
4
3
-3
4
3
30
7
74
5
Threat C
400
125
68
18
10
4
8
6
--3
8
5
30
7
68
2.5
40
140
85
18
I0
0.1
0.2
- 10
-3
1
0
60
9
79
9
80
130
80
18
10
0.2
0.1
-7
0
2
2
60
9
76
6
Broadband
We have
IPL = SL - N
Threat
SL
N
DI
5 log d
A
re
5 log BT,
n
5 log n
PL
R (km)
Threat A
140
54
18
6
1000
2
16
30
7
122
>lo0
~~~~
Threat R
1 10
54
18
6
I000
2
16
30
7
92
30
Threat C
100
54
18
6
1000
10
20
60
9
88
20
146
Passive Sonar
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.
147
Triangulation
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
R=
Rmax =
k i n =
148
Passive Sonar
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 available 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.
149
__
--
I R = (223 - p - t)/tan61
150
Passive Sonar
SB (not shown)
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
151
and for this example 6T = 67 ms. 6T is now used as an initial value for the offset 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:
For the downward-looking beam
= 2100/
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
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
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.
153
.+ f
/
Target
Axis of arrays
R2 = p2 + ( R - dr)
where d r is the additional path length due to curvature. We have p
therefore
R2 = (dsin8)2 + ( R - (rr)2
and when R is large compared with d, we have
= dsine,
154
Passive Sonar
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
Figure 8.12
HDPR example
IS5
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
Delay between arrays A and B = 0.0067 s
Delay between arrays A and C = 0.0133 s
Referring to Figure 8.1 1, we have
BE
= 0.0067
c = 10.05 m
CF = 0.0133 c = 20 m
= 10.05 -
10 = 0.05 m
0.05) = 1000 m
I56
Passive Sonar
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 hydrophones 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 programming) 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.
y1
Length (m)
40
40
40
100
50
25
~~
200-400
400-800
800-1600
300
600
1200
Figure 8.14
2.5
1.25
0.625
Nested array
158
Passive Sonar
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
No
Twin urrays
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
In a triplet array, each element now comprises three hydrophones - a triplet - in
the vertical plane (Figure 8.16). Because time delays are measured between all
three pairs, the lefthght ambiguity can be resolved regardless of any rotation of
the array. The method needs a fairly large array diameter to house the triplets and
to produce measurable time delays.
Figure 8.16
Hydrophone triplet
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. Propagation is spherical spreading plus absorption.
9
Active Sonar
162
Active
Sonar
The first and fourth components are normally removed by own doppler nuZZiJication (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,
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.
The replica correlation process used to detect an FM pulse is effectively a matched
filter. Independent outputs from the correlator occur every 1/ B seconds, i.e., at the
resolved pulse length.
9.3 CW Processing
CW processing is similar to narrowband passive processing where, after beamforming, the data is analysed using a fast Fourier transform (FFT) process and
then displayed to an operator andor input to an automatic detection and classification 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.
Active Sonar
164
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
power - 1 (mean noise power
after correlation
BT before correlation
and the DT equation from this process is
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 knowledge of the signal) whereas broadband processes are incoherent (nothing known
of the signal).
A coherent gain is equal to the number of independent samples = BT
0
samples = (BT>'/*
165
= SL
+ TS
2PL - (N- DI
+ lolog B )
(5logd - lolog BT
510g n )
I2PL = SL + TS - N + DI + lolog T
- 5 logd
+ 5 log n I
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
167
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 a CWpulse, B = 1/T, where T is the pulse length.
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
T
CW pulse:
RI = 1010g(10/10) = 0 dB
Example 9.2
FMpulse:
= 300 HZ
4 = 10"
23
+ TS
5 logd
+ 510g nJ
(i:
I$)
168
Active Sonar
9.7 FMpulses
In the strict radio communications sense, sonar FM pulses are not frequency
modulated at all. They are simply wideband pulses where the frequency changes
throughout the duration, T, of the pulse (Figure 9.2).
Reference
Frequency
Frequency
Time
Time
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
At = (T/B)Af
where A j is the frequency shift due to target doppler, and
cT
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
250
2
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
8
Frequency (Hz)
12
16
Active Sonar
172
for I tl d T / 2
otherwise
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
8
Frequency (Hz)
12
I '
16
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
173
CW Pzilses
Af3 =-
dB
Of
-70,
-800
8
Frequency ( Hz)
12
16
(J
Active Sonar
174
Rj (a)
- 20
I s, 10 kH2 pulses
-40
-60
I
-1.2
-0.8
-0.4
Figure 9.6
0
0.4
Target doppler (knots)
0.8
1.2
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.
175
Target at 45"
Target on beam
f O
h0.6
i0.9
f1.2
Zt0
-1-0.9
A0
$10.04
ct0.06
-1-0.08
0
10
il.4
15
-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.
Active Sonar
176
Table 9.2 Doppler spreadderrors due to ship motion
Sea state
Pulse length
(ms)
Target on bow
Target at 45"
Target on beam
1000
1000
0
0
2.7
250
250
0
0
1.8
0.9
0.4
0.2
1.3
0.6
0.3
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
Doppler spread (knots)
2
1.2
4
1.6
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
Their relative importance is highly dependent on the sonar platform. Beamwidth
177
f 0 . 5 knot
f 0 . 7 knot
f 1.O knot
f1.5 knots
178
Active Sonar
v knots
t
Doppler
'y;
Ahead
Beam
Astern
Platform
-v+
-+
-v
-2v
Figure 9.7
179
12
--
-4
0
4
Target doppler (knots)
12
180
Active Sonar
Table 9.4
Sonar parameters
~~~
Parameter
Frequency
Pulse length, T
Source level, SL
FM bandwidth
CW bandwidth
Beamwidth, full aperture
Beamwidth, half-aperture
Phase bin width (2 bins)
Reverberation index, RI
FM, full beams
FM, half-beams
cw
DI (receive)
Integrated target strength, TS
5 log d
5 log n
Background noise level
Value
Units
10
kHz
2 10
dB
400
HZ
2
8
16
2
Hz
degrees
degrees
degrees
17
23
-6
12
10
I0
3
dB
dB
dB
dB
dB
dB
dB
dB
50
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 halfbeam processing. The last two items are applicable to both FM and CW pulses.
Enhanced FM detection performance for all doppler values
Noise-limited CW detection performance once the doppler spectrum of the
reverberation within the main lobe of a beam is exceeded
Beamforming
181
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.
Active Sonar
182
half-beams
\
120"
60"
-0"
- 60"
- 120"
- 16"
- 8"
8"
16"
Full beam
Left half-beam
Right half-beam
Figure 9.9
I83
r2
f - - - - A 7
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
PD gain = 5 log(BTinty1) = 5 log(400 X 20
= 4.5
dB
Active Sonar
184
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
Bulk steer and add to form two full beams
A
Left
channel
Phase
compare
Doppler
channel 1
Doppler
channel 2
Left
half-beam
Doppler
channel
(filter) N
Phase
compare
Right
half-beam
--+
'
channel
I (filter) N 1
1 Doppler
ichannel .V
185
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 determining 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
9.17
Figure 9.12 shows the beam pattern of one of a pair of codirectional half-beams
formed using half the aperture of a 32A array of 64 elements, together with the
phase difference plot resulting from the complex cross-correlation of the pair of
half-beams. Hamming shading is used. The f 2 " beamwidth of each half-beam
corresponds to a phase difference of about .f120". By using the bulk steer and
phase binning process described earlier, phase differences greater than this can
eliminate all returns from not only the sidelobes but also the skirts of the main
lobe. Note that the phase difference remains greater than f120" for all values of
azimuth beyond the 3 dB points of the main lobe. In practice, imperfections in the
array, and unwanted curvature in the case of a towed array, may cause re-entrants
to occur in the phase difference plot beyond the i120" limits which will result in
returns through sidelobes. To minimize this, the sidelobe levels should be reduced
by beam shading - Hamming shading in this example.
186
Active Sonar
-33-36-39-42-45-48-
-80
-60
-40
-20
20
40
(a)
60
80
Azimuth (degrees)
0
-30
-60
-90
- 120
- 150
-180
-80
-60
(b)
Figure 9.12
-40
-20
20
40
60
80
Azimuth (degrees)
Detection Performance
+ DI + lolog T
5 10gd + 5 log TZ
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 .
Table 9.5 LPFM reverberation-limited range
Reverberation-limitedrange (km)
Sb
Full-beam processing
Half-beam processing
= -30
Sb = -37
Sb = -40
8*
40
160
80
320
32
"Note that, in high reverberation conditions, the noise-limited detection range, 9 km. is only achieved by
using half-beam processing
sb =
Sb
Sb =
Sb
sb
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.
ClassiJicationand display, against any background, are enhanced by the availability of fine bearing (phase bins) providing an extra dimension for echo shape
recognition and improved bearing accuracy for tracking.
+ TS
5 log d
+ 5 log t~ - Rjl
Note the additional term, Rj; this is the reverberation rejection and it determines
how the performance of the CW pulse varies with doppler. ( R j is zero for FM
pulses.) We are effectively calculating the detection range of a target with zero
doppler and then multiplying it (adding dB) by some factor determined by Rj.
For the CW transmissions, we are less interested in the reverberation-limited
range and pose the question, How much doppler before the detection is noise
limited? To answer this, assume the reverberation-limited range is 20 km (i.e.,
always longer than the noise-limited range) giving 10 log R = 13 dB. Then solve
the reverberation-limited equation for R j, again using the parameters of Table 9.4.
The target dopplers required to achieve the R j values (Table 9.6) are then obtained
from Figures 9.6 and 9.8, but adding another 0.7 knot as suggested in Section 9.5.
Note that the doppler magnitudes required using full-beams are for a moderate
platform speed, 10 knots and assume quite low sidelobes (20 to 30 dB down). In
practice, noise-limited detections will frequently only be achieved for target
Detection Pevfornzance
Table 9.6 Target doppler
-30
-37
- 40
--30
-37
-40
Full-beam processing
-27
-20
-17
Half-beam processsing
-27
-20
- 17
8
3
2
1.1
1 .o
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.
Because of the target doppler requirement, it is operationally unsafe to use only
a CW pulse for any active sonar - unless there is a high probability of sufficient
target doppler - and the preferred mode of operation is to transmit both FM and
CW pulses every ping.
190
Active Sonar
ttxi'(s.S+ a )
- 10
Figure 9.13
+ 10
+20
Target strength (dB)
191
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
Figure 9.14
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 diagrams for CW, FM and PRN pulses and concentrates on results and practical
examples.
193
A rnbigui? Diagrums
CW pulses
By equating the ambiguity function for a rectangular CW pulse to 0.5 (half-power)
and first putting t = 0 and then cffo = 0, we obtain dfo = 0.88/T and f = 0.6T,
and the resulting ambiguity diagrams for long and short CW pulses are sketched
in Figure 9.15. The area of the elliptical ambiguity diagram - for any length of
pulse - is approximately 0.6T(O.88/T)(n/4) = 0.4. The values for dfo and t can
be used to determine the resolution (extent of the ambiguity) of the pulse and the
number of references (in this case matched filters) required to cover the expected
target doppler.
Figure 9.15
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
T
z ___
Therefore
0 . 8 8 ~ 0.88 X 3000
= 0.7 knot
2Tfo 2 X 0.5 X 4000
7Jz-=-
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
R=
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).
f = 0.88/B
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:
Sfo
=-2vfo = 0.6B
C
therefore
900 B
v = - (knots)
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
t=--
therefore
660
R = __
B
(metres)
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
R=-=
B
1.32m
900 B
u=------= 45 knots
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
Figure 9.17
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?
197
A mhiguity Diagrams
We have
= 6 6 0 / B = " ~ 0 . 6m
6
ZI
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?
Suppose we begin like this:
0
= 5200/BT
knots for BT
500
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
An alternative to a broadband FM pulse is the pseudorandom noise (PRN) pulse.
Figure 9.18
Ambiguity Diagrums
f 99
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.
= 660/B =
1.65 m
200
Active Sonar
Table 9.7
Pulse type
cw
450T
LPM
LFM (BT < 100)
LFM (BT > 500)
PRN
750/ B
750/B
750/ B
750/ B
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.
9.22
20 I
202
Active Sonar
.4 ctive Displa~vs
203
Table 9.8 Parameters of a typical sonar
Parameter
Frequency
Pulse length, T
FM bandwidth
CW bandwidth
Beamwidth, half-aperture
Beam spacing
No. of beams, 360" cover
Doppler channel width
Doppler cover
FM increments
classification (3 ms)
surveillance (20 ms)
CW increment ( I 33 ms)
Phase bin width
Range scales
Value
Units
10
300
400
kHz
ms
Hz
Hz
degrees
degrees
-7
16
11.25
32
1
k30
2
15
100
1
4, 8, 16
knots
knots
m
tn
m
degrees
km
RANGE
5 12 pixels
BEARING
5 12 pixels
32 X 16 phase bins
Range scale
Range increments
Pixels per increment
4
267
2
8
533
1
16
1067
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 equivalent 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
BEARING, DOPPLER
5 12 pixels
Range scale
Range increments
Pixels per increment
4
40
12
8
80
6
16
160
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.
205
Active Displa-ys
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.
Noise
Range
Figure 9.20
Ping n
Pingn
Ping n
Pingn
_______+
206
Active Sonar
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 classification, 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:
Other modes of the sonar such as passive or intercept
Other sonars on the same platform
Sonars and other sensors on consort vessels or aircraft
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.
Amplitude Prqfiles
213
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).
9.31
Amplitude Profiles
The amplitude profiles approach imitates the way an experienced operator assesses
the contact echo. If the resolution is sufficiently fine (1 m or so) then the echo
from a contact will display some detail of the reflecting structure of the contact.
The amplitude of the returns along the length of the contact will be proportional to
the reflectivity of the contacts structure. A submarine may well return strong
echoes from the bow and the fin, and weaker echoes from the hull and the casing.
These reflections yield a set of strongly aspect-dependent amplitudehange profiles
2 14
Active Sonar
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 processing, 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.
215
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.
Figure 9.22
Multipath example
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 displacements will now be estimated for the most relevant propagation modes, i.e., surface
duct in deep water and in shallow water.
216
Active Sonar
[8(i)
hd]
121
where
hd = duct depth
c = speed of sound
= velocity gradient
(
E)
X 1001
0.017
112
= 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.
217
Duct
10 008
5 Olh
10 032
5 064
8
16
32
64
Active Sonar
218
Figure 9.25
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 representative 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.
Table 9.10 Path length differences at different angles
e
1"
4"
9"
1"
4"
9"
Direct path, R
Indirect path
Difference
(m)
(m)
(m)
11 458
11 449
11 365
22 916
22 898
22 730
1 1 460
11 464
11 506
22 920
22 928
23 012
2
15
141
4
30
282
Bottom
reflections
Loss per
reflection
(dB)
Total loss
(dB)
1
2
5
1
3
9
1
2
3
1
2
3
1
4
15
1
6
27
219
Problems
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?
10
Echo Sounding and Side Scan
Sonars
222
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.
223
500 nt
40 kHz
5"
59 dB at 500 m
41 dB at 100 m
26 dB
133 ps
-39 dB
35 dB
+10 dB
-50 dB
187 dB
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).
PL = 20log Y
Y
= 500 m,
Y =
+ ar X lou3
a
=9
(dB)
d B / h therefore PL = 59 dB
A solid angle beamwidth of 5" will result from a square array of side L, or a
224
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
12PL = SL
Therefore
SL = 2PL - TS + N
+ 50 + 35 - 26 + 39 + 10 - 3 = 187 dB
22 s
Therefore
lOlOgP=187-171 - 2 6 z - l O d B
giving P = 100 mW.
The acoustic power intensity is 10-1/0.382 = 0.7W/m2 and, clearly, cavitation
will not be a problem. The array design is therefore practical and the remaining
parameters may be inserted in Table 10.1.
The echo sounder or fish finder would in practice have a selection of source
levels and pulse lengths under operator control. It would clearly be practical to
have a source level much greater than 187 dB - a value of 207 dB would still only
need I0 W of acoustic power.
The broadband, FM, techniques essential for the detection of many larger
targets at longer ranges, and offering excellent performance against a background
of reverberation, are not essential for these systems. The necessary range resolution has been shown to be obtainable using a very short CW pulse, and fish
detection is normally against a noise background only; the major source of
reverberation is the sea bed, which is beyond the range of most fish species of
interest. An exception is the volume reverberation due to dense concentrations
of plankton, which can limit the detection of individual fish.
226
Figure 10.2
An echo sounder transducer pointing directly downwards with a solid beamwidth of 5".
227
250 tK
I
I
I
I
1
50m
4----'1
I
Range (&)
Beamwidth 2"
'Beamwidth'
4.5 m
'Beamwidth'
9.0 m
228
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
Example 10.1
Suppose the display has 250 lines, each one displaying a ping's worth of data, then
V = - l5Oo = 3 m/s (6 knots)
2 X 250
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.
229
200 rn
100 kHz
1"
45"
53 dB at 200 rn
21 dB
133 ps
-39 dB
30 dB
-40 dB
-50 dB
204 dB
The propagation loss will be approximated by spherical spreading plus absorption (not quite so accurately in this instance, because the sea bed presents one
boundary after the wavefront reaches it).
PL
Y
= 2010g Y
+ ur X lo-'
(dB)
PL = 53 dB
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&. Therefore, the width is 76/(100 X 2) = 0.38 m and the height is 76/(100 X 90)
230
= 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)
= 10 log [(4 X 0.38 X 0.0084)/0.0152] = 18 dB
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
T=-
2 x m - 2XO.l
- -=
C
1500
133
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 noiselimited 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
cT
= 0.34 m2
+ 1010gO.34 = -40
dB
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
+ N - DI
1010gT
+ 5logd - 510g H
+ 39 + 10 - 3 = 201 dB
SL = 106 + 50 + 30 - 21 + 39
SL,,,
+ 10 - 3 = 211 dB
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?
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 antisubmarine 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
235
The Threats
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.
0.3 m, TS = -30 dB
236
Maximum range
Frequency
Transmit beamwidths
vertical
horizontal
Receive beamwidths
vertical
horizontal
PL
Directivity index, D1
transmit
receive
Source level, SL
Pulse length, T FM
300 rn
200 kHz
5"
90"
1"
1"
65 dB at 300 m
cw
Bandwidth, B
FM
cw
10logT FM
cw
16 dB
43 dB
192 dB
10 ms
125 p
80 kHz
8kHz
-20 dB
-29 dB
30 dB
-20 dB
-20 dB
10 dB
3 dB
237
Design Example
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
m2. To avoid cavitation the maxiradiating surface is about 75 X 4 = 300 X
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
SL = 10 log P
+ 171 + 16 = 192 dB
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
DI,
=3
+ lOlog(100 X 100) = 43 dB
238
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 background 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.)
CW pulse: a comparable bandwidth would imply a very short pulse
(1 /80 000 = 12.5 p).Such a short pulse would have negligible noise-limited
performance and a compromise is necessary. Let the CW pulse have a duration
of 125 ps. The bandwidth is then 8 kHz.
239
11.8 Performance
The noise-limited active sonar equation is
I2PL = SL + TS - N
+ DI + 1010gT
5logd
+ 5 10gn I
1 lolog R = lOlog(B/&) -
S b --
11
+ TS - 5 10gd + 5 log TZ I
* 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.
CW pulse: 1OlogX = 3 9 + 30 - 11 - 30 - 1 0 + 3 = 21 dB therefore R =
139 m. The reverberation-limited range requirement is not 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
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
Figure 11.2
Shadowgraph effect
Classlfication
24 1
242
243
sb
Value
10
100
210
5000
8
16
1
37
12
-30
10
3
50
-30
Units
kHz
ms
dB
Hz
degrees
degrees
degrees
dB
dB
dB
dB
dB
dB
dB
Detection performance
The noise-limited performance is given by
=I
210 - 30 - 5 0 + 12 - 10 - 1 0 + 3 = 125 dB
244
If we assume spherical spreading and absorption, the range (from Figure 3.3) is
1000 m.
The reverberation-limited performance is given by
11 - 3 0 - 1 0 + 3 =
19dB
11.11.2
Dedicated design
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, however, be high enough to provide sufficient bandwidth for a good reverberationlimited 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).
245
100 m
HO kHz
40 kHz
20"
90"
30
2"
83 dB at 1000 m
43 dB
10 dB
45 dB
193 dB
100 ms
-10 dB
35 dB
-30 dB
10 dB
3 dB
(dB)
= 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
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. Therefore 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:
0.55 X 802) = 42 dB
3
Ahead beam uses the 60 staves in segments BandC
247
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-limited performance is given by
I2PL = SL
+ TS
2PL
193 - 30 - 35
+ 42
10 - 10
+ 3 = 153 dB
-30-
10+3=25dB
248
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
250
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
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
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:
253
Worked Examples
Assume that DI
=6
PL = 2010g Y + 0.9~
X
138 dB
Intercept range = 50 km
How does this compare with the detection range from the ship?
2PL = SL + TS - N
2PL=220+10-60+20+0PL = 20 log Y + 0 . 9 ~
X
10+3=183dB
= 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
PL=200+3-30-
PL = 20 log r + 16r X
I18 dB
Intercept range = 3 km
What is the probable detection range of the submarine by the torpedo?
2PL = SL + TS - N
25s
Reducing SL
Increasing frequency to limit range by absorption losses
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?
Assume DI
=6
= 40
dB:
256
Communications Sonars
257
258
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 correspondingly 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
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?
Let DI = 30 dB, 5 logd = 10 dB, N = 35 dB (SS4 at 50 kHz), B, = 50 kHz and
assume spherical spreading and absorption:
At50kHz
PL
= 2Olog4000
+ (16 X 4) = 136 dB
S L z 136-3O-t 3 5 + 1 0 + 4 7 = 1 9 8 d B
260
At 1OOkHz
+ (35 X 4) = 212 dB
SL = 212 - 30+35 + 10 + 4 7 = 274 dB
PL = 2010g4000
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 successfully. 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
+ (16 X 1) = 76 dB
SL = 76 - 30 + 35 + 10 + 47 = 138 dB
PL = 2010g 1000
At 1OOkHz
+ (35 X 1) = 95 dB
30 + 35 + 10 + 47 = 157 dB
PL = 2010g 1000
SL = 95 -
Now the difference in SL between the ends of the band is only 19 dB and easy to
equalize.
Fortunately, a modest range is often acceptable: it is easy for a quiet submarine to
detect and locate a relatively noisy surface vessel and then approach to, say, 500 m
before establishing communication. Both vessels should reduce their source levels to
the minima required for reliable communication (perhaps 10 dB greater than the usual
values determined on detection criteria alone).
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 interactions between the initial requirements and the limitations imposed by platform
size and the environment will be demonstrated by considering practical systems.
262
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.
263
8h
= 8"
Vertical beamwidth 8,
12"
10log(2h/l) = 9 dB
+ 171 + DIt
= 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 transmitting 25 000/600 = 40 W of acoustic power. Suppose the radiating surface of each
) ~ , the total radiating surface is 600r~(0.05)~
= 4.7 m2.
element to be ~ ( 0 . 0 5 then
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
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.
cw
DI (receive)
Target strength, TS
5 logd
5 log n
Background noise
Value
Units
5.5 kHz
1
S
224
dB
750
Hz
2
Hz
8 degrees
20
-6
25
10
12
3
55
-33
dB
dB
dB
dB
dB
dB
dB
dB
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.
265
FM bandwidth
The FM bandwidth is chosen to resolve the structure of targets and non-targets, as
discussed in Chapter 9.
Reverberation indices
For FM, 10 log(750/8) = 20 dB
0
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
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
The noise-limited performance of the sonar is given by the following equation:
I2PL = SL + TS - N + DI
+ 10 - 55 + 25 + 0
12
+ 3 = 195 dB
12+3=13dB
267
268
Table 13.2
Parameter
Frequency
Pulse length, T
Source level, SL
FM bandwidth
CW bandwidth
Horizontal beamwidth,
Reverberation index, RI
FM
cw
DI (receive)
Target strength, TS
5 log d
5 log n
Background noise
sb
Value
Units
kHz
2
224
300
dB
Hz
Hz
3* degrees
1
20
-2
18
10
12
3
70
-33
dB
dB
dB
dB
dB
dB
dB
dB
*Maximum
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
The FM bandwidth is now limited to 300 Hz by the available bandwidth for an
acceptable projector transmit efficiency. This is less than optimum to resolve the
structure of both targets and non-targets, but it is acceptable.
Reverberation indices
= 20 dB
Target strength
The target strength (TS) is an average value for a typical submarine.
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
Perjormance
The noise-limited performance of the sonar is given by the following equation:
2PL = 224
+ 10
70 + 18
+ 3 - 12 + 3 = 176 dB
+ 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
= +3
Increase SL
= +3 dB
dB
$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:
+ TS - 2PL
Reverb level = RL = SL + (Sb + 10 log A)
Echo level
= EL =
SL
2PL
DI
Using the parameters fiom Table 13.2 and assuming spherical spreading and
negligible absorption, here is what we obtain:
EL
= 234 - 40 logR
RLCW 224 - 33
R L ~ ~ = 2 2 4 - 3 3 +lOlog(O.13R)-40l0gR=
30 log R
182-3OlogR
272
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 -
----I
10
20
..
--
Noise level
30
50
100
Range (km)
Figure 13.2 Echo, noise and reverberation levels plotted against range
273
274
cw
DI (receive)
Target strength, TS
5 log d
5 log n
Background noise
sb
Value
Units
5.5
100
kHz
ms
224
750
20
8
dB
Hz
20
-6
25
-15
12
3
55
-33
Hz
degrees
dB
dB
dB
dB
dB
dB
dB
dB
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
The noise-limited performance of the sonar is given by the following equation:
I2PL = SL + TS - N
2PL = 224 - 15 - 55 + 25 - 10 - 12
+ 3 = 160 dB
1 10log R = 10 log(B/&) -
S b - 41
+ TS - 5 log d + 5 log nl
276
+ 3 = - 12 dB and the
+ 3 = -38
dB and the
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
278
Conclusion
Conclusion
279
Solutions to Problems
Problem 1.1
I
= p 2 / p c = (10-4)2/(1.5 x
lo6) = 0.67 x
W/m2
+ 170.8 + DIt
= IOl0g40000 + 171 + 15
SL = 101og P
= 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
5000/75 = 67 elements. At 8 kHz we have
DI = 3
= 27
rows of
Alternatively
DI = 3
The trivial difference is due to the element numbers being the nearest integer values.
Problem 2.4
The DI of a baffled cylinder is given by
DI = 10 log 5hdfi = 10 log(2 X 3
X 5 2 ) = 22
dB
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
At 20 kHz 80 = 20log r + ( 3 r X
therefore
= 3.3 km
+ (0.3r X
At 20 kHz 80 = 10 log r + (3r X
At 5 kHz 80 = lolog
therefore r = 100 km
therefore r = 13 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):
At 10 kHz TS = IOlog[(0.5 X 22)/(2
0.15)] = -2 dB
284
Solutions to Problems
Problem 4.3
The target strength of any plate at an angle 8 to normal is given by
TS = 10 log(A/A)2
+ 20 log(&
sin x ) + 20 log(cos 8)
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
20 log v = sh
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
Problem 6.1
Start with
285
A = (1500 X 10-'/2)
TSR
-40
4000
0.17 = 51 000 m2
+ 47 = 7 dB
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,:
For 5logd = 6 dB, Pfa= 2
For 5 log d
r=
5 dB, I>f, 1
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:
3000 HZ
400 HZ
200 HZ
80 HZ
At 40 HZ
At
At
At
At
DI = 3
DI = 3
DI = 3
DI = 3
DI = 3
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
PL, a = 0.2 (dB)
128
R (km)
>loo
98
35
93
25
286
Solutions to Problems
Narrowband
PL(dB)
R (km)
89
30
86 59 56 64 61
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
PL = SL - N
Therefore
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.
Af = 2
The total shift is 25 Hz, therefore the shift due to the target is 6.2 Hz:
6.2 = 2
(S cos cp)
(5000/1500)
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
(SCOSX
~ ~()f / ~ )
120 = 2
(S cos cp)
(4000/1500)
(S cos cp)
(4000/1500)
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
1.9 = 2
288
Solutions to Problems
At 1000 Hz
1.9 = 2
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
For noise-limited detection,
2PL = SL + TS - N
289
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
PL = 115 = 20 log r
+7r X
(dB)
Problem 11.1
Noise limited
The noise-limited active sonar equation is
2PL = SL + TS - N
2PL = 192 - 25 - 30
+ 43 - 20 - 10 + 3 == 153 dB
therefore PL = 77 dB:
PL = 77 = 2010g r + ur
( a = 50 dB/km)
giving r = 470 m.
2PL = 192 - 25
30 + 43 - 39 - 10 + 3 = 134 dB
therefore PL = 67 dB:
PL=67=201ogr+arX
giving Y = 330 m.
(a=50dB/km)
290
Solutions to Problems
Reverberation limited
The reverberation-limited active sonar equation is
10 log R = 10 log(B/&)
- 11
+ TS - 5 log d + 5 log TZ
10+3=36dB
lOlogR=39+30-
11 - 2 5 -
10+3=26dB
Problem 11.2
Use the approximate formula for absorption attenuation, a
spherical spreading.
Frequency (kHz)
a (dB/km)
PL (dB)
SL (dB)
60
15
70
205
80
23
76
205
= 0.05f'.4,
and
100
32
84
205
assume
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
100 dB
At 15 kHz, a = 1.6and
PL = 2 0 1 0 g ~ +1 . 6 X
~ lo-'
=
Y =
100 dB
12000 m
= 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 = SL + TS - N
5000 m
Solutions to Problems
292
+ 1.6r X
104 dB
r = 13000m
As always, this (noise-limited) performance will only be achieved if the target
submarine shows sufficient doppler. Note that the 40 dB (100 times) increase in
transmitted power only increases the range by a factor of 2.6.
Index
Absorption of sound
mechanisms, 46
Active detection and classification
computer assisted detection and
classification
amplitude profiles, 2 13- 14
levels of, 2 11- 12
pulse features relevant to, 2 12- 13
statistical analysis, 2 13
effects of multipath on, 2 14- 18
unified approach to, 206-9
Active displays
CW surveillance, 204-5
FM surveillance, 203-4
general format, 202-3
history, range recorder, 205
Active sonar, xxiii
CWpulse, 163
Hamming pulse, 172-3
rectangular pulse, 171
reverb. rejection by, 174
doppler, 161-2
equations,
detection threshold, 165
noise limited, 165
reverb. limited, 166
FM pulse, 163
linear period, 169
linear frequency, 169
full beam processing, 180-9
half beam processing, 180-9
matched filter processing, 163-4
ping, 161
PRN pulse, 168
pulse types, 163
Index
Attenuation coefficient (continued)
table of values, 47
variation with:
frequency, 47
salinity, 47
temperature, 47
Azimuth, xvii
Baffles, 40- 1,92
Band level, 10
Bathythermograph, 5 1
errors, 64
Beamformers, 39
Beam pattern, 14
continuous line, 2 1-3
sidelobe levels and positions, 2 1-2
dipole, 17
effects of steering, 18- 19
effects of shading, 24-7
line array, 20- 1
Beamwidth
continuous line, 23
cylinder, 35
dipole, 17
line array, 19
Bearings only analysis, 148
Bottom bounce, 57,278
Bottom scattering strength, 1 14- 16
Boundary reverberation, 110- 1 1
Boundary roughness, 104-5
Cavitation, 6-7
threshold, 6
variation with depth, 6
variation with frequency, 6
variation with pulse length, 64
Cladding, 78-9
Communications sonar
bandwidth requirements, 257-8
examples of, 259-60
sonar equation, 258
Convergence gain, 56
Convergence zone, 56
Critical bandwidth of ear, 133-4
Decibel
defined, xviii
use in sonar equations, xix
Index
benefits, 180-1,208-9,278
comparison with FB processing, 180- 1.
209
detection performance, 187-9
signal processing
FM, 183-4
CW, 184-5
refutation of criticism, 278-9
Hearing
threshold frequency, I34
threshold of, 1 1
Helicopter sonar, 37, 177
Hodgeson model, 59-60
example, 6 1-3
Horizontal direct passive ranging, 153-5
Hull mounted ship sonars, 261 -7
performance, 266-7
representative design, 263 -6
Hydrophone, xvii
sensitivity, 9
lntegration factor, 124, 131, 139
Intensity of sound, 3
Interaction effects, 8
Intercept sonar
effectiveness, 257
function, 250- 1
overview, 249
reducing probability of exploitation, 256-7
reducing probability of intercept, 255-6
by sector transmissions, 255-6
sonar equation, 252
Isotropic spectrum level, 90
Large aperture array, half beam plots, 185-6
Layer scattering strength, 1 11-2
LOFARGRAM, 136
Longitudinal wave, I
Low frequency active sonar,
bandwidth limitations, 273
doppler requirements, 273
performance, 270-2
typical design, 267-9
295
dedicated design, 244-8
overview, 242
Mine hunting sonar
acoustic shadows, 240- 1
arrays for, 236-7
backgrounds to detection, 234-5
classes of mines, 234
classification, 240- 1
overview, 233-4
performance, 239
pulses for, 238
range and bearing resolutions, 235
Minimum discernable signal, 12I
Motional resistance, 84
Multipath propagation, 66,2 14-2 18
in deep water, 2 16- 17
in shallow water, 218
Netsonde, 22 1
Noise
addition of levels, 93-4
ambient, xvii, 86-8
biological, 87
directivity of, 93
electrical, 90- 1
factor, 94-8
flow, 89, 91,93
isotropic spectrum level, 90
machinery, 89, 91
propeller, 89, 92
radiated, 93, 125-8
rain, 87--8
self, xvii, 83, 89-93
shipping, 87-8
thermal. 84-5
vessel, 89-93
Noise factor
calculations, 97 - 100
of receiver, 94
of receiver & hydrophone, 95
Noise levels, practical values, 100
Normalisation, 131
Operational degradation, 20 1,277
Machinery noise, 9 1
Maximum response angle, 147, 154
Mine avoidance sonar,
adaption of A/'S sonar, 243-4
Particle velocity, 1, 1 1
Passive arrays, 132
surface ship hull mounted, I32
296
Passive arrays (continued)
submarine hull mounted, 132
towed arrays, 132, 156-60, 143-6
Passive displays, 135-6
broadband, 135
narrowband / LOFAR, 136
visual integration on, 135
Passive ranging, xviii, 146- 155
BOA, 148
HDPR, 153-5
triangulation, 147-8
VDPR, 149-52
Passive sonar, xvii, xxii, 125-60
aural, 133-4
broadband, 129-30,135, 137-8, 141-6
complete system, 130
DEMON, 129
electronic integration, 135, 137
equations
broadband cross-correlator, 138
broadband energy, 137
narrowband amplitude, 140
narrowband power, 139-40,
narrowband, 129-31, 136, 139-40, 141-6
vernier processing, 129
Performance prediction, 1 19
accuracy, 64
Ping, 161
Ping interval, 161
Ping/pulse repetition rate, 161
Plankton, 10.1, 10.3
Pressure of sound, 1
Probability of detection, 120-2
Probability of false alarms, 120-2
Projector, xviii
interaction effects, 8
sensitivity, 9
Propagation loss, 43-47
allowable, xxii
definition, 43
models, 58-60
two way, xxiii
Propagation modes
bottom bounce, 57
convergence zone, 56
deep sound channel, 52-53
reliable acoustic path, 53-54
surface duct, 54-55
Index
Propeller noise, 92
Plane wave, I
Radar
cross section, 70
equation, xix
detection of submarines, xxi
Radiated noise, 93, 125-8
continuum, 126
line structure, 125
nature of, 126
representativevalues, 127-8
source level, 125
tonals, 126
Radiated power, 5
Radiation resistance, 84
Rain noise, 87-8
Range recorder, 205
Ranging
Active, xviii
Passive, xviii, 146-55
Ray theory, 58
Hodgeson model, 59-60
examples, 6 1-63
Receiver operating characteristics, 122-3
recommended values for 5 log d, 124
Reference intensity of sound, 3
Reliable acoustic path, 53-54
Remotely operated vehicle, 1 1.1
Reverberation, xviii
backscattering strength, 103, 106, 1 1 1- 17
boundary, 1 10- 1 1
in shallow water, 1 14- 16
index, 166, 168
reflection, 103-4
sea-surface, 113
sea-bed, 114-16
target strength of, 106, 166
under ice, 1 17
variation with frequency, 116
volume, 107-9, 167
Reverberating area, 1 10
Reverberating volume, 107-8
Reverberation spreading
beamwidth, 175
environment, 176
practical envelopes, 179
297
ship motion, 175-6
sidelobes, 178
Target strength
definition, 67
dependence on pulse type and length, 69
integrated TS, 68
measurement, 68
near resonance, 72
of mines, 72-3
of reverberation, 106, 166
of simple shapes, 70-1
of small targets, 72
of sphere, 69
of submarines, 74-6
of torpedoes, 73 -4
of towed arrays, 77
peak TS, 68
practical values, 80
reduction, 78-9
ThreshoId frequency of hearing, 134
Threshold of hearing, 1 1
Thermal noise, 84- 5
Torpedo active detection
optimum pulse, 276
using modified A/S sonar, 274-5
Torpedo passive detection, 141-5
Towed arrays, 132,143-6,156-60
bearing ambiguity, 158-9
cable, 156
curvature of, 156
heading sensors, I56
hose, 156
hydrophone triplets, 159
nested, 157
resolving bearing ambiguity, 158-9
self noise of, 159-60
tail, 156
target strength, 77
twin arrays, 159
VIM, 156
Transducer, xviii
298
Transverse wave, 1
Triangulation, 147-8
Under ice scattering strength, 1 17
Underwater telephone, 12.8
Vertical direct passive ranging, 149-52
Very long pulses, 200
Vessel noise, 89-93
Vibration isolation module, 156
Visual integration, 135, 137
Index
Volume reverberation, 107-9, 167
Volume scattering strength, 1 12
Wavefiont, f
Wave height, 87
Wavelength, xx
Wave theory, 58
Wind speed, 87
Yaw, 175