CN112146744A - Method for accurately measuring reciprocal constant of reverberation pool - Google Patents

Method for accurately measuring reciprocal constant of reverberation pool Download PDF

Info

Publication number
CN112146744A
CN112146744A CN202010940592.0A CN202010940592A CN112146744A CN 112146744 A CN112146744 A CN 112146744A CN 202010940592 A CN202010940592 A CN 202010940592A CN 112146744 A CN112146744 A CN 112146744A
Authority
CN
China
Prior art keywords
reverberation
pool
reciprocal
signal
constant
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.)
Granted
Application number
CN202010940592.0A
Other languages
Chinese (zh)
Other versions
CN112146744B (en
Inventor
杨柳青
李琪
黄勇军
陈毅
尚大晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
715th Research Institute of CSIC
Original Assignee
715th Research Institute of CSIC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 715th Research Institute of CSIC filed Critical 715th Research Institute of CSIC
Priority to CN202010940592.0A priority Critical patent/CN112146744B/en
Publication of CN112146744A publication Critical patent/CN112146744A/en
Application granted granted Critical
Publication of CN112146744B publication Critical patent/CN112146744B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses an accurate measurement method for a reciprocal constant of a reverberation pool, which mainly comprises a hydrophone and a reciprocal transducer, wherein the hydrophone and the reciprocal transducer respectively move on two sides of the pool, collect sound field information of the reverberation pool, and calculate the reciprocal constant of the reverberation pool by utilizing the reciprocity of the reciprocal transducer and combining the collected sound field information. The invention has more used space information, lower requirements on the sound source level of the transmitting transducer and higher precision of the test result; the method is used for measuring the reciprocal constant of the reverberation water pool in the research subject, the measurement result is basically consistent with the result of the reverberation time method, the effectiveness of the method is proved, and the method can be gradually applied to the parameter calibration of each reverberation field.

Description

Method for accurately measuring reciprocal constant of reverberation pool
Technical Field
The invention relates to the field of acoustic metering test technology, in particular to an accurate measurement method for a reciprocal constant of a reverberation water pool.
Background
The underwater acoustic measurement and test usually adopts free field condition, i.e. the condition of isotropy and no interface reflection, and usually needs to combine with impulse sound technology to reach the ideal free field. Therefore, a large anechoic pool needs to be built for free field calibration at lower frequencies. The method has higher accuracy, but has high construction cost and low-frequency test lower limit. The reverberant sound field is relatively easy to obtain, and the lower limit of the test frequency is far lower than that of the free field method.
The reverberation sound field is also called diffuse reflection sound field, diffusion sound field or dispersion sound field, and the sound field is characterized in that the sound energy density is uniform, the incident sound energy in all directions is uniform, and the sound waves transmitted to a certain point in all directions are irregularly distributed. In 1961, the method for calibrating the acoustic sensor in the reverberant field was proposed by h.g. diesel, and then applied to absolute (reciprocal) calibration of the hydrophone by domestic wuweinshoo, shenpaul, etc., with good results. In the reciprocal calibration of reverberant fields, it is most important to measure the reciprocal constant of the sound field. Common methods are reverberation time method and reverberation radius method. The reverberation time method is a method for measuring the reverberation time of a reverberation water pool and calculating by using a formula to obtain a reciprocal constant of the reverberation water pool, and is the most common method at present. However, the measurement of the reverberation time has a spatial error and a repetitive error, so that a method of averaging the measurements of multiple sampling points in the space needs to be adopted, and a stable result can be obtained by averaging the measurements of more than 7 sampling points for multiple times. The reciprocal constant of the reverberation pool is obtained by dividing twice the reverberation radius by the density of water and the transmission frequency. The reverberation radius is the distance between the direct sound and the reflected sound, so the reverberation radius can be obtained by respectively measuring the direct wave sound pressure and the reverberation field sound pressure generated by the same sound source in the anechoic water pool and the reverberation water pool, but the anechoic effect of the anechoic water pool at low frequency is poor, so the method has great limitation, and is not common.
In summary, both methods need to go through multiple measurement tests, and the accuracy of the measurement result is proportional to the measurement times, and a large number of repeated tests at different positions are needed to realize the accurate measurement of the reciprocity constant.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a method for accurately measuring the reciprocity constant of a reverberation water pool.
The purpose of the invention is achieved by the following technical scheme: the method for accurately measuring the reciprocity constant of the reverberation water pool mainly comprises the following steps:
1) the reciprocal transducer randomly moves in a region which is at one side of the reverberation pool and is more than a half wavelength away from each pool surface and water surface, and transmits a broadband continuous signal T;
2) the hydrophones are arranged on the other side of the reverberation pool and also randomly move in a region which is more than a half wavelength away from each pool surface and the water surface;
3) the signal-to-noise ratio of an open-circuit voltage signal E1 at the output end of the hydrophone is not lower than 30dB, A/D sampling is carried out on the open-circuit voltage signal at the output end of the hydrophone and a transmitting current signal E2 of the reciprocal transducer, and the signal acquisition length is 2 minutes;
4) e1 and E2 are processed to obtain the power P at the frequency point of interestUAnd PIAnd calculating to obtain a reciprocity constant J of the reverberation pool, wherein the calculation formula is as follows:
Figure BDA0002673514740000021
in the formula, MJSensitivity of the hydrophone; mHIs the sensitivity of the reciprocal transducer.
The reciprocal transducer is a broadband transducer, the working bandwidth of the wideband transducer is larger than the bandwidth B of the broadband continuous signal T, the sending voltage response in the bandwidth is not lower than 160dB, the reciprocity deviation is lower than 0.3dB, and the sensitivity of the reciprocal transducer needs to be known.
The hydrophone has a sensitivity deviation within the bandwidth B of the band continuous signal T of better than +/-0.5 dB.
The broadband continuous signal T has a power spectrum curve similar to the transmission response curve of the reciprocal transducer.
The broadband continuous signal T adopts a signal with the power spectral density being approximately constant.
The E1 and E2 signal processing method mainly comprises the following steps:
1) e1 and E2 are divided into 120 segments,is marked as Us1，Us2，…，US120And Is1，Is2，…，IS120Each segment is 1s in length, a Hanning window with the same time length is added to each segment of signal, and then power spectrum estimation is carried out on the Hanning window, wherein the calculation formula is
Figure BDA0002673514740000022
In the formula (I), the compound is shown in the specification,
Figure BDA0002673514740000023
represents UsiAs a function of the power spectrum of (a),
Figure BDA0002673514740000024
is represented bysiPT represents the power spectrum transformation of the signal, L represents the length of the window signal, and hann represents the hanning window;
2) for the frequency point of interest f0The power of the point is averaged to obtain PUAnd PIThe calculation method is as follows
Figure BDA0002673514740000025
The invention has the beneficial effects that: the invention breaks through the original measuring method of the reciprocal constant of the reverberation pool, the hydrophone and the reciprocal transducer respectively move at the two sides of the pool, the sound field information of the reverberation pool is collected, the reciprocal constant of the reverberation pool is calculated by combining the reciprocity of the reciprocal transducer with the collected sound field information, which is equivalent to that the hydrophone collects numerous space points in a reverberation control area, the used space information is more, the requirement on the sound source level of the transmitting transducer is lower, and the test result has higher precision; the method is used for measuring the reciprocal constant of the reverberation water pool in the research subject, the measurement result is basically consistent with the result of the reverberation time method, the effectiveness of the method is proved, and the method can be gradually applied to the parameter calibration of each reverberation field.
Drawings
Fig. 1 is a schematic view of the sound field arrangement of the present invention.
FIG. 2 is a graph comparing the reciprocity constants measured by the reverberation time method of the present invention.
Description of reference numerals: a reverberation pool 1, a space region 2, a reciprocal transducer 3, a space region 4 and a hydrophone 5.
Detailed Description
The invention will be described in detail below with reference to the following drawings:
example 1: as shown in the attached drawings, the method for accurately measuring the reciprocity constant of the reverberation water pool mainly comprises the following steps:
1) the reciprocal transducer 3 moves randomly in a region which is more than a half wavelength away from each pool surface and water surface at one side of the reverberation pool 1 and transmits a broadband continuous signal T;
2) the hydrophones 5 are arranged on the other side of the reverberation pool 1 and also move randomly in a region which is more than a half wavelength away from each pool surface and water surface;
3) the signal-to-noise ratio of the open-circuit voltage signal E1 at the output end of the hydrophone 5 is not lower than 30dB, A/D sampling is carried out on the open-circuit voltage signal at the output end of the hydrophone 5 and the emission current signal E2 of the reciprocal transducer 3, and the signal acquisition length is 2 minutes;
4) e1 and E2 are processed to obtain the power P at the frequency point of interestUAnd PIAnd calculating to obtain a reciprocity constant J of the reverberation pool 1, wherein the calculation formula is as follows:
Figure BDA0002673514740000031
in the formula, MJThe sensitivity of the hydrophone 5; mHThe sensitivity of the reciprocal transducer 3.
The reciprocal transducer 3 is a broadband transducer, the working bandwidth of which is larger than the bandwidth B of the broadband continuous signal T, the sending voltage response in the bandwidth is not lower than 160dB, the reciprocity deviation is lower than 0.3dB, and the sensitivity of the reciprocal transducer 3 needs to be known. The hydrophone 5 has a sensitivity deviation of better than + -0.5 dB within the bandwidth B of the band continuous signal T. The broadband continuous signal T has a power spectrum curve similar to the transmit response curve of the reciprocal transducer 3. The broadband continuous signal T adopts a signal with approximately constant power spectral density, such as a chirp signal, white noise and the like.
The E1 and E2 signal processing method mainly comprises the following steps:
1) e1 and E2 are divided into 120 segments, which are marked as Us1，Us2，…，US120And Is1，Is2，…，IS120Each segment is 1s in length, a Hanning window with the same time length is added to each segment of signal, and then power spectrum estimation is carried out on the Hanning window, wherein the calculation formula is
Figure BDA0002673514740000032
In the formula (I), the compound is shown in the specification,
Figure BDA0002673514740000033
represents UsiAs a function of the power spectrum of (a),
Figure BDA0002673514740000034
is represented bysiPT represents the power spectrum transformation of the signal, L represents the length of the window signal, and hann represents the hanning window;
2) for the frequency point of interest f0The power of the point is averaged to obtain PUAnd PIThe calculation method is as follows
Figure BDA0002673514740000035
Example 2: this example measured the reciprocity factor in the frequency range of 2kHz to 10kHz in a 15 m.times.9 m.times.6 m reverberant pool 1. The measuring system consists of a reciprocal transducer 3, a hydrophone 5, a signal source, a current sampler, a power amplifier, a digital oscilloscope and the like.
A reciprocal transducer 3 and a hydrophone 5 are arranged in the reverberation pool 1, the reciprocal transducer 3 is positioned in the space region 2, and the hydrophone 5 is positioned in the spaceAnd (4) a region. The distance between the space areas 2 and 4 and the pool wall, the pool bottom and the water surface is 1m respectively, and the distance is 2.5m from the center line, namely the space areas 2 and 4 are symmetrical about the center line of the reverberation pool 1. The reciprocal transducer 3 is excited to generate a noise signal with a bandwidth of 1.6kHz to 12.5kHz, and the reciprocal transducer 3 is moved slowly and at a uniform speed in the spatial region 2, while the hydrophone 5 is moved slowly and at a uniform speed in the spatial region 4. The sampling rate of the digital oscilloscope is set to be 50kHz, and an open-circuit voltage signal E1 at the output end of the hydrophone 5 and a transmitting current signal E2 of the reciprocal transducer 3 are respectively acquired. After 2 minutes, the reciprocal transducer 3 and hydrophone 5 were stopped from moving, and the acquisition of E1 and E2 were stopped. Signals E1 and E2 are equally divided into 120 segments, denoted as Us1，Us2，…，US120And Is1，Is2，…，IS120Adding Hanning window with same time length to each signal segment, and estimating power spectrum with the formula
Figure BDA0002673514740000041
In the formula (I), the compound is shown in the specification,
Figure BDA0002673514740000042
represents UsiA power spectrum function of;
Figure BDA0002673514740000043
is represented bysiA power spectrum function of; PT represents the power spectrum transformation of the signal; an indication of a dot product; l represents the length of the window signal; hann denotes hanning window. Finally, averaging the power of 1/3 octave center frequency in the frequency band of 2 kHz-10 kHz to obtain PUAnd PIThe calculation method is as follows
Figure BDA0002673514740000044
And then calculating a reciprocity constant J of the reverberation pool, wherein the calculation formula is as follows:
Figure BDA0002673514740000045
fig. 2 is a comparison of the reciprocity constants measured by the reverberation time method and the present invention. As can be seen, the reciprocity constants measured by the two methods are substantially the same. The reverberation time method is to obtain the reciprocity constant of the reverberation pool 1 by measuring the reverberation time of 7 space points, and theoretically, the precision of the reciprocity constant can be effectively improved by increasing the number of the measurement points. However, the invention adopts dynamic acquisition, the whole acquisition process contains the reverberation sound field information of an infinite number of space points, and compared with a reverberation time method, the method has higher precision.
It should be understood that equivalent substitutions and changes to the technical solution and the inventive concept of the present invention should be made by those skilled in the art to the protection scope of the appended claims.

