Connectivity on Underwater MI-Assisted Acoustic Cooperative MIMO Networks.

In traditional underwater wireless sensor networks (UWSNs), it is difficult to establish reliable communication links as the acoustic wave experiences severe multipath effect, channel fading, and ambient noise. Recently, with the assistance of magnetic induction (MI) technique, cooperative multi-input-multi-output (MIMO) is utilized in UWSNs to enable the reliable long range underwater communication. Compared with the acoustic-based UWSNs, the UWSNs adopting MI-assisted acoustic cooperative MIMO are referred to as heterogeneous UWSNs, which are able to significantly improve the effective cover space and network throughput. Due to the complex channel characteristics and the heterogeneous architecture, the connectivity of underwater MI-assisted acoustic cooperative MIMO networks is much more complicated than that of acoustic-based UWSNs. In this paper, a mathematical model is proposed to analyze the connectivity of the networks, which considers the effects of channel characteristics, system parameters, and synchronization errors. The lower and upper bounds of the connectivity probability are also derived, which provide guidelines for the design and deployment of underwater MI-assisted acoustic cooperative MIMO networks. Monte Carlo simulations were performed, and the results validate the accuracy of the proposed model.

1. Introduction

Underwater wireless sensor networks (UWSNs) have been applied in various oceanography missions such as underwater pollution detection, offshore oil extraction, and surveillance where human operation is impossible [1,2,3,4,5]. Since the main function of UWSNs is to process extracted data and transmit them to remote locations, network connectivity is one of the most fundamental problems in UWSNs [6]. Due to unstable underwater environment, limited availability, and non-rechargeability of sensor nodes, the design of UWSNs would reduce cost, be robust to node failures, reduce communication overhead, and guarantee a high level of network connectivity [7,8,9].

Connectivity is affected by communication range of the sensor node, which ensures that every sensor node is connected to the surface BS directly or through a multi-hop path [10]. Compared with electromagnetic wave (EM), acoustic wave is able to propagate much longer distances in the ocean, which is often used as the media for information exchange in UWSNs. However, the acoustic wave during propagation suffers from multipath effect, channel fading, and ambient noise; as a result, it is a great challenge to establish reliable communication links among different nodes of UWSNs [11,12,13]. Furthermore, the acoustic spectrum becomes more and more scarce with the increasing usage of underwater acoustic systems [14]. The underwater acoustic multi-input-multi-output (MIMO) system may be used to satisfy the long-range and high-throughput requirements [15,16]. However, many wireless devices cannot support multiple antennas due to size, cost, and/or hardware limitations in practice. Hence, cooperative MIMO is adopted, the basic idea of which is to group multiple devices into virtual antenna arrays to emulate MIMO communications. Multiple sensor nodes cooperatively transmit the same signals to a receiver and the receiving power gain is enhanced. Consequently, it can sufficiently increase the range of wireless links without amplifying the transmission power, while saving the energy of each transmitting agent [17]. The virtual antenna array composed of multiple devices in cooperative MIMO relies on wireless communication among the sensor nodes. It requires that all the participating nodes share the same transmitted signals and carrier frequency, and they have accurate synchronous clocks.

However, underwater acoustic channels are not proper for the synchronization as the low propagation speed takes long synchronization time and has large synchronization errors. As a popular near-field communication technique, magnetic induction (MI) has great potential in UWSNs because it can provide a constant channel, high data rate, and negligible propagation delay with the underwater propagation speed of m/s [18,19,20]. Such high speed signal propagation is able to significantly facilitate the synchronization of distributed nodes, which enables the long distance underwater communication and high throughput. Therefore, heterogeneous wireless underwater sensor networks adopting MI-assisted acoustic cooperative MIMO are proposed [21]. In such networks, the distributed nodes are divided into multiple clusters by region, as shown in Figure 1. Each cluster consists of one master sensor node and multiple slave sensor nodes. The information exchange between different clusters are realized by the cooperative MIMO, while MI communication among intra-cluster nodes guarantee the accurate synchronization. The system architecture and process diagram of underwater MI-assisted acoustic cooperative MIMO networks are introduced in Section 2, and much more details can be also found in [21].

Figure 1

The system architecture of underwater MI-assisted acoustic cooperative MIMO networks.

