Next Article in Journal
Satellite Edge Computing for the Internet of Things in Aerospace
Previous Article in Journal
Traversability Assessment and Trajectory Planning of Unmanned Ground Vehicles with Suspension Systems on Rough Terrain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Ultra-Short Baseline Positioning Model Based on Rotating Array & Reusing Elements and Its Error Analysis

Key Laboratory of Micro-Inertial Instruments and Advanced Navigation Technology, Ministry of Education, School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(20), 4373; https://doi.org/10.3390/s19204373
Submission received: 13 August 2019 / Revised: 27 September 2019 / Accepted: 4 October 2019 / Published: 10 October 2019
(This article belongs to the Section Sensor Networks)

Abstract

:
The USBL (Ultra-Short Base Line) positioning system is widely used in underwater acoustic positioning systems due to its small size and ease of use. The traditional USBL positioning system is based on ‘slant range and azimuth’. The positioning error is an increasing function with the increase in distance and the positioning accuracy depends on the ranging accuracy of the underwater target. This method is not suitable for long-distance underwater positioning operations. This paper proposes a USBL positioning calculation model based on depth information for ‘rotating array and reusing elements’. This method does not need to measure the distance between the USBL acoustic array and target, so it can completely eliminate the influence of long-distance ranging errors in USBL positioning. The theoretical analysis and simulation experiments show that the new USBL positioning model based on ‘rotating array and reusing elements’ can completely eliminate the influence of the wavelength error and spacing error of underwater acoustic signals on the positioning accuracy of USBL. The positioning accuracy can be improved by approximately 90%, and the horizontal positioning error within a positioning distance of 1000 m is less than 1.2 m. The positioning method has high precision performance in the long distance, and provides a new idea for the engineering design of a USBL underwater positioning system.

1. Introduction

Due to the severe attenuation of electromagnetic waves under water, the application of high precision satellite navigation systems underwater is limited. As the positioning error of the INS (Inertial Navigation System) accumulates with time, one single INS cannot complete the high precision navigation and positioning task in the long-range alone [1]. Acoustic waves are the only information carriers that have been proven to be able to travel long distances under water [2]. Underwater acoustic positioning is the main means of underwater target positioning and tracking at present [3]. Pan-Mook Lee and others proposed the range aiding method to integrate with IMU(Inertial Measurement Unit) /DVL (Doppler Velocity Log) [4]. According to the length of the baseline, an underwater acoustic positioning system can be divided into LBL (Long Base Line), SBL (Short Base Line), and USBL, where USBL is a variant of SBL [5]. By virtue of the small size and convenient use, USBL is widely used in underwater positioning. At present, there are mature USBL products from Kongsberg in Norway; IXBlue in France; LinkQuest in the United States; and AAE, Nautronix, and Sonardyne in the UK. Among them, Kongsberg is one of the earliest research institutes for USBL acoustic positioning systems. Kongsberg launched the first generation USBL in 1996 and their new HiPAP700 is one of a few long-range USBL positioning systems where its working water depths can reach 10,000 m. However, a single SINS (Strap-down Inertial Navigation System) or USBL cannot complete underwater positioning tasks with high precision performance in long distance submerged travel. The SINS/USBL integrated navigation mode is widely adopted in existing underwater positioning products [6]. GAPS-USBL, as a high precision USBL positioning system launched by Ixblue in France, adopts a combined design of USBL/INS /GNSS (Global Navigation Satellite System). This integrated structure eliminates the installation error calibration between a transceiver array of USBL equipment and INS devices before use. Therefore, the influence of installation error angle on the positioning accuracy of USBL is reduced [7]. In [8,9,10,11], the loose combination and tight combination were adopted in the SINS/USBL integrated navigation positioning algorithm. Bingbing Gao and other researchers from Northwestern Polytechnical University proposed the SINS/USBL/DVL integrated navigation scheme [12]. Jianhua Cheng, Ping Dong and others from Harbin Engineering University proposed an INS/USBL nonlinear combination method based on relative measurement information for nonlinear problems [13]. The filtering theory of SINS/USBL integrated navigation is based on the mature filtering algorithms of SINS and other observation information such as SINS/GNSS [14] and SINS/DVL [15]. Since there is an installation error angle between the base surface of the acoustic array and the mounting base surface. The calibration is indispensable before used. Professor Hsin-Hung Chen from the National Sun Yat-Sen University of Taiwan proposed a multi-circle sailing and straight-line sailing method, which was verified experimentally that to be feasible [16,17]. A USBL installation error angle calibration method based on double transponders was proposed in 2018 [2], but the calibration process was more complicated. In 2019, a USBL installation error angle calibration single transponder USBL based on dual vector reconstruction was proposed, where high-precision and a fast calibration of the USBL installation error angle can be achieved [18].
The traditional USBL system adopts the narrow-band underwater acoustic signal design, which uses the phase difference or phase ratio method for the positioning calculation. In general, the larger the spacing of acoustic array elements, the higher the positioning accuracy of the USBL. Unfortunately, the spacing of the acoustic array elements cannot be designed to be infinite for two reasons. The main reason is that the large-scale array structure of the USBL loses its compaction, which is inconvenient to use. The other is that when the spacing of the acoustic array elements is greater than half a wavelength, the phase difference ambiguity problem will occur, which is a very troublesome problem. On this issue, Min Yu, a scholar from Wuhan University of Technology, proposed a method of expanding the spacing of array elements by adding auxiliary elements in her PhD thesis (Harbin Engineering University, China, 2006) [19]. This method not only improves the positioning accuracy of USBL, but also solves the problem of phase difference ambiguity. However, this method creates new problems that with the increase of the number of acoustic array elements, a large amount of work on the USBL element spacing error calibration is needed. In addition, as the number of array elements increases, the complexity of the manufacturing of USBL acoustic array devices increases, which will reduce the reliability of the USBL positioning system [20]. Therefore, a method to improve the accuracy of the USBL positioning calculation has been a research hotspot in this field. From the existing reference, the main method to improve the accuracy of the USBL positioning solution is to calibrate the installation error angle [16,21,22] and the spacing error between the acoustic array elements [23]. The existing technology can reduce the spacing error of the USBL acoustic array elements to less than 1 mm or even less than 0.12 mm [23], but the calibration process is extremely complex. With the increase of the usage time, the spacing of the elements will be changed. In particular, after replacing the acoustic array elements, this parameter needs to be corrected again. Yunhui Wang, a graduate student of the research team of Guiqing Sun from the Oceanography School of Zhejiang University proposed a USBL positioning model based on depth information [20]. In the long-distance situation, the model has better positioning accuracy than the traditional method base on slant range and the azimuth, and it is not needed to measure the distance, which improves the positioning efficiency, so it is especially suitable for the passive positioning of cooperative targets. However, the spacing error between different acoustic array elements is neglected, so this method fails to eliminate the influence of the element spacing error on the accuracy of USBL positioning. How to eliminate the influence of the spacing error array of elements in USBL positioning calculation is of great significance to improve the accuracy of the USBL positioning system.
To overcome the shortcoming where the positioning error is correlated with the positioning distance in the traditional USBL positioning algorithm, a new model for USBL positioning calculation has been designed. The passive positioning of the cooperative target is considered as the positioning object of this model researched in this paper. That is, there is no acoustic transmitter in the USBL acoustic array. The USBL positioning system does not emit underwater acoustic signals. The USBL positioning calculation depends on the depth information of the target and the phase difference between the acoustic array elements. This model can perform a USBL positioning calculation without slant distance measurement and can completely eliminate the influence of the spacing error on positioning accuracy, thus achieving a high-precision USBL positioning calculation.
The rest of this paper is structured as follows. Section 2 illustrates the basic working principle of the traditional USBL positioning system and the characteristics of the spacing error of the acoustic array elements, then analyses the horizontal positioning error of two positioning models based on the slant range and azimuth and phase difference ratio. Section 3 gives the design and error analysis of a USBL positioning model based on the rotating array and array element reusing in detail. This is also the main contribution of this paper. Section 4 presents the simulation and verification, and the numerical simulations of the proposed model are carried out and compared with the traditional positioning model. Section 5 is the conclusion and outlook.