Claims (6)

1. A method for accurately measuring a reciprocal constant of a reverberation pool is characterized by comprising the following steps: the method comprises the following steps:
1) the reciprocal transducer (3) moves randomly in a region which is on one side of the reverberation water pool (1) and is more than half wavelength away from each pool surface and water surface, and emits a broadband continuous signal T;
2) the hydrophones (5) are arranged on the other side of the reverberation pool (1) and also randomly move in a region which is more than a half wavelength away from each pool surface and water surface;
3) the signal-to-noise ratio of an open-circuit voltage signal E1 at the output end of the hydrophone (5) is not lower than 30dB, A/D sampling is carried out on the open-circuit voltage signal at the output end of the hydrophone (5) and a transmitting current signal E2 of the reciprocal transducer (3), and the signal acquisition length is 2 minutes;
4) e1 and E2 are processed to obtain the power P at the frequency point concernedUAnd PIAnd calculating to obtain a reciprocity constant J of the reverberation pool (1), wherein the calculation formula is as follows:
Figure FDA0002673514730000011
in the formula, MJIs the sensitivity of the hydrophone (5); mHIs the sensitivity of the reciprocal transducer (3).
2. The method for accurately measuring the reciprocal constant of the reverberation pool of claim 1, wherein: the reciprocal transducer (3) is a broadband transducer, the working bandwidth of the transducer is larger than the bandwidth B of the broadband continuous signal T, the voltage response sent in the bandwidth is not lower than 160dB, and the reciprocity deviation is lower than 0.3 dB.
3. The method for accurately measuring the reciprocal constant of the reverberation pool of claim 1, wherein: the hydrophone (5) has a sensitivity deviation within the bandwidth B of the band-continuous signal T of better than +/-0.5 dB.
4. The method for accurately measuring the reciprocal constant of the reverberation pool of claim 1, wherein: the broadband continuous signal T has a power spectrum curve similar to the transmission response curve of the reciprocal transducer (3).
5. The method for accurately measuring the reciprocal constant of the reverberation pool of claim 1, wherein: the broadband continuous signal T adopts a signal with the power spectral density being approximately constant.
6. The method for accurately measuring the reciprocal constant of the reverberation pool of claim 1, wherein: the E1 and E2 signal processing mode comprises the following steps:
1) e1 and E2 are divided into 120 segments, which are marked as Us1，Us2，…，US120And Is1，Is2，…，IS120Each segment is 1s in length, a Hanning window with the same time length is added to each segment of signal, and then power spectrum estimation is carried out on the Hanning window, wherein the calculation formula is
Figure FDA0002673514730000012
In the formula (I), the compound is shown in the specification,
Figure FDA0002673514730000013
represents UsiAs a function of the power spectrum of (a),
Figure FDA0002673514730000014
is represented bysiPT represents the power spectrum transformation of the signal, L represents the length of the window signal, and hann represents the hanning window;
2) for the frequency point of interest f0The power of the point is averaged to obtain PUAnd PIThe calculation method is as follows
Figure FDA0002673514730000015
CN202010940592.0A 2020-09-09 2020-09-09 Method for accurately measuring reciprocal constant of reverberation pool Active CN112146744B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010940592.0A CN112146744B (en) 2020-09-09 2020-09-09 Method for accurately measuring reciprocal constant of reverberation pool

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010940592.0A CN112146744B (en) 2020-09-09 2020-09-09 Method for accurately measuring reciprocal constant of reverberation pool

Publications (2)

Publication Number Publication Date
CN112146744A true CN112146744A (en) 2020-12-29
CN112146744B CN112146744B (en) 2022-10-21

Family

ID=73889301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010940592.0A Active CN112146744B (en) 2020-09-09 2020-09-09 Method for accurately measuring reciprocal constant of reverberation pool

Country Status (1)

Country Link
CN (1) CN112146744B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4434648A (en) * 1981-02-26 1984-03-06 Cornell Research Foundation, Inc. Electroacoustic transducer calibration method and apparatus
US4453238A (en) * 1982-04-15 1984-06-05 The United States Of America As Represented By The Secretary Of The Navy Apparatus and method for determining the phase sensitivity of hydrophones
CN104501939A (en) * 2014-11-19 2015-04-08 哈尔滨工程大学 Inversion method utilizing single hydrophone to measure reverberation time of non-anechoic pool
US20160091361A1 (en) * 2014-09-25 2016-03-31 Steven E. Crocker System and method for the calibration of a hydrophone line array
CN106501795A (en) * 2016-11-09 2017-03-15 哈尔滨工程大学 A kind of method that utilization reverberation tank carries out underwater acoustic transducer reciprocity calbration
CN108037496A (en) * 2017-10-11 2018-05-15 中国船舶重工集团公司第七〇五研究所 A kind of free field hydrophone plural number sensitivity accurate measurement method
CN109302667A (en) * 2018-09-01 2019-02-01 哈尔滨工程大学 A kind of underwater sound transmitting transducer broadband sends the method for fast measuring and device of response

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4434648A (en) * 1981-02-26 1984-03-06 Cornell Research Foundation, Inc. Electroacoustic transducer calibration method and apparatus
US4453238A (en) * 1982-04-15 1984-06-05 The United States Of America As Represented By The Secretary Of The Navy Apparatus and method for determining the phase sensitivity of hydrophones
US20160091361A1 (en) * 2014-09-25 2016-03-31 Steven E. Crocker System and method for the calibration of a hydrophone line array
CN104501939A (en) * 2014-11-19 2015-04-08 哈尔滨工程大学 Inversion method utilizing single hydrophone to measure reverberation time of non-anechoic pool
CN106501795A (en) * 2016-11-09 2017-03-15 哈尔滨工程大学 A kind of method that utilization reverberation tank carries out underwater acoustic transducer reciprocity calbration
CN108037496A (en) * 2017-10-11 2018-05-15 中国船舶重工集团公司第七〇五研究所 A kind of free field hydrophone plural number sensitivity accurate measurement method
CN109302667A (en) * 2018-09-01 2019-02-01 哈尔滨工程大学 A kind of underwater sound transmitting transducer broadband sends the method for fast measuring and device of response

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑士杰: "《水声计量测试技术》", 31 January 1995, 哈尔滨工程大学出版社 *

Also Published As

Publication number Publication date
CN112146744B (en) 2022-10-21

Similar Documents

Publication Publication Date Title
CN106501795B (en) A method of underwater acoustic transducer reciprocity calbration is carried out using reverberation tank
CN109302667B (en) Method and device for rapidly measuring broadband sending response of underwater acoustic emission transducer
US11391863B2 (en) Method of free-field broadband calibration of hydrophone sensitivity based on pink noise
US4391124A (en) Electroacoustic transducer calibration method and apparatus
US4434648A (en) Electroacoustic transducer calibration method and apparatus
CN109238436B (en) Method for measuring transient sound source characteristics in reverberation pool
CN109991590B (en) System and method for testing low-frequency emission characteristic of transducer in pressure tank in limited space
JPS5847026B2 (en) How to calibrate acoustic radiation transducers
CN110530510B (en) Method for measuring sound source radiation sound power by utilizing linear sound array beam forming
CN109324320A (en) A method of hydrophone batch, which is carried out, using reverberation tank calibrates
CN111586546B (en) Method and system for measuring resonance point transmission response of low-frequency transducer
CN112146744B (en) Method for accurately measuring reciprocal constant of reverberation pool
CN110440896B (en) Ultrasonic measurement system and measurement method
Vorländer et al. Novel broad-band reciprocity technique for simultaneous free-field and diffuse-field microphone calibration
CN111397721A (en) Method and system for absolute calibration of co-vibrating vector hydrophone based on water surface boundary vibration measurement technology
CN112146745B (en) Method for accurately measuring reverberation time of pool
CN111207820A (en) Method for calibrating array elements of buoy hydrophone array in reverberation pool
CN113758617B (en) Efficient nondestructive testing system and method for stress gradient based on broadband frequency-sweep signal frequency domain calculation
CN110703186A (en) Bi-hydrophone phase difference calibration and measurement system and method for near field of pool
CN110208778B (en) Transducer broadband complex response measuring device and method based on logarithm variable window function
RU2655049C1 (en) Method of calibration of the hydrophone on the field at low frequencies
Pedersen et al. Application of time delay spectrometry for rough surface characterization
Sawada et al. Precision calibration of echo sounder by integration of standard sphere echoes
Liu et al. Acoustic method for obtaining the pressure reflection coefficient using a half-wave layer
Behler et al. Reciprocal measurements on condenser microphones for quality control and absolute calibration

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant