Method for optimization of a parametric array shallow profile difference frequency conversion performance

Abstract

A method for optimization of a parametric array shallow profile difference frequency conversion performance, comprising the following steps: building a three-layer parametric array shallow profile emission acoustic wave propagation path model, and establishing a detection horizontal position DL equation capable of detecting a detection performance of the parametric array shallow profile; and building a model for optimization of the parametric array shallow profile difference frequency conversion performance. A parametric receiving array includes the receiving hydrophone and a pump transducer, comprising: a silicone rubber column is added between the receiving hydrophone and the pump transducer to improve the conversion efficiency of the parametric receiving array. The present disclosure improves the conversion efficiency of the parametric receiving array by increasing the high-frequency power and acoustic pressure level and reducing the ratio to improve the conversion efficiency and the parameter receiving array uses silicone rubber columns.

Claims

1. A method for optimization of a parametric array shallow profile difference frequency conversion performance, comprising the following steps: S1: building a three-layer parametric array shallow profile emission acoustic wave propagation path model, the three-layer parametric array shallow profile emission acoustic wave propagation path model includes a water layer, a sedimentary layer and a rock layer; the sedimentary layer is a pore fluid formed by particles, and establishing a detection horizontal position DL equation capable of detecting a detection performance of the parametric array shallow profile and estimating an echo level at a ground profile and a top of a rock formation:
DL=SL−(PG+DL.sub.W.Math.S+AL.sub.W)−TL.sub.W.Math.S−AL.sub.S−RL.sub.S.Math.R−TL.sub.S.Math.W+RVS+RG+SPG: the SL is a horizontal position of an acoustic wave emission source, the PG is a generation efficiency from an initial wave frequency to a difference frequency, the DL.sub.W.Math.S is a diffusion loss of an initial wave in a water and the sedimentary layer, and the AL.sub.W is an absorption loss of the initial wave in a water, and a unidirectional wave propagation between a transducer and a surface of a sediment is regarded as a nonlinearity and a wave propagation of the transducer reflected from the surface of the sediment is regarded as a linearity; the TL.sub.W.Math.S is a transmission loss of the initial wave at a boundary between the water layer and the sedimentary layer and the TL.sub.S.Math.W is the transmission loss of an echo at the boundary between the water layer and the sedimentary layer; the AL.sub.S is the absorption loss of the initial wave in the sediment; the RL.sub.S.Math.R is a reflection loss at the boundary between the sediment and the rock layer; the RVS is a voltage sensitivity of a receiving hydrophone for an acoustic echo; the RG is a gain of the receiving hydrophone; and the SPG is a signal processing gain of the receiving hydrophone; S2: building a model for optimization of the parametric array shallow profile difference frequency conversion performance and calculating a low-frequency wave radiation power W.sub.ed: $W_{ed} = \frac{ρ_{f} e^{α_{S} r}}{2 ρ_{S} c_{S}} 4 π r;$ $= \frac{w^{4} β^{2} w_{1} w_{2}}{2 {πρ}_{S} c_{s}^{7} γ};$ a propagation distance of an acoustic wave in the water is r, the propagation distance in the sediment is d; and a propagation speed of the acoustic wave in the water is c.sub.R, the propagation speed of the acoustic wave in the sediment is c.sub.S, and the propagation speed of the acoustic wave in a rock is c.sub.W; a water density is ρ.sub.W, a sediment density is ρ.sub.S, and a rock density is ρ.sub.R; an attenuation coefficient of the acoustic wave absorbed in the sediment is α.sub.S; the β is a porosity of the sediment, the e is an expansion rate of a member attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance, and the γ is a structure factor of a transmitted acoustic wave propagation path model; wherein, the w.sub.1 and the w.sub.2 are frequencies of two high-frequency frequencies respectively, since the frequencies of the two high-frequency frequencies are similar, $W_{1} = W_{2} = \frac{W}{2}$ and the w is a total radiated wave frequency; the ρ.sub.f is a pore fluid density in the sediment, and ρ.sub.s=βρ.sub.f+(1−β)ρ.sub.g and the ρ.sub.g is a particulate matter density in the sediment; S3: calculating a low-frequency wave conversion efficiency, and the low-frequency wave conversion efficiency is a ratio of a low-frequency wave frequency W.sub.ed to a total radiated wave frequency W: $\frac{W_{ed}}{W} = 2 π^{3} β^{2} f_{d}^{4} W / ρ_{S} c_{s}^{7} γ^{2};$ the f.sub.d is a differential frequency; S4: adjusting and increasing a high-frequency power and a high-frequency acoustic pressure level and reducing the ratio $\frac{W_{e d}}{W} = 2 π^{3} β^{2} f_{d}^{4} W / ρ_{S} c_{s}^{7} γ^{2}$ to improve a conversion efficiency. $\frac{W_{e d}}{W} .$

2. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 1, wherein: under the condition that the acoustic wave emission source is 1 m away from the transducer and a bidirectional distance diffusive loss is −6 db/DD, the ALS=2 α.sub.Sd DLW.Math.S=20 log 2(r+d), and ALS=2α.sub.Sd.

3. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 1, wherein: an acoustic wave loss caused by going back and forth between a water surface and the sediment is a sum of TL.sub.W.Math.S and TL.sub.S.Math.W:
TL.sub.W.Math.S+TL.sub.S.Math.W=−20 log 2(1+R.sub.W.Math.S)(1−R.sub.W.Math.S), and $R_{W .Math. S} = \frac{ρ_{S} c_{S} - ρ_{W} c_{W}}{ρ_{S} c_{S} + ρ_{W} c_{W}};$ and the R.sub.W.Math.S is a reflection coefficient at the water surface and the sediment.

4. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 1, wherein: RL.sub.S.Math.R=−20 log(R.sub.S.Math.R), $R_{S .Math. R} = \frac{ρ_{R} c_{R} - ρ_{S} c_{S}}{ρ_{R} c_{R} + ρ_{S} c_{S}};$ the R.sub.S.Math.R is the reflection coefficient at the sediment and the rock.

5. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 1, wherein: the echo has a transmission loss and the reflection loss in the sediment, the transmission loss and the reflection loss in an echo sediment are related to the ρ.sub.S, c.sub.S and α.sub.S, and establishing a Biot-Stoll model to calculate ρ.sub.S, c.sub.S and α.sub.S: $\nabla^{2} (H e - C ξ) = \frac{\partial^{2}}{\partial t^{2}} (ρ e - ρ_{f} ξ);$ $\nabla^{2} (H e - M ξ) = \frac{\partial^{2}}{\partial t^{2}} (ρ_{f} e - m ξ) - \frac{F η}{k} \frac{\partial ξ}{\partial t};$ wherein, the expansion rate of a member attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance is e=div(u, U), and a fluid volume flowing into or out of a volume unit attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance is ξ=β div(u-U,U), the β is the porosity of the sediment, the u is a position after the parametric array shallow profile is moved, the U is a position after the pore fluid is moved and the η is a flow velocity of the pore fluid; the div function is a result of finding a quotient and a remainder of two parameters; the H is a calculation weighting parameter of the expansion rate e, the C is the calculation weighting parameter of the fluid volume ξ flowing into or out of a volume unit attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance, the M is the calculation weighting parameter of a volume variation of the fluid volume ξ, and the m is the structure factor of the transmitted acoustic wave propagation path model.

6. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 5, wherein: a H formula is as following: $H = \frac{{(K_{g} - K_{b})}^{2}}{D - K_{b}} + K_{b} + \frac{4}{3} u;$ wherein: the K.sub.b a frame bulk modulus, and the K.sub.g is a particle bulk modulus; a C formula is as following: $C = \frac{K_{g} (K_{g} - K_{b})}{D - K_{b}};$ a M formula is as following: $M = \frac{K_{g}^{2}}{D - K_{b}};$ a m formula is as following: $m = γ \frac{ρ_{f}}{β};$ and the m is the structure factor of the transmitted acoustic wave propagation path model; a D formula is as following: $D = K_{g} [1 + β (\frac{K_{g}}{K_{f}} - 1)];$ the K.sub.f is a bulk modulus of the pore fluid.

7. The method for optimization of the parametric array shallow profile difference frequency conversion performance according to claim 1, wherein: a γ formula is as following: $γ = 1 + g (\frac{1 - β}{β});$ and the g is a radius of the particles in the sediment.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a schematic diagram of the optimization method;

(2) FIG. 2 is a schematic diagram of the acoustic wave propagation path of the parametric array shallow profile in the water body;

(3) FIG. 3 is a schematic diagram of pore fluid with particles in sediment;

(4) FIG. 4 is a schematic diagram of the structure of the parametric receiving array.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(5) The following content will clearly and completely describe the technical solutions in the embodiments of the present disclosure with reference to the accompanying drawings in the embodiments. Obviously, the described embodiments are only a part of the embodiments of the present disclosure, rather than all the embodiments. Based on the embodiments of the present disclosure, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of the present disclosure.

(6) The parametric array transmits two signals by a slightly different high frequency at a high acoustic pressure. Under high acoustic pressure, the non-linear characteristics of the acoustic propagation of the two signals, interactions and new frequencies (sum frequency, difference frequency, and harmonics) increase; and the difference in the emission frequency is called the secondary frequency. The parametric array is a base array that uses two high-frequency acoustic waves with similar frequencies to produce a nonlinear effect in the water medium, forming a difference frequency narrow beam. Usually two close transmitting transducers are used as primitives to respectively emit high-frequency sinusoidal acoustic waves with close frequencies, and when the acoustic amplitude is large enough, the water medium produces a nonlinear effect, and the acoustic velocity is no longer constant, but a function of pressure. When the acoustic amplitude increases positively, the acoustic velocity increases, and when the acoustic amplitude increases negatively, the acoustic velocity decreases, which distort the sine wave into a sawtooth wave, and form an interaction zone at a certain distance from the sound axis of the transmitter array, generating two high-frequency sum frequency, difference frequency and other combined frequencies. But the frequencies of the two high-frequency sum frequency, difference frequency and other combined frequencies are similar, the difference frequency wave is lower than the acoustic source frequency and the sum frequency of the two. During propagation, the high-frequency acoustic waves decay rapidly, leaving only difference frequency waves.

(7) The difference frequency conversion efficiency of the parametric array is related to different media. When the transmitted wave frequency, total radiated wave frequency and difference frequency are fixed, the conversion efficiency in fresh water is generally higher than that in sea water, and the conversion efficiency in air is higher than that in fresh water. The conversion efficiency can be improved by increasing the power and acoustic pressure level of the high-frequency and reducing the ratio, achieving the goal of optimization. However, with the increasing power, sawtooth waveforms will appear, and the absorption of acoustic waves by the medium will increase significantly, therefore, the equivalent length will be limited, so there is a certain limit to improve the efficiency by increasing the power. But the conversion efficiency is related to the characteristics of the medium. By changing the medium, the conversion efficiency can be improved by using a medium with a low acoustic velocity, a low density and a larger nonlinear coefficient.

(8) Based on the above principles, as shown in FIG. 1-2, the method for optimization of a parametric array shallow profile difference frequency conversion performance considering the wave propagation path between the sensor and the top of the rock formation on the ground profile, comprising the following steps:

(9) S1: building a three-layer parametric array shallow profile emission acoustic wave propagation path model, the three-layer parametric array shallow profile emission acoustic wave propagation path model includes a water layer, a sedimentary layer and a rock layer; as shown in FIG. 3, the sedimentary layer is a pore fluid formed by particles, and establishing a detection horizontal position DL equation capable of detecting a detection performance of the parametric array shallow profile and estimating an echo level at a ground profile and a top of a rock formation:
DL=SL−(PG+DL.sub.W.Math.S+AL.sub.W)−TL.sub.W.Math.S−AL.sub.S−RL.sub.S.Math.R−TL.sub.S.Math.W+RVS+RG+SPG;

(10) since the loss in water is less than the loss in sediments and rocks, the initial wave absorption loss in water is ignored in the general equation;

(11) the SL is a horizontal position of an acoustic wave emission source, the PG is a generation efficiency from an initial wave frequency to the difference frequency, the DL.sub.W.Math.S is a diffusion loss of an initial wave in a water and the sedimentary layer, and the AL.sub.W, is an absorption loss of the initial wave in the water, a unidirectional wave propagation between a transducer and a surface of a sediment is regarded as a nonlinearity and a wave propagation of the transducer reflected from the surface of the sediment is regarded as a linearity; the TL.sub.W.Math.S is a transmission loss of the initial wave at a boundary between the water layer and the sedimentary layer and the TL.sub.S.Math.W is the transmission loss of an echo at the boundary between the water layer and the sedimentary layer; the AL.sub.S is the absorption loss of the initial wave in the sediment; the RL.sub.S.Math.R is a reflection loss at the boundary between the sediment and the rock layer; the RVS is a voltage sensitivity of an acoustic echo receiving hydrophone; the RG is a gain of the receiving hydrophone; and the SPG is a signal processing gain of the receiving hydrophone;

(12) S2: building a model for optimization of the parametric array shallow profile difference frequency conversion performance and calculating a low-frequency wave radiation power W.sub.ed:

(13) $\begin{matrix} W_{ed} = \frac{ρ_{f} e^{α_{S} r}}{2 ρ_{S} c_{S}} 4 π r \\ = \frac{w^{4} β^{2} w_{1} w_{2}}{2 {πρ}_{S} c_{S}^{7} γ}; \end{matrix}$

(14) a propagation distance of the acoustic wave in water is r, the propagation distance in the sediment is d; and a propagation speed of the acoustic wave in water is c.sub.R, the propagation speed of the acoustic wave in the sediment is c.sub.S, and the propagation speed of the acoustic wave in a rock is c.sub.W, the unit is m/s, a water density is ρ.sub.W, a sediment density is ρ.sub.S, and a rock density is ρ.sub.R and the unit is kg/m.sup.3; an attenuation coefficient of the acoustic wave absorbed in the sediment is α.sub.S, and unit is dB/m; the β is a porosity of the sediment, the e is an expansion rate of a member attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance, and the γ is a structure factor of a transmitted acoustic wave propagation path model;

(15) wherein, the w.sub.1 and the w.sub.2 are frequencies of two high-frequency frequencies respectively, since the frequencies of two high-frequency frequencies are similar,

(16) $W_{1} = W_{2} = \frac{W}{2}$
and the w is a total radiated wave frequency; the ρ.sub.f is a pore fluid density in the sediment, and ρ.sub.s=βρ.sub.f+(1−β)ρ.sub.g and the ρ.sub.g is a particulate matter density in the sediment;

(17) S3: calculating a low-frequency wave conversion efficiency, and the low-frequency wave conversion efficiency is a ration of a low-frequency wave frequency W.sub.ed to the total radiated wave frequency W:

(18) $\frac{W_{ed}}{W} = 2 π^{3} β^{2} f_{d}^{4} W / ρ_{S} c_{S}^{7} γ^{2}$

(19) the f.sub.d is a differential frequency;

(20) S4: adjusting and increasing a high-frequency power and a high-frequency acoustic pressure level and reducing the ratio

(21) 0 $\frac{W_{ed}}{W} = 2 π^{3} β^{2} f_{d}^{4} W / ρ_{S} c_{S}^{7} γ^{2}$
to improve a conversion efficiency

(22) $\frac{W_{ed}}{W} .$

(23) Under the condition that the acoustic wave emission source is 1 m away from the transducer and a bidirectional distance diffusive loss is −6 db/DD, the ALS=2 α.sub.Sd DLW.Math.S=20 log 2(r+d), and ALS=2α.sub.Sd.

(24) An acoustic wave loss caused by going back and forth between a water surface and the sediment is a sum of TL.sub.W.Math.S and TL.sub.S.Math.W:
TL.sub.W.Math.S+TL.sub.S.Math.W=−20 log 2(1+R.sub.W.Math.S)(1−R.sub.W.Math.S), and

(25) $R_{W .Math. S} = \frac{ρ_{S} c_{S} - ρ_{W} c_{W}}{ρ_{S} c_{S} + ρ_{W} c_{W}};$

(26) and the R.sub.W.Math.S is a reflection coefficient at the water surface and the sediment.

(27) The RL.sub.S.Math.R=−20 log(R.sub.S.Math.R),

(28) $R_{S .Math. R} = \frac{ρ_{R} c_{R} - ρ_{S} c_{S}}{ρ_{R} c_{R} + ρ_{S} c_{S}}$

(29) the R.sub.S.Math.R is the reflection coefficient at the sediment and the rock.

(30) In order to adapt to different types of sediments and avoid the inaccuracy of the final applied nonlinear parameters, which results in a poor optimization of conversion efficiency, it is necessary to calculate the transmission loss and reflection losses of the echo in the sediment. The transmission loss and the reflection loss in an echo sediment are related to the ρ.sub.S, c.sub.S and α.sub.S, and calculating ρ.sub.S, c.sub.S and α.sub.S by establishing Biot-Stoll model, and the Biot-Stoll model is based on the consideration that the three-layer model constructed will carry out the overall frame movement in the water body and will absorb water and expand in the water body:

(31) $\nabla^{2} (He - C ξ) = \frac{\partial^{2}}{\partial t^{2}} (ρ e - ρ_{f} ξ);$

(32) $\nabla^{2} (He - M ξ) = \frac{\partial^{2}}{\partial t^{2}} (ρ_{f} e - m ξ) - \frac{F η}{k} \frac{\partial ξ}{\partial t};$

(33) wherein, the expansion rate of a member attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance is e=div(u, U), and a fluid volume flowing into or out of a volume unit attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance is ξ=β div(u-U,U), the f is the porosity of the sediment, the u is a position after the parametric array shallow profile is moved, the U is the position after the pore fluid is moved and the η is a flow velocity of the pore fluid; the div function is a result of finding a quotient and a remainder of two parameters; the H is a calculation weighting parameter of the expansion rate e, the C is the calculation weighting parameter of the fluid volume ξ flowing into or out of a volume unit attached to the model for optimization of the parametric array shallow profile difference frequency conversion performance, the M is the calculation weighting parameter of a volume variation of the fluid volume ξ, and the m is the structure factor of the transmitted acoustic wave propagation path model.

(34) The H formula is as following:

(35) $H = \frac{{(K_{ℊ} - K_{b})}^{2}}{D - K_{b}} + K_{b} + \frac{4}{3} u;$

(36) the K.sub.b is frame bulk modulus, and the K.sub.g is particle bulk modulus;

(37) the C formula is as following:

(38) $C = \frac{K_{ℊ} (K_{ℊ} - K_{b})}{D - K_{b}};$

(39) the M formula is as following:

(40) $M = \frac{K_{ℊ}^{2}}{D - K_{b}};$

(41) the m formula is as following:

(42) $m = γ \frac{ρ_{f}}{β};$

(43) the D formula is as following:

(44) 0 $D = K_{ℊ} [1 + β (\frac{K_{ℊ}}{K_{f}} - 1)];$

(45) the K.sub.f is the bulk modulus of the pore fluid.

(46) The structure factor of the transmitted acoustic wave propagation path model γ formula comprising:

(47) $γ = 1 + ℊ (\frac{1 - β}{β});$

(48) and the g is a radius of the particles in the sediment.

(49) Establishing a detection horizontal position DL equation capable of detecting a detection performance of the parametric array shallow profile and estimating an echo level at a ground profile and a top of a rock formation and a Biot-Stoll model to evaluate the performance of PSBP system and to improve the measuring system, especially in terms of frequency setting.

(50) As shown in FIG. 4, the present disclosure provides a parametric receiving array adopting the method for optimization of a parametric array shallow profile difference frequency conversion performance, and the parametric receiving array includes the receiving hydrophone and a pump transducer, comprising: a silicone rubber column is added between the receiving hydrophone and the pump transducer to improve the conversion efficiency of the parametric receiving array. Furthermore, the optimized parametric array shallow profile and narrow beam without sidelobe are further used, reducing the sound virtual image and reverberation of other reflectors. Therefore, the parametric receiving array which is suitable for detection in slopes where submarine hydrothermal deposits can be found on and in the irregular surface areas is able to emit a highly directional acoustic beam, thereby improving the positioning accuracy of the detection target. Meanwhile, the parametric receiving array can optimize the shallow profile of some parameter arrays to meet the low-frequency power consumption mode. And the parametric receiving array can be installed on autonomous underwater robots to detect buried sediments in deep sea areas and reduce the power consumption of autonomous underwater robots.

(51) Due to the increase of the non-linear coefficient of the silicone rubber medium and the decrease of the acoustic velocity, the beam width has a sharper reducing directivity, and the conversion efficiency is improved. Therefore, in the parameter receiving array composed of pump transducer and receiving hydrophone, the conversion efficiency of the parametric array receiving array can be improved by adding silicone rubber columns to the transducer and the receiving hydrophone.

(52) Although the present disclosure has been described with reference to the preferred embodiments, without departing from the scope of the present disclosure, various modifications can be made to the disclosure and the components therein can be moved with equivalents. In particular, as long as there is no structural conflict, the various technical features mentioned in the various embodiments can be combined in any manner. The present disclosure is not limited to the specific embodiments disclosed in the text, but includes all technical solutions falling within the scope of the claims.

Method for optimization of a parametric array shallow profile difference frequency conversion performance

Assignee

Inventors

Cpc classification

Classification Explorer

G01N29/348

PHYSICS

Classification Explorer

G01N2291/02491

PHYSICS

Classification Explorer

G01S7/539

PHYSICS

Classification Explorer

G01S7/521

PHYSICS

Classification Explorer

G01S15/88

PHYSICS

International classification

Classification Explorer

G01S7/521

PHYSICS

Classification Explorer

G01N29/34

PHYSICS

Classification Explorer

G01S15/88

PHYSICS

Classification Explorer

G01S7/539

PHYSICS

Abstract

Claims

Description