The performance of the underwater MI-assisted acoustic cooperative MIMO networks is obtained in [21], in which the performance gain of underwater MI-assisted acoustic cooperative MIMO networks is validated in terms of synchronization errors, signal-to-noise ratio (SNR), effective communication time, and the upper bound of the throughput in physical layer. In this paper, we focus on the system level performance of underwater MI-assisted acoustic cooperative MIMO networks in terms connectivity and analyze the impact of clock synchronization on connectivity. Because of the complex channel characteristics and the heterogeneous network architecture, the connectivity analysis of MI-assisted acoustic cooperative MIMO networks is much more complicated than that of UWSNs. In particular, the acoustic cooperative MIMO and MI communications are performed in inter-cluster and intra-cluster manners, respectively. Therefore, the underwater acoustic channel and MI channel should both be considered, which increases the complexity of connectivity analysis of the networks. Besides the complex channel characteristics, the heterogeneous architecture of MI-assisted acoustic cooperative MIMO networks also poses a great challenge for connectivity analysis. Compared with the traditional UWSNs [7,22], the inter-cluster and intra-cluster connectivity need to be calculated. Meanwhile, the synchronization error is inevitable in the connectivity analysis of the networks.

Without loss of generality, the networks are supposed to be deployed in the bounded three-dimensional space . In this underwater MI-assisted acoustic cooperative MIMO networks, the clusters are distributed according to a homogeneous Poisson point process and the nodes in a cluster are distributed according to another homogeneous Poisson point process. When considering the probability that a certain cluster can be directly connected with surface BS, it is not only related to its position in space , but also related to the maximum transmission range of the cluster, which depends on the number of nodes in the cluster [23]. When considering the multi-hop fashion, some work has been well explored to analyze the connectivity of wireless sensor networks assuming two-dimensional coordinates [24]. We extend it to the three-dimensional scene in this network model and give the upper and lower bounds of connectivity.

The major contributions of this paper can be summarized as follows:

A mathematical model is developed to analyze the connectivity of underwater MI-assisted acoustic cooperative MIMO networks, which considers the three-dimensional underwater bounded environment. Connectivity probabilities of direct and multi-hop fashions are derived.

Due to the characteristics of heterogeneous networks, inter-cluster and intra-cluster connectivities need to be calculated, and the influence of synchronization errors on connectivity is also considered. Monte Carlo simulations were carried out, and the results agree well with the analytical model that validate the accuracy of the model.

Based on the framework, we present numerical and simulation results that quantitatively analyze the effects of various system parameters on the connectivity in this networks, considering the effects of channel characteristics, operating frequency, operating range, and density of nodes. The results provide the guideline for the deployment of underwater MI-assisted acoustic cooperative MIMO networks according to the practical requirement.

The remainder of this paper is organized as follows. In Section 2, a brief introduction to the network architecture of underwater MI-assisted acoustic cooperative MIMO networks is given. Based on the system architecture, the connectivity problem is formulated in Section 3. In Section 4, the transmission ranges of node and cluster are derived. Next, connectivity probabilities of direct and multi-hop fashions are further analyzed in Section 5. In Section 6, the connectivity performance evaluation by Monte Carlo simulations is presented. Finally, the conclusions of this paper are drawn in Section 7.

2. System Model

The complex system architecture constitutes one of the major challenges in the connectivity analysis. In this section, we describe the system architecture and channel model, and the synchronization accuracy is analyzed.

2.1. System Architecture

The system architecture of underwater MI-assisted acoustic cooperative MIMO networks is illustrated in Figure 1. It consists a surface base station (BS) and a great number of underwater sensor nodes. The surface BS aggregates the data from underwater nodes for environment sensing, which includes water velocity, temperature, salinity, etc. Underwater sensor nodes are managed in the form of clusters with assigned IDs. In each cluster, the sensor selected as the master node collects the monitoring data from slave nodes and coordinatea the intra-cluster synchronization. The acoustic cooperative MIMO is utilized between the surface BS and clusters as well as different clusters for information exchange and command interaction. As the cooperative MIMO technique collects all the intra-cluster nodes power, the long range underwater acoustic communication is enabled, which expands the effective cover space. Compared with traditional UWSNs [1], the cluster of the networks can be considered as an enhanced node due to the adoption of MI-assisted acoustic cooperative MIMO.

The process diagram of MI-assisted acoustic cooperative MIMO networks is shown in Figure 2. It includes five phases: Data request, Aggregation, Broadcast, Synchronization, and Communication. During the Data request period, the BS sends the data request to clusters. As the cluster receives the request, the monitoring data are transmitted from slave nodes to the master node by MI communication during the Aggregation period. The master node aggregates the monitoring data from different slave nodes and then broadcasts them to all the slave nodes during the Broadcast period. As a result, all the intra-cluster slave nodes share the same monitoring data. Timing-sync Protocol for Sensor Networks (TPSN) is adopted to adjust the clocks of slave nodes to that of the master node, which is operated in the Synchronization period [25]. After the Synchronization period, all the nodes within a cluster have the same monitoring data as well as the same clock. Finally, the monitoring information converged to the surface BS via the direct link or the multi-hop fashion in the Communication period.

Figure 2

The process diagram of underwater MI-assisted acoustic cooperative MIMO networks.

2.2. Channel Characteristics

There are two types of channel in the underwater MI-assisted acoustic cooperative MIMO networks: the underwater acoustic channel and the magnetic induction channel.

In underwater MI-assisted acoustic cooperative MIMO networks, underwater acoustic communication is adopted for communication between clusters in terms of cooperative MIMO. In this paper, we consider Thorp’s underwater acoustic channel model [26], which is based on empirical data. Thorp’s model provides the attenuation for the underwater acoustic channel, which is influenced by the frequency and the communication distance between the sender node and the receiver node. The transmission power model is formulated as [26] where is received power, is the transmitted power, and is the channel attenuation depending on the distance d and the frequency f (kHz) where is the normalizing constant and is the spreading factor. The spreading factor has a value between 1 and 2 depending on the depth. is the absorption coefficient and is expressed in decibels per kilometer

In underwater MI-assisted acoustic cooperative MIMO networks, MI is adopted for information aggregation and synchronization within a cluster. As a near-filed communication technique, the transmission power model is formulated as [27] where is received power, is the transmitted power, and r is the distance of the master node and one slave node; where and are quality factors of transmitter and receiver, respectively; and are the inductances of the transmitting and receiving coils; and are the transmitting and receiving coil resistances; and are efficiencies of the transmitter and receiver coils; and are source and load resistances; and are the radius of the transmitting and the receiving coil; is resonant frequency; is a parameter related to the channel state; and is conductivity. For sea water, fresh water, and dry oil, is S/m, S/m, and S/m, respectively.

Highly conductive sea water induces significant eddy current that incurs very high path loss. Therefore, in seawater, the path loss of MI communication also needs to take into account the influence of medium loss [28]. The attenuation coefficient is where is vacuum permeability.

Therefore, the path loss caused by the influence of seawater medium (dB) is

Finally, the path loss of MI communication in seawater is

2.3. Synchronization Accuracy Analysis

For the TPSN synchronization algorithm [25], the clock drift of two nodes is assumed to be unchanged in this small span of time. As shown in Figure 2, k − 1 slave nodes adjust their clocks to that of the master node in turn. The node z is taken as an example to illustrate the TPSN synchronization process. The node z sends synchronization-pulse packet to the master node at , and this package contains time . The master node receives the packet at , where is equal to . Here, and represent the clock drift between the two nodes and propagation delay, respectively. The master node returns an acknowledgement packet to the node z at , which contains the values of , , and , as well as beacon signal generated by the master node. The node z receives the packet at . The clock drift and propagation delay of the node z can be calculated as

As the drift is indicated in Equations (8) and (9), the node z is able to correct its clock aligned with that of the master node.

However, the node’s clock is driven by crystal oscillation. Doe to the existence of impure crystals, environmental changes, and component aging, frequency and synchronization errors are unavoidable. First, the frequency error of the synchronization is analyzed herein [21]. The average relative clock drift in of the node z to the master node can be expressed as where is the oscillator frequency of the node z and is the oscillator frequency of the master node.

For node z, its working frequency is proportional to its crystal oscillator frequency, thus where K is the amplification factor of the frequency multiplier.

After receiving the acknowledgement packet at the node z, the frequency is estimated according to the local oscillator of the node where is the working frequency of the master node and is the error of the frequency estimation.

The frequency offset at the node z for the frequency synchronization can be determined as

According to Equation (13), the optimal frequency offset is . Due to the estimated error , the frequency cannot be perfectly synchronized and the error can be defined as

For the node z, the duration of the time slot is , which mainly includes the sending time and transmission time of synchronization-pulse packet and acknowledgement packet where and are the size of synchronization-pulse packet and acknowledgement packet, B is the bandwidth, is the distance between master node and node z, and v is the propagation speed of signal.

Due to the frequency error of crystal oscillator at each node, the total drifting time of each node will gradually increase along with time. For the node z, the time synchronization error to the beginning of beamforming can be calculated by

3. Problem Formulation

Based on the system architecture and operation mechanism of underwater MI-assisted acoustic cooperative MIMO networks given in Section 2, the problem of connectivity analysis is formulated in this section.

Without loss of generality, the network is supposed to be deployed in a bounded 3D space . In this area, the clusters are regarded as independent units, which are distributed according to a homogeneous Poisson point process with the density of . The nodes in one certain cluster are subject to another homogeneous Poisson distribution with the density of .

In MI-assisted acoustic cooperative MIMO networks, a cluster with k nodes and the surface BS can successfully establish a connected link if and only if the following two conditions are satisfied: (1) slave nodes can connect topologically with the master node; and (2) the cluster with k nodes participating in cooperative MIMO and the surface BS can connect topologically via a direct link or a multi-hop fashion. Specifically, the definitions of node connected and cluster connected are provided as follow [24]:

When there is a master node within the maximum transmission range of the slave node, it is considered that it can connect to the master node, and the slave node is defined as connected.

When the surface BS is within the maximum transmission range of the cluster, it is considered that it can connect to the surface BS directly, and the cluster is defined as connected.

When there are relays that can connect topologically with each other until the surface BS, it is considered that it can connect to the surface BS in multi-hop fashion, and the cluster is defined as connected.

From the above definitions, the connectivity of the jth cluster in the networks can be defined as where is probability that the jth cluster including k nodes connects to the surface BS, denoted as Equation (18), and is the probability that there are k nodes participating in the cooperative MIMO in the jth cluster, indicated as Equation (19). Specifically, where is the position of the jth cluster, and is the probability that the jth cluster at is connected to the surface BS, which depends on the distribution and maximum transmission range of the cluster. As the maximum transmission range of the cluster depends on the number of nodes participating in the cooperative MIMO, the maximum transmission range of the cluster with k nodes is recorded as .

According to Definition 1, we have where and is the maximum transmission range of the node. Substituting Equations (18) and (19) into Equation (17), we have

From Equation (20), to derive the connectivity probability , the following terms are requested: the maximum transmission range of the node , the maximum transmission range of the cluster , and the probability that the jth cluster at connects to the surface BS . These are analyzed in the following sections.

4. Transmission Range

4.1. Inter-Cluster Acoustic Transmission Range

It is assumed that there are k nodes distributed in the jth cluster. Let denote the time delay of the ith sensor, where . The signal from the k nodes arrives at the surface BS can be written as where is the amplitude of the ith transmitted signal, is the underwater acoustic channel envelope [26], f is the signal frequency, and is the noise. Similar to the assumption in [21,29], as k nodes share the same transmitted signal, the amplitude can be recorded as H. Considering the surface BS deployed in the far-field, all the channel envelops are the same.

Therefore, the signal at the surface BS can be rewritten as

We consider the perfect CSI information, the channel envelope and the delay are known. From Equation (22), the received signals by surface BS cannot be enhanced as the phase difference caused by the time delay. To maintain the reliable communication, the clusters should be close to the surface BS that leads to a limited effective cover space. To tackle this problem, pre-compensation phases are added to the nodes to eliminate the time delays, and this technique is referred as cooperative MIMO. After the phase compensation, the received signal at the surface BS can be expressed as where is the pre-compensation phase and is the wave length. Define , , and . The received signal can be expressed in a compact form [29]:

Then, the SNR of received signal at the surface BS can be calculated as

The optimal compensation vector can be obtained by solving the following optimization problem:

The optimization problem in Equation (26) can be solved by the Lagrangian multiplier method, and the optimal compensation vector is given as below [29]

The above analysis does not consider the synchronization error. Due to the existence of the impure crystals, component aging, and time drifting, frequency and time estimation errors are unavoidable. Taking account of the frequency and time estimation errors, the signals received at the surface BS can be expressed as where , ∘ is Hadamard product. is frequency estimation error of the ith node and is the time estimation error of the ith node, which are explained in Section 2.

Therefore, the SNR of signals received at the surface BS is given as [21]

Based on Equation (29) the maximum transmission range of a cluster is determined as where indicates the minimum SNR that the surface BS is able to receive and recover the transmitted signal from the cluster.

4.2. Intra-Cluster MI Transmission Range

In underwater MI-assisted acoustic cooperative MIMO networks, within a cluster MI is adopted for information aggregation and synchronization. To obtain the intra-cluster transmission range, the SNR of MI is needed. From the transmission power model in Equation (4) in Section 2, the SNR of the signal received by the master node can be expressed as

Furthermore, the maximum transmission range of intra-cluster communication can be provided as where indicates the minimum SNR that the master node is able to receive and recover the transmitted signal from slave nodes.

5. Connectivity Analysis

From the problem formulation in Section 3, the probability participates in the connectivity analysis. This section provides the derivations of in direct and multi-hop fashions. As the analysis form of in multi-hop is very complicated, the lower and upper bounds are developed, which can provide the guidelines for the networks design in practice.

5.1. Directly Connected

Without loss of generality, it is supposed that the is a cube. There is a surface BS and m clusters () distributed in the . According to the homogeneous distribution of the clusters, is the probability that the jth cluster falls into the hemisphere with surface BS as the center of the sphere and the radius of . We define the area where clusters can establish a direct link with surface BS. Therefore, can be derived as where . However, when considering the boundary effect, the transmission range gets larger with the increase of the number of intra-cluster nodes. Based on the values of and R, can be calculated as:

Case 1: When , is a hemisphere, as shown in Figure 3a. We have

Figure 3

The geometries of underwater MI-assisted acoustic cooperative MIMO networks when directly connected.

Case 2: When , the hemisphere partially extends beyond the sides of the cube to form four hemispheric crowns, as shown in Figure 3b. We have

Case 3: When , the hemisphere extends beyond the sides of the cube and does not reach the bottom of the cube, as shown in Figure 3c. is calculated by calculus. We have

Case 4: When , the hemisphere extends beyond the sides of the cube and extends beyond the bottom of the cube to form a spherical crown, as shown in Figure 3d. We have

Case 5: When , continues to increase, the hemisphere does not completely contain the entire cube, as shown in Figure 3e. We have

Case 6: When , the cluster at any point in can be connected to BS, and then .

Based on the above analysis, can be calculated as Equation (39).

5.2. Lower Bound of Connectivity Probability in Multi-Hop Fashion

Assuming there is a connected link with a relay cluster B between the cluster A and surface BS, there are a nodes in cluster A and b nodes in relay cluster B. Then, the connectivity probability can be calculated as where is the probability that the cluster A has a nodes, is the probability that the relay B has b nodes, is the probability that the cluster A is connected with relay B, and is the probability that the relay B is connected with surface BS. Furthermore,

From the above analysis, the lower bound of the connectivity probability can be obtained by assuming that all clusters have the same number of nodes, which means that the maximum communication range of all clusters is the same.Therefore, we can obtain where is the connectivity probability that the cluster with k nodes is connected with surface BS in multi-hop fashion, and it is worth noting that the number of nodes in the cluster as a relay is also k.

To derive the lower bound of connectivity probability of jth cluster with k nodes , we first map the networks on a discrete cube, as shown in Figure 4 The straight line cb connecting the jth cluster and the surface BS is set as the horizontal line. The jth cluster is located at the center of the cube and sides of the cube are parallel to cb. The length of each side is , which means that clusters in two adjacent cubes can be connected with each other. The clusters are distributed according to a homogeneous Poisson point process; as a result, they can be randomly divided into different cubes.

Figure 4

The geometry of underwater MI-assisted acoustic cooperative MIMO networks in multi-hop fashion.

If the jth cluster can connect with the surface BS through multi-hop, there could be one or more paths that are distinguished by different relay choices. The set of these paths is denoted as , where denote all possible paths. When the cluster density is high, we use the maximum probability that a certain path exists as the lower bound of , i.e., . When the cluster density is low, we calculate the lower bound of the probability that there is at least one open path, i.e., 1 − . The lower bound is , where is the number of the existing paths and means that there is no paths between the jth cluster and surface BS.

To get the valid lower bound of connectivity probability , the larger of the two bounds is selected as the lower bound of . Therefore,

For the different paths that the jth cluster can connect with surface BS, the shorter path can yield the higher connectivity probability. Therefore, the shortest path of the relay nodes that exist in the cubes for which a straight line passes can yield the maximum existing probability. The number of cubes W can be derived as where D is the distance between the jth cluster and the surface BS and means rounding a to the nearest integer no less than a.

When the cluster density is high, the lower bound of connectivity probability is where is the probability that there is no cluster in a cube.

When the cluster density is low, we need to calculate the probability . When the jth cluster cannot connect with the surface BS, there must be a hull C that contains the jth cluster and excludes the surface BS, as shown in Figure 4. The cubes inside the hull contain at least one cluster, while those cubes outside the hull have no cluster. Hence, can be evaluated by counting the number of hulls and the hulls must be passed by the straight line .

Firstly, the number of faces of a hull is denoted as n, the number of hulls that begin at a face and has n faces is denoted as . Obviously, is less than the number of self-avoiding walks beginning at a face, i.e., . The total number of such closed hulls is denoted as

There are relays on one side of the hull and no relays on the other side. Hence, the probability that the jth cluster cannot connect with the surface BS is

Substituting Equations (45) and (47) into Equation (43), we have

5.3. Upper Bound of Connectivity Probability in Multi-Hop Fashion

To get the upper bound of connectivity probability, we consider the probability in terms of hops. Therefore, the connectivity probability of the jth cluster is where is probability that the jth cluster is connected with surface BS directly, is probability that the jth cluster is connected with surface BS through one relay, and is probability that the jth cluster is connected with surface BS through two relays. From Equation (49), we can derive that the more paths there are, the greater is the probability of connectivity. The connectivity probability increases with the number of interconnected nodes until there are no isolated nodes. The probability that there are no isolated clusters is the upper bound for the connectivity probability:

These clusters are distributed according to a homogeneous Poisson point process with density , where the probability of the cluster is isolated is

Substituting Equation (52) into Equation (51), we have

6. Simulation Results

To evaluate the accuracy of developed mathematical models, some simulations were carried out. The simulation conditions were set as follows. The was a 20 km × 20 km × 20 km cube. The clusters were deployed according to a homogeneous Poisson point process with the density of , while the intra-cluster nodes were distributed according to another homogeneous Poisson point process with the density of . The transmitting power of each node was 10 dBm, and the underwater noise level was dBm. The MI communication was operated at the frequency of 10 MHz with the bandwidth of 20 kHz. The underwater acoustic communication was performed at the frequency of 10 kHz with the bandwidth of 10 kHz. The frequency of crystal oscillator used for nodes’ clock synchronization was 100 kHz. The packet lengths of and were set to 100 bits. The propagation speeds of MI and sound were m/s and 1500 m/s, respectively.

6.1. Transmission Range

The maximum transmission ranges of the cluster are plotted in Figure 5. With the increase of SNR threshold, the maximum transmission range decreases subsequently, which is consistent with the analysis. At the same time, connectivity is affected by communication range of the sensor node, which is enlarged by cooperative MIMO compared with single node. We show that the cooperative MIMO using MI-assisted synchronization and pure acoustic-based synchronization both have a great improvement in maximum transmission range compared with a single node in Figure 5. Therefore, the underwater MI-assisted acoustic cooperative MIMO networks can leverage the gains of MIMO in a distributed fashion to enhance connectivity of UWSNs. Meanwhile, MI communication works at ultra-high frequency and the signal transmission delay can be neglected. Thus, the influence of frequency and time estimation errors on its SNR is almost 0, and the maximum transmission range is almost the same as with perfect synchronization. However, in a pure acoustic-based synchronization system, due to the existence of error, the maximum transmission range is obviously less than the perfect synchronization.

Figure 5

The inter-cluster maximum transmission range.

To explore the transmission range versus the variety of coil radius and SNR threshold, a simulation experiment was conducted, and the results are indicated in Figure 6. As shown in Figure 6, the range of the cluster can be extended by increasing the coil radius. Compared with underwater acoustic communication, MI communication has a smaller transmission range, so it is suitable for short distance inter-cluster communication. At the same time, it is worth noting that there is great difference in the transmission range of the MI communication between sea water and lake water. Highly conductive sea water induces significant Eddy current that incurs very high path loss, which results in the inevitable decrease of the transmission range. However, the transmission range of underwater MI communications can be increased by using the optimal operating frequency and larger coil antennas.

Figure 6

The intra-cluster maximum transmission range.

6.2. The Directly Connected Probability of the jth Cluster Consisting of k Nodes

To get the directly connected probability of the jth cluster with consisting of k nodes, i.e., , it is necessary to calculate the probability that the actual distance between the surface BS and the jth cluster is no more than the maximum transmission range , i.e., . To obtain the statistical performance, Monte Carlo simulations were performed. In a certain simulation, an underwater MI-assisted acoustic cooperative MIMO network was generated. If , the result was marked as 1; if , the result was marked as 0. Finally, the probability was obtained by the 10,000 simulation results. As shown in Figure 7, the simulation results verify the correctness of the theoretical analysis in Equation (39).

Figure 7

The directly connected probability of the jth cluster with consisting of k sensors .

6.3. Directly Connected Probability of the jth Cluster with Different SNR Threshold

The directly connected probability of the jth cluster with different SNR thresholds is investigated in Figure 8. To evaluate the performances of the proposed approach, 5000 Monte Carlo simulations were carried out in each case. Figure 8 shows that the higher SNR threshold leads to the smaller connectivity probability under the same node density. As increases from 5 to 20 dB, the transmission range of sensor gets shorter; as a result, fewer sensors participate in the cooperative MIMO, which means the smaller transmission range. Therefore, the probability of connectivity is lower. Note that, when node density is low, there is notable gap between analysis and simulation results. That is because we can gat that the transmission range of a cluster has a great influence on the theoretical analysis of its connectivity from Equation (39). While the cluster communicates with the BS in the form of cooperative MIMO, the number of nodes participating in the cooperative MIMO determines its transmission range. When the connectivity of the jth cluster is theoretically analyzed, the number of nodes may be non-integer, which is calculated by . However, the number of nodes in the cluster is determined in simulation. When the node density is low, there are few nodes participating in cooperative MIMO, and the non-integer number of nodes causes notable gap between analysis and simulation results, as shown in Figure 8.

Figure 8

The directly connected probability of the jth cluster with different SNR thresholds.

6.4. Theoretical Bounds of the jth Cluster

In Figure 9, the theoretical upper and lower bounds of underwater MI-assisted acoustic cooperative MIMO networks are compared with the simulation results with a variety of sensor densities. Each simulated connectivity probability was calculated based on 500 simulation iterations. The SNR threshold of the received signal was set to 20 dB. The connectivity probability increases as the cluster density increases. This result is consistent with the connectivity analysis and provide the guideline for the deployment of the networks according to the practical requirement.

Figure 9

The theoretical bounds of the jth cluster with dB.

6.5. Directly Connected Probability of the jth Cluster

The jth cluster’s directly connected probabilities in the cases of perfect synchronization, pure acoustic-based synchronization, and MI-assisted synchronization were investigated under the same simulation conditions, and the results are indicated in Figure 10. The connectivity probabilities increase as the node density increases. As the node density in the cluster becomes higher, it demonstrates that more sensors within a cluster participate in the cooperative MIMO, resulting in the long-range transmission. Considering the frequency and time estimation errors in the synchronization process, the maximum transmission range of the cluster is smaller than that of the perfect synchronization. Therefore, for the cluster with the same sensor density, the connectivity probability becomes smaller.

Figure 10

The directly connected probability of the jth cluster with dB.

7. Conclusions

Recently, underwater MI-assisted acoustic cooperative MIMO networks have attracted much attention as they significantly improve the effective cover space and network throughput. However, the complex channel characteristics and the heterogeneous architecture make the connectivity more complex than that of traditional UWSNs. In this paper, we investigate the connectivity of the networks considering the influences of channel characteristics, system parameters, and synchronization error. To provide the theoretical guide for the design of underwater MI-assisted acoustic cooperative MIMO networks, the upper and lower bounds are both derived. Finally, Monte Carlo simulations were performed to address the correctness of the proposed mathematical models.