CN104502627B - Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP - Google Patents
Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP Download PDFInfo
- Publication number
- CN104502627B CN104502627B CN201410836271.0A CN201410836271A CN104502627B CN 104502627 B CN104502627 B CN 104502627B CN 201410836271 A CN201410836271 A CN 201410836271A CN 104502627 B CN104502627 B CN 104502627B
- Authority
- CN
- China
- Prior art keywords
- pulse
- delay
- msub
- mrow
- doppler
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses solve method based on transmitting Design of Signal and the phase ambiguity of processing in a kind of ADCP, it is characterised in that：It is delayed based on different transmitting signal, mainly solves Doppler and be concerned with phase fuzzy problem present in flow measurement mode, when carrying out flow velocity measurement using relevant flow measurement mode, change the delay between pulse, even the delay τ between pulse 1 and pulse 21Delay τ between pulse 2 and pulse 32Difference, if the Doppler difference measured twice is farther out, there are phase ambiguity, carries out Doppler frequency shift modification by frequency search at this time, and then obtain accurate radial flow speed.The present invention solves phase fuzzy problem, so as to improve the maximum detection flow velocity of relevant flow measurement mode and broadband flow measurement mode.
Description
Technical field
The present invention relates to a kind of phase ambiguity solution for being used for the relevant flow measurement mode of acoustic Doppler fluid velocity profile instrument (ADCP)
Certainly method, belongs to underwater sound signal field of measuring technique.
Background technology
Acoustic Doppler fluid velocity profile instrument (ADCP) is a kind of sonar that tests the speed to be worked according to underwater sound doppler principle, mainly
For remote measurement interior current speed in a big way.Experiment proves in seawater that there are substantial amounts of scattering object such as fine particle, swim
Biology and bubble etc., they are flowed with seawater.These scattering objects and seawater combine together, their speed represents sea
The speed of stream.Under this premise, launch sound wave into seawater as a kind of ADCP of active sonar, by above-mentioned scattering object
Sound wave is scattered, ADCP is received and handled to echo-signal again.According to doppler principle, due to ADCP and scattering object
Between there are relative motion, launch and there is a Doppler frequency between sound wave and scatter echo frequency.By measuring this
Doppler frequency can directly calculation go out the relative velocity of ADCP and scattering object.If if ocean current is divided into dried layer, pass through survey
Amount can be obtained by fluid velocity profile of each layer of ocean current with respect to ADCP.
Echo signal processing algorithm is an importance of ADCP flow measurements research.For returning for different transmitting signal forms
Ripple sampled data, also has different flow measurement signal processing algorithms.Arrowband ADCP (Kuroda, Y.Development of a
shipboard acoustic Doppler current profiler,OCEANS'88.A Partnership of Marine
Interests.Proceedings) general transmitting single-frequency wide pulse signal, make use of non-between pulse when carrying out flow velocity measurement
Relevant flow measurement mode, its main signal Processing Algorithm is complex self-correlation algorithm.Then, in order to improve spatial resolution without damage
Lose measurement accuracy, someone make use of the flow measurement mode that is concerned between pulse, the echoes of two or more coherent pulses that pass a test it
Between phase difference measure flow velocity.This mode improves spatial resolution compared with non-coherent fashion, but which is larger in flow velocity
When there are phase ambiguity phenomenon to need to solve, the signal processing algorithm of pulse coherence flow measurement mode is also complex self-correlation algorithm.By
When complex self-correlation algorithm carries out Doppler-frequency estimation, Doppler frequencyWherein, τ is signal time delay, and φ (τ) is back
The phase of ripple signal complex autocorrelation, works as fdBe wide enough so that φ (τ) beyond (- π, π] scope when, phase ambiguity will be produced and showed
As.Equally, in the wideband encoding Doppler current measured technology proposed later, there is also phase ambiguity phenomenon to need to solve.
The content of the invention
Goal of the invention：For problems of the prior art and deficiency, the present invention provides one kind and is used for acoustic Doppler
The phase ambiguity of the relevant flow measurement mode of fluid velocity profile instrument (ADCP) solves method, and this method can effectively solve coherent manner and broadband
Phase ambiguity phenomenon existing for mode.
Technical solution：A kind of phase ambiguity based on transmitting Design of Signal and processing in ADCP solves method (with relevant side
Exemplified by formula), it is delayed based on different transmitting signal, mainly solves Doppler's phase ambiguity present in flow measurement mode that is concerned with and ask
Topic, when carrying out flow velocity measurement using relevant flow measurement mode, changes the delay between pulse, even prolonging between pulse 1 and pulse 2
When τ1Delay τ between pulse 2 and pulse 32Difference, if the Doppler difference measured twice is farther out, there are phase mode
Paste, carries out Doppler frequency shift modification by frequency search at this time, and then obtains accurate radial flow speed；Mainly include following step
Suddenly：
(1) pulse coherence flow measurement mode is used, the one-shot measurement cycle includes 3 pulses, and between pulse 1 and pulse 2
It is delayed as τ1=n1T, is determined by maximum detection depth H：C is 1500m/s in formula, and θ is transducer beams direction
With the angle of vertical direction, T is pulse width, n1For integer；Delay between pulse 2 and pulse 3 is τ2=n2T>τ1, n2To be whole
Number.
(2) by maximum detection Doppler frequency fmaxIt is calculated：
(3) it is f that the echo produced by pulse 1 and pulse 2, which calculates Doppler frequency,d1, returned by what pulse 2 and pulse 3 produced
It is f that ripple, which calculates Doppler frequency,d2.Search for m1And m2FindMinimum value so that obtain it is accurate
Doppler frequency shiftPhase fuzzy problem, wherein m can so be eliminated1∈[0,n1]、m2∈[0,n2], and
m1、m2For integer.
(4) by Doppler frequency shift fdAnd applyRadial flow speed v is calculated, wherein, λ is transmitting impulse wave
Long, f is transmitting pulse frequency.
Beneficial effect：The present invention using relevant flow measurement mode and broadband flow measurement mode when carrying out flow velocity measurement, by changing
Become the delay between pulse, obtain 2 Doppler frequency shifts there are phase ambiguity, obtained using frequency search method really more
General Le frequency displacement, solves phase fuzzy problem, so as to improve the maximum detection flow velocity of relevant flow measurement mode and broadband flow measurement mode.
Brief description of the drawings
Fig. 1 is for nonperiodic signal as the transmitting signal of ADCP to solve the schematic diagram of phase fuzzy problem；
Fig. 2 is the specific implementation process of present example.
Embodiment
With reference to specific embodiment, the present invention is furture elucidated, it should be understood that these embodiments are merely to illustrate the present invention
Rather than limit the scope of the invention, after the present invention has been read, various equivalences of the those skilled in the art to the present invention
The modification of form falls within the application appended claims limited range.
As shown in Figure 1, the one-shot measurement cycle is 3 pulses, transmitting signal pulsewidth is T, prolonging between pulse 1 and pulse 2
When be n1T, the delay between pulse 2 and pulse 3 are n2T, respectively believes the echo of pulse 1 and pulse 2 using complex self-correlation algorithm
Number and pulse 2 and the echo-signal of pulse 3 handled, obtain 2 Doppler frequency fd1And fd2, search again for m1And m2So thatF at this timedAs correct Doppler frequency.Finally by Doppler frequency shift fdCalculate footpath
To flow velocity v, using the following formulaWherein λ is transmitting pulse wavelength, and f is transmitting pulse frequency.
Fig. 2 is the specific implementation process based on transmitting Design of Signal with the bilevel Linear programming of processing.
Measurement result under the different radial flow speeds of table 1
Table 1 is a specific embodiment according to Fig. 2.F in table 1maxFor maximum detection Doppler frequency, vmaxCan for maximum
Survey the maximum detection flow velocity corresponding to Doppler frequency, fdTo need the standard Doppler frequency measured, v fdCorresponding radial direction
Flow velocity.Specific measuring condition is as follows：
A) flow depth is 6.5m, and current thickness is 0.065m, signal-to-noise ratio 30dB, velocity of sound 1500m/s, estimates performance
It is by 100 statistical results averagely obtained；
B) transmitting carrier frequency is 600kHz, and the sample frequency of receives echo-signal is 1.5MHz, and beams incident angle is 30 °；
C it is single carrier frequency rectangle CW pulse signals) to launch signal, signal length be it is corresponding with the current thickness of 0.065m,
Specific to should be 0.1ms, signal delay is determined by flow depth, the delay τ between pulse 1 and pulse 21For 100T, pulse 2 and arteries and veins
Delay τ between punching 32By τ1And fmaxDetermine.
If without using Ambiguity Technique, maximum detection Doppler frequency is only 50Hz.In table 1, how general maximum detection is
Strangle frequency to have clear improvement, 5kHz, 10kHz, 24kHz are got in embodiment.Therefore the invention can effectively solve pulse coherence survey
The phase fuzzy problem of stream mode, and then obtain correct radial flow speed information.
Claims (1)
1. method is solved based on transmitting Design of Signal and the phase ambiguity of processing in a kind of ADCP, it is characterised in that：Based on difference
The delay of transmitting signal, solve Doppler and be concerned with phase fuzzy problem present in flow measurement mode, using relevant flow measurement mode into
When row flow velocity measures, change the delay between pulse, specially make the delay τ between pulse 1 and pulse 21With pulse 2 and pulse
Delay τ between 32Difference, if the Doppler difference measured twice is farther out, there are phase ambiguity, passes through frequency at this time
Search carries out Doppler frequency shift modification, and then obtains accurate radial flow speed；
(1) pulse coherence flow measurement mode is used, the one-shot measurement cycle includes 3 pulses, and the delay between pulse 1 and pulse 2
For τ1=n1T, is determined by maximum detection depth H：C is 1500m/s in formula, and θ is transducer beams direction with hanging down
Nogata to angle, T is pulse width, n1For integer；Delay between pulse 2 and pulse 3 is τ2=n2T>τ1, n2For integer；
(2) by maximum detection Doppler frequency fmaxIt is calculated：
<mrow>
<msub>
<mi>n</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msubsup>
<mi>n</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>Tf</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>n</mi>
<mn>1</mn>
</msub>
<msub>
<mi>Tf</mi>
<mi>max</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
(3) it is f that the echo produced by pulse 1 and pulse 2, which calculates Doppler frequency,d1, the echo meter that is produced by pulse 2 and pulse 3
Calculation Doppler frequency is fd2；Search for m1And m2FindMinimum value so that obtain accurately it is more
General Le frequency displacementPhase fuzzy problem, wherein m can so be eliminated1∈[0,n1]、m2∈[0,n2], and m1、m2
For integer；
(4) by Doppler frequency shift fdAnd applyRadial flow speed v is calculated, wherein, λ is transmitting pulse wavelength, and f is
Launch pulse frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410836271.0A CN104502627B (en) | 2014-12-29 | 2014-12-29 | Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410836271.0A CN104502627B (en) | 2014-12-29 | 2014-12-29 | Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104502627A CN104502627A (en) | 2015-04-08 |
| CN104502627B true CN104502627B (en) | 2018-04-27 |
Family
ID=52944039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410836271.0A Active CN104502627B (en) | 2014-12-29 | 2014-12-29 | Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104502627B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105049045B (en) * | 2015-08-12 | 2018-04-17 | 江苏中海达海洋信息技术有限公司 | Low noise high-precision signal processing system based on ADCP |
| CN106248991B (en) * | 2016-07-15 | 2018-12-21 | 南京世海声学科技有限公司 | A kind of ADCP solution under low signal-to-noise ratio tests the speed blur method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0573431A1 (en) * | 1990-09-26 | 1993-12-15 | Rowe Deines Instr Inc | Acoustic doppler current profiler. |
| CN102226839A (en) * | 2011-06-07 | 2011-10-26 | 北京理工大学 | Estimation method for time delay of line scanning pulse with low sampling rate |
| CN103630706A (en) * | 2013-11-12 | 2014-03-12 | 方世良 | Method for acquiring radial direction water velocity in acoustic Doppler current profiler |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7007555B2 (en) * | 2004-03-15 | 2006-03-07 | Teledyne Rd Instruments, Inc. | System and method of horizontal wave measurement |
| US7539082B2 (en) * | 2006-09-28 | 2009-05-26 | Teledyne Rd Instruments, Inc. | System and method for acoustic Doppler velocity processing with a phased array transducer including using a wide bandwidth pulse transmission to resolve ambiguity in a narrow bandwidth velocity estimate |
| JP5980113B2 (en) * | 2009-05-27 | 2016-08-31 | テレダイン アールディー インスツルメンツ，インコーポレイテッド | System and method for determining wave characteristics from a mobile platform |
-
2014
- 2014-12-29 CN CN201410836271.0A patent/CN104502627B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0573431A1 (en) * | 1990-09-26 | 1993-12-15 | Rowe Deines Instr Inc | Acoustic doppler current profiler. |
| CN102226839A (en) * | 2011-06-07 | 2011-10-26 | 北京理工大学 | Estimation method for time delay of line scanning pulse with low sampling rate |
| CN103630706A (en) * | 2013-11-12 | 2014-03-12 | 方世良 | Method for acquiring radial direction water velocity in acoustic Doppler current profiler |
Non-Patent Citations (2)
| Title |
|---|
| ADCP数据采集和信号处理的分析与设计;殷俊琳;《中国优秀硕士学位论文全文数据库信息科技辑》;20060115(第1期);5-13 * |
| 宽带声学多普勒流速剖面仪的中频正交采样算法;彭东立 等;《声学技术》;20130228;第32卷(第1期);15-18 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104502627A (en) | 2015-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100195576B1 (en) | Apparatus for measuring the velocity of moving body | |
| CN103630706B (en) | A kind of method obtaining radial water Flow Velocity in acoustic Doppler fluid velocity profile instrument | |
| CN104126126A (en) | Ship speed meter and ship speed measurement method | |
| CN102183741B (en) | Method for realizing high frame rate unambiguous positioning in long baseline asynchronous acoustic positioning system | |
| Murakawa et al. | A dealiasing method for use with ultrasonic pulsed Doppler in measuring velocity profiles and flow rates in pipes | |
| CN105785067A (en) | Seawater flow velocity measurement method based on sonic propagation time difference | |
| CN104502627B (en) | Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP | |
| US8437222B2 (en) | System and method of range estimation | |
| CN107783137B (en) | Sound Doppler and acoustic correlation velocity measurement method based on five-beam configuration | |
| Zedel et al. | Resolving velocity ambiguity in multifrequency, pulse-to-pulse coherent Doppler sonar | |
| RU2626380C1 (en) | Selection system of moving targets with measurement of range, radial velocity and direction of motion | |
| Cui et al. | Strong scattering targets separation based on fractional Fourier transformation in pulse-to-pulse coherent acoustical Doppler current profilers | |
| CN109164441A (en) | A kind of method of radar range finding | |
| Jaafar et al. | Velocity and turbulence measurements by ultrasound pulse Doppler velocimetry | |
| KR20060135715A (en) | Device for avoiding obstacles for high-speed multi-hulled watercraft | |
| Boltryk et al. | An ultrasonic transducer array for velocity measurement in underwater vehicles | |
| KR101480170B1 (en) | Method for estimaing velocity of underwater vehicles by acoustic doppler measurement and apparatus for estimaing velocity of underwater vehicles using thereof | |
| Park et al. | Array tilt effect induced by tidal currents in the northeastern East China Sea | |
| Mosquera et al. | Salinity estimation from acoustic Doppler velocimeter measurements | |
| LIU et al. | Combined method of conventional and coherent Doppler sonar to avoid velocity ambiguity | |
| CN106248991A (en) | ADCP solution under a kind of low signal-to-noise ratio tests the speed blur method | |
| Hui et al. | Reverberation Channel | |
| CN111337881B (en) | Underwater target detection method utilizing propeller noise | |
| JP7429448B2 (en) | underwater detection device | |
| Goncharov et al. | Local acoustic tomography on shelf of the Black Sea |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |