Method for determining effective sound velocity in deep sea

The present invention relates to the sound velocity propagation field, and in particular, to a method for determining an effective sound velocity in the deep sea.

An underwater acoustic positioning system measures a propagation time from a transmission point to a receiving point, then converts the propagation time into a slant distance between the transmission point and the receiving point, and finally implements underwater positioning. However, due to uneven distribution of a sound velocity in a vertical direction and resulting sound ray bending, the propagation distance of a sound signal from the transmission point to the receiving point is longer than the slant distance between the transmission point and the receiving point, which has great influence on the ranging accuracy and the positioning accuracy of the underwater acoustic positioning system. Therefore, sound ray bending needs to be corrected.

One idea is to use a sound ray tracing technology. In the sound ray tracing technology, a hierarchical approximation method and a layer-by-layer calculation method are used to compensate sound ray bending, so as to improve the positioning accuracy of the underwater acoustic positioning system. However, the sound ray tracing technology needs to balance between the positioning accuracy and the calculation amount, greatly restricting actual applications.

Another idea is to use an effective sound velocity technology. An effective sound velocity v e is a ratio of a slant distance R between a transmission point and a receiving point to a propagation time T t from the transmission point and the receiving point. Correspondingly, the slant distance R between the transmission point and the receiving point is equal to a product of the propagation time T t and effective sound velocity v e . The measurement accuracy of an underwater propagation time is very high, and therefore the accuracy of the slant distance R depends on the accuracy of effective sound velocity v e .

Only a direct wave area (both an ultra short baseline positioning system and a long baseline positioning system work in this area) is considered. The sound velocity v e satisfies a function between a transmission point depth z A , a receiving point depth z B , a pitch angle β between the transmission point depth z A and the receiving point depth z B , and a sound velocity profile c(z) between the transmission point depth z A and the receiving point depth z B , which is referred to a Taylor series analysis method:

v e= v (1 −C 2 J 2 −C 3 J 3 −C 4 /J 4 )   (1)

where v is an arithmetic mean sound velocity; w is a mean square-arithmetic mean sound velocity (used for rapid calculation); C 2 , C 3 , and C 4 are a second-order sound velocity coefficient, a third-order sound velocity coefficient, and a fourth-order sound velocity coefficient; J 2 , J 3 , and J 4 are a second-order J function, a third-order J function, and a fourth-order J function:

Therefore, information required for calculating the effective sound velocity v e by using the Taylor series analysis method is the transmission point depth z A , the receiving point depth z B , the pitch angle β between the transmission point depth z A and the receiving point depth z B , and the sound velocity profile c(z) between the transmission point depth z A and the receiving point depth z B .

A main calculation amount of calculating the effective sound velocity v e by using the Taylor series analysis method is reflected in two integral operations, that is, Formula (2) and Formula (4). As long as the transmission point depth or the receiving point depth changes, it is necessary to recalculate corresponding v e . In the deep sea, a difference between the transmission point depth and the receiving point depth is usually very large, and an actual sound velocity profile between the transmission point depth and the receiving point depth is very complex, leading to a very large main integral calculation amount of Formula (2) and Formula (4). As a result, it is difficult to satisfy an actual application requirement.

An objective of the present invention is to provide a method for determining an effective sound velocity in the deep sea, so as to rapidly and accurately determine an effective sound velocity.

To achieve the above purpose, the present invention provides the following technical solution.

The present invention provides a method for determining an effective sound velocity in the deep sea, where the method is applied to an apparatus for determining an effective sound velocity in the deep sea. The apparatus includes a transmission point, a receiving point, and an underwater mobile carrier. The transmission point is installed on the sea surface and the depth of the transmission point is unchanged. The receiving point is installed on the underwater mobile carrier, and the depth of the receiving point changes with movement of the underwater mobile carrier. The underwater mobile carrier can measure a sound velocity profile between the transmission point and the receiving point and a horizontal distance between the transmission point and the receiving point.

The method includes obtaining a transmission point depth at a moment k and a receiving point depth at the moment k and obtaining a receiving point depth at a moment k+1 and the horizontal distance between the transmission point and the receiving point. A pitch angle is determined between the transmission point and the receiving point according to the transmission point depth, the receiving point depth at the moment k+1, and the horizontal distance between the transmission point and the receiving point. A sound velocity profile is obtained between the transmission point and the receiving point at the moment k+1, and a difference between the receiving point depth at the moment k and the receiving point depth at the k+1 moment is also obtained. An arithmetic average sound velocity variable quantity and a mean square arithmetic average sound velocity variable quantity are determined according to the receiving point depth at the moment k+1, the receiving point depth at the moment k, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, and the sound velocity profile between the transmission point and the receiving point at the moment k+1.

A depth difference between the receiving point and the transmission point at the moment k is obtained and then an arithmetic average sound velocity at the moment k+1 and a mean square average sound velocity at the moment k+1 is determined according to a depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, an arithmetic average sound velocity at the moment k, a mean square average sound velocity at the moment k, the arithmetic average sound velocity variable quantity, and the mean square arithmetic average sound velocity variable quantity.

An n-order sound velocity function variable quantity (n=2,3,4) is determined according to the sound velocity profile between the transmission point and the receiving point at the moment k+1, the receiving point depth at the moment k, the receiving point depth at the moment k+1, and the arithmetic average sound velocity at the moment k+1.

An n-order sound velocity function (n=2,3,4) at the moment k+1 is then determined according to the n -order sound velocity function variable quantity, the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the sound velocity profile between the transmission point and the receiving point at the moment k+1, and the arithmetic average sound velocity;

An n-order sound velocity coefficient (n=2,3,4) at the moment k+1 is calculated according to a pitch angle at the moment k+1, and an effective sound velocity at the moment k+1 is determined according to the n-order sound velocity coefficient at the moment k+1, the n-order sound velocity function corresponding to the n-order sound velocity coefficient at the moment k+1, and the arithmetic average sound velocity at the moment k+1.

Optionally, the obtaining of the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1 includes determining the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1 according to the receiving point depth at the moment k and the receiving point depth at the moment k+1 by using a formula Δz k =z B,k+1 −z B,k , where

Δz k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; z B,k is the receiving point depth at the moment k; and z B+1 is the receiving point depth at the moment k+1.

Optionally, the obtaining of the depth difference between the receiving point and the transmission point at the moment k includes determining the depth difference between the receiving point and the transmission point at the moment k according to the receiving point depth at the moment k and the transmission point depth at the moment k by using a formula z ABk =z B,k −z A , where

z ABk is the depth difference between the receiving point and the transmission point at the moment k; z A is the transmission point depth; and z B,k is the receiving point depth at the moment k.

Optionally, the determining an arithmetic average sound velocity variable quantity and a mean square arithmetic average sound velocity variable quantity according to the receiving point depth at the moment k+1, the receiving point depth at the moment k, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, and the sound velocity profile between the transmission point and the receiving point at the moment k+1 includes respectively determining the arithmetic average sound velocity variable quantity and the mean square arithmetic average sound velocity variable quantity according to the receiving point depth at the moment k+1, the receiving point depth at the moment k, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, and the sound velocity profile between the transmission point and the receiving point at the moment k by using a formula ΔV k =∫ z B,k z B,k +Δz k c(z)dz and a formula ΔW k =∫ z B,k z B,k +Δz k c 2 (z)dz, where

ΔV k is the arithmetic average sound velocity variable quantity; ΔW k is the mean square arithmetic average sound velocity variable quantity; z B,k +Δz k is the receiving point depth at the moment k+1; z B,k is the receiving point depth at the moment k; Δz k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; and c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1.

Optionally, the determining an arithmetic average sound velocity at the moment k+1 and a mean square average sound velocity at the moment k+1 according to a depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, an arithmetic average sound velocity at the moment k, a mean square average sound velocity at the moment k, the arithmetic average sound velocity variable quantity, and the mean square arithmetic average sound velocity variable quantity includes obtaining a sound velocity profile at the moment k, determining the arithmetic average sound velocity at the moment k and the mean square average sound velocity at the moment k according to the sound velocity profile c(z) at the moment k, the transmission point depth z A , and the receiving point depth z B,k ; and

determining the arithmetic average sound velocity at the moment k+1 and the mean square average sound velocity at the moment k+1 according to the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the arithmetic average sound velocity at the moment k, the mean square average sound velocity at the moment k, the arithmetic average sound velocity variable quantity, and the mean square arithmetic average sound velocity variable quantity by using a formula

v _ k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁢ z [ ABk ) ⁢ vk + Δ ⁢ ⁢ Vk

and a formula

w _ k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁡ [ z ABk ⁢ w _ k + Δ ⁢ ⁢ W k ] ,

where

v k+1 is the arithmetic average sound velocity at the moment k+1; w k+1 is the mean square average sound velocity at the moment k+1; z ABk is the depth difference between the receiving point and the transmission point at the moment k+1; v k is the arithmetic average sound velocity at the moment k; w k is the mean square average sound velocity at the moment k; ΔV k is the arithmetic average sound velocity variable quantity; and ΔW k is the mean square arithmetic average sound velocity variable quantity.

Optionally, the determining an n-order sound velocity function variable quantity (n=2,3,4) according to the sound velocity profile between the transmission point and the receiving point at the moment k+1, the receiving point depth at the moment k, the receiving point depth at the moment k+1, and the arithmetic average sound velocity at the moment k+1 specifically includes determining the n-order sound velocity function variable quantity (n=2,3,4) according to the sound velocity profile between the transmission point and the receiving point at the moment k+1, the receiving point depth at the moment k, the receiving point depth at the moment k+1, and the arithmetic average sound velocity at the moment k+1 by using a formula

Q k = ∫ z B , k z B , k + 1 ⁢ ( c ⁡ ( z ) - v _ k + 1 v _ k + 1 ) n ⁢ dz ,

where

Q k is the n-order sound velocity function variable quantity; c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1; z B,k is the receiving point depth at the moment k; z B,k+1 is the receiving point depth at the moment k+1; and v k+1 is the arithmetic average sound velocity at the moment k+1.

Optionally, the determining an n -order sound velocity function (n=2,3,4) at the moment k+1 according to the n-order sound velocity function variable quantity, the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the sound velocity profile between the transmission point and the receiving point at the moment k+1, and the arithmetic average sound velocity specifically includes determining the n-order sound velocity function (n=2,3,4) at the moment k+1 according to the n -order sound velocity function variable quantity, the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the sound velocity profile between the transmission point and the receiving point at the moment k+1, and the arithmetic average sound velocity by using a formula

J n k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁡ [ ∫ z A z B , k ⁢ ( c ⁡ ( z ) - v _ k + 1 v _ k + 1 ) n ⁢ dz + Q k ] ,

where

J n k+1 is the n-order sound velocity function corresponding to the n-order sound velocity coefficient at the moment k+1; Q k is the n-order sound velocity function variable quantity; z ABk is the depth difference between the receiving point and the transmission point at the moment k+1; Δz k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1; and v k+1 is the arithmetic average sound velocity at the moment k+1.

Optionally, the calculating an n-order sound velocity coefficient (n=2,3,4) at the moment k+1 according to a pitch angle at the moment k+1 includes calculating a second-order sound velocity coefficient, a third-order sound velocity coefficient, and a fourth-order sound velocity coefficient according to the pitch angle at the moment k+1 by using formulas: C 2 k+1 =1−½cot 2 β k+1 ; C 3 k+1 =−C 2 k+1 −1/2cot 4 β k+1 ; and C 4 k+1 =1/2cot 4 β k+1 (3−5cot 2 β k+1 )+C 2 k+1 , where C 4 k+1 J 4 k+1 ), where

β k+1 is the pitch angle at the moment k+1, and C 2 k+1 , C 3 k+1 , and C 4 k+1 are respectively the second-order sound velocity coefficient, the third-order sound velocity coefficient, and the fourth-order sound velocity coefficient.

Optionally, the determining an effective sound velocity at the moment k+1 according to the n-order sound velocity coefficient at the moment k+1, the n-order sound velocity function corresponding to the n-order sound velocity coefficient at the moment k+1, and the arithmetic average sound velocity at the moment k+1 includes determining the effective sound velocity at the moment k+1 according to a formula v e k+1 = v k+1 (1−C 2 k+1 J 2 k+1 −C 3 k+1 J 3 k+1 −C 4 k+1 J 4 k+1 ), where

v e k+1 is the effective sound velocity at the moment k+1, v k+1 is the arithmetic average sound velocity at the moment k+1, J 2 k+1 , J 3 k+1 , and J 4 k+1 and respectively a second-order sound velocity function, a third-order sound velocity function, and a fourth-order sound velocity function, and C 2 k+1 , C 3 k+1 , and C 4 k+1 are respectively the second-order sound velocity coefficient, the third-order sound velocity coefficient, and the fourth-order sound velocity coefficient.

According to specific embodiments provided in the present invention, the present invention discloses the following technical effects: The present invention provides a method for determining an effective sound velocity in the deep sea, which can convert an integral operation from original whole integral from the transmission point to the receiving point to local integral corresponding to a change of the receiving point depth, greatly reducing an integral calculation amount, so as to implement rapid calculation of t effective sound velocity in the deep sea.

To describe the technical solutions in the embodiments of the present invention or in the prior art more clearly, the following briefly introduces the accompanying drawing required for describing the embodiments. Apparently, the accompanying drawings in the following description show merely some embodiments of the present invention, and a person of ordinary skill in the art may still derive other drawings from these accompanying drawings without creative efforts.

FIG. 1 is a flowchart of a method for determining an effective sound velocity in the deep sea according to an embodiment of the present invention.

FIG. 2 is a schematic diagram of an apparatus for determining an effective sound velocity in the deep sea according to an embodiment of the present invention;

The following clearly and completely describes the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are merely a part rather than all of the embodiments of the present invention. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

An objective of the present invention is to provide a method for determining an effective sound velocity in the deep sea, so as to rapidly and accurately determine an effective sound velocity.

To make the foregoing objective, features, and advantages of the present invention clearer and more comprehensible, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

FIG. 1 is a flowchart of a method for determining an effective sound velocity in the deep sea according to an embodiment of the present invention. The method for determining an effective sound velocity in the deep sea is provided, where the method is applied to an apparatus for determining an effective sound velocity in the deep sea. As shown in FIG. 2 , the apparatus includes a transmission point 1 , a receiving point 2 , and an underwater mobile carrier 3 . The transmission point 1 is installed on the sea surface and the depth of the transmission point is unchanged. The receiving point 2 is installed on the underwater mobile carrier 3 and the depth of the receiving point changes with movement of the underwater mobile carrier 3 . The underwater mobile carrier 3 can measure a sound velocity profile between the transmission point 1 and the receiving point 2 and a horizontal distance between the transmission point and the receiving point 2 . As shown in FIG. 1 , the method includes the following steps:

Step 101 . Obtain a transmission point depth at a moment k and a receiving point depth at the k moment.

Step 102 . Obtain a receiving point depth at a moment k+1 and the horizontal distance between the transmission point and the receiving point.

Step 103 . Determine a pitch angle between the transmission point and the receiving point according to the transmission point depth, the receiving point depth at the moment k+1, and the horizontal distance between the transmission point and the receiving point.

Step 104 . Obtain a sound velocity profile between the transmission point and the receiving point at the moment k+1.

Step 105 . Obtain a difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1.

Step 106 . Determine an arithmetic average sound velocity variable quantity and a mean square arithmetic average sound velocity variable quantity according to the receiving point depth at the moment k+1, the receiving point depth at the k moment, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, and the sound velocity profile between the transmission point and the receiving point at the moment k+1.

Step 107 . Obtain a depth difference between the receiving point and the transmission point at the moment k.

Step 108 . Determine an arithmetic average sound velocity at the moment k+1 and a mean square average sound velocity at the moment k+1 according to a depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, an arithmetic average sound velocity at the moment k, a mean square average sound velocity at the moment k, the arithmetic average sound velocity variable quantity, and the mean square arithmetic average sound velocity variable quantity.

Step 109 . Determine an n-order sound velocity function variable quantity (n=2,3,4) according to the sound velocity profile between the transmission point and the receiving point at the moment k+1, the receiving point depth at the moment k, the receiving point depth at the moment k+1, and the arithmetic average sound velocity at the moment k+1.

Step 110 . Determine an n-order sound velocity function (n=2,3,4) at the moment k+1 according to the n -order sound velocity function variable quantity, the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the sound velocity profile between the transmission point and the receiving point at the moment k+1, and the arithmetic average sound velocity.

Step 111 . Calculate an n-order sound velocity coefficient (n=2,3,4) at the moment k+1 according to a pitch angle at the moment k+1.

Step 112 . Determine an effective sound velocity at the moment k+1 according to the n-order sound velocity coefficient at the moment k+1, the n-order sound velocity function corresponding to the n-order sound velocity coefficient at the moment k+1, and the arithmetic average sound velocity at the moment k+1.

determining the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1 according to the receiving point depth at the moment k and the receiving point depth at the moment k+1 by using a formula Δz k =z B,k+1 −z B,k , where

66 k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; z B,k is the receiving point depth at the moment k; and z B,k+1 is the receiving point depth at the moment k+1.

determining the depth difference between the receiving point and the transmission point at the moment k according to the receiving point depth at the moment k and the transmission point depth at the moment k by using a formula Z ABk =z B,k −Z A , where

z ABk is the depth difference between the receiving point and the transmission point at the moment k; z A is the transmission point depth; and z B,k is the receiving point depth at the moment k.

respectively determining the arithmetic average sound velocity variable quantity and the mean square arithmetic average sound velocity variable quantity according to the receiving point depth at the moment k+1, the receiving point depth at the moment k, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, and the sound velocity profile between the transmission point and the receiving point at the moment k by using a formula ΔV k =∫ z B,k z B,k +Δz k c(z)dz and a formula ΔW k =∫ z B,k z B,k +Δz k c 2 (z)dz, where

ΔV k is the arithmetic average sound velocity variable quantity; ΔW k is the mean square arithmetic average sound velocity variable quantity; z B,k +Δz k is the receiving point depth at the moment k+1; z B,k is the receiving point depth at the moment k; Δz k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; and c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1.

obtaining a sound velocity profile at the moment k;

determining the arithmetic average sound velocity at the moment k and the mean square average sound velocity at the moment k according to the sound velocity profile c(z) at the moment k, the transmission point depth z A , and the receiving point depth z B,k ; and

determining the arithmetic average sound velocity at the moment k+1 and the mean square average sound velocity at the moment k+1 according to the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the arithmetic average sound velocity at the moment k, the mean square average sound velocity at the moment k, the arithmetic average sound velocity variable quantity, and the mean square arithmetic average sound velocity variable quantity by using a formula

v _ k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁢ z [ ABk ) ⁢ vk + Δ ⁢ ⁢ Vk

and a formula

w _ k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁡ [ z ABk ⁢ w _ k + Δ ⁢ ⁢ W k ] ,

where

v k+ 1 is the arithmetic average sound velocity at the moment k+1; w k+1 is the mean square average sound velocity at the moment k+1; z ABk is the depth difference between the receiving point and the transmission point at the moment k+1; v k is the arithmetic average sound velocity at the moment k; w k is the mean square average sound velocity at the moment k; ΔV k is the arithmetic average sound velocity variable quantity; and ΔW k is the mean square arithmetic average sound velocity variable quantity.

determining the n-order sound velocity function variable quantity (n=2,3,4) according to the sound velocity profile between the transmission point and the receiving point at the moment k+1, the receiving point depth at the moment k, the receiving point depth at the moment k+1, and the arithmetic average sound velocity at the moment k+1 by using a formula

Q k = ∫ z B , k z B , k + 1 ⁢ ( c ⁡ ( z ) - v _ k + 1 v _ k + 1 ) n ⁢ dz ,

where

Q k is the n-order sound velocity function variable quantity; c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1; z B,k is the receiving point depth at the moment k; z B,k+1 is the receiving point depth at the moment k+1; and v k+1 is the arithmetic average sound velocity at the moment k+1.

determining the n-order sound velocity function (n=2,3,4) at the moment k+1 according to the n-order sound velocity function variable quantity, the depth difference between the receiving point and the transmission point at the moment k+1, the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1, the sound velocity profile between the transmission point and the receiving point at the moment k+1, and the arithmetic average sound velocity by using a formula

J n k + 1 = 1 z ABk + Δ ⁢ ⁢ z k ⁡ [ ∫ z A z B , k ⁢ ( c ⁡ ( z ) - v _ k + 1 v _ k + 1 ) n ⁢ dz + Q k ] ,

where

∫ n k+1 is the n-order sound velocity function corresponding to the n-order sound velocity coefficient at the moment k+1; Q k is the n-order sound velocity function variable quantity; z ΔBk is the depth difference between the receiving point and the transmission point at the moment k+1; Δz k is the difference between the receiving point depth at the moment k and the receiving point depth at the moment k+1; c(z) is the sound velocity profile between the transmission point and the receiving point at the moment k+1; and v k+1 is the arithmetic average sound velocity at the moment k+1.

calculating a second-order velocity coefficient, a third-order sound velocity coefficient, and a fourth-order sound velocity coefficient according to the pitch angle at the moment k+1 by using formulas: C 2 k+1 =1−½cot 2 β k+1 ; C 3 k+1 =−C 2 k+1 =½cot 4 β k+1 ; and C 4 k+1 =⅛cot 4 β k+1 (3−5cot 2 β k+1 )+C 2 k+1 , where

β k+1 is the pitch angle at the moment k+1, C 2 k+1 , C 3 k+1 , and C 4 k+1 , are respectively the second-order sound velocity coefficient, the third-order sound velocity coefficient, and the fourth-order sound velocity coefficient.

determining the effective sound velocity at the moment k+1 according to a formula v e k+1 = v k+1 (1−C 2 k+1 J 2 k+1 −C 3 k+1 −C 4 k+1 J 4 k+1 ), where

v e k+1 is the effective sound velocity at the moment k+1; v k+1 , is arithmetic average sound velocity at the moment k+1; J 2 k+1 , J 3 k+1 , and J 4 k+1 are respectively a second-order sound velocity function, a third-order sound velocity function, and a fourth-order sound velocity function, and C 2 k+1 , C 3 k+1 , and C 4 k+1 are respectively the second-order sound velocity coefficient, the third-order sound velocity coefficient, and the fourth-order sound velocity coefficient.

Integral operations of Formula (2) and Formula (4) are converted from whole integral to local integral in a depth change area, to reduce an integral calculation amount of Formula (2) and Formula (4), so as to reduce a calculation amount of the effective sound velocity v e . Only in this way, in an underwater acoustic positioning system, when a depth of the transmission point or the receiving point changes, a corresponding effective sound velocity v e can be calculated very rapidly, and then a measured propagation time is converted into a slant distance between the transmission point and the receiving point, and finally implements real-time underwater positioning.

In the present invention, it is assumed that a position of a sound source is near a sea surface and is stationary; a position of a receiving point is below the position of the sound source and the position can be moved; when both the sound source and the receiving point change, a specific method of the present invention is obtained through the following derivation method.

(1) Calculate an arithmetic average sound velocity v k+1 at the moment k+1 through the arithmetic average sound velocity v k at the moment k.

From the moment k to the moment k+1, the corresponding depth difference Δz k of the receiving point is as follows:

Δ z k =z B,k+1 −z B,k   (5)

The corresponding depth difference z ABk between the receiving point and the transmission point at the moment k is as follows:

z ABk =z B,k −z A   (6)

The following formulas can be obtained through Formula (2a), Formula (2b), and Formula (6):

∫ z A z B,k c ( z ) dz=z ABk v k   (7a):

∫ z A z B,k C 2 ( z ) dz=z ABk w k   (7b)

At the moment k+1, the arithmetic average sound velocity v k+1, and the mean square arithmetic average sound velocity Δ w k+1 are respectively as follows:

For the receiving point, from the moment k to the moment k+1, let the arithmetic average sound velocity variable quantity ΔV k and the mean square arithmetic average sound velocity variable quantity ΔW k respectively satisfy the following:

Δ V k =∫ z B,k z B,k +Δz k c ( z ) dz   (9a)

Δ W k =∫ z B,k z B,k +Δz k c 2 ( z ) dz   (9b)

In this way, Formula (8) changes to

In this case, an integral operation of Formula (2) is converted from original whole integral of the transmission point to the receiving point to local integral corresponding to a change of the receiving point depth, greatly reducing an integral calculation amount.

(2) Calculate J n k (n=2,3,4) at the moment k+1 through J n k (n=2,3,4) at the moment k.

The following formula can be obtained through Formula (4) and Formula (7):

At the moment k+1, an n-order J function J n k+1 (n=2,3,4) is as follows:

For the receiving point, from the moment k to the moment k+1, let an n-order J function variable quantity Q k satisfy the following:

In this case, Formula (12) changes to:

A function ƒ( v k ) is defined as follows:

A integral of the function ƒ( v k ) is as follows:

According to the Taylor expansion, a function ƒ( v k+1 ) at the moment k+1 is as follows:

ƒ( v k+1 )≈ƒ( v k )+ƒ′( v k )·Δ v +½ƒ″( v k )·Δ v 2 +⅙ƒ′″( v k )·Δ v 3 + 1/24 ƒ (4) ) v k )·Δ v v 4   (17)

Δ v k = v k+1 − v k   (18)

Substitute Formula (16a-d) is substituted into Formula (17) to conduct a series of algebraic operations, and the function ƒ( v k+1 ) at the moment k+1 may be simplified to:

To further simplify Formula (19), let

In this case, a function ƒ( v k+1 ) at the moment k+1 is as follows:

ƒ( v k+1 )≈ƒ( v k )+Φ k ·c ( z )+Ψ k ·c 2 ( z )   (22)

Formula (22) is substituted into Formula (12), and the n-order J function J n k+1 (n=2,3,4) at the moment k+1 is as follows:

J n k+1 Formula (7) and Formula (11) are substituted into Formula (23), and the n-order J function J n k+1 (n=2,3,4) at the moment k+1 is as follows:

In this case, an integral operation of Formula (4) is converted from original whole integral of the transmission point to the receiving point to local integral of Formula (13) corresponding to a change of the receiving point depth, greatly reducing an integral calculation amount.

Each embodiment of the present specification is described in a progressive manner, each embodiment focuses on the difference from other embodiments, and the same and similar parts between the embodiments may refer to each other.

Several examples are used for illustration of the principles and implementation methods of the present invention. The description of the embodiments is used to help illustrate the method and its core principles of the present invention. In addition, those skilled in the art can make various modifications in terms of specific embodiments and scope of application in accordance with the teachings of the present invention. In conclusion, the content of this specification shall not be construed as a limitation to the invention.