CC BY
11
V ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ
UDC 519.87 Articles
doi:10.31799/1684-8853-2018-4-15-23
APPLICATION OF DIFFERENTIAL MODE FOR AUV LOCATION
L. A. Martynovaa'b, Dr. Sc., Tech., Associate Professor, martynowa999@bk.ru G. G. Bezruka, PhD, Tech., geoint99@mail.ru A. A. Myslivyia, PhD, Tech., aam-07@mail.ru
aPrincipal Research and Development Test Centre of Robotic Technology Ministry of Defence of the Russian Federation (PRDCR), 5, Seryogin St., 125167, Moskow, Russian Federation
bConcern CSRI Elektropribor, JSC State Research Center of Russia, 30, Malaya Posadskaya St., 197046, Saint-Petersburg, Russian Federation
Introduction: Signals from satellite radio navigation systems cannot penetrate seawater. Therefore, autonomous underwater vehicles are located with the use of one or more support objects. Because of the layered nature of a marine environment, a sound beam goes not in a straight line as radio signals do, but along an arcuate trajectory. Therefore, the correction of the vehicle location by the differential method has a large error. Purpose: Assessing the requirements and conditions necessary for using the differential mode of correcting a vehicle location. Results: We have simulated the propagation of sound beams in a vertical plane according to Snell's law, under the condition that sound beam propagation follows the Rice distribution. We have determined the distance covered by a sound beam, evaluated its dependence on the emission direction, specified the procedure of locating a vehicle with the use of the differential mode for correcting the location, and formulated the conditions for using the differential method which are determined by the sound emission direction, by the sound speed profile and by the mutual position of the emitting buoys, the correcting base station and the underwater vehicle. Practical relevance: The obtained results can be used for locating vehicles in areas where you cannot use the traditional methods of precise location.
Keywords — Autonomous Underwater Vehicle, Route Trajectory, Current, Efficiency Evaluation.
Citation: Martynova L. A., Bezruk G. G., Myslivyi A. A. Application of Differential Mode for AUV Location. Informatsionno-upravliaiushchie sistemy [Information and Control Systems], 2018, no. 4, pp. 15-23. doi:10.31799/1684-8853-2018-4-15-23
Introduction
Increasing the accuracy of locating autonomous underwater vehicles (AUV), especially in the case of using them in groups, for example, when conducting seismic survey or searching for anthropogenic objects [1-10], is of paramount importance. By now, a large number of methods for positioning an underwater object have been developed [1116], including AUV positioning, a comprehensive review of which is given in [17]. To locate an AUV under water, range-measurement hydroacoustic navigation systems [18] or systems with a single mobile beacon [19] are used. Some recent attempts have been made to conduct underwater positioning in accordance with the principles of satellite radio navigation systems [20-22].
At the same time, the peculiarity of locating an AUV in a group is that:
— the location should be conducted aboard the AUV;
— there can be no emitting devices or hy-droacoustic modems aboard the AUV;
— there can be no inertial navigation system or auxiliary equipment like lags or echosounders aboard the AUV due to the limitations on weight-size parameters and power consumption.
A way to increase the accuracy of locating an object in satellite radio navigation is the use of differ-
ential corrections to pseudoranges or pseudospeeds of radio signal propagation [23]. The transfer of the differential method of object location correction into a marine environment requires taking into account the laws of sound signal propagation in water and the dependence of the propagation on the mutual position of the AUV and reference objects (other underwater or surface vessels with known coordinates, bottom beacons, surface buoys, escort vessels, etc.), as well as the dependencies on the state and parameters of the marine environment determined by the location, weather conditions, time of year and time of day.
Marine environment has the following features:
— significant attenuation of sound in water;
— sound propagation along multiple directions;
— refraction of beams causing a curvature of the trajectory along which a sound beam propagates;
— layered nature of the medium and the dependence of the sound speed on depth.
The above-mentioned studies on navigational location of AUV do not discuss the ways to apply the differential method to correcting underwater object location.
Differential Mode Description
The idea of the differential mode of correcting an object location is based on developing corrections
to pseudoranges from the source of a hydroacoustic signal to its consumer. As signal sources, we will consider transmitter buoys (further referenced to as buoys). The consumer is an AUV. The prefix "pseudo-" is added because the time scales of the radiation source and the AUV are not synchronized with one another and hence the measured range value differs from the actual one. To maintain the differential mode of pseudorange correction, in a specified seabed area a base station (BS) is established in advance, whose coordinates are being determined over a long period of time by means of precise measurement of its position [24].
We will assume that the buoys on the sea surface are located using the signals of satellite radio navigation systems. The same signals are used to synchronize the time scales of the buoys with high accuracy.
The formation of differential corrections of the BS with their transfer to consumers, including AUV, goes as follows. The buoys periodically emit signals into the seawater. The period of the signals is maintained with high accuracy. The frequency of the signals is known to the BS and the AUV. Aboard the AUV, using the signal received from the buoy, the signal arrival delay is determined, counting from the moment of its emission. According to the delay, distance to the buoy is calculated. In this manner, distances to at least four buoys are found. The distances to the buoys determine the AUV position by the difference-ranging method. This method is used because of the mismatch between the time scales of the buoys and the AUV: the beginning of the AUV time scale does not coincide with the beginnings of the time scales of the buoys mutually coordinated by the signals of satellite radio navigation systems.
With this measurement, the evaluation of the pseudoranges between the AUV and the buoys has a certain error similar to the errors in satellite radio navigation systems caused by the features of both signal propagation environment and signal delay measurement hardware. AUV onboard delay calculation is so sophisticated because marine environment is very different from aerial: the signals of satellite radio navigation systems propagate in space rectilinearly and almost without distortion, while acoustic signals propagating through water are influenced by a host of unfavorable factors. The world ocean is an inhomogeneous medium; there are warm and cold currents, whirlpools, areas with different concentrations of dissolved salts, etc. Acoustic signals propagating with a low speed relative to radio signals will physically move, shift in phase, be reflected and distorted while passing through the inhomogeneities of a marine environment. Partially, their influence can be compensated by choosing appropriate mathematical processing methods, as all the above-mentioned factors are
constant [25, 26]. However, not everything can be compensated only mathematically; there are errors in the measurements as such, which also contribute to the error of pseudorange calculation.
The base station, similarly to AUV, also estimates the distance to the buoys by the delays of the received signals. Simultaneously, the BS evaluates the calculated distance to the buoy, knowing its exact and constant location and the position of the buoy at the signal emission moment. It is assumed that the emitted signal contains the coordinates of the buoy and the emission time. The BS compares the obtained distance values with the calculated ones; on the base of this comparison it produces a correction, and emits it into the water for the consumers, including the AUV. Receiving these corrections and adding them up, the AUV updates its estimation of the distance to the buoy. The updated distances to four buoys are then used in the difference-ranging method.
Mathematically, it can be described as follows. Let it be given that:
xBS, yBS, zBS are BS coordinates known with the required accuracy;
xi, yi, zi are the coordinates of the ith buoy with a high accuracy; i = 1, ..., M, where M is the number of the buoys, M = 4.
At some instant of time t, all the M buoys emit sonar signals carrying the information about the time t and each ith buoy coordinates xi, yi, zi.
Let the time in which a signal from the ith buoy reaches the BS be xBg i and the time in which it reaches the AUV be xAUV i.
Then the actual distance between the ith buoy and the BS is defined by the following expression:
-^BS i = TBS
where C is the speed of sound in water in this particular area.
From the ith buoy coordinates xi, yi, zi and the base station coordinates xBS, yBS, zBS precisely known at the base station, the distance between the buoy and the base station is calculated:
^calcBsi =\j(xi - *BS )2 + (Vi - VBS )2 + (zi - ZBS )2 •
It determines the differential correction to the distance to each ith buoy as the difference between the actual and calculated values:
ARi — -RactBS; _ ^oaioBS; •
The distance from the ith buoy to the AUV is defined by the following expression:
-^AUV i = TAUV
Using the correction ARi produced at the BS and transferred to the AUV, the distance between the buoy and the AUV is calculated:
R.
calc
AUV1
= R,
act
AUV1
~AR:.
Having these values for the distances from the AUV to four or more buoys, the AUV locates itself by the difference-ranging method through solving a system of equations
tJ(x2 - X )2 + (y2 - y )2 + (z2 - zf -
-■\(xi - *)2 + (í/i - y)2 + - zf
- *)2 +(ys - y)2 +(z3 - 2)2
-\¡(xi - x)2 + (z/i - y)2 + {zi - zf
^j(x4 - x )2 +(y4 - y )2 +(z4 - zf - X )2 + (?/! - y )2 + (z-y - zf
= An
21
= An
31
= Ar<
41
where Ari1 is the difference between the distances from a couple of reference points {Oi, O1} to the point in question; x, y, z are the AUV position point coordinates; xi, yi, zi are the coordinates of the ith buoy; i = 1, ..., 4.
Fig 1 explains the difference-ranging method; the symbol pbi designates the distance from the buoy reference points to the AUV position point.
The calculation of ARi corrections and their transfer to the consumers occur with a certain constant periodicity, which is why this is called a differential mode of pseudorange correction. The differential mode is successfully used in satellite navigation, improving object location accuracy.
Applying the differential mode in a marine environment requires that we consider how the marine environment affects the propagation of sound and the determination of distances between objects.
Let us assume that:
— the beam propagation follows the Rice distribution [27] with the dominance of one beam;
— a BS knows the sound speed profile (SSP) which is periodically measured with a pulsed SSP meter [28];
— the AUV also knows the SSP which can either coincide with the BS SSP or differ from it;
— signals formed by a buoy correlate with one another.
Estimating the Influence of Seawater Stratification on the Sound Beam Trajectory
The curvature of a sound beam propagation trajectory caused by seawater stratification is schematically shown in Fig. 2.
The distance to which a beam deviates from the vertical because it was emitted at an angle to the horizontal is determined by the expression [29]
ARj = pj (sin a j - sin a ),
(1)
where j is a layer number j = 1, ..., n; n is the number of the layers; pj is the radius of the circle whose arc is the beam trajectory in the jth layer; aj is the angle of the beam entry into the jth layer; a^ is the angle of the beam entry into the (j - 1)th layer; when j = 1 the direction of the signal emitted by the buoy is a—1.
The radius of the circle is determined by the expression
Pi =-
a is a relative gradient:
a =-
Gc_ c0
■ Fig. 1. Diagram explaining the difference-ranging method
Fig. 2. Schematic representation of beam curvature
where GC is a gradient of the speed of sound in water determined by the expression
Gc =-
cj+1 ~ cj zi+1 - V
Let us assume that GC = const. Then the full horizontal distance of the beam deviation from the vertical will be determined by the following expression:
Ri = ■
7=1
(2)
From triangle A00102 we can express the distances between the point of the beam entry into the /h layer and the point where the beam leaves the jth layer:
di =
ej+i - xi) + (zj+i - zi) ■ (3)
Besides, from triangle A00102 we can derive an expression for finding the angle Pj at a vertex as a half-angle between the directions to the point of the beam entry into the jth layer and the point where the beam leaves it (see Fig. 2). This expression looks like
di • ß; — = P; Sin — ,
2 1 2
from which it follows that
ßy = 2a
2py
(4)
Then, knowing the vertex angle Pj and radius p, we can determine the length of the arc l, along which the beam propagated:
lj = PyPy ■
(5)
The angle Pj is expressed in radians. The above expressions (1)-(5) allow us to estimate the distance covered by the beam in the 7th layer.
Then, depending on the AUV immersion depth, the beam path trajectory length is estimated as a sum of the beam trajectories through the layers:
n
£=1 h ■ 7=1
If the AUV does not know the speed of sound, the following averaging formulas can be used: — weighted average speed of sound
Z(c7 - C7h - Zi-1)
izl_
2(zauv _ zi)
where zi is the ith buoy antenna immersion depth; in the problem under discussion zi = 0;
— harmonic mean of the speed of sound
Ch =
_ (ZAUV " zi)
2 £ (Ci - C7-l)' 7=1 {Zi ~ Z1 "I )
To estimate the change in the calculated range when the curvilinear trajectory of the sound propagation is taken into account, we estimated the deviations of an arcuate beam trajectory from the rectilinear path connecting these points. As initial data, a real SSP was considered (Fig. 3, a and b). In the calculations, the width (thickness) of a step layer was equal to 50 m. The size of the chosen sampling interval was determined by the data error [30]. The beam pattern representation is shown in Fig. 4.
The results of the calculations (Fig. 5) showed that the deviation of an arcuate trajectory from the rectilinear one is about 9 m. This value was obtained for near-surface emission of an acoustic beam at an angle of 10° to the horizontal. When the emission direction changed to 80°, the deviation decreased from 9 down to 1.5 m.
a) Sr3_1 67.053N 79.923E 0 50 100 150 200 t 250 Q 300 350 400 450
28.07.2014 04:40:00
532 m
500
b)
1435 1440 1445 1450 1455 1460 1465 Sound speed, m/s
0
100
200
300
400
500
1435 1440 1445 1450 1455 1460 1465
■ Fig. 3. The sound speed profile for the Kara Sea: a — data of the source [30]; b — digital representation for use in modeling
Fig. 4. Beam pattern representation
10
50 100 150 200 250 300 350 400 450 500 Depht, m
-0-10 -»-20 -A-30 —40 —50 —60 —70 --80
■ Fig. 5. Deviation of an arcuate beam propagation trajectory from the straight line connecting the points of the buoy position on the sea surface and the position of the receiver (AUV)
a)
0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0
* I
'! \ 'i-
V
V
\
---f-y-1
10 20 30 40 50 60 70 80 Angle of slope, deg
b)
700 600 500 400 300 200 100 0
t ~ ~ \
i A \ \
\\
\
i r--1-
10 20 30 40 50 60
Angle of slope, deg
--c ■ cw
70
80
■ Fig. 6. An estimate of signal travel time (a) and covered distance (b) along the beam, at a known sound speed C (dashed line) and a weighted average Cw (solid line)
Besides, we estimated the duration of the beam travel and the distance it covers for the case when the speed of sound is considered weighted average (Fig. 6, a and b).
The obtained results allow us to state that when the emission angle relative to the water surface is 20° or more, we can use the weighted average value of the sound speed, rather than counting it literally by layers, because the results for the time and distance at a> 20° nearly coincide.
Determining a Distance Covered
by a Beam Along an Arcuate Trajectory
by Layers
As discussed above, a beam propagation trajectory can be sophisticated. Therefore, when you calculate the length covered by a beam using the measured delay between the emission of a signal and its reception in a certain point, you have to take into account the complex curvilinearity of its trajectory. To use the formulas (1)-(5), you need to know the direction in which the buoy emits the beam, i. e. the beam angle a in the top layer. Since the trajectory is curvilinear, the emission direction is not obvious, requiring a special approach to its determination.
Mathematically, the problem of determining the emission direction a can be formulated as follows:
Let the beam arrive to an AUV point whose coordinates are unknown and need to be found. We know the time necessary for the signal to cover the curvilinear trajectory, and we know the SSP.
We need to find the coordinates of the AUV point to which the beam arrives.
This means we need to find a triplet of coordinates which would meet the following criterion:
(XAUV, ^AUV, ZAUV ) -(x, y, z) e K : AT[(*, y, z); (xAUV, yAVY, zAUV )] =
7=1
where K is the set of all points of possible AUV position in the given area, and AT [(1); (2)] is the time of sound travel between points 1 and 2.
In order to determine the emission direction, let us consider two cases. In the first case, the AUV has the information about its immersion depth; in the second case, it has not.
1. The AUV has the information about its immersion depth.
Let us assume that the AUV measures the depth by a special pressure gauge. Aboard the AUV, the delay is determined between the moments of emission and reception of the sound. To determine the
AUV position, you need, for various values of a, to go through all the variants for the time of signal travel through the layers and to calculate the total signal travel time. The value closest to the measured time will be the problem solution. It will allow you to uniquely determine the emission direction for the signal accepted by the AUV.
Mathematically, this approach is called a net method: we build a net and find a solution in its knots. A net is characterized by its cell size.
The result of modeling the beam travel time for different emission directions from 10 to 80° is shown in Fig. 7.
Knowing the depth of the AUV position and the duration of the sound signal travel from a buoy to the AUV, we estimate the distance and, according to the chart (see Fig. 7), find the beam angle at the moment of its emission.
For example, the time of signal travel along the arc from the surface to the depth of 270 m is 0.5 s. For this conditions, the following variants are possible:
— the emitted signal inclination angle is 10°, and the depth of the AUV is H = 105 m;
— the emitted signal inclination angle is 20°, and the depth of the AUV is H = 270 m.
Out of these two variants, the more suitable one is the second (a = 20°), because the actual depth coincides with the calculated one. This variant is chosen as a solution. To find the AUV coordinates, the coordinates of the signal reception point are used.
However, as shown above, at emission angles bigger than 20° we can use weighted average value of the sound speed. In this case, we estimate the distance covered by the beam immediately by the delay. Knowing the depth H (cathetus) and inclined range of the beam travel L (hypotenuse), we can find the emission direction angle
a = arccos
50 100 150 200 250 300 350 400 450 500 Depht, m
— 10 —20—30—40 — 50 —60 —70 —80
■ Fig. 7. Calculated time of signal travel along the arc for different emission angles
2. The AUV has no information about its immersion depth.
In this case, we need to simultaneously use two measurements from two buoys, because one measurement can only provide a multivalued solution.
If the depth is unknown, we will take into account the fact that everything happens in one vertical plane: the two buoys and the AUV are three points through which one and only one plane passes. Since the buoys emit sound in all directions, and the arc length does not depend on the direction in the horizontal plane, depending only on the inclination angle in the vertical plane, the solution will be the intersection point of two arcs whose length depends on the inclination angle and the AUV depth.
In the same way as described above, the emission direction is determined for the signal from the first buoy. For the buoy signal arrival delay measured aboard the AUV, several variants come out which differ from each other in the signal arrival depths and emission directions. These variants form a set of solutions for the first buoy.
After that, the same procedure is applied to the delay of the second buoy signal. The solution variants form a set of solutions for the second buoy.
Out of the two solution sets related to the first and second buoys, an element is found which is an intersection of the two. This element is the solution. It is uniquely associated with the emission direction of the first and second buoys and the AUV position coordinates which correspond to the arrival point of the signal from the first and second buoys.
Similarly to the known depth variant, we get expressions to determine the inclination angles aj and a2:
ai = arccos
a9 = arccos
Li L2
H\ - H2.
The given expressions can be used for inclination angles a1 and a2 bigger than 20°.
Finding AUV Coordinates
with the Difference-Ranging Method
The arcuate curvilinear trajectory of the beam travel from two spatially separated buoys gets into a certain point at a certain depth where the AUV is positioned which has received a sound signal from both the buoys. In the same way as if they were rectilinear inclined trajectories, after determining
the point to which the beams from both the buoys come, inclined straight line segments are calculated. The set of possible solutions for two buoys form a hyperbola. By the intersection of three hyperbolas formed by four buoys using the difference-ranging method, the AUV determines its coordinates.
Differential Method Description
The mutual position of the AUV, BS and one of the four buoys in the marine environment is shown in Fig. 8. The AUV location is adjusted using the differential mode as follows.
A buoy emits a navigational hydroacoustic signal into the seawater. The BS receives it. Knowing its precise location, SSP and the precise location of the buoy, the BS estimates the signal travel time: calculated Tcalc and actual Tact. The calculated time is determined taking into account the refraction which depends on the current SSP in this area according to Snell's law, so that the beam comes to the BS point. For that, all possible positions of the beam arrival to the BS area are considered and the best one is chosen for various emission directions a by the criterion how close the beam arrival point is to the BS location. Finally, the actual beam travel time from the buoy to the BS is found, and the distance Ract BS is determined which differs from the calculated distance Rcalc BS by AR:
Ract BS = Rcalc BS +
On this basis, the correction AR is determined on the BS as the difference between the calculated
Buoy
beam travel distance Rcalc BS and the actual distance Ract BS which is precisely known at the BS:
AR = Ract BS - Rcalc BS.
When estimating the rectilinear range of the beam travel between the emission point and the BS position, you can use not only the average value of the sound speed according to the SSP, but also any other value, provided that this value is used by the AUV. This should be agreed upon in advance when you perform a differential correction of the AUV location. Besides, it is desirable that the processing principles of the navigation receivers at the BS and the AUV are similar; ideally, the receivers should be identical. The signal travel correction AR calculated at the BS is then emitted from the buoy to the consumers in the marine environment, including the AUV. Such corrections are emitted to the signals of all the four buoys.
The autonomous underwater vehicle near the BS measures the delays of the signals from the four buoys. The obtained delays are converted into the ranges to the buoys as discussed above. Based on these ranges, the beam trajectory projections are determined. An obtained range is a sum of the calculated range and the correction:
Ract AUV = Rcalc AUV +
Simultaneously with the measurement of the actual range, the AUV receives from the BS the information about the pseudorange corrections AR for each of the four buoys, and then specifies the pseudoranges to the buoys:
Rcalc AUV = Ract AUV - AR,
using them for locating itself by the difference-ranging method.
Conclusion
The work is aimed at finding the ways of using the differential method to locate an AUV in a marine environment. Our research was focused on the parameters which considerably affect the usage of difference-ranging method for AUV location: sound emission direction, vertical distribution of the sound speed, sound travel trajectory, and the distance covered by a sound beam. The paper discusses an approach to applying the differential method of correcting an AUV location in a marine environment. The obtained results can be used for a more precise location of AUV or other underwater objects.
The research was supported by RFBR grant 17-07-00125.
1. Gizitdinova M. R., Kuz'mitski M. A. Mobile Underwater Robots in Modern Oceanography and Hydro-physics. Fundamental'naya iprikladnaya gidrofizika [Fundamental and Applied Hydrophysics], 2010, vol. 3, no. 1 (7), pp. 4-13 (In Russian).
2. Bozhenov Yu. A. Using Autonomous Underwater Vehicles for Research in Arctic and Antarctic. Fundamental'naya i prikladnaya gidrofizika [Fundamental and Applied Hydrophysics], 2011, vol. 4, no. 1, pp. 4-68 (In Russian).
3. Ilarionov G. Yu., Sidenko K. S., Bocharov L. Yu. Ugroza iz glubiny: XXI vek [Threat from the Depth: 21st Century]. Khabarovsk, Khabarovskaya krae-vaya tipografiya Publ., 2011. 304 p. (In Russian).
4. Bezruk G. G., Martynova L. A., Saenko I. B. Dynamic Method of Searching Anthropogenic Objects in the Seabed with use of Autonomous Underwater Vehicles. Trudy SPIIRAN [SPIIRAS Proceedings], 2018, no. 4, pp. 203-226 (In Russian). doi:10.15622/sp.58.9
5. Martynova L. A., Karsaev O. V. А Method of Coordinating the Behavior of Autonomous Underwater Vehicles Group on a Multi-Agent Basis in the Conduct of the Seismic Survey. Izvestiia YUFU. Tekhnicheskie nauki [Izvestiya SFedU. Engineering Sciences], 2018, no. 1(195), pp. 52-67 (In Russian).
6. Martynova L. A., Grinenkov A. V., Pronin A. O., Kulikovskikh J. V. Research of the Functioning of a Multi-Agent Control System for an Autonomous Underwater Vehicle using Simulation. Naukoemkie tekh-nologii v kosmicheskikh issledovaniyakh Zemli [High Technologies in Earth Space Research], 2017, vol. 9, no. 5, pp. 52-65 (In Russian).
7. Martynova L. A. Concerted Action of a Radiator and Autonomous Uninhabited Submersibles for Effective Seismic Exploration. Informatsionno-upravliaiush-chie sistemy [Information and Control Systems], 2017, no. 1, pp. 83-92 (In Russian). doi:10.15217/issn1684-8853.2017.1.83
8. Martynova L. A. Differential Method of Positioning a Standalone Unmanned Submersible in Seismic Exploration. Informatsionno-upravliaiushchie sistemy [Information and Control Systems], 2017, no. 4, pp. 7785 (In Russian). doi:10.15217/issn1684-8853.2017.4.77
9. Martynova L. A. Underwater Observation under Intellectual Interference. Informatsionno-upravliai-ushchie sistemy [Information and Control Systems], 2018, no. 1, pp. 31-41 (In Russian). doi:10.15217/ issn1684-8853.2018.1.31
10. Martynova L. A., Konyukhov G. V., Pashkevich I. V., Rukhlov N. N. Peculiarities of Group Control of AUV in Seismic Exploration. Izvestiia YUFU. Tekhnicheskie nauki [Izvestiya SFedU. Engineering Sciences], 2017, no. 9(194), pp. 53-63 (In Russian).
11. Kebkal K. G., Kebkal V. K., Kebkal A. G., Petroc-cia R. Experimental Estimation of Delivery Success of Navigation Data Packages Transmitted via Digital
Hydroacoustic Communication Channel. Giroskopi-ya i navigatsiya [Gyroscopy and Navigation], 2016, no. 2 (93), pp. 107-122 (In Russian). doi:10.17285/ 0869-7035.2016.24.2107-122
12. Kebkal K. G., Banasch R. Swep-Spread Carier for Underwater Communication over Acoustic Chanels with Strong Multipath Propagation. Journal of the Acoustical Society of America, 2002, vol. 12(5), pp. 2043-2052.
13. Kebkal K. G., Kebkal V. K., Kebkal A. G. Digital Hy-droacoustic Networks for Communication in Conditions of Long Delays and Breaks the Connection: an Experimental Research. Podvodnye issledovaniia i robototekhnika [Underwater Research and Robotics], 2001, no. 2(20), pp. 12-19 (In Russian).
14. Dubrovin F. S., Scherbatyuk A. F. Study of the Algorithms for the Single Beacon Mobile Navigation of Unmanned Underwater Vehicles: Results of Simulation and Sea Trials. Giroskopiya i navigatsiya [Gyros-copy and Navigation], 2015, no. 4 (91), pp. 160-172 (In Russian). doi:10.17285/0869-7035.2015.23.4.160-172
15. Kebkal K. G., Kebkal A. G., Kebkal O. G. Synchronization Tools of Acoustic Communication Devices in Control of Underwater Sensors, Distributed Antennas, and Autonomous Underwater Vehicles. Gi-roskopiya i navigatsiya [Gyroscopy and Navigation], 2014, no. 2, pp. 70-85 (In Russian).
16. Kebkal K. G., Kebkal O. G., Glushko E., Kebkal V. K., Sebastiao L., Pascoal A., Gomes J., Ribeiro J., Silva H., Ribeiro M., Indivery G. Underwater Acoustic Modems with Integrated Atomic Clocks for OneWay Travel-Time Underwater Vehicle Positioning. UACE2017 4th Underwater Acoustics Conference and Exhibition, At Skiathos, Greece, 2017, pp. 315-323.
17. Kebkal K. G., Mashoshin A. I. Hydro-Acoustic Methods of Positioning Autonomous Unmanned Underwater Vehicles. Giroskopiya i navigatsiya [Gyroscopy and Navigation], 2016, vol. 24, no. 3 (94), pp. 115-130 (In Russian).
18. Scherbatyuk A. Ph., Dubrovin F. S., Rodionov A. Yu.,
Unru P. P. Group Navigation and Control for Marine Autonomous Robotic Complex Based on Hydroacous-tic Communication. Proc. of the IROS 2016 IEEE/ RSJ Conf., October 9-14, 2016, Daejeon, Korea, 2016. Available at: https://www.semanticscholar.org/pa-per/Group-navigation-and-control-for-marine-au-tonomous-Scherbatyuk-Dubrovin/4853b452a93a-52c8288823a83b583d1ee18eeaef (accessed 5 March 2018).
19. Scherbatyuk A. Ph., Dubrovin F. S. About Accuracy Estimation of AUV Single-Beacon Mobile Navigation using ASV, Equipped with DGPS. Proc. of the OCEANS 2016 MTS/IEEE Conf., Shanghai, China, 2016. Available at: https://www.semanticscholar. org/paper/About-accuracy-estimation-of-AUV-sin-gle-beacon-ASV%2C-Dubrovin-Scherbatyuk/ 41763212d8011a4d06f6014b63a6d4ee8fb84373 (accessed 5 March 2018).
20. Podvodnyi GPS s nulia za god [Underwater GPS from Scratch for the Year]. Available at: https://geektimes. ru/post/284902/ (accessed 5 March 2018).
21. Podvodnyi GPS: prodolzhenie. Available at: https:// geektimes.ru/post/289159/ (accessed 5 March 2018).
22. POSYDON - proekt podvodnoi navigatsionnoi siste-my, kotoraia pozvolit opredeliat' mestopolozhenie v glubinakh morei i okeanov [POSYDON is a Project of an Underwater Navigation System that will allow to Determine the Location in the Depths of the Seas and Oceans]. Available at: https://dailytechinfo.org/news/ 8169-posydon-proekt-podvodnoy-naviga-cionnoy-sistemy-kotoraya-pozvolit-oprede-lyat-mestopolozhenie-v-glubinah-morey-i-okeanov.html (accessed 5 March 2018).
23. Iatsenkov B. B. Osnovy sputnikovoi navigatsii [Fundamentals of Satellite Navigation]. Moskow, Katalog Publ., 2002. 106 p. (In Russian).
24. Buia M., Flores P. E., Hill D., Palmer E., Ross R., Walker R., Houbiers M., Thompson M., Laura S., Menlikli C., Moldoveanu N., Snyder E. Shooting Seismic Surveys in Circles. Available at: https://www.slb. com/~/media/Files/resources/oilfield_review/ors08/ aut08/shooting_seismic_surveys_in_circles.pdf (accessed 5 March 2018).
25. Raskita M. A. Sound Velocity Profiles Reconstruction Error Estimate at Water Masses Remote Ecological Monitoring. Izvestiia YUFU. Tekhnicheskie nauki [Izvestiya SFedU. Engineering Sciences], 2009, no. 6(95), pp. 91-95 (In Russian).
УДК 519.87
doi:10.31799/1684-8853-2018-4-15-23 Исследование возможности применения дифференциального режима уточнения местоположения АНПА под водой
Мартынова Л. А.а б, доктор техн. наук, старший научный сотрудник, martynowa999@bk.ru Безрук Г. Г.а, канд. техн. наук, geoint99@mail.ru Мысливый А. А. а, канд. техн. наук, aam-07@mail.ru
аГлавный научно-исследовательский испытательный центр робототехники Министерства обороны Российской Федерации, Серегина ул., 5, Москва, 125167, РФ
бАО «Концерн «ЦНИИ «Электроприбор», Малая Посадская ул., 30, Санкт-Петербург, 197046, РФ
Постановка проблемы: ввиду отсутствия под водой сигналов спутниковых радионавигационных систем позиционирование автономных необитаемых подводных аппаратов происходит с использованием одного или нескольких опорных объектов. Слоистость морской среды является причиной распространения звукового луча не по прямой, как у радиосигналов, а по дугообразной траектории, из-за чего корректировка местоположения аппарата дифференциальным методом происходит с большой погрешностью. Цель: оценка требований и условий, необходимых для использования дифференциального режима уточнения местоположения аппарата. Результаты: проведено моделирование распространения звуковых лучей в вертикальной плоскости согласно закону Снеллиуса при условии подчинения распространения звуковых лучей распределению Райса, определен пройденный звуковым лучом путь, оценена его зависимость от направления излучения звукового сигнала, а также определен порядок оценки местоположения аппарата с использованием дифференциального режима уточнения его местоположения. Установлены условия использования дифференциального метода корректировки местоположения аппарата, определяемые направлением излучения; профилем распределения скорости звука; взаимным положением излучающих буев, корректирующей базовой станции и подводного аппарата. Практическая значимость: полученные результаты могут быть использованы для уточнения местоположения аппарата в районах с отсутствием возможности применять традиционные методы точного позиционирования аппарата.
Ключевые слова — автономный необитаемый подводный аппарат, маршрутная траектория, течение, оценка эффективности.
Цитирование: Martynova L. A., Bezruk G. G., Myslivyi A. A. Application of Differential Mode for AUV Location // Информационно-управляющие системы. 2018. № 4. С. 15-23. doi:10.31799/1684-8853-2018-4-15-23
Citation: Martynova L. A., Bezruk G. G., Myslivyi A. A. Application of Differential Mode for AUV Location. Informatsionno-upravliaiushchie sistemy [Information and Control Systems], 2018, no. 4, pp. 15-23. doi:10.31799/1684-8853-2018-4-15-23
26. Gurbatov S. N., Kurin V. V., Kustov L. M., Proncha-tov-Rubtsov N. V. Physical Modeling of Nonlinear Sound Wave Propagation in Oceanic Waveguides of Variable Depth. Acoustical Physics, 2005, vol. 51, no. 2, pp. 152-159.
27. Karabanov I. V., Burdinskiy I. N., Mironov A. S., Linnik M. A. Model of Spread Spectrum Phase-Shift Keyed Signal Propagation in Underwater Acoustic Channel. Izvestiia YUFU. Tekhnicheskie nauki [Izvestiya SFedU. Engineering Sciences], 2013, no. 3 (140), pp. 102-112 (In Russian).
28. Kolomiets S. M. The Determination of the Vertical Profile of the Speed of Sound in the Ocean with using the Doppler Locator. Issledovaniia v oblasti est-estvennykh nauk [Researches in Science], 2014, no. 9, pp. 16-25. Available at: http://science.snauka.ru/ en/2014/09/8305 (accessed 5 March 2018) (In Russian).
29. Zaraiskii V. A., Tiurin A. M. Teoriia hydrolokacii [Theory of Sonar]. Leningrad, Voenno-morskaya akademiya Publ., 1975. 604 p. (In Russian).
30. Balakin R. A., Golavskii V. E., Timets V. M. Applied Research of Hydroacoustic Fields in the Arctic Seas. Trudy Gosudarstvennogo okeanograficheskogo in-stituta im. N. N. Zubova [Proc. of N. N. Zubov State Oceanographic Institute], 2016, no. 217, pp. 312-322 (In Russian).