2. The Principle and Error Analysis of Traditional USBL Positioning System Based on Distance and Angle of Incidence

2.1. Coordinate Frame Definition

Multiple frames were involved in this paper: a USBL array frame O U - X U Y U Z U (abbreviated as U -frame) with the Z axis pointing upward; the three coordinate axes were set in the right-hand frame, and the coordinate axes were orthogonal to each other. Receiving array elements 1, 2, 3, and 4 were located on the x-axis and y-axis, respectively, as shown in Figure 1. It also can be seen that the acoustic emitter was located at the origin of the U -frame, which is the acoustic center of the acoustic array (there is no acoustic emitter in the USBL acoustic array with passive location mode). The USBL position calculation of the target was carried out in the U -frame.
The body frame (or base frame O B X B Y B Z B , abbreviated as the B -frame) refers to the frame of the carrier where the USBL transceiver array equipment is installed. The positive direction of the y-axis points to the front of the bow line. The x-axis points to the starboard side. The z-axis points upward, which is perpendicular to the plane x o z . Ideally, the B -frame coincides with the origin of the U -frame, and the axes are also completely parallel. However, in engineering application, the origins of the frames are not completely coincident, and there is an angular deviation between their corresponding axes, that is, the angle misalignment error of the USBL, which is shown as Figure 2.
In the global frame ( O G X G Y G Z G , abbreviated as the G -frame), the x-axis points to the east, the y-axis points to the north, and the z-axis points upward. The WGS-84 frame, where the absolute position of the target can be obtained, was adopted in this paper.

2.2. The Basic Principle of the USBL Positioning System

It can be seen from Figure 1 that the basic equations of the USBL positioning algorithm based on the slant range and azimuth are as follows.
x u = R cos θ x
y u = R cos θ y
z u = R 2 X u 2 Y u 2 = R 1 cos 2 θ x cos 2 θ y
R = C T
where x u , y u , and z u are the calculated values of the USBL positioning in the U -frame. θ x is the angle between the slant range and the x-axis. θ y is the angle between the slant range and the y-axis. θ z is the angle between the slant range and the z-axis. R is the distance between the center of the acoustic array and target. C is the equivalent underwater acoustic velocity and T is the one-way propagation time of the acoustic signal.
Figure 3 is the diagram of the USBL phase difference.
φ = 2 π d λ cos θ
λ = C / f
Equation (5) represents the phase difference between the two array elements, where d is the distance between the two array elements; θ is the direction angle of the acoustic signal; λ is the wavelength of the underwater acoustic signal; f is the frequency of the underwater acoustic signal; and C is the underwater acoustic velocity. From Equation (5), it can be obtained that
θ x = arccos ( λ φ 1 , 2 2 π d 1 , 2 )
θ y = arccos ( λ φ 3 , 4 2 π d 3 , 4 )
Substituting Equations (7) and (8) into Equations (1) and (2)
x u = R λ φ 1 , 2 2 π d 1 , 2
y u = R λ φ 3 , 4 2 π d 3 , 4
where λ is the wavelength; φ 1 , 2 , φ 3 , 4 are the phase difference of 1–2, and 3–4, respectively; and d 1 , 2 , d 3 , 4 are the spacing of array elements 1–2 and 3–4, respectively. Equations (9), (10), and (3) are the basic formulas of the USBL positioning system.

2.3. Array Element Spacing Error Modeling

The array element spacing error is one of the main error sources of the USBL positioning calculation. The high-precision calibration and compensation for the element spacing error can effectively suppress the USBL positioning error. High-precision calibration methods can be used in most underwater positioning applications; however, under the existing technical conditions, due to the limited calibration accuracy of the spacing error, the existing calibration methods in high-precision and ultra-high-precision underwater applications still cannot meet the requirements of positioning accuracy. How to suppress or completely eliminate the influence of elements spacing error on USBL positioning accuracy is still a hot topic in this field. This section establishes a USBL array elements spacing error model to study the characteristics of the USBL array elements spacing error.
The USBL array elements spacing error can be divided into the measurement error δ d m and the device drift error δ d d . The measurement error δ d m is the main error, which has a direct impact on the accuracy of the USBL positioning calculation. The error is constant in one measurement. Device drift error δ d d refers to the spacing error caused by the acoustic center drift of the element and the change in device geometry caused by rust and aging, which is very small and has little change over a long time. It usually has little influence on the positioning accuracy in engineering applications and therefore was neglected in this paper.
d ˜ = d + δ d
δ d = δ d m + δ d d
δ d ˙ m = 0
where d ˜ is the measured value of the element spacing; d is the truth-value of the element spacing; δ d is the error of the element spacing; δ d m is the measurement error which is the constant error; and δ d d is the drift error. Figure 4 shows the schematic diagram of the USBL receiving array based on six receiving array elements.
The array elements 1, 2, 3, 4, 5, and 6 are all located on the coordinate axis. The intersection point of the connection of array elements 1–2, 3–4, and 5–6 is the coordinate origin. The spacing of array elements 1–2, 3–4, and 5–6 are as follows:
d ˜ 1 , 2 = d 1 , 2 + δ d 1 , 2
d ˜ 3 , 4 = d 3 , 4 + δ d 3 , 4
d ˜ 5 , 6 = d 5 , 6 + δ d 5 , 6
where d ˜ i , j is the measurement values of the spacing between array elements i and j ; d i , j is the true value of the spacing between array elements i and j ; δ i , j is the spacing error between array elements i and j .

2.4. The Error Analysis of the USBL Positioning Solution Based on the Slant Range and Azimuth Method

Using the complete differential method in Equation (9), the positioning error equation of the USBL positioning solution in the x-axis can be obtained.
δ x u = R λ 2 π d 1 , 2 δ φ 1 , 2 + R φ 1 , 2 2 π d 1 , 2 δ λ + λ φ 1 , 2 2 π d 1 , 2 δ R R λ φ 1 , 2 2 π δ d 1 , 2 d 1 , 2 2
Considering that each error item is independent of each other, the mean square error of the target location solution value in the x-axis direction is as follows:
δ d x 2 = ( R λ φ 1 , 2 2 π d 1 , 2 ) 2 ( ( δ φ 1 , 2 φ 1 , 2 ) 2 + ( δ λ λ ) 2 + ( δ R R ) 2 + ( δ d 1 , 2 d 1 , 2 ) 2 )
where δ φ 1 , 2 is the phase difference error; δ λ is the wavelength error; and δ R is the ranging error and δ d 1 , 2 is the spacing error between the array elements 1 and 2. δ R / R is the relative distance error factor; δ λ / λ is the relative wavelength error factor; and δ d 1 , 2 / d 1 , 2 is the relative spacing error factor.
In order to analyze the influence of errors, Equation (17) is rewritten in the form of relative errors.
δ x u R = λ φ 1 , 2 2 π d 1 , 2 ( δ R R + δ λ λ δ d 1 , 2 d 1 , 2 ) + λ δ φ 1 , 2 2 π d 1 , 2
According to Equation (5) rewriting Equation (19), it can be obtained that
δ x u R = cos θ x ( δ R R + δ λ λ δ d 1 , 2 d 1 , 2 ) + λ δ φ 1 , 2 2 π d 1 , 2
With Equations (4) and (6), Equation (20) can be written as follows:
δ x u R = cos θ x ( δ T T + 2 δ C C δ d 1 , 2 d 1 , 2 ) + λ δ φ 1 , 2 2 π d 1 , 2
where δ T is the time delay measurement error; δ C is the equivalent acoustic velocity estimation error; δ T / T is the relative time delay measurement error factor; and δ C / C is the relative equivalent acoustic velocity estimation error factor. From Equation (20), it can be seen that the USBL positioning error is caused by the time delay measurement error δ T , the equivalent acoustic velocity estimation error δ C , and the spacing error of array elements δ d 1 , 2 is related to the incident angle θ x . When the target is located right below the acoustic array, that is, the value of the incidence angle θ x is close to 90°, the positioning error of the USBL is less affected by these three errors and is mainly affected by the phase difference error. Therefore, the positioning accuracy of the USBL has the best performance in the range of cone angle directly below it.
When the target is not directly under the USBL transceiver array, in addition to the phase difference error δ φ , the USBL positioning error caused by the time delay measurement error δ T , the equivalent acoustic velocity error δ C , and the spacing error of the array elements δ d should not be neglected. Based on the performance parameters of the current products, the effects of these three errors were analyzed and the details are shown in Table 1.
Considering the location environment with a location distance of 1000 m, an equivalent underwater acoustic velocity of 1500 m/s, and a distance of the USBL array elements of 250 mm, which is shown in Table 2.
The relative time delay measurement error factor is as follows:
δ T T 0.0015
where δ T is the time delay error term, and its relative error decreases with the increase in the positioning distance without considering the acoustic line bending. The relative equivalent acoustic velocity estimation error factor is
2 δ C C =   0.0000666
The relative spacing error factor is
δ d d 1 , 2 = 0.008
The equivalent acoustic velocity estimation error term and the spacing error term are relatively stable values. Comparing Equations (22)–(24), it can be seen that the relative spacing error factor was much larger than the relative time delay measurement error factor and the relative equivalent acoustic velocity estimation error factor. Therefore, it is indispensable to calibrate the spacing error of the array elements δ d prior to using it in a high-precision underwater acoustic positioning environment.

2.5. The Error Analysis of USBL Positioning Calculation Based on the Phase Ratio Method

In Figure 4, taking the array elements spacing error into account, the phase differences between two array elements on the same coordinate axis are as follows:
φ 1 , 2 = 2 π d ˜ 1 , 2 cos θ x λ
φ 3 , 4 = 2 π d ˜ 3 , 4 cos θ y λ
φ 5 , 6 = 2 π d ˜ 5 , 6 cos θ z λ
According to the calculation method of the USBL based on the slant range and azimuth, it can be obtained that
cos θ x = x R
cos θ z = z R = h R
where h is the depth of the target. From Equations (25), (27), and (28), it can be derived that
x ˜ = R cos θ x = h cos θ z cos θ x = h φ 1 , 2 φ 5 , 6 d ˜ 5 , 6 d ˜ 1 , 2
x ˜ = x d ˜ 5 , 6 d ˜ 1 , 2
where x ˜ is the coordinate value of the target in the U -frame, only considering the influence of the element spacing error on the USBL positioning accuracy and x is the coordinate value of the target in the USBL U -frame without considering each error term. In [20], the array elements spacing error was ignored. In fact, not only is there a spacing error between the elements, but the spacing error between the different elements is also different, that is d 1 , 2 d 3 , 4 d 5 , 6 . From Equations (14) and (16), it can be known that
d ˜ 5 , 6 d ˜ 1 , 2 = d + δ d 5 , 6 d + δ d 1 , 2 = d d + δ d 1 , 2 + δ d 5 , 6 d + δ d 1 , 2 1 + δ d 5 , 6 d + δ d 1 , 2
Substituting Equation (32) into Equation (31), it can be obtained that
x ˜ x ( 1 + δ d 5 , 6 d + δ d 1 , 2 )
It can be seen from Equation (33) that there are two effective methods to improve the horizontal positioning accuracy, reduce the element spacing error δ d 5 , 6 , and expand the element spacing δ d . Numerical simulations based on these two methods are performed below.

2.6. The Simulation and Analysis of the Spacing Error of Array Elements

In order to study the influence of the spacing error of array elements on the horizontal positioning accuracy, the influence of the spacing error of array elements on the horizontal positioning accuracy in the horizontal positioning range of 4000 m was studied, regardless of other error terms. Research was conducted on the horizontal positioning error (x-axis direction) caused by the spacing error of array elements of 250 mm with errors of 2 mm and 1 mm. The calibration method studied in [23] could reduce the spacing error of array elements less than 0.12 mm after compensation. This paper also used the numerical simulations and analyses of the spacing error of array elements of 12 mm.
It can be seen from Table 3. that by using the current ideal element calibration method [23] in a positioning environment where the horizontal distance is 4000 m and the element spacing error is less than 0.12 mm, the horizontal positioning error caused by the elements spacing error was 1.9191 m, a small positioning error that is hard to achieve in mainstream products. With a distance of 1000 m, the 1 mm elements spacing error could bring about a 4 m horizontal positioning error, and the 2 mm element spacing error could bring about a 7.9 m horizontal positioning error. This is an intolerable error term in high precision underwater positioning environments. For example, in order to reduce the positioning accuracy to 2 m with the slanting distance of 1000 m, the depth of 100 m, and the horizontal distance of about 995 m by only considering the elements spacing error, the element spacing error should be less than 0.5 mm. However, when considering the influence of other error items in the engineering project, the elements spacing error value should be smaller to meet the design requirements.
It is mentioned in the literature that increasing the spacing of the array elements can suppress the influence of the array elements spacing error on the positioning of the USBL. In the case of d = 500 mm, a simulation was also carried out. The details are shown as Table 4.
It can be seen from Table 3 and Table 4 that when the spacing of the array elements is doubled, the USBL horizontal positioning error term caused by the equivalent element spacing error was about half of the original. Therefore, increasing the spacing of the array elements can indeed suppress the USBL horizontal positioning error caused by the spacing error of the array elements. However, the consequence of the large spacing of the array elements is that the size of the USBL array becomes larger, so the portability of the USBL will be affected. In addition, it also brings about a phase difference ambiguity problem that needs to be solved. Furthermore, as the spacing of the array elements increases, the spacing error of the array elements will also increase. Therefore, increasing the spacing of the array elements is not an ideal means by which to effectively suppress the influence of the spacing error of array elements on the USBL positioning.

3. A USBL Positioning Model and Its Error Analysis Based on Rotating Array and Reusing Elements

3.1. The Design of the Algorithmic Model

The method based on phase ratio in [20] neglected the influence of the spacing error of the array element. In this paper, we considered that there was a spacing error between the array elements, and that there were differences in the spacing errors between different groups of array elements, as seen in Figure 4.
δ d 1 , 2 δ d 3 , 4 δ d 5 , 6
Since the truth-value of the spacing of array elements d i , j is unknown, and the spacing error of array elements δ d i , j is also difficult to calibrate accurately, this paper designed a USBL positioning calculation model based on the rotating array and reusing elements method, which is used to eliminate the influence of the spacing error of array elements on the USBL positioning system.
According to Equations (14)–(16) and (34), the purpose of this model is to equalize the spacing errors between the three elements groups, that is,
δ d 1 , 2 = δ d 3 , 4 = δ d 5 , 6
The spacing error between array elements is difficult to calibrate accurately. To make the three spacing errors completely equal, the most appropriate method is to use the elements reusing method. Thus, array elements 3, 4, 5, and 6 were removed, leaving only array elements 1 and 2, as shown in Figure 5. When it is necessary to calculate the phase difference φ 3 , 4 and the phase difference φ 5 , 6 , the coordinate system is rotated so that array elements 1 and 2 on the coordinate axis are respectively rotated to the positions of elements 3–4 and 5–6, and the phase differences φ 3 , 4 and φ 5 , 6 are respectively calculated. In this process, only two array elements connected rigidly are used. In the whole positioning process, the spacing error between the two array elements is always δ d 1 , 2 .
Figure 5 shows the USBL structure based on virtual six receiving elements array. The array was composed of array elements 1 and 2, which were connected rigidly. The positioning principle of the USBL acoustic receiving array was the same as that of the 6-element USBL receiving array shown in Figure 4, and the error principle also was the same.
Figure 6 shows the rotating schematic diagram of the USBL array based on the rotating array and reusing elements. Figure 6a is the schematic diagram of the U -frame rotating 90° clockwise around the y-axis, so that array elements 1 and 2 of the original coordinate x-axis can rotate to the position of array elements 5 and 6 on the z-axis. Similarly, Figure 6b is the schematic diagram of the U -frame rotating 90° clockwise around the x -axis, so that array elements 5 and 6 on the original z-axis can rotate to the position of array elements 3 and 4 on the y-axis. Figure 6c is the schematic diagram of the USBL array U -frame rotating 90° clockwise around the z-axis, so that array elements 3 and 4 on the original y-axis can rotate back to the position of the array elements 1 and 3 on the x-axis.
The calculation process is shown as follows:
Figure 7 depicts the flow chart of the USBL positioning calculation model based on the rotating array and reusing elements. First, the USBL positioning system acquires the depth information of the target h t . Then, the USBL array obtains the phase difference φ 1 , 2 between array elements 1 and 2 by receiving the underwater acoustic signal from the target.
After acquiring the phase difference φ 1 , 2 , the U -frame is rotated around the y-axis to rotate the x-axis to the position of the z-axis, and the phase difference φ 5 , 6 is measured. Similarly, the U -frame is rotated from the z-axis around the x-axis to the position of the y-axis, and the phase difference φ 3 , 4 is measured. After the three phase differences φ 1 , 2 , φ 5 , 6 , and φ 3 , 4 are collected, the i-th USBL horizontal positioning calculation is performed. Then, the U -frame is rotated around the z-axis and the y-axis is rotated back to the x-axis. At this time, if the positioning task is not finished, the program enters the next round of the target positioning solution, otherwise the positioning task is completed and ends this procedure.
In Figure 6a,
φ 1 , 2 = 2 π d ˜ 1 , 2 cos θ x λ
where φ 1 , 2 is the phase difference between array elements 1 and 2; d ˜ 1 , 2 is the spacing value between array elements 1 and 2, which contains the constant spacing error δ d 1 , 2 ; and θ x is the angle between the underwater acoustic signal line and the x-axis.
In Figure 6b,
φ 5 , 6 = 2 π d ˜ 5 , 6 cos θ z λ
where φ 5 , 6 is the phase difference between array elements 5 and 6; d ˜ 5 , 6 is the spacing value between array elements 5 and 6, which contains the constant spacing error δ d 5 , 6 ; and θ z is the angle between the underwater acoustic signal line and the z-axis.
It can be derived from Equation (30) that
x ˜ = h φ 1 , 2 φ 5 , 6 d ˜ 5 , 6 d ˜ 1 , 2
Since array elements 5 and 6 are obtained by rotating array elements 1 and 2 of the x-axis on the USBL acoustic array around the y-axis to the z-axis, and the two elements are rigidly connected before and after the rotation, it can be seen that d ˜ 1 , 2 = d ˜ 5 , 6 . Thus, Equation (38) can be written as
x = x ˜ = h φ 1 , 2 φ 5 , 6
The calculation of the y-axis coordinate value of the target T in the USBL U -frame is performed as follows. First, the USBL matrix U -frame is rotated around the z-axis, so that array elements 1 and 2 on the x-axis are rotated to the positions of array elements 3 and 4 on the y-axis, and the phase difference φ 3 , 4 is measured. Then, the U -frame is rotated around the x-axis, so that array elements 3 and 4 on the y-axis are rotated to the positions of array elements 5 and 6 on the z-axis, and the phase difference φ 5 , 6 is measured. Finally, the y-axis coordinate value of the target T in the U -frame is calculated according to the depth value h of the target T and the phase difference values φ 3 , 4 and φ 5 , 6 . The USBL acoustic array structure and its rotation sequence are shown in Figure 8.
In order to save the positioning calculation time and adapt to the dynamic positioning environment, the horizontal position of the target T can be solved in the order shown in Figure 7. After the USBL acoustic positioning system completes the x-axis positioning calculation, the phase difference φ 3 , 4 of array elements 3 and 4 on the y-axis can be obtained immediately. As shown in Figure 6b, after acquiring the phase difference φ 5 , 6 of array elements 5 and 6 on the z-axis, the U -frame rotates around x-axis, and rotates array elements 5 and 6 on the z-axis to the position of elements 3 and 4 on the y-axis to measure the phase difference φ 3 , 4 . Then, combined with the target depth value h , complete the positioning calculation of the target (y-axis direction) in the U -frame.
y = R cos θ y = h cos θ z cos θ y = h φ 3 , 4 φ 5 , 6
z = | h T h A |
where h T is the depth value of the target and h A is the depth value of the USBL acoustic array. Equations (39)–(41) are the basic equations for the USBL positioning calculation based on rotating array and reusing elements method. This positioning method is suitable for USBL positioning calculations where the target depth h T is known.
It can be seen from Equation (39) that the positioning accuracy of the USBL is only related to the depth value h T of the target and the error of the two phase differences φ 1 , 2 and φ 5 , 6 . The positioning accuracy is independent of the spacing of the array elements d i , j . Therefore, this model can completely eliminate the influence of the spacing error on the positioning accuracy of the USBL.

3.2. Error Analysis

In this section, the error source analysis of the USBL positioning model based on rotating array and reusing elements is performed. Apply the complete differential to Equation (39), and the positioning error (x-axis) of the target in the U -frame can be obtained.
δ x = δ h φ 1 , 2 φ 3 , 4 + h δ φ 1 , 2 φ 3 , 4 φ 1 , 2 δ φ 3 , 4 φ 3 , 4 2 = h φ 1 , 2 φ 3 , 4 ( δ h h + δ φ 1 , 2 φ 3 , 4 δ φ 3 , 4 φ 3 , 4 )
Under the condition that each error term is independent of each other, the mean square error of the target in the x-axis of the U -frame is
δ d x 2 = ( h φ 1 , 2 φ 3 , 4 ) 2 ( ( δ h h ) 2 + ( δ φ 1 , 2 φ 3 , 4 ) 2 + ( δ φ 3 , 4 φ 3 , 4 ) 2 )
Contrasting Equation (43) and Equation (18), it can be found that when compared with the USBL positioning calculation method based on the slant range and azimuth method, the proposed USBL positioning calculation method based on the rotating array and reusing elements method eliminates the horizontal positioning error caused by the signal wavelength error δ λ and the USBL element spacing error δ d . The slant range measurement error δ R is replaced by the depth measurement error δ h . In the case of long distance, the depth measurement error is much smaller than the slant measurement error.
In the case of high-precision underwater positioning, the USBL positioning model based on the rotating and array element reuse should consider the USBL horizontal positioning error caused by the rotary angle error of the rotating device. During the process of the rotary device rotation, it is difficult to rotate the rotary device at exactly 90°. When the line connecting the two array elements after rotation does not completely coincide with the target coordinate axis, the distance between the two array elements projected onto the coordinate axis will be less than the spacing value of the actual array elements. As shown in Figure 6a, the acoustic array coordinate system U is rotated about the y -axis, rotating the array element 1–2 on the x -axis to the positions of the array elements 5–6 on the z -axis. However, in the actual rotation process, when the x -axis and the z -axis do not completely coincide, there is a small error angle δ θ x z between them, as shown in Figure 9, where the projection of the array elements 1–2 on the x -axis are array elements 1 2 . The distance between them is defined as the actual array spacing d ˜ 1 , 2 .
The space error δ d r of the array elements caused by the rotation error of the rotary device is
δ d r = d 5 , 6 d 1 , 2
In Figure 9
δ d r 2 = d o , 5 d o , 1 = d o , 5 d o , 1 × cos δ θ x z = d 2 d 2 × cos δ θ x z = d 2 ( 1 cos δ θ x z )
So,
δ d r = d ( 1 cos δ θ x z )
δ x = | x ˜ x | = | ( h φ 1 , 2 φ 5 , 6 ) ( 1 d ˜ 5 , 6 d ˜ 1 , 2 ) | = | ( h φ 1 , 2 φ 5 , 6 ) ( 1 d + δ d r d ) | = h φ 1 , 2 φ 5 , 6 δ d r d
In Equation (47), δ d r / d is the rotation angle error factor, where it can be seen that the rotation angle error of the USBL rotating device will produce the spacing error of the array elements, and the array element spacing error is one of the important parameters affecting the positioning accuracy of the USBL. The influence of the rotation angle error on the positioning accuracy of the USBL will be discussed by the numerical simulation method.
It can be seen from Table 5 that a 1° rotation angle error can produce a positioning error of about 0.15‰ in the horizontal distance in the rotating and array element reuse model. That is, a rotation error of 1° can cause about a 1.5 m horizontal positioning error in the horizontal positioning distance of 10,000 m. Usually, the rotation error will be less than 0.1°, and the horizontal rotation position error of 10,000 m at this level will cause a horizontal positioning error of 0.015 m, which is a very small error value. This error can be neglected even in high-precision, full-depth positioning operations. It can be seen that the rotation angle error of the USBL rotating device was less than 0.1°, which has little effect on the horizontal positioning accuracy in the USBL positioning model based on the rotary and array reusing elements method.

4. Simulation and Verification

In order to analyze the influence of the element spacing error δ d on the positioning accuracy of the USBL, considering the wavelength error δ λ , ranging error δ R , and element spacing error δ d of the underwater acoustical signal, this paper used the numerical simulations to set the elements spacing error of different orders of magnitude and then evaluated their impact on the USBL positioning accuracy. In addition, in order to verify the effectiveness of the proposed positioning calculation model in this paper, it was also compared with the existing USBL positioning calculation model.
Since all USBL positioning calculation methods have a slightly lower accuracy in the vertical direction than other methods, this paper only examined the positioning accuracy of the USBL positioning calculation model in the horizontal direction. The USBL solution model proposed in this paper was mainly used to eliminate the spacing error of the array elements δ d and the wavelength error δ λ of the underwater acoustic signal. In the simulation, the influence of the installation error angle on the positioning accuracy was not considered, that is, the installation error was assumed to be fully compensated and that there was no error in the phase difference measurement. The sailing schematic diagram of the mother ship based on the traditional model of the slant range and azimuth is shown in Figure 10.
The target was mounted with a transponder that received the underwater acoustic signal from the acoustic transmitter of the acoustic array on the mother ship, and responded to the acoustic array with an acoustic signal. In this process, the distance between the target and the USBL acoustic array was measured, and the phase difference of the acoustic array was measured when the USBL acoustic array received the underwater acoustic signal of the transponder.
Figure 11 is the sailing schematic diagram of the mother ship based on the rotating array and reusing elements method.
The underwater targets can adopt the beacon mode. The USBL acoustic array realizes the measurement of the phase difference by passively receiving the underwater acoustic signal from the beacon, thereby achieving an accurate measurement of the underwater target orientation. At the same time, the depth value h T of the underwater target is measured by a high-precision depth gauge and is sent to the USBL acoustic array through the underwater acoustic signal, and the depth value h A of the acoustic coordinate origin of the USBL acoustic array is recorded at the moment of the phase difference measurement (also measured by a high-precision depth meter). The distance difference h = | h T h A | between the origin of the U-frame and the underwater beacon is considered as the z-axis coordinate value of the target in the U-frame.
The simulation parameters were set as Table 6,
In the circumnavigation sailing process of the mother ship, the circumnavigation radius was designed as 100 m, 500 m, and 1000 m, respectively, and the USBL positioning accuracy under different positioning distances was investigated. Figure 12 shows the navigational trajectory diagram of the mother ship with the circumnavigation radius of 100 m.

4.1. Simulation Experiment 1

In order to study the influence of the spacing error of the acoustic array elements of different orders of magnitude on the USBL positioning accuracy of the traditional USBL positioning model based on the slant range and azimuth method, in this simulation, the USBL acoustic array was designed with a 4-elements acoustic array structure, and the positioning calculation was based on the slant range and azimuth method. Assuming that the spacing errors of the two elements were equal, that is δ d = δ d 1 , 2 = δ d 3 , 4 , and the error values were set to 3 mm, 2 mm, 1 mm, and 0.1 mm, respectively, we investigated the influence of the element spacing error of different orders of magnitude on the USBL horizontal positioning accuracy, the result show in Figure 13.
It can be seen from Figure 13 that in the case of δ d , at a horizontal distance of 1000, when δ d = 3 mm, the average positioning error of USBL in the x-axis direction was 13.62 m, and when it was 2 mm, the average positioning error of the USBL in the x-axis direction was 9.46 m. When δ d = 1 mm, the average positioning error of the USBL in the x-axis direction was 5.30 m. When δ d = 0.1 mm, the average positioning error of the USBL in the x-axis direction was 1.56 m. The mean square error was 0.43 m, 0.30 m, 0.17 m, and 0.049 m, respectively. It can be seen that the element spacing error is one of the main error sources of the USBL positioning system. The larger the element spacing error, the larger the positioning error. A 1 mm element spacing error can cause about a 6 m horizontal positioning error at the positioning distance of 1000 m. It can be seen that the positioning error is a steady state error from the theoretical analysis and simulation.

4.2. Simulation Experiment 2

In order to verify the effectiveness of the USBL positioning model based on the rotating array and reusing elements method proposed in this paper, we analyzed the mother ship’s horizontal positioning accuracy under different positioning distances. For this reason, the traditional positioning error based on the slant range and azimuth method and Sun’s method in the same environment were simulated and compared.
Figure 14 shows the positioning error of the target in the x-axis of the U-frame in the case where the circumnavigation radius of the mother ship was 100 m, 500 m, and 1000 m, respectively, with the element spacing error of δ d = 2 mm and it was not calibrated. It can be seen from Figure 14 that when the circumnavigation radius of the mother ship was 100 m, the positioning error of the USBL in the U-frame (x-direction) was about 1 m. When the circumnavigation radius of the mother ship was 500 m, the positioning error of the USBL in the USBL frame (x-direction) was about 5.2 m. When the circumnavigation radius of the mother ship was 1000 m, the positioning error of the USBL in the USBL frame (x-direction) was about 10.4 m. This positioning error is related to the distance of the target transponder, and the positioning error was about 1% of the slant distance. This is a relatively large error in a long distance and high-precision positioning environment.
In [23], the geometrical element spacing of the USBL could be calibrated within 0.12 m. Figure 15 is the result of the simulation of the USBL horizontal positioning according to Sun’s method. It can be seen from Figure 15 that this method could greatly improve the positioning accuracy, especially in the case of a horizontal distance of 1000 m, where the horizontal positioning accuracy (x-axis direction) could reach 2.6 m.
Figure 16 is the horizontal positioning error diagram of the USBL positioning model based on the rotating array and reusing elements method proposed in this paper. It can be seen from Figure 16 that the positioning error was 0.12 m when the horizontal bypass distance was 100 m. When the horizontal bypass distance was 500 m, the positioning error was 0.57 m, and when the horizontal bypass distance was 1000 m, the positioning error was 1.14 m.

4.3. Simulation Experiment 3

In order to study the influence of the rotation angle error on the USBL positioning accuracy based on the rotating and array elements reuse method proposed in this paper, this experiment simulated the USBL rotation angle error of different magnitudes. The USBL acoustic array was designed with a 2-element structure, as shown in Figure 5, using the positioning calculation method based on depth and phase ratio information. The depth of the transponder h was 100 m, the depth measurement error δ h was 0.001*h+0.01 m. The mother ship carrying the USBL acoustic array slowly sailed around the underwater transponder with the radius of the bypass r set as 1000 m. The rotation angle error values of the USBL rotating device were set to 5°, 3°, 2°, 1°, 0.1°, and 0.01°, respectively. We investigated the effect of the USBL rotation angle error δ θ at different magnitudes on the USBL horizontal positioning accuracy. Figure 17 presents the horizontal positioning error that considers the influence of different rotation angle errors at the horizontal distance of 1000 m.
It can be seen from Figure 17 that when the USBL rotation angle error was 5°, the x -direction positioning error of the USBL based on the rotary and element reuse model was about 5.3 m. The horizontal positioning errors were 2.6 m, 1.8 m, and 1.3 m when the USBL rotation angle error was set to 3°, 2°, and 1°, respectively. When the rotation angle error was 0.1° and 0.01°, the USBL horizontal positioning error was less than 1.2 m.

5. Conclusions

The theoretical analysis and simulation experiments showed that the element spacing error was the main error source of the USBL positioning system. Under the condition of a positioning range of 1000 m, the array element spacing of 250 mm, and regardless of other error terms, a 1 mm element spacing error can cause a horizontal positioning error of about 4 m in the USBL positioning model based on the slant range and azimuth method. Additionally, there was a relationship between the error and the distance of the positioning distance where the larger the distance, the larger the positioning error. The expansion of the spacing of the array elements can suppress its influence on the USBL positioning error to a certain extent, but it cannot substantially eliminate the influence of the spacing error of the array elements. Moreover, the expansion of the spacing of the array elements will cause phase difference ambiguity, and the USBL receiving array will become huge, which will destroy the portability of the USBL.
A new USBL positioning model based on the rotating array and reusing element method was proposed in this paper. Only two receiving array elements are required in this method. The simulation showed that for the USBL positioning model based on the method of the rotating array and reusing elements method in the range of 100 m, 500 m, and 1000 m in the USBL positioning system with a spacing of 250 mm, regardless of other errors, the positioning error can theoretically be reduced to 0.12 m, 0.57 m, and 1.14 m, respectively. In the case of considering the rotation angle error of the USBL rotating device and the depth measurement error, in the horizontal positioning environment of 1000 m, when the rotation angle error was less than 0.1°, the underwater positioning error was less than 1.2 m. The rotation angle error of the USBL rotating device was less than 0.1°, which is technically easy to implement. The model proposed in this paper can completely eliminate the influence of the wavelength error of the underwater acoustic signal and the error of the spacing of array elements on the positioning accuracy of the USBL. The positioning method has the performance of high precision in the long distance, which has important engineering application value, and provides a new idea for the engineering design of the USBL underwater positioning system.

Author Contributions

Conceptualization, J.T.; methodology, J.T.; software, J.T.; validation, L.Z.; formal analysis, J.W.; investigation, L.H.; writing—original draft preparation, J.T.; writing—review and editing, J.T., L.H., Y.L.; supervision, X.X.; project administration, X.X.; funding acquisition, X.X.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under grants 51775110, 61473085, and 51375088; the Open Fund of Key Laboratory of Inertial Technology under grant 614250607011709; the Fundamental Research Funds for the Central Universities under grants 2242018K40065, 2242018K40066, 2242019K40035, and 2242019K40040; the Foundation of Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology (Ministry of Education, China) under grant SEU-MIAN-201803; and the Foundation of Shanghai Key Laboratory of Navigation and Location Based Services.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jin, B.N.; Xu, X.S. USBL Technology and Its Applications in Ocean Engineering. Navig. Position. Timing 2018, 5, 8–20. [Google Scholar]
  2. Tong, J.W.; Xu, X.S.; Zhang, T.; Zhang, L.; Li, Y. Study on Installation Error Analysis and Calibration of Acoustic Transceiver Array Based on SINS/USBL Integrated System. IEEE Access 2018, 6, 66923–66939. [Google Scholar]
  3. Sun, D.J.; Zheng, C.E.; Zhang, J.C.; Han, Y.F.; Cui, H.Y. Development and Prospect for Underwater Acoustic Positioning and Navigation Technology. J. Chin. Acad. Sci. 2019, 3, 331–338. [Google Scholar]
  4. Lee, P.M.; Jun, B.H.; Kim, K.; Lee, J.; Aoki, T.; Hyakudome, T. Simulation of an Inertial Acoustic Navigation System With Range Aiding for an Autonomous Underwater Vehicle. IEEE J. Oceanic Eng. 2007, 32, 2. [Google Scholar] [CrossRef]
  5. Tian, T. Underwater Positioning and Navigation Technology; Liu, P., Ed.; National Defense Industry Press: Beijing, China, 2007; pp. 1–150. [Google Scholar]
  6. Hu, H.Q. Research on AUV navigation technology based on USBL aided SINS. Master’s Thesis, Southeast University, Nanjing China, 2017. [Google Scholar]
  7. Ixblue. Gaps. Available online: https://www.ixblue.com/products/gaps-pre-calibrated-insusbl-system (accessed on 6 October 2019).
  8. Morgado, M.; Oliveira, P.; Silvestre, C.; Vasconcelos, J.F. USBL/INS Tightly-Coupled Integration Technique for Underwater Vehicles. In Proceedings of the 9th International Conference on Information Fusion, Florence, Italy, 10–13 July 2006; pp. 1–8. [Google Scholar]
  9. Morgado, M.; Oliveira, P.; Silvestre, C.; Vasconcelos, J.F. Embedded Vehicle Dynamics Aiding for USBL/INS Underwater Navigation System. IEEE Trans. Control Syst. Technol. 2014, 22, 322–330. [Google Scholar] [CrossRef]
  10. Morgado, M.; Oliveira, P.; Silvestre, C. Tightly coupled ultrashort baseline and inertial navigation system for underwater vehicles: An experimental validation. J. Field Robotics 2013, 30, 142–170. [Google Scholar] [CrossRef]
  11. Zhang, Y.W. An algorithm used in underwater SINS/USBL integrated navigation. Navig. Position. Timing 2016, 7–13. [Google Scholar]
  12. Gao, B.B. Research on SINS/USBL/DVL integrated navigation system and high performance algorithm for underwater vehicles. Master’s Thesis, Northwestern Polytechnical University, Xi’an, China, 2015. [Google Scholar]
  13. Dong, P.; Cheng, J.H.; Liu, L.Q.; Mou, H.J. INS/USBL nonlinear integrated navigation method based on observations of relative information. J. Syst. Eng. Electron. 2019, 41, 179–185. [Google Scholar]
  14. Hu, G.G.; Gao, S.S. A derivative UKF for tightly coupled INS/GPS integrated navigation. ISA Trans. 2015, 56, 135–144. [Google Scholar] [CrossRef] [PubMed]
  15. Chang, L.B.; Hu, B.Q. Robust Initial Attitude Alignment for SINS/DVL. IEEE/ASME Trans. Mechatron. 2018, 23, 2016–2021. [Google Scholar] [CrossRef]
  16. Chen, H.H. The Estimation of Angular Misalignments for Ultra Short Baseline Navigation Systems. Part I: Numerical Simulations. J. Navig. 2013, 66, 561–578. [Google Scholar] [CrossRef]
  17. Chen, H.H. In-situ alignment calibration of attitude and ultra short baseline sensors for precision underwater positioning. Ocean Eng. 2008, 35, 1448–1462. [Google Scholar] [CrossRef]
  18. Tong, J.W.; Xu, X.S.; Zhang, T.; Li, Y.; Yao, Y.Q.; Weng, C.C.; Hou, L.H.; Zhang, L. A misalignment angle error calibration method of underwater acoustic array in strapdown inertial navigation system/ultrashort baseline integrated navigation system based on single transponder mode. Rev. Sci. Instrum. 2019, 90, 085001. [Google Scholar] [CrossRef] [PubMed]
  19. Yu, M. Research on Long Range Ultra Short Baseline System. Ph.D. Thesis, Harbin University of Engineering, Harbin, China, 2005; pp. 1–143. [Google Scholar]
  20. Wang, Y.H.; Sun, G.Q.; Lai, X.Q.; Liu, R.B. Research of ultrashort baseline positioning method based on depth information. In Proceedings of the IEEE 2nd Information Technology, Networking, Electronic and Automation Control Conference, Chengdu, China, 15–17 December 2017; pp. 18–22. [Google Scholar]
  21. Chen, H.H. The Estimation of Angular Misalignments for Ultra Short Baseline Navigation Systems. Part II: Experimental Results. J. Navig. 2013, 66, 773–787. [Google Scholar] [CrossRef] [Green Version]
  22. Yu, M. In-situ calibration of transceiver alignment for a high-precision USBL system, In Proceedings of the International Conference on Computer Application and System Modeling, Taiyuan, China. 22–24 October 2010. pp. V11-84–V11-87.
  23. Sun, D.J.; Ding, J.; Zheng, C.E. Array geometry calibration for underwater compact arrays. Appl. Acoust. 2019, 145, 374–384. [Google Scholar] [CrossRef]
Figure 1. The schematic diagram of the Ultra-Short Base Line array structure based on four receiving array elements.
Figure 1. The schematic diagram of the Ultra-Short Base Line array structure based on four receiving array elements.
Sensors 19 04373 g001
Figure 2. The diagram of the receiving and sending array with the angle misalignment error.
Figure 2. The diagram of the receiving and sending array with the angle misalignment error.
Sensors 19 04373 g002
Figure 3. Diagram of the USBL phase difference.
Figure 3. Diagram of the USBL phase difference.
Sensors 19 04373 g003
Figure 4. The schematic diagram of the USBL array structure based on six receiving array elements.
Figure 4. The schematic diagram of the USBL array structure based on six receiving array elements.
Sensors 19 04373 g004
Figure 5. The schematic diagram of the USBL array structure based on the ‘rotating array and reusing elements’.
Figure 5. The schematic diagram of the USBL array structure based on the ‘rotating array and reusing elements’.
Sensors 19 04373 g005
Figure 6. The rotating schematic diagram of the USBL array based on the rotating array and reusing elements. (a) The schematic diagram of the U -frame rotating 90° clockwise around the y-axis; (b) the schematic diagram of the U -frame rotating 90° clockwise around the x -axis; and (c) the schematic diagram of the U -frame rotating 90° clockwise around the z-axis.
Figure 6. The rotating schematic diagram of the USBL array based on the rotating array and reusing elements. (a) The schematic diagram of the U -frame rotating 90° clockwise around the y-axis; (b) the schematic diagram of the U -frame rotating 90° clockwise around the x -axis; and (c) the schematic diagram of the U -frame rotating 90° clockwise around the z-axis.
Sensors 19 04373 g006aSensors 19 04373 g006b
Figure 7. The flow chart of the USBL positioning calculation model based on the rotating array and reusing elements.
Figure 7. The flow chart of the USBL positioning calculation model based on the rotating array and reusing elements.
Sensors 19 04373 g007
Figure 8. The rotation diagram of the USBL acoustic array.
Figure 8. The rotation diagram of the USBL acoustic array.
Sensors 19 04373 g008
Figure 9. Schematic diagram of the USBL rotation error angle.
Figure 9. Schematic diagram of the USBL rotation error angle.
Sensors 19 04373 g009
Figure 10. The sailing schematic diagram of the mother ship in the traditional model based on slant range and azimuth method.
Figure 10. The sailing schematic diagram of the mother ship in the traditional model based on slant range and azimuth method.
Sensors 19 04373 g010
Figure 11. The sailing schematic diagram of the mother ship based on the rotating array and reusing elements method.
Figure 11. The sailing schematic diagram of the mother ship based on the rotating array and reusing elements method.
Sensors 19 04373 g011
Figure 12. Navigational trajectory diagram of the mother ship with the circumnavigation radius of 100 m.
Figure 12. Navigational trajectory diagram of the mother ship with the circumnavigation radius of 100 m.
Sensors 19 04373 g012
Figure 13. Schematic diagram of the influence of array element spacing error on positioning accuracy (based on the slant range and azimuth method).
Figure 13. Schematic diagram of the influence of array element spacing error on positioning accuracy (based on the slant range and azimuth method).
Sensors 19 04373 g013
Figure 14. Horizontal positioning error at different positioning distances (based on the slant range and azimuth method).
Figure 14. Horizontal positioning error at different positioning distances (based on the slant range and azimuth method).
Sensors 19 04373 g014
Figure 15. Horizontal positioning error at different positioning distances (by Sun’s method).
Figure 15. Horizontal positioning error at different positioning distances (by Sun’s method).
Sensors 19 04373 g015
Figure 16. Horizontal positioning error at different positioning distances (based on the rotating array and reusing elements method).
Figure 16. Horizontal positioning error at different positioning distances (based on the rotating array and reusing elements method).
Sensors 19 04373 g016
Figure 17. The influence of rotation angle error on USBL horizontal positioning accuracy.
Figure 17. The influence of rotation angle error on USBL horizontal positioning accuracy.
Sensors 19 04373 g017
Table 1. Error parameters.
Table 1. Error parameters.
Error Term δ T   ( ms ) δ C   ( m / s ) δ d   ( mm )
Error value10.052
Table 2. Location environment and equipment parameters.
Table 2. Location environment and equipment parameters.
C (m/s) R (m) d (mm)
15001000250
Table 3. Horizontal positioning error of USBL caused by the spacing error of the array elements. ( d = 250   m m , x-axis direction).
Table 3. Horizontal positioning error of USBL caused by the spacing error of the array elements. ( d = 250   m m , x-axis direction).
d = 250   m m
x (m) 10501002005001000200030004000
δ x (m) δ d = 2   m m 0.07940.39680.79371.58733.96837.936515.873023.809531.7460
δ x (m) δ d = 1   m m 0.03980.19920.39840.79681.99203.98407.968111.952215.9363
δ x (m) δ d = 0.12   m m 0.00480.02400.04800.09600.23990.47980.95951.43931.9191
Table 4. Horizontal positioning error of the USBL caused by the spacing error of the array elements. ( d = 500   m m , x-axis direction).
Table 4. Horizontal positioning error of the USBL caused by the spacing error of the array elements. ( d = 500   m m , x-axis direction).
d = 500   m m
x (m) 10501002005001000200030004000
δ x (m) δ d = 2   m m 0.03980.19920.39840.79681.99203.98417.968111.952215.9363
δ x (m) δ d = 1   m m 0.01990.09980.19960.39920.99801.99603.99205.98807.9840
δ x (m) δ d = 0.12   m m 0.00240.01200.02390.04790.11990.23990.47990.71980.9598
Table 5. Analysis of the influence of the USBL rotation angle error on the positioning error.
Table 5. Analysis of the influence of the USBL rotation angle error on the positioning error.
NO. δ θ x z   ( ° ) δ θ x z   ( rad ) δ d r   ( mm ) δ d r / d   ( )
150.0872660.9513253.805302
230.0523600.3426161.370465
320.0349070.1522930.609173
410.0174533.807621 × 10−20.152305
50.10.0017453.807717 × 10−41.523087 × 10−3
60.010.0001753.807718 × 10−51.523087 × 10−5
Table 6. Simulation parameters.
Table 6. Simulation parameters.
ParameterValue
Depth of transponder (beacon) (m)100 ± 2
Circumnavigate sailing radius of the mother ship (m)100, 500, 1000
Water acoustic velocity (m/s)1500
Water acoustic velocity error (m/s)±0.3%v + 0.005
Time delay estimation error (s)0.001
Depth measurement error (m)±0.1%h + 0.01
Distance of array elements d (mm)250
Spacing error of array elements δ d 1 , 2 , δ d 3 , 4 (mm)±3, ±2, ±0.1
Hydro acoustic signal frequency (kHz)25
Hydro acoustic signal wavelength (m)0.06

Share and Cite

MDPI and ACS Style

Tong, J.; Xu, X.; Hou, L.; Li, Y.; Wang, J.; Zhang, L. An Ultra-Short Baseline Positioning Model Based on Rotating Array & Reusing Elements and Its Error Analysis. Sensors 2019, 19, 4373. https://doi.org/10.3390/s19204373

AMA Style

Tong J, Xu X, Hou L, Li Y, Wang J, Zhang L. An Ultra-Short Baseline Positioning Model Based on Rotating Array & Reusing Elements and Its Error Analysis. Sensors. 2019; 19(20):4373. https://doi.org/10.3390/s19204373

Chicago/Turabian Style

Tong, Jinwu, Xiaosu Xu, Lanhua Hou, Yao Li, Jian Wang, and Liang Zhang. 2019. "An Ultra-Short Baseline Positioning Model Based on Rotating Array & Reusing Elements and Its Error Analysis" Sensors 19, no. 20: 4373. https://doi.org/10.3390/s19204373

APA Style

Tong, J., Xu, X., Hou, L., Li, Y., Wang, J., & Zhang, L. (2019). An Ultra-Short Baseline Positioning Model Based on Rotating Array & Reusing Elements and Its Error Analysis. Sensors, 19(20), 4373. https://doi.org/10.3390/s19204373

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop