1. Introduction
Underwater wireless optical communication (UWOC) is a useful way to transmit data for short distance applications in the ocean, thanks to its high transmission speed [
1,
2,
3]. Recently, many researchers have focused on this topic. They are mainly investigating increasing the communication rate with high-order modulation [
4,
5,
6,
7] and high-performance devices [
8], or improving the communication distance with a high sensitivity detector [
9] and photon counting detection [
10,
11]. Meanwhile, in the area of ocean engineering, several sea trials have demonstrated its potential application prospects. From 2008 to 2014, Woods Hole Oceanographic Institution (WHOI) demonstrated the value of this technology in deep sea with up to 10M bps speed, more than 100 m distance, and more than 2000 m depth [
12,
13,
14,
15,
16]. In these trials, UWOC devices were usually carried on mobile platforms, such as human occupied vehicles (HOV), remote operated vehicles (ROV) and autonomous underwater vehicles (AUV). They demonstrate the potential value of using such technology to harvest data from subsea nodes and seafloor observatories by AUV and optical modem [
12,
13,
14], or transmitting data between ROV, HOV, and other underwater platforms [
15,
16]. Using the same technical scheme, commercial UWOC products have been applied in ocean engineering with up to 10M bps rate. It has been used to stream ocean exploration missions live [
17]. Furthermore, in 2017, the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) demonstrated 20M bps speed with 120 m distance in dark clear water. In the system, the UWOC device was carried on an ROV [
18].
These sea trials show that carrying the devices on mobile platform is a practical application of UWOC. In this way, UWOC devices can move underwater to establish optical communication links with cooperative targets, which is helpful to exploit its advantages. However, while UWOC devices and such platforms work in the ocean, they usually need a camera to observe the surroundings and guide their action. In this situation, because the majority of ocean is always dark [
19], active illumination is necessary for imaging. When UWOC works in such an environment, the illumination light affects its performance, and can even make it unworkable. Therefore, it is necessary to study the influence of light illumination on UWOC. As an underwater light emitting diode (LED) is usually applied to illumination, in this work, we analyze the influence of such light on bidirectional UWOC with the Monte Carlo method.
The organization of this paper is as follows:
Section 2 sets up a theoretical model to study light noise with the Monte Carlo method.
Section 3 demonstrates the simulation results.
Section 4 makes a conclusion.
2. Theoretical Model
The simulation model is shown in
Figure 1. It has two identical terminals, the distance between which is
S. In each terminal, there is a UWOC transmitter, a UWOC receiver, and an underwater LED device. In the two terminals, the diameters of the luminous surface for LEDs are both
D. The diameters of the UWOC receivers are both
d. The separated distances between the UWOC receiver and LED are both
L. When the LEDs and UWOC system work simultaneously, the illumination noise would arrive to both the two receivers, because there is scattering in the water. For example, light from “LED 1#” in
Figure 1 affects “Receiver 1#” and “Receiver 2#” simultaneously.
We simulated such influence with the Monte Carlo method. This method is quite popular to study underwater optical channels [
20]. It is based on the fact that light is made up of many photons. Supposing the process of light propagating underwater is a linear time-invariant system, the transmitting characteristic of every photon follows the same statistical rule. Thus, by tracing the photons’ paths underwater one by one and counting their characteristics, the light feature can be obtained. For example, the time-domain feature of light can be analyzed by counting the photons’ optical path underwater, while the characteristics of the light field can be studied by counting the photons’ positions. In this paper, in order to study the influence of LED illumination noise on UOWC receiver, we count the photons from LED into receivers with different positions. There are three steps.
(1) The process of photons from LED into water is simulated one by one. The coordinates at which that the photons arrived at the plane of the two UWOC receivers are recorded.
(2) The photons whose coordinates are in the receiving surface of UWOC receivers and arriving angle is within the field of view (FOV) of receiver are counted. Then, the number of photons arrived at the receivers is acquired, which includes the opposite receiver in the cooperative terminal and the adjacent receiver in the same terminal.
(3) The numbers of photons arrived at the receivers are divided by the total number we simulated. Then, the relative intensity in the two receivers is obtained. They can be used to evaluate the forward noise from LED illumination to the opposite receiver and backscattering noise to the adjacent receiver, respectively.
As shown in
Figure 1, for simplicity, we study the influence of light from “LED 1#” on “Receiver 1#” and “Receiver 2#”. As shown in
Figure 2, this is supposing the luminous surface of “LED 1#” is uniform and round with diameter
D. We establish coordinate system with the center of “LED 1#”. Photons’ position and direction cosine are expressed as follows.
Here, (ρ, α) is the polar coordinate for the photon position in the plane z = 0. θ0 and φ0 are theelevation angle and azimuth angle for the starting direction. randr, randJ, and randφ are random numbers, which follow uniform distribution in range of 0 and 1. Random numbers are produced by a computer, which will appear in several places below in the model.
The light intensity of LED usually follows Lambert distribution. Its radiation intensity
I(
θ0) with an elevation angle
θ0 can be expressed as follows,
where
I0 is the intensity at
θ0 = 0. Furthermore, it can be understood that the photon probability density at
θ0 is cos
θ0. Then, a random number
randθ can be produced with an arbitrary
θ0, as shown in Equation (4).
In Equation (4),
randθ is a random number, which follows uniform distribution in the range of 0 and 1. Therefore,
θ0 can be expressed as follow.
Once a photon leaves “LED 1#”, it will transmit straight in water until arriving at the scattering particle. The free path
s is
Here
a,
b, and
c are the absorption coefficient, scattering coefficient, and attenuation coefficient of water, respectively.
rands is a random number. Then, the new position (
x,
y,
z) for photon is
Once
z ≤
S and
z ≥ 0, the photon will be absorbed or scattered in water, which is determined with the method of “Russian roulette”. That means we firstly should define a threshold
WH as follows:
Then, we produce a random number
randW from computer. If
randW >
WH, the photon is absorbed. We stop the tracing process. Otherwise, it is scattered, the new direction (
μx,
μy,
μz) of which is as follows:
φ and
θ are the azimuth and scattering angle, respectively. We suppose
φ obeys uniform distribution as shown in Equation (10), in which
randφ1 is a random number.
θ is determined by the volume scattering function (VSF). There are several models for it, such as the Petzold average particle phase function, Fournier- Forand (FF) function, Henyey–Greenstein (HG) function and its related approximations, etc. Only the HG function can give an analytical expression for scattering angle in a Monte Carlo simulation [
20,
21]. Although some research indicates that the function value in this model is less than the real water with the scattering angle close to zero or in the backward direction, the simulation error is acceptable in most cases. Hence, we will use the HG function to generate the scattering angle as shown in Equation (11), not only for simplicity, but also because the analytical expression is useful to reduce the computing error in simulation.
g is the asymmetry parameter that depends on the characteristics of medium, which equals to the average cosine of the scattering angle over all scattering directions. In this paper, we set
g = 0.924, because it is considered as a good approximation for most practical situations [
22].
randθ1 is a random number. Then, the process of a photon’s motion in water can be described with Equations (1)–(11).
Once
z ≤ 0, we modify the coordinate to the plane
z = 0 by Equation (12). We can use these photons to study the backscattering noise from “LED 1#” on “Receiver 1#”. If
z ≥
S we modify the coordinate to the plane
z =
S by Equation (13). We can use these photons to study the noise from “LED 1#” on “Receiver 2#”. Then, illumination from the forward light and backscattering noise on receivers are both acquired.
The light distribution from “LED1#” is usually symmetrical to the center. We set the center of “Receiver 1#” and “Receiver 2#” at coordinates (0,
L, 0) and (0,
L,
S) for simplicity. If a photon’s coordinate (
X1,
Y1,
Z1) satisfy Equation (14), and its arriving angle to the receiver is within the field of view (FOV), we record it and stop the tracing.
Here,
d is the diameter of receiver. We denote the total number of tracing photons as
N. If the number of photons arrived at “Receiver 2#” and “Receiver 1#” are
NF and
NB, respectively. The relative intensity
ηF and
ηB are denoted as Equations (15) and (16), which can be used to evaluate the forward noise from LED illumination to the opposite receiver and backscattering noise to the adjacent receiver.
3. Simulation Result
In order to study the influence of illumination noise, we set the diameter of the luminous surface D as 0.1 m for “LED 1#”, the diameter d of “Receiver 1#” and “Receiver 2#” as 0.1 m, and the total number N of tracing photons as 109. Then, we analyze the variation of relative intensity ηB and ηF with the separated distance L from UWOC receiver to the optical axis of LED light, by setting a different absorption coefficient a, scattering coefficient b, transmitting distance S, and FOV of the receiver.
3.1. Analysis of the Influence of Illumination Noise with Different Absorption Coefficient
In the real ocean, the absorption coefficient and scattering coefficient (
a,
b) are different with water [
20], which can be close to zero and larger than one. Then, the communication distance
S is inconstant. Usually, the UWOC system is suitable for relatively clean water. Therefore, we set
b = 0.10 m
−1,
S = 20 m, while
a = 0.05 m
−1, 0.10 m
−1, and 0.15 m
−1, respectively. Then, we count the number of photons arriving at “Receiver 2#” and “Receiver 1#” to calculate the relative intensity
ηF and
ηB, when FOV is 60°, 120°, and 180°, respectively. The results
ηF and
ηB vary with
L and are shown in
Figure 3 and
Figure 4, respectively.
In
Figure 3 and
Figure 4, the general tendency of
ηF and
ηB versus
L is decreasing and the tendency of
ηB is more evident. When absorption coefficient
a is greater, or the FOV is smaller, the relative intensity
ηF and
ηB would be smaller. Furthermore, the absorption coefficient has a greater effect on
ηF, while the FOV has more influence on
ηB.
It needs to be noted here that, in
Figure 4, when
L is larger or
a is smaller, the monotonicity of the curves is not perfect. This is because the number of tracing photons is 10
9 in the simulation. When
L is larger or
a is smaller, the number of photons arrived at the receiver is very small. Then, the error is larger. This is caused by the Monte Carlo method itself because it uses a stochastic process to study a physical problem. The accuracy is determined by the sample size. We tried to increase the number of tracing photons to 10
10 to reduce this error. The computing time increased tremendously. The monotonicity of the curves became better, but the error still existed. If the number became larger, the monotonicity would be more perfect, but the error could not be removed thoroughly. As the monotonicity of curves in
Figure 4 is modified with the number of tracing photons, the tendency we observe here is credible. This kind of situation would exist in the following simulation.
Therefore, we conclude that the forward noise from LED illumination to the opposite receiver and backscattering noise to the adjacent receiver are both decreasing with the absorption coefficient and the change of forward noise is more evident. Narrowing the receiver’s FOV would reduce such noises. Enlarging the separated distance between receivers and the optical axis of LED would reduce such noises too, while it is more evident to backscattering noise.
3.2. Analyze the Influence of Illumination Noise with Different Scattering Coefficient
We set
a = 0.10 m
−1,
S = 20 m, while
b = 0.05 m
−1, 0.10 m
−1 and 0.15 m
−1, respectively. Then, we count the number of photons arriving at “Receiver 2#” and “Receiver 1#” to calculate the relative intensity
ηF and
ηB, when the FOV is 60°, 120°, and 180°, respectively. The results
ηF and
ηB vary with
L and are shown in
Figure 5 and
Figure 6, respectively.
In
Figure 5, with the increase of the separated distance
L, the relative intensity
ηF decreases. When the scattering coefficient
b is greater, or the FOV is smaller, the relative intensity
ηF is smaller. At a smaller FOV and at larger distances, the relative intensities converge for different scattering coefficients. The difference is quite small. In
Figure 6, with the increase of the separated distance
L, the relative intensity
ηB decreases. Its tendency is more evident than
ηF, especially at a smaller FOV and at larger distances. When the scattering coefficient
b and the FOV are smaller, the relative intensity
ηF is smaller.
Therefore, we conclude that the forward noise from LED illumination to the opposite receiver decreases with the scattering coefficient, while the backscattering noise to the adjacent receiver increases with it. Narrowing the receiver’s FOV would reduce such noises, while the effect is more evident with backscattering noise. Enlarging the separated distance between receivers and the optical axis of LED would reduce such noises too, while it is more evident to backscattering noise.
3.3. Analyze the Influence of Illumination Noise with Different Transmitting Distance
We set
a = 0.10 m
−1,
b = 0.10 m
−1, and S = 20 m, 25 m, and 30 m, respectively. Then, we count the number of photons arriving at “Receiver 2#” and “Receiver 1#” to calculate the relative intensity
ηF and
ηB, when FOV is 60°, 120°, and 180°, respectively. The results
ηF and
ηB vary with
L and are shown in
Figure 7 and
Figure 8 respectively.
In
Figure 7, with the increase of separated distance
L, the relative intensity
ηF decreases. When the transmitting distance
S is greater, or the FOV is smaller, the relative intensity
ηF is smaller. In
Figure 8, with the increase of separated distance
L, the relative intensity
ηB decreases and its tendency is more evident than
ηF. When the FOV is smaller, the relative intensity
ηF is smaller too. When the transmitting distance
S varies,
ηF nearly has no change.
Therefore, we conclude that the forward noise from LED illumination to the opposite receiver decreases with the transmitting distance, while the backscattering noise to the adjacent nearly has no relation with it. Narrowing the receiver’s FOV would reduce such noises, while the effect is more evident with backscattering noise. Enlarging the separated distance between receiver and the optical axis of LED would reduce such noises too, while it is more evident with backscattering noise.
3.4. Discussion
From the results above, it is clear that forward noise from the LED illumination affects the opposite receiver, while the backscattering noise impacts the adjacent receiver in the same terminal. With the increase of separated distance between receiver and the optical axis of LED, both the noises decrease. The tendency is more evident for the backscattering one. With the increase of absorption coefficient, the two noises are decreased, while the change of the forward noise is more evident. With the increase of scattering coefficient, the forward noise to the opposite receiver decreases, while the backscattering noise to the adjacent receiver increases. With the increase of the transmitting distance, the forward noise to the opposite receiver decreases, while the backscattering noise to the adjacent nearly has no change. With the increase of FOV, both the noises increase. What needs illustration is that when we change other simulation parameters, such results still stand.
The simulation results above can be explained theoretically. For the opposite receiver, the forward noise is determined by the loss in the transmitting channel. When increasing the absorption coefficient, scattering coefficient, transmitting distance, and the separated distance between receiver and the center of LED illumination, or narrowing the FOV of receiver, the loss becomes greater. Then, the forward noise is less.
For the adjacent receiver in the same terminal, the backscattering noise is mainly determined by the backscattering light near the receiver. When the absorption coefficient enlarges, the noise light into receiver reduces a little. When the scattering coefficient enlarges, the backscattering becomes more evident and the noise is increased. When the transmitting distance increases, the backscattering has little change. When FOV reduces, or the separated distance between receiver and the center of LED increases, less light comes into the receiver.
Furthermore, from the results above, it is clear that backscattering noise from LED illumination to the adjacent receiver in the same terminal cannot be neglected in application, in particular when there is strong scattering in the water.
In order to reduce such light noise from LED illumination, besides inserting an optical filter in the receivers, the FOV should be smaller and the separated distance between receiver and the center of LED illumination should be longer. These two methods are more useful to the backscattering noise of the adjacent receiver.