estimating underwater acoustic propagation acoustic...
TRANSCRIPT
![Page 1: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/1.jpg)
2/23/20062/23/2006 11
Estimating UnderwaterEstimating Underwater Acoustic PropagationAcoustic Propagation
Ethem Mutlu SEthem Mutlu SöözerzerResearch EngineerResearch Engineer
MIT Sea Grant College ProgramMIT Sea Grant College Program
![Page 2: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/2.jpg)
OutlineOutline
�� What is decibel?What is decibel?�� Transducers and hydrophonesTransducers and hydrophones�� Underwater acoustic propagationUnderwater acoustic propagation�� Ray tracingRay tracing�� Delay and signal strength calculationsDelay and signal strength calculations�� Channel impulse responseChannel impulse response
�� Estimating the range of a sourceEstimating the range of a source�� Estimating the direction of a sourceEstimating the direction of a source
2/23/20062/23/2006 2 2
![Page 3: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/3.jpg)
DeciDeci--belbel (dB)(dB)
�� A decibel (dB) is a unit for measuring theA decibel (dB) is a unit for measuring the relative strength of a signal in logarithmic scalerelative strength of a signal in logarithmic scale
�� P(dBP(dB) = 10 log10(P/Pr)) = 10 log10(P/Pr) = 20 log10(V/Vr)= 20 log10(V/Vr)
�� V(dBV(dB) = 10 log10(V/Vr)) = 10 log10(V/Vr) = 10 (log10(V)= 10 (log10(V) –– log10(Vr))log10(Vr))
2/23/20062/23/2006 3 3
![Page 4: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/4.jpg)
Transducer and HydrophoneTransducer and Hydrophone SpecificationsSpecifications
�� Open Circuit Receiving Response (OCRR)Open Circuit Receiving Response (OCRR)�� Transmitting Voltage Response (TVR)Transmitting Voltage Response (TVR)�� Directionality PatternDirectionality Pattern�� PrePre--AmplifierAmplifier
2/23/20062/23/2006 4 4
![Page 5: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/5.jpg)
Transmitting Voltage ResponseTransmitting Voltage Response (TVR)(TVR)
�� output SIL generated per 1 Voutput SIL generated per 1 V of input Voltage at 1m rangeof input Voltage at 1m range as a function of frequencyas a function of frequency
�� dB redB re µµPaPa / V/ V
V = 200 V @@ fcfc=22 kHz=22 kHz
SIL (µµPaPa) = V (V) x TVR (µµPaPa / V)/ V)
SIL = 10 log10(V) + TVR(fc)
ITC 1001 Transducer VTR = 26 + 144 = 170 dB re µµPaPa
2/23/20062/23/2006 5 5
110
2 10
120
130
140
150
18Frequency in kHz
dB re
µPa
/V @
1m
26 34
Transmitting Voltage Response
Figure by MIT OCW.
![Page 6: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/6.jpg)
Open Circuit Receiving ResponseOpen Circuit Receiving Response (OCRR)(OCRR)
�� output voltage (V) generatedoutput voltage (V) generated by the transducer perby the transducer per µµPaPa ofof sound pressure as a functionsound pressure as a function of frequencyof frequency
�� dB re 1V /dB re 1V / µµPaPa
SIL = 190 dB re µµPaPa @@ fcfc=22 kHz=22 kHz
V = SIL (V/µµPaPa) x OCRR (µµPaPa))
VdB = SIL + OCRR(fc)
ITC 1001 Transducer OCRR = 190 + (-190) = 0 dB re 1V
VdB re 1V = 10 log10 ( V / 1 )
V = 10( VdB / 10 ) = 1 V
2/23/20062/23/2006 6 6
102
-220
-210
-200
-190
-180
18 26 34Frequency in kHz
dB re
1V
/µPa
Open Circuit Receiving Response
Figure by MIT OCW.
![Page 7: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/7.jpg)
Directionality PatternDirectionality Pattern
�� ITC 1001 spherical transducerITC 1001 spherical transducer �� ITC 2010ITC 2010 toroidaltoroidal transdtransducerucer
�� Uniform resUniform respponse oonse ovver aller all �� More gain over the sidesMore gain over the sides angles ( 0 to 2angles ( 0 to 2ππ) on both) on both (horizontal plane) than the(horizontal plane) than the horizontal and vertical planehorizontal and vertical plane over the top and bottomover the top and bottom
(vertical plane)(vertical plane)
2/23/20062/23/2006 7 7
Directivity Pattern at 2.0 kHz
240
10dB/div
270
300
120
90
60
210 180 150
330 0 30
Directivity Pattern at 18.0 kHz
240
10dB/div
270
300
120
90
60
210 180 150
330 0 30
Figures by MIT OCW.
![Page 8: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/8.jpg)
2/23/20062/23/2006 88
PrePre--AmplifierAmplifier
Amplifies the signals generated at the receiving element Amplifies the signals generated at the receiving element (hydrophone or transducer)(hydrophone or transducer)Gain defined in dBGain defined in dB
mV
mV
amplifier
V
pre-amplifier
Vin = 1 mV, Gain = 20 dB
Vout(dB) = Vin(dB) + Gain(dB)
Vin(dB) = 10 log10(Vin) = 10 log(10e-3) = -30 dB
Vout(dB) = -30 + 20 = 0 dBVout = 10(-10/10) = 0.1 Volt
![Page 9: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/9.jpg)
Shallow Water PropagationShallow Water Propagation
�� Assumptions:Assumptions:�� CoConstant sound speed (c = 1500nstant sound speed (c = 1500 m/sm/s))�� Surface and bottom are smoothSurface and bottom are smooth
r=100m
d=20m
h=80m
source destination
θ2 θ2
θ1 θ1
2/23/20062/23/2006 9 9
![Page 10: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/10.jpg)
Length of propagation pathsLength of propagation paths Direct path => d0 = 100mDirect path => d0 = 100mSurface reflection => d1 = 2d/cos(Surface reflection => d1 = 2d/cos(θθ1) = 107.7 m1) = 107.7 m
θθ1 = atan(r/2d)1 = atan(r/2d)Bottom reflection => d2 = 2h/cos(q2) = 188.7 mBottom reflection => d2 = 2h/cos(q2) = 188.7 m
θθ2 = atan(r/2h)2 = atan(r/2h)SBS reflection => d3 = 2(2d/cos(SBS reflection => d3 = 2(2d/cos(θθ3)+ h/cos(3)+ h/cos(θθ3)) = 260 m3)) = 260 mBSB reflection => d4 = (2d/cos(BSB reflection => d4 = (2d/cos(θθ4)+ 2(h/cos(4)+ 2(h/cos(θθ4))) = 399.5 m4))) = 399.5 m
2/23/20062/23/2006 1010
r=100m
d=20m
h=80m
source destination
θ2 θ2
θ1 θ1
![Page 11: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/11.jpg)
2/23/20062/23/2006 1111
Time of arrival to the receiverTime of arrival to the receiver
�� ττii == ddii/c/c�� ττ0 = 66.70 = 66.7 msecmsec�� ττ1 = 71.81 = 71.8 msecmsec�� ττ2 = 125.82 = 125.8 msecmsec�� ττ3 = 173.33 = 173.3 msecmsec�� ττ4 = 266.34 = 266.3 msecmsec
![Page 12: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/12.jpg)
Transmission lossTransmission loss �� fcfc=22kHz=22kHz�� T = 15T = 15 ˚̊C (fmC (fm==100100 kcycleskcycles/sec, A=6e/sec, A=6e--4, B=2.4e4, B=2.4e--7)7)�� a = (A fm f2)/(f2+fm2)+Bf2 dB/ma = (A fm f2)/(f2+fm2)+Bf2 dB/m�� Surface reflection loss (Surface reflection loss (RLsRLs) = 1 dB) = 1 dB�� Bottom reflection loss (Bottom reflection loss (RLbRLb) = 3 dB) = 3 dB
TL = TLs + TLa + RLs + RLb TLi = 20log(di) + adi + RLs + RLb
TL0 = 40.3 dB TL1 = 42.0 dB TL2 = 49.1 dB TL3 = 54.0 dB TL4 = 60.2 dB
2/23/20062/23/2006 1212
![Page 13: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/13.jpg)
2/23/20062/23/2006 1313
Impulse ResponseImpulse Response
�� output of a system to anoutput of a system to an excitation with a unit impulseexcitation with a unit impulse
⎩ ⎨ ⎧ =
=)( otherwise n
n 0
11δ
r=100m
d=20m
h=80m
source destination
θ1 θ1
θ2 θ2
![Page 14: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/14.jpg)
Transmitted and Received SILTransmitted and Received SIL
�� ITC1001,ITC1001, fcfc=22kHz, V=100 V=22kHz, V=100 V�� SILsSILs = 144 + 10 log(100) = 170 dB re= 144 + 10 log(100) = 170 dB re mPamPa at 1mat 1m
�� SILriSILri == SILsSILs –– TLiTLi�� SILr0 = 129.7 dBSILr0 = 129.7 dB�� SILr1 = 128.0 dBSILr1 = 128.0 dB�� SILr2 = 120.9 dBSILr2 = 120.9 dB�� SILr3 = 116.0 dBSILr3 = 116.0 dB�� SILr4 = 109.8 dBSILr4 = 109.8 dB
2/23/20062/23/2006 1414
![Page 15: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/15.jpg)
Received VoltageReceived Voltage
�� OCRR =OCRR = --162 dB re V /162 dB re V / µµPAPA�� PrePre--amplifier gain, G = 40 dBamplifier gain, G = 40 dB�� VdBVdB == SILriSILri --162 + 40162 + 40�� V = 10(VdB/10)V = 10(VdB/10)
�� V0 = 5.9 VV0 = 5.9 V�� V1 = 4.0 VV1 = 4.0 V�� V2 = 0.8 VV2 = 0.8 V�� V3 = 0.2 VV3 = 0.2 V�� V4 = 0.1 VV4 = 0.1 V
2/23/20062/23/2006 1515
![Page 16: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/16.jpg)
2/23/20062/23/2006 1616
Results for 1000mResults for 1000m
![Page 17: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/17.jpg)
Determining the Range of a SourceDetermining the Range of a Source
�� Tracker sends a pulse,Tracker sends a pulse, p(tp(t) = A sin(2) = A sin(2ππffc1c1t), 0<t<Tst), 0<t<Ts�� Target replies, p1(t) = A sin(2Target replies, p1(t) = A sin(2ππffc2c2(t(t--ττpp--ττtt))))�� Tracker receives, p2(t) = A sin(2Tracker receives, p2(t) = A sin(2ππffc2c2(t(t--ττpp--ττtt--ττpp))))�� How can we measureHow can we measure ττpp++ττtt++ττ ??pp
2/23/20062/23/2006 1717
![Page 18: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/18.jpg)
2/23/20062/23/2006 1818
CorrelationCorrelation
�� Definition:Definition:
�� Autocorrelation:Autocorrelation:
( ) ∫ ∞
∞−
−= λ)dta(t)b(tλR b)(a,
( ) ∫ ∞
∞−
−= λ)dta(t)a(tλR a)(a,
autocorrelation
![Page 19: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/19.jpg)
Estimating the RangeEstimating the Range
�� CorrelateCorrelate p(tp(t) with) with p(tp(t--ττpp--ττtt--ττpp))�� Find the peak of the correlation,Find the peak of the correlation, λλ�� λλ = 2= 2ττpp--ττtt
�� ττpp is the propagation delayis the propagation delay�� Range is, d=cRange is, d=c ττpp = 1500= 1500 ττpp
2/23/20062/23/2006 1919
![Page 20: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/20.jpg)
2/23/20062/23/2006 2020
Correlation ResultsCorrelation Results
100m, delay estimate is 66.7 msec 100m, delay estimate is 666.7 msec
![Page 21: Estimating Underwater Acoustic Propagation Acoustic Propagationdspace.mit.edu/bitstream/handle/1721.1/74140/2-017j... · 2019. 9. 13. · Underwater acoustic propagation ... V = 200](https://reader035.vdocument.in/reader035/viewer/2022071611/614a2d6412c9616cbc693f2e/html5/thumbnails/21.jpg)
Determining the Direction of theDetermining the Direction of the TargetTarget
θθ r1-r4
d
Quadrant 1
Quadrant 2Quadrant 3
Quadrant 4
Η1
Η2Η3
Η4
�� Four hydrophonesFour hydrophones�� Measure delay atMeasure delay at
each hydrophoneeach hydrophone�� Compare delay pairsCompare delay pairs
((ττ11,, ττ22), (), (ττ22,, ττ33),), ((ττ33,, ττ44), (), (ττ44,, ττ11) to find) to find which quadrantwhich quadrant
�� Estimate the angleEstimate the angleθ = sign(r1-r4)acos( |r1-r4| / d)
2/23/20062/23/2006 2121