This paper proposes a new imaging geometry model for multi-receiver synthetic aperture sonar (SAS). The model considers the change of the speed of sound in seawater, the effect of platform movement on the acoustic velocity vector (AVV), and the Doppler effect. Based on the proposed model, a solution to determine the phase distribution was generated to improve the SAS image quality. The simulation results demonstrate the merits of proposed model compared to the traditional models that consider the speed of sound in seawater as a fixed value, ignore the change of AVV during transmission, and suppress the Doppler effect.

※ The user interface design of www.jsts.org has been recently revised and updated. Please contact inter@theieie.org for any inquiries regarding paper submission.

### Journal Search

## 1. Introduction

Synthetic aperture sonar (SAS) is based on the coherent combination of received
signals in successive pings (pulse repetition interval) when a sonar platform moves
along a straight path. With this principle, SAS can generate a virtual aperture that
is many times larger than the physical size to increase the along-track (azimuth)
resolution ^{[1}, ^{2]}. Thanks to the achievement of high resolution independent of both the range and frequency,
SAS is widely used for many applications, such as searching for small objects, imaging
of wrecks, underwater archaeology, and pipeline inspection ^{[2}, ^{3]}.

The working principle of SAS is similar to that of synthetic aperture sonar (SAR).
However, the area mapping rate of SAS is lower than that of SAR because the speed
of sound in seawater is much lower than the propagation speed of electromagnetic waves
in air. To increase the area coverage rate, multi-receiver SASs configured with a
transmitter (projector) and a receiver (hydrophone) array have been used commonly
in SAS applications ^{[3}, ^{4]}.

It is challenging to construct a geometry model that completely describes the
physical processes for SAS imaging reconstruction due to the complicated configuration
and challenges of navigation, sound velocity errors, topographic errors, vehicle stability,
and shallow waters ^{[5]}. A conventional model uses phase center approximation (PCA), which considers a projector/hydrophone
pair as an element at a central point of the two elements, which is used to simplify
the calculation ^{[6]}. Due to the difference of propagation distance in these two cases, PCA error is generated
when calculating the propagation time from the transmitter to a target and then back
to the receivers, which causes degradation of the SAS image quality ^{[7]}.

To reduce the PCA error, a geometry model was made to consider the platform motion
during reception ^{[7]}. However, this model does not consider the change of the speed of sound in seawater
according to depth, it ignores the change of the acoustic velocity vector (AVV) during
transmission, and it suppresses the Doppler effect. With these limited conditions,
the physical processes affecting the received signals are not considered completely.

This paper proposes an imaging geometry model for multi-receiver SAS that takes
into account the change of the speed of sound according to depth, the change of AVV
during transmission, and the Doppler shift. It is based on more complete consideration
of the physical processes compared with conventional models. With these complete descriptions,
the proposed model can provide simulation data for researching and developing SAS
image reconstruction algorithms. Based on the back projection algorithm (BPA) ^{[8]}, the SAS image quality obtained with the proposed model can be improved compared
with that derived from conventional models.

## 2. Proposal of Imaging Geometry Model for Multi-receiver SAS

Two-dimensional (2D) models including range dimension (depth) and azimuth dimension
have been used widely in research to improve the along-track resolution for multi-receiver
SAS ^{[3}, ^{7}, ^{8]}. In the depth dimension, the speed of sound in seawater is a nonlinear function of
temperature, salinity, and depth ^{[9]}:

##### (1)

$c=1402.5+5T-5.44\times 10^{-2}T^{2}+2.1\times 10^{-4}T^{3}+1.33S\\ -1.23\times 10^{-2}ST+8.7\times 10^{-5}ST^{2}+1.56\times 10^{-2}Z\\ +2.55\times 10^{-7}Z^{2}-7.3\times 10^{-12}Z^{3}+1.2\times 10^{-6}Z\left(Φ -45\right)\\ -9.5\times 10^{-13}TZ^{3}+3\times 10^{-7}T^{2}Z+1.43\times 10^{-5}SZ$where $\textit{T}$ is the temperature in $^{\circ}$C, $\textit{Z}$ is the depth in m, $\textit{S}$ is the salinity in ‰, and ${\Phi}$ is the geographic latitude in $^{\circ}$.

The dependence of the speed of sound in seawater on depth is also described by
the sound velocity profile (SVP) ^{[10]} obtained by using sound velocity profilers. The variation of the speed of sound in
seawater can be up to 2% over a typical depth profile ^{[1}, ^{5]}, which can cause defocus and degrade SAS image quality. To deal with these variations,
phase gradient autofocus (PGA) or an average value of the speed of sound can be used
for reconstructing SAS images ^{[5}, ^{11]}.

PGA is carried out based on acquisition data at receivers and phase adjustment according to neighbor points to achieve a sharp image. Therefore, it is difficult to construct mathematical models describing received signals and phase distribution using PGA. The proposed 2D imaging geometry model for multi-receiver SAS is based on averaging the speed of sound along the depth combined with the consideration of the Doppler effect and the change of AVV during transmission, as shown in Fig. 1.

In Fig. 1, $\textit{c}$$_{av}$ is the average speed of sound in seawater, and $\textit{d}$
is the distance between two adjacent elements in an array with $\textit{N}$ uniformly
spaced receivers. The baseline distance between the transmitter and $\textit{i}$th
receiver is $\textit{d}$$_{i}$. The SAS is moved along axis $\textit{x}$ with velocity
$\textit{v}$ and is located at O(0, 0) at $\textit{T}$ = 0. Axis $\textit{r}$ expresses
the depth $\textit{Z}$ in (1) (it is also the range dimension ^{[3}, ^{7}, ^{8}, ^{13]}).

The platform velocity vector can change due to rotational movements (pitch, roll,
and yaw) and translational movements (heave, sway, and surge). These motions are accurately
determined by high-grade aided inertial navigation systems (INSs), including an inertial
measurement unit, Doppler velocity logger, and error state Kalman filter ^{[5]}. To obtain images with high quality, subwavelength navigation accuracy is required
for reconstructing SAS images ^{[14]}. The requirement is satisfied by using the sonar itself as a navigation sensor combined
with INS ^{[5}, ^{14]}. Based on the navigation, the platform can be positioned, which is considered as
a uniform velocity motion (speed $\textit{v)}$ to simplify SAS imaging reconstruction
^{[3}, ^{7}, ^{8}, ^{15}, ^{17]}.

In Fig. 1, the total velocity vector $\overrightarrow{v_{\Sigma T}}$ is the sum of vectors
$\overrightarrow{v}$ and$\overrightarrow{c_{a v}}$ . The angle between $\overrightarrow{v_{\Sigma
T}}$ and $\overrightarrow{v}$ is $\alpha_{1}$ ^{[13]}:

The module of the total velocity vector is expressed as:

The propagation time of the signal from the transmitter to point P($\textit{x, r}$) is given by:

##### (4)

$\tau _{1}=\frac{\sqrt{\left(x-vt\right)^{2}+r^{2}}}{v_{\Sigma T}}=\frac{\sqrt{\left(x-vt\right)^{2}+r^{2}}}{v\cos \alpha _{1}+\sqrt{{c_{av}}^{2}-v^{2}\sin ^{2}\alpha _{1}}}$After transmitting to point P, the signal is obtained at the $\textit{i}$th receiver according to the direction from point P to receiver R$_{\mathrm{i}}$, which is determined by angle $\alpha _{2i}$:

##### (5)

$\alpha _{2i}=\arccos \left(\frac{x-\left[d_{i}+v\left(\tau _{1}+\tau _{2i}\right)+vt\right]}{\sqrt{\left[x-\left(d_{i}+v\left(\tau _{1}+\tau _{2i}\right)+vt\right)\right]^{2}+r^{2}}}\right)$
The propagation time of the signal from P to R$_{\mathrm{i}}$ is $\tau _{2i}$
^{[13]}:

##### (6)

$\tau _{2i}=\frac{v\left(\sqrt{\left(x-vt\right)^{2}+r^{2}}-\left(d_{i}+v\tau _{1}\right)\cos \alpha _{1}\right)+\sqrt{\Delta _{i}}}{{c_{av}}^{2}-v^{2}}$
where $\Delta _{i}$ is expressed as ^{[13]}

##### (7)

Types of signals used in modern SAS include gated continuous wave (CW), linear
frequency modulation (LFM), and hyperbolic frequency modulation (HFM) pulses ^{[11]}. The CW pulse is widely used in classical sonar ^{[1]}, whereas the LFM pulse is commonly used in modern SAS ^{[3}, ^{5}, ^{7}, ^{8}, ^{11}, ^{15]}. Therefore, this section presents signal model for SAS using the CW pulse and LFM
pulse. When using the CW pulse, the transmitted signal with carrier frequency $\textit{f}$$_{c}$
is expressed by:

where $w\left(\tau \right)$ is the windowing function of the transmitted signal
(a Gaussian function or rectangle function). A rectangle function with pulse duration
$\textit{T}$$_{p}$ is defined as ^{[15]}:

The received signal at point P is determined as:

##### (10)

$s_{P}\left(\tau ,t\right)=\omega _{TP}\left(t\right)\mathrm{w}\left[\eta _{1}\left(\tau -\tau _{1}\right)\right]\\ \times \exp \left\{j2\pi f_{c}\left[\eta _{1}\left(\tau -\tau _{1}\right)\right]\right\}$
where $\omega _{TP}\left(t\right)$ is the beam pattern of the transmitter. $\eta
_{1}$ is the time-stretching factor at point P due to the Doppler effect ^{[16}, ^{17]}:

The received signal at the $\textit{i}$th receiver due to the scattering from P is determined as:

##### (12)

$s_{i}\left(\tau ,t\right)=\omega _{TP}\left(t\right)\omega _{Ri}\left(t\right)\mathrm{w}\left[\eta _{1}\eta _{2i}\left(\tau -\tau _{2i}\right)-\eta _{1}\tau _{1}\right]\\ \times \exp \left\{j2\pi f_{c}\eta _{1}\eta _{2i}\left(\tau -\tau _{2i}\right)-\eta _{1}\tau _{1}\right\}$
where $\omega _{Ri}\left(t\right)$ is the beam pattern of the $\textit{i}$th receiver.
$\eta _{2i}$ is the time-stretching factor at the $\textit{i}$th receiver due to the
Doppler effect, which is calculated as ^{[16}, ^{17]}:

With the combination of received signals in the overlap of the main beams of the transmitter and each receiver, $\omega _{TP}\left(t\right)$ and $\omega _{Ri}\left(t\right)$ can be suppressed to focus on the coherent processing of the phase. After ignoring the scattering from the sea surface, the received signal at the $\textit{i}$th receiver is determined as:

##### (14)

$s_{i}\left(\tau ,t\right)=\mathrm{w}\left[\eta _{1}\eta _{2i}\left(\tau -\tau _{2i}\right)-\eta _{1}\tau _{1}\right]\\ \times \exp \left\{j2\pi f_{c}\eta _{1}\eta _{2i}\left(\tau -\tau _{2i}\right)-\eta _{1}\tau _{1}\right\}$When LFM pulses are used, the transmitted signal is expressed as:

##### (15)

$s_{e\_ LFM}\left(\tau ,t\right)=w\left(\tau \right)\exp \left(j2\pi f_{c}\tau +j\pi \gamma \tau ^{2}\right)$where ${\gamma}$ is the chirp rate [Hz/s]. Similarly to (14), the received signal at the $\textit{i}$th receiver in SAS using the LFM pulses is given by:

##### (16)

Expressions (14) and (16) represent echo signals when considering the change of AVV during transmission, the Doppler effect, and the change of the speed of sound in seawater according to depth. Based on these equations, received data can be generated for researching and developing SAS imaging reconstruction algorithms.

## 3. Reconstructing SAS Image in Multi-receiver SAS using BPA

SAS image reconstruction estimates the reflectivity function from backscattered
echoes when an SAS platform moves along a known track ^{[1}, ^{2]}. Due to its high resolution, BPA is implemented for SAS applications requiring image
quality gain with arbitrary platform histories and a nonlinear SVP ^{[4]}. With the BPA, the image reconstruction in the range dimension and in the azimuth
can be employed by a matched filter and beamforming, which is also called the delay-and-sum
method ^{[1}, ^{2]}. To evaluate the effectiveness of the proposed model, this section discusses SAS
using the CW pulse for simplicity.

After match filtering, the pulse width at the output considered at half the level of the peak is equal to the pulse width of the CW pulse at the input. With wide CW pulses, the beamforming by the delay-and-sum method can be converted to the phase compensation after removing the carrier frequency. Based on the proposed model, the phase distribution is calculated as below when the main beam is steered to the target at point P$_{0}$($\textit{x}$$_{0}$\textit{, r}$_{0}$)

##### (17)

$\psi _{i}\left(t,x_{0},r_{0}\right)=2\pi f_{c}\left(\frac{c_{av}}{c_{av}-v\cos \alpha _{01}}\tau _{01}+\frac{c_{av}+v\cos \alpha _{02i}}{c_{av}-v\cos \alpha _{01}}\tau _{02i}\right)$where ${\alpha}$$_{01}$ and ${\alpha}$$_{02}$ are determined similarly to ${\alpha}$$_{1}$ and ${\alpha}$$_{2}$ in (4) and (5), but $\textit{x}$ and $\textit{r}$ are replaced by $\textit{x}$$_{0}$ and $\textit{r}$$_{0}$, respectively. After compensating the phase, the backscattered signals from $\textit{N}$ receivers in $\textit{M}$ pings are coherently combined as:

##### (18)

$ff\left(x,r,x_{0},r_{0}\right)=\sum _{m=1}^{M}\sum _{i=1}^{N}ss_{eni}\left(\tau ,t\right)\exp \left(j\psi _{i}\left(t,x_{0},r_{0}\right)\right)$where $ss_{eni}\left(\tau ,t\right)$ is the envelope signal obtained from (14) after demodulation.

The beam pattern in Eq. (18) is a function of the variables $\textit{x}$ and $\textit{r}$, which are known as
the azimuth slice when the range variable $\textit{r}$ is fixed ^{[3}, ^{7}, ^{8}, ^{13]}. To reduce the side-lobe level (SLL) in the beam pattern, the elements in the receiver
array must be arranged densely without gaps between two adjacent elements ^{[1}, ^{4]}. With this condition, the main beam of each element is given by ^{[18]}

For convenient calculation, this beam width can be approximately calculated as
${\lambda}$/d. The length of the synthetic aperture at range $\textit{r}$ was determined
by ^{[18]}:

where $\eta $ is a programmable parameter controlling the process beamwidth ($\eta
\leq 1$) ^{[1]}. The number of pings for coherent processing according to the synthetic aperture
length $L_{sa}$ and the length of receiver array $\textit{L}$ is determined as:

Table 1. The parameters of multi-receiver SAS.

## 4. Simulation Results

To emphasize the effectiveness of the proposed model, we considered an example
of multi-receiver SAS with the parameters shown in Table 1. These parameters were chosen to decrease SLL (the along-track sampling criterion)
^{[1]}.

SVP can be determined using accurate sound velocity profilers (e.g., $\textit{miniSVP}$
(www.valeport.co.uk), which has with an accuracy of 0.02 m/s). With this accuracy,
SAS imaging reconstruction is implemented without correcting for sound velocity errors,
in which case there are no accumulated errors from navigation inaccuracies ^{[11]}. When navigation errors are not compensated completely, the average sound velocity
(ASV) for processing is chosen around the measured value until obtaining sharp images
^{[11]}.

To focus on reconstructing SAS images, we used SSP from measured data ^{[5}, ^{11]}, as shown in Fig. 2, and we assumed that there are no navigation errors. Based on this SSP, the ASV between
SAS and the seafloor is equal to 1492 m/s when the SAS is located at a depth of 4
m and the target is located on the seafloor with a depth of 17.3 m ^{[11]}. This ASV was used to generate SAS data. When the difference between the depth of
the target and the depth of the SAS is calculated as 13.3 m, the number of pings $\textit{M}$
is 12 with $\eta =0.6047$.

To highlight the merits of the proposed model, we considered an ideal target located
at (13.3 m; 2 m) in coordinate system O$\textit{rx}$. Azimuth slices were generated
with MATLAB software, as shown in Fig. 3. There were three cases of beamforming based on the proposed model, the conventional
model ^{[8]} ($\textit{c}$ = 1500 m/s, ignoring the change AVV and the Doppler effect), and PCA
($\textit{c}$ = 1500~m/s). In Fig. 3, the azimuth slices obtained from the proposed solution, the conventional solution
^{[8]}, and the PCA are denoted by a solid red curve, dashed blue curve, and dotted brown
curve, respectively.

Fig. 3 shows that the proposed solution has no azimuth deviation. In cases of beamforming
based on the conventional model ^{[8]} and PCA, the azimuth deviations are 0.105 m and 0.235 m, respectively. The deviations
can cause defocus and degrade the image quality with distributed targets.

SLL in the azimuth slice derived from the proposed solution was also lower than
that derived from the two conventional models. Compared with PCA, the proposed solution
generates a significantly lower SLL (11.383 dB). The improvement of SLL by the proposed
solution in comparison with the conventional solution ^{[8]} is 1.67 dB. The improvements of SLL can increase the contrast of the SAS image and
reduce reflections from the sea surface due to side lobes.

Fig. 4 shows the relative normalized azimuth slices obtained when the three azimuth slices
in Fig. 3 are normalized according to the same maximum value (relative normalization). The
proposed solution increased the peak level of the main beam (or signal to noise ratio
(SNR)) by 0.35 dB and 5.582 dB compared to the conventional solution and PCA, respectively.
The improvement of SNR can increase the SAS imaging quality ^{[19]}. With the consideration of severe conditions, the computation time of the phase distribution
using the proposed solution is larger than that using the conventional solution ^{[8]} and PCA. Table 2 shows the computation times based on MATLAB 2015A and a laptop with an i5-7200U Intel
processor @ 2.5 GHz and 8 GB of RAM for the three cases with the target at (13.3 m,
2 m).

According to Table 2, the proposed solution uses more than 5 times and 4 times the computation time in
comparison with the solution based on PCA and the conventional solution^{[8]}, respectively. To reduce the computation time, designers can use more powerful computing
devices than ours. In addition to accurate sound velocity profilers, the proposed
model required bathymetry equipment to map the scene prior to reconstructing SAS images,
as in conventional SAS processing ^{[1}, ^{5}, ^{11]}. This requirement can be satisfied by accurate bathymetry equipment such as the MIDAS
Surveyor (www.valeport.co.uk) with accuracy being greater than ${\pm}$0.01 m or ${\pm}$0.02\%.

For detection in the presence of multiple reflections (multipath), vertical SAS arrays are used for minimizing the seafloor-scattered echoes. With the configuration of the SAS array along the horizontal dimension used in this paper, multipath rejection was not implemented so as to concentrate on processing according to the azimuth dimension. This issue will be considered in more detail in the future.

## 5. Conclusion

This paper proposed a new imaging geometry model for multi-receiver SAS considering the variation of the speed of sound in seawater, the change of AVV during transmission, and the Doppler effect. From the proposed model, the SAS data could be generated to research and develop SAS image reconstruction algorithms. This paper also proposed a beamforming solution for multi-receiver SAS considering physical processes more completely than conventional models. With the proposed solution, the SAS imaging quality was improved by determining the target position more accurately, reducing SLL, and increasing SNR. The simulation results demonstrated the effectiveness of the proposed model and the merits of the proposed solution.

### REFERENCES

## Author

Nguyen Dinh Tinh received his B.E. degree in Electronics -Telecommuni-cations and an M.E. degree in Radar Navigation Engineering from Le Quy Don Technical University, Vietnam, in 2008 and 2012, respectively. He has been working as a lecturer at Le Quy Don Technical University since 2009. His research interests include antennas, signal processing, synthetic aperture sonar, and sonar engineering.

Trinh Dang Khanh received his B.E. degree in Electronics -Telecommuni-cations and Ph.D. degree in Electronic Engineering from Le Quy Don Technical University, Vietnam, in 1983 and 2000, respectively. He has been working as an Associate Professor since 2011 at Le Quy Don Technical University. His research interests include microwave circuit design, signal processing, inertial navigation system, and sonar engineering.