分层纵向水下无线光信道建模研究

孙若凡; 刘学蕊; 娄毅; 乔钢; SUN Ruofan; LIU Xuerui; LOU Yi; QIAO Gang

引 en

分层纵向水下无线光信道建模研究

孙若凡^1,2,3

机构：

1. 哈尔滨工程大学水声技术重点实验室，黑龙江哈尔滨 150001

2. 海洋信息获取与安全工信部重点实验室（哈尔滨工程大学）工业和信息化部，黑龙江哈尔滨 150001

3. 哈尔滨工程大学水声工程学院，黑龙江哈尔滨 150001

×
，刘学蕊^1,2,3

机构：

1. 哈尔滨工程大学水声技术重点实验室，黑龙江哈尔滨 150001

2. 海洋信息获取与安全工信部重点实验室（哈尔滨工程大学）工业和信息化部，黑龙江哈尔滨 150001

3. 哈尔滨工程大学水声工程学院，黑龙江哈尔滨 150001

×
，娄毅^1,2,3,

机构：

1. 哈尔滨工程大学水声技术重点实验室，黑龙江哈尔滨 150001

2. 海洋信息获取与安全工信部重点实验室（哈尔滨工程大学）工业和信息化部，黑龙江哈尔滨 150001

3. 哈尔滨工程大学水声工程学院，黑龙江哈尔滨 150001

×
，乔钢^1,2,3

机构：

1. 哈尔滨工程大学水声技术重点实验室，黑龙江哈尔滨 150001

2. 海洋信息获取与安全工信部重点实验室（哈尔滨工程大学）工业和信息化部，黑龙江哈尔滨 150001

3. 哈尔滨工程大学水声工程学院，黑龙江哈尔滨 150001

×

1. 哈尔滨工程大学水声技术重点实验室，黑龙江哈尔滨 150001；
2. 海洋信息获取与安全工信部重点实验室（哈尔滨工程大学）工业和信息化部，黑龙江哈尔滨 150001；
3. 哈尔滨工程大学水声工程学院，黑龙江哈尔滨 150001；

Modeling of Hierarchical Vertical Underwater Wireless Optical Channel

SUN Ruofan^1,2,3

Affiliation：

1. Acoustic Science and Technology Laboratory，Harbin Engineering University，Harbin 150001 ，China

2. Key Laboratory of Marine Information Acquisition and Security（Department of Industry and Information Technology），Harbin Engineering University，Harbin 150001 ，China

3. College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001 ，China

×
， LIU Xuerui^1,2,3

Affiliation：

1. Acoustic Science and Technology Laboratory，Harbin Engineering University，Harbin 150001 ，China

2. Key Laboratory of Marine Information Acquisition and Security（Department of Industry and Information Technology），Harbin Engineering University，Harbin 150001 ，China

3. College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001 ，China

×
， LOU Yi^1,2,3,

Affiliation：

1. Acoustic Science and Technology Laboratory，Harbin Engineering University，Harbin 150001 ，China

2. Key Laboratory of Marine Information Acquisition and Security（Department of Industry and Information Technology），Harbin Engineering University，Harbin 150001 ，China

3. College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001 ，China

×
， QIAO Gang^1,2,3

Affiliation：

1. Acoustic Science and Technology Laboratory，Harbin Engineering University，Harbin 150001 ，China

2. Key Laboratory of Marine Information Acquisition and Security（Department of Industry and Information Technology），Harbin Engineering University，Harbin 150001 ，China

3. College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001 ，China

×

1. Acoustic Science and Technology Laboratory，Harbin Engineering University，Harbin 150001 ，China；
2. Key Laboratory of Marine Information Acquisition and Security（Department of Industry and Information Technology），Harbin Engineering University，Harbin 150001 ，China；
3. College of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001 ，China；

作者简介:

孙若凡(1996-),男,博士生,主要从事水下无线光通信系统的性能分析研究。

通讯作者:

娄毅(1987-),男,博士,副教授,主要从事高速水下光通信和全双工水声通信的研究。

中图分类号:TN929.1

文献标识码:A

文章编号:2096-5753(2022)04-0322-07

DOI:10.19838/j.issn.2096-5753.2022.04.007

全文
评论
参考文献
出版信息

参考文献 1

殷敬伟,惠俊英,郭龙祥.点对点移动水声通信技术研究[J].物理学报,2008(3):1753-1758.

查找原文

参考文献 2

YANG M,LI X K,MENG X X.Characteristic analysis of underwater acoustic scattering echoes in the wavelet transform domain[J].船舶与海洋工程学报(英文版)2017(1):93-101.

查找原文

参考文献 3

ANOUS N,ABDALLAH M,UYSAL M,et al.Performance evaluation of LOS and NLOS vertical inhomogeneous links in underwater visible light communications[J].IEEE Access,2018,6:22408-22420.

查找原文

参考文献 4

ELAMASSIE M,UYSAL M.Performance characterization of vertical underwater VLC links in the presence of turbulence[C]//2018 11^th International Symposium on Communication Systems,Networks & Digital Signal Processing(CSNDSP).Budapest:IEEE,2018.

查找原文

参考文献 5

ELAMASSIE M,SAIT S M,UYSAL M.Underwater visible light communications in cascaded gammagamma turbulence[C]//2018 IEEE Globecom Workshops(GC Wkshps).Abu Dhabi:IEEE,2018.

查找原文

参考文献 6

LOU Y,CHENG J L,NIE D H,et al.Performance of vertical underwater wireless optical communications with cascaded layered modeling[J].IEEE Transactions on Vehicular Technology,2022,71(5):5651-5655.

查找原文

参考文献 7

YILMAZ A,ELAMASSIE M,UYSAL M.Diversity gain analysis of underwater vertical MIMO VLC links in the presence of turbulence[C]//2019 IEEE International Black Sea Conference on Communications and Networking(BlackSeaCom).Sochi:IEEE,2019.

查找原文

参考文献 8

CHRISTOPOULOU C,SANDALIDIS H G,VAVOULAS A.Optimal trade-off between depth and coverage for a vertical UWOC link[J].IEEE Photonics Technology Letters,2020,32(4):208-211.

查找原文

参考文献 9

CHRISTOPOULOU C,SANDALIDIS H G,VAIOPOULOS N.Performance of an underwater optical wireless link with a randomly placed or moving receiver[J].IEEE Journal of Oceanic Engineering,2020,46(3):1068-1079.

查找原文

参考文献 10

SAIT M,GUO Y,ALKHAZRAGI O,et al.The impact of vertical salinity gradient on non-line-of-sight underwater optical wireless communication[J].IEEE Photonics Journal,2021,13(6):1-9.

查找原文

参考文献 11

JI X Y,YIN H X,JING L Y,et al.Modeling and performance analysis of oblique underwater optical communication links considering turbulence effects based on seawater depth layering[J].Optics Express,2022,30(11):18874-18888.

查找原文

参考文献 12

HALTRIN V I.Chlorophyll-based model of seawater optical properties[J].Applied Optics,1999,38(33):6826-6832.

查找原文

参考文献 13

YAP Y X,JASMAN F,MARCUS T C E.Impact of Chlorophyll concentration on underwater optical wireless communications[C]//2018 7^th International Conference on Computer and Communication Engineering(ICCCE).Kuala Lumper:IEEE,2018.

查找原文

参考文献 14

JOHNSON L J,GREEN R J,LEESON M S.Underwater optical wireless communications:depth dependent variations in attenuation[J].Applied Optics,2013,52(33):7867-7873.

查找原文

参考文献 15

HALTRIN V I.Theoretical and empirical phase functions for Monte Carlo calculations of light scattering in seawater[C]//Fourth International Conference on Remote Sensing for Marine and Coastal Environments.Orlando:Environmental Research Institute of Michigan,1997.

查找原文

参考文献 16

COX W C.Simulation,modeling,and design of underwater optical communication systems[M].Los Angeles:North Carolina State University,2012.

查找原文

目录contents

摘要Abstract
关键词Keywords
0 引言
1 系统和信道模型
1.1 系统模型
1.2 吸收与散射模型
1.2.1 吸收系数
1.2.2 散射系数
1.2.3 衰减系数
1.3 散射相位函数
2 蒙特卡洛仿真
2.1 光子的初始状态
2.2 光子的传播
2.3 光子的接收
3 仿真结果与讨论
4 结束语
参考文献

摘要

水下无线光通信（Underwater Wireless Optical Communication，UWOC）是水下物联网感知层的新兴技术之一。在感知层中，UWOC 技术常被用于支持数据的纵向上行传输。针对几何光学理论无法准确描述 UWOC 链路衰减的问题，研究了基于 MIE 散射理论的纵向 UWOC 信道模型；提出了分层的纵向海水介质模型，在此基础上设计了一种利用蒙特卡洛仿真方法对 UWOC 纵向链路进行评估的方案。该方案考虑了收发装置空间位置的非对齐以及出现位置失配的情况，通过该方案可求得 UWOC 链路冲激响应的数值结果，研究成果对于 UWOC 设计具有重要意义。

Abstract

Underwater Wireless Optical Communication（UWOC）is one of the emerging technologies for the underwater IoT sensing layer. In the sensing layer，UWOC technology is often used to support the vertical uplink transmission of data. This paper investigates a vertical UWOC channel model based on MIE scattering theory to address the problem that geometric optical theory cannot accurately describe the UWOC link fading. A hierarchical vertical seawater medium model is proposed，based on which this paper designs a scheme for evaluating UWOC vertical links using Monte Carlo simulation methods. The scheme considers the non-alignment of the spatial position of the transceiver and the occurrence of position mismatch. This scheme can obtain the numerical results of the impulse response of the UWOC link. The results of the study are of great significance for UWOC design.

关键词

水下无线光通信；蒙特卡洛方法；信道建模

Keywords

underwater wireless optical communication ； Monte Carlo method ； channel modeling

0 引言
水下物联网（Internet of Underwater Things， IoUT）是物联网的水下延伸，被定义为通过数字设备使水下物体进行智能互联的水下设备网络。 IoUT为人类探索海洋环境提供了崭新的方式，其应用场景包括海洋资源勘探、海洋环境检测、水下设备维护、军事防御等。IoUT依网络构架可被划分为网络层、应用层以及感知层。其中感知层被应用于识别物体以及收集信息，例如海底传感器、水下航行器、水下中继节点等。随着人类水下活动的日益频繁以及复杂化，为了达到智能互联的目标， IoUT对其感知层收集信息的效率提出了更高的要求，对水下高性能通信技术的需求日益提升。
IoUT感知层使用的水下无线通信技术主要包含3种：水声通信、射频通信以及水下无线光通信（UWOC）技术。研究证明，在短距离传输中，与水声通信技术相比，UWOC可以提供更高的数据速率，而且功耗显著降低，设备使用寿命更长。然而与水声通信类似，当环境发生改变时，被设计用于原水下环境的UWOC设备可能失效^[1-2]。因此在设计环节对UWOC进行预估是推动系统实际布放的关键。准确的UWOC信道模型是进行性能分析的基础。
目前UWOC建模的研究集中在水平及固定距离的链路上。在水平链路上，光传播常被视作具有恒定衰减特性，其信道衰减模型通常具有较为简单的数学形式。而在垂直链路中，由于温度、压力、溶解物成分及浓度出现深度依赖性变化，使得垂直链路的衰减特性并不恒定。近年来，只有较少的纵向UWOC信道模型被提出。LOS通信方面，2018年，ANOUS首次对纵向UWOC信道进行了建模以及性能分析。ANOUS将纵向信道近似为若干层的均匀信道，应用几何光学理论，将传播损失建模为指数衰减与分界面透射率乘积的形式^[3]。其后ELAMASSIE考虑了湍流对纵向链路的影响。在^[4-5] 中，ELAMASSIE分别计算了湍流衰落服从DoubleGamma、对数正态分布情况下纵向链路传播损失表达式。此外，在文献^[6]中，娄考虑了湍流衰落服从对数-Gamma分布时纵向链路的传播损失。以Elamassie的模型为基础，2019年，YILMAZ首次对受对数正态分布湍流影响的多入多出（MIMO） UWOC系统的传播损失进行了建模^[7]。非视距通信（NLOS）方面，CHIRSTOPOULOU对纵向链路中光束的覆盖面积进行了建模，并考虑了接收机处于不同位置时UWOC信道传播损失^[8-9]。SAIT实验验证了存在温度梯度和盐度波动时，纵向NLOSUWOC信道传播损失模型^[10]。最近，为符合实际应用情况，Ji提出了一种考虑湍流影响的斜光链路模型^[11]。尽管近年来已经有若干纵向链路的信道损失模型被提出，但据调研，大部分纵向链路研究均以几何光学理论为基础，从波动光学角度解释光传播中的衰减。尚未有基于MC方法的纵向LOS链路模型被提出。与波动理论相比，MC方法能够获得更多的信道信息。
本文的贡献包括：
1）本文提出了基于MC方法的LOS-UWOC纵向链路传播损失模型。其中海洋介质被近似为多层水平的均匀介质，假定介质中的粒子服从MIE散射理论。通过MC方法能够求得信道冲激响应，为进一步研究高速通信情况下码间串扰的影响提供了基础。
2）为使模型更符合实际情况，本文考虑了收发装置非对准对纵向UWOC传播损失的影响。
1 系统和信道模型
1.1 系统模型
在实际水下应用中，光通信常被用于数据的短距离上行传输，例如水下节点与海面基站之间的通信。如图1所示，水下存在一条纵向的上行UWOC链路。发射机是安装了高斯光源的水下传感器节点U。接收机是被吊放在水面舰船底部的光电二极管S。发射机及接收机的位置相对固定，不发生主动位移。在图1所示UWOC系统中，海水介质被近似为两层水平的介质，其中每层介质均具有不随深度变化的恒定衰减特性。介质I的纵向距离为 L，该层介质的吸收、散射及衰减系数分别为 $a (λ) b (λ) c (λ)$ ，且 $c (λ) = a (λ) + b (λ)$ ，其中 $λ$ 是光的波长。介质II的纵向距离为 L，其吸收散射及衰减系数分别为 $a (λ), b (λ), c (λ)$ 。尽管水下光传播的衰减具有波长选择性，但是在本文中波长并不是重点讨论的内容。因此下文论述中将忽视波长变化对衰减的影响，仅考虑波长为532nm时信道的衰减特性，即 $λ =$ 532nm。
图1 纵向UWOC链路的模型
Fig.1 Model of vertical UWOC link
1.2 吸收与散射模型
1.2.1 吸收系数
海洋介质的吸收主要由海水中无机物的固有吸收以及溶解/悬浮有机物的吸收组成。HALTRIN将海水的吸收建模为4种物质的共同作用^[12]，分别是纯水、黄腐酸、腐殖酸以及叶绿素的吸收作用。海水的吸收系数表示为

\begin{matrix} a (L) = a_{w} + a_{f} C_{f} (L) e^{- p_{f} λ_{0}} + \\ a_{h} C_{h} (L) e^{- p_{h} λ_{0}} + a_{c} C_{c} (L)^{0.602} \end{matrix}

(1)

式中： $a_{w}$ 是纯水的吸收系数； $a_{f}$ 是黄腐酸的吸收系数； $a_{h}$ 是腐殖酸的吸收系数； $a_{c}$ 是叶绿素的吸收系数； $p_{f}$ 与 $p_{h}$ 均是模型中的常数因子；L 是深度； $C_{f} (L)$ ， $C_{h} (L)$ ， $C_{c} (L)^{}$ 分别表示黄腐酸、腐殖酸与叶绿素的浓度，其中 $C_{f} (L)$ ， $C_{h} (L)$ 同时受深度及叶绿素浓度的影响。其关系为

\begin{matrix} C_{f} (L) = 1.74098 C_{c} (L) e x p [0.12327 (C_{c} (L) / C_{c}^{0})] \\ C_{h} (L) = 0.19334 C_{c} (L) e x p [0.12343 (C_{c} (L) / C_{c}^{0})] \end{matrix}

(2)

式中， $C_{c}^{0}$ 是用于确定浓度量级的标准浓度， $C_{c}^{0} =$ 1mg/m³。
在文献^[13]中，YAP考虑了纵向链路中叶绿素浓度对传播损失的影响，并提出了叶绿素浓度表达式。

\begin{matrix} C_{c} (L) = C_{b g} + S L + \\ \frac{C_{t o p}}{σ \sqrt{2 π}} e x p (- \frac{1}{2} {(\frac{L - L_{m a x}}{σ})}^{2}) \end{matrix}

(3)

式中：C_bg表示表面的背景叶绿素浓度；S 表示叶绿素浓度的垂直梯度；C_top 指高于背景水平的叶绿素总量； L_max 表示不同类型介质中叶绿素浓度出现最大值的深度，此外 $σ$ 是叶绿素浓度的标准差。不同水体的叶绿素分布特征已经在文献^[14] 中被研究。
1.2.2 散射系数
Haltrin认为，海水介质中的散射主要由3种物质共同作用，分别是纯水、海水中的小微粒以及较大的微粒。折射率是区分微粒大小的参考，通常认为折射率小于1.03的是小微粒，大于1.15的是大微粒。海水的散射系数被表示为

b (L) = b_{w} + b_{s} C_{s} (L) + b_{1} C_{1} (L)

(4)

式中： $b_{w}$ 是纯水的散射系数； $b_{s}$ 和 $b_{1}$ 分别指小微粒和大微粒的吸收系数； $C_{s} (L)$ 和 $C_{1} (L)$ 分别表示小微粒与大微粒的浓度，其浓度与深度及叶绿素浓度有关。

\begin{matrix} C_{s} (L) = 0.01739 C_{c} (L) e x p (0.11631 [C_{c} (L) / C_{c}^{0}]) \\ C_{1} (L) = 0.76284 C_{c} (L) e x p (0.03092 [C_{c} (L) / C_{c}^{0}]) \end{matrix}

(5)

式（5）中，3种物质的散射系数由光的波长决定，其表达式为

\begin{matrix} b_{w} (λ) = 0.005826 {(400 / λ_{0})}^{4.322} \\ b_{s} (λ) = 1.151302 {(400 / λ_{0})}^{1.7} \\ b_{1} (λ) = 0.341074 {(400 / λ_{0})}^{0.3} \end{matrix}

(6)

1.2.3 衰减系数
如式（1）及式（4）所示，海水的吸收与散射特性具有深度依赖性，使得纵向链路的衰减特性并不恒定。为简化蒙特卡洛仿真的流程，在本文中介质被近似成了2层均匀介质。因此为充分表示该层介质的衰减特性，尽可能地降低图1所示假设引起的误差，本文使用对应深度范围内衰减因子的平均值表示该层均匀介质的衰减特性。以图1中的介质I为例，其衰减特性为

\begin{matrix} a_{1} = \frac{1}{L_{1}} \int_{L_{2}}^{L_{1} + L_{2}} a (x) d x \\ b_{1} = \frac{1}{L_{1}} \int_{L_{2}}^{L_{1} + L_{2}} b (x) d x \end{matrix}

(7)

且 $a_{1} + b_{1} = c_{1}$ 。
同理，介质II的衰减特性也可用相同的积分方法得到。
1.3 散射相位函数
散射相位函数描述了光子与粒子碰撞后向不同角度散射的概率。由于Henyey-Greenstein散射函数（HG）能够与水中微粒的散射取得良好的近似，因而被广泛应用于UWOC领域。HG函数表示为

p_{H G} (μ) = \frac{1 - g^{2}}{{(1 + g^{2} - 2 g c o s (μ))}^{3 / 2}}

(8)

式中： $p_{H G} (μ)$ 是与散射角有关的散射相位函数；g 是影响前向散射和后向散射程度的因子； $μ$ 表示散射的方向，取值范围为[0,2 $π$ ]。文献^[16]中列举了因子 g、衰减系数 c 以及渐进辐射传递方程本征值 $γ$ 的参考数据。其中渐进辐射传递方程本征值 $γ$ 与衰减系数 c 的关系为

γ = \frac{k_{w}}{c}

(9)

式中： $k_{w}$ 是介质的扩散衰减系数，该值表示了距离随辐照度的变化。
2 蒙特卡洛仿真
MC方法是一种经典的概率方法。该方法能够模拟介质中每个光子运动，记录光子的运动轨迹，最后可通过统计学方法求得信道的冲激响应函数，并计算出信道传播损失。但从微观领域解释光通信时，MC方法存在一定缺陷，即不能解释一部分波动现象，例如全内反射等几何光学理论。尽管近年来已有若干关于光子在折射面处运动的研究，但相关理论仍存在一定争议。
在本文中MC方法以MIE散射理论作为基础。将光子的散射描述为光子与介质中粒子的弹性碰撞。由于光子的碰撞是弹性的，不会改变光子的能量，因此本文假设光束的波长保持恒定。MC方法将光子的吸收表示为光子以一定概率被粒子“截停”。此外，由于偏振光在具有散射性质的介质中传播时会失去其偏振态，因此本文忽视了光束传播中的偏振现象，并假定介质中的光是完全非偏振的。
文献^[19]使用MC方法对完全对齐的横向LOS-UWOC进行了仿真。UWOC的MC仿真可被分为3个模块：光子的初始状态、光子的传播、光子的接收。
2.1 光子的初始状态
如图1所示，在发射平面处建立三维笛卡尔坐标系。x，y 轴均位于水平的初始发射面上；z 轴垂直于初始发射面，指向海面。当UWOC中发射机采用高斯光源时，在初始发射面上，光子的初位置分布是以镜头中心为中点的高斯分布。光子初位置距中心的距离 $r_{0}$ 为^[16]

r_{0} = w_{s} \sqrt{- l n (1 - ς)}

(10)

式中： $w_{s}$ 是光束的束腰； $ς$ 是服从[0，1]之间均匀分布分随机变量。光子初始方向的方位角 $φ_{0}$ 服从[0，2 $π$ ] 的均匀分布。光子初始运动方向的天顶角为

θ_{0} = \frac{r_{0}}{f_{s}}

(11)

式中： $f_{s}$ 是发射机透镜的焦距， $f_{s} = w_{s} / t a n (θ_{s})$ ； $θ_{s}$ 是光束的扩散角的半角。
考虑UWOC中收发机非对齐的情况。初始发射面的法向量平行于 z 轴正向，因此只需将初始发射面转动天顶角 $θ_{t}$ 以及方位角 $φ_{t}$ ，即可使实际发射面的法向量指向空间中任意位置。以上转动过程相当于将发射机先绕 z 轴转动 $φ_{t}$ ，再绕 y 轴为转动 $θ_{t}$ 。
2.2 光子的传播
在文献^[16]中，光子的传播流程已被讨论。光子新的运动方向可以通过体积散射函数描述。

ζ = \int_{0}^{θ^{'}} \tilde{β} (θ) s i n θ d θ

(12)

式中： $ζ$ 是服从[0，1]均匀分布的随机数； $θ^{'}$ 是散射的最大散射角。这个方程依赖于体积散射函数的积分来形成一个累积分布函数（CDF）。由于海水的体积散射函数没有解析表达式，因此只能对方程进行数值计算。针对这一问题，Haltrin提出了几种近似于模拟天然水介质中散射函数的经验解。
光子在与介质粒子碰撞前的光路长度 $l$ ，可以通过累积分布函数求得。

P_{l} (l) = 1 - e^{- l}

(13)

式中， $P_{l} (l)$ 是[0，1]均匀分布的随机数，描述了光子在2次散射事件间中继续运动光路 $l$ 的概率；光路长度 $l$ 是衰减系数与光子单次实际运动距离的乘积 $l = c r$ 。由以上关系求得光子单次实际运动距离 r 为

r = - \frac{1}{c} l o g P_{l} (l)

(14)

上文提到，吸收主要表现形式为光子为粒子截停。然而若对每个光子的每一步运动均判决是否被吸收，则会极大地增加MC仿真的计算量。针对这一问题，COX提出了光子权重的概念，COX将MC仿真中的光子视作是初始值为1的光子团，在光子团与介质例子发生碰撞后，光子团中有 a/c的光子被吸收，有 b/c。因此在发生碰撞后，剩余光子团将乘以一个权值，以表示光子团中剩余光子数，该权值被称为反照度 $ω = b / c$ ^[16]。
基于MIE散射理论，当光子在介质中传播时，在未发生碰撞时，光子仍保持匀速直线运动。因此如图2所示，尽管光子已经穿越了2种介质宏观上的分界面，但由于光子在穿越界面过程中未与介质中的粒子发生碰撞，因此光子将保持直线运动，直至其与介质II中的粒子再度发生碰撞。
图2 光子在分界面处的传播
Fig.2 Propagation of photons at interface
由式（14）可知，光子在2次散射之间的随机运动距离与介质衰减系数有关。当衰减系数c增加时，光子运动距离的期望随之降低。如图2所示，一个光子在介质I内A点处与粒子发生碰撞，算得直至下次碰撞前的直线运动距离为 $r_{1} = - l o g (P_{l}^{1}) / c_{1}$ ，其中 $P_{l}^{1}$ 是[0，1]之间的常数。在经过 $r_{1}$ 距离后，光子将运动到介质II中的B位置。然而，由于光子在2种介质中发生碰撞的概率不同，光子可能运动到C位置即发生碰撞。为补偿光子穿介质后运动路径的变化，本文将MC过程中光子穿界面时的随机运动距离设置为

r_{0} = κ r_{1} + (1 - κ) r_{1} \frac{c_{1}}{c_{2}}

(15)

式中， $κ = d_{1} / (d_{1} + d_{2})$ ， $d_{1}$ 与 $d_{2}$ 分别是A点与B点到分界面的距离，该值可在MC仿真过程中对光子位置进行追踪得到。
2.3 光子的接收
MC仿真跟踪了每一个光子的运动状态，能够记录光子的路径、运动时间、入射位置、入射角度、权重等数据。当光子到达接收面后，接收机会根据光子的运动状态进行2项判决，分别是光子入射角时候小于接收机视场角，以及光子入射位置是否在接收机上。

F O V / 2 ⩾ θ_{i}

(16)

D_{R} ⩾ \sqrt{{(x^{e} - x^{i})}^{2} + {(y^{e} - y^{i})}^{2} + {(z^{e} - z^{i})}^{2}}

(17)

式中： $θ_{i}$ 是光子入射方向与接受面法线的夹角。 FOV是视场角的全角； D_R 表示接收机镜头的半径； $(x^{i}, y^{i}, z^{i})$ 是光子到达接收平面的位置， $(x^{e}, y^{e}, z^{e})$ 是接收机镜头中心的位置。文献^[15]详细介绍了MC仿真中信道冲激响应的获取方法。
3 仿真结果与讨论
本节展示了MC仿真的结果，通过统计学方法计算了UWOC链路的传播损失。表1展示了MC仿真涉及的部分参数。
图3 展示了在S2以及S5水体中，散射、吸收、衰减系数以及反照度随深度的变化。由图可知，在2种水体中，随着深度的增加，衰减系数总体上都会逐渐降低。但在某一深度均会出现衰减系数的峰值。例如在S2水体中，95m处衰减系数达到峰值。而在S5水体中，峰值位于45m处。这是因为在HALTRIN提出的衰减系数经验公式中，衰减正比于水体中的叶绿素含量。因此在叶绿素峰值附近，衰减系数最大。
表1 MC仿真的参数
Table1 Parameters of MC simulation
图3 叶绿素浓度呈S2与S5分布时，衰减系数随深度的变化
Fig.3 When concentration of chlorophyll is S2and S5， attenuation coefficient changes with depth
根据图3所示的衰减系数分布对介质进行分层。在S2分布的水体中，令发射机位于水下100m深度，接收机被吊放在40m深度，收发装置完全对齐。在图4展示了通过MC方法求得的UWOC冲激响应数据。在图4所示S2介质中。均假设介质I距离 L=20m，介质II距离 L=40m。每层介质衰减系数的平均值由式（7）计算得到。考虑假设介质I与介质II中溶解物成分相同以及成分不同的情况。分别使2层介质的HG因子均为0.93以及0.93&0.8。图4S5水体中的额UWOC中 L=40m，L=20m。由图可知，当介质中的微粒的前向散射减弱，后向散射增强时，UWOC信道的冲激响应将出现时间扩展。
图4 叶绿素浓度呈S2与S5分布时，双层UWOC信道的冲激响应
Fig.4 Impulse response of two-layer UWOC channel with S2and S5distribution of chlorophyll concentration
图5 描述了在穿越2层介质后，接收机接收到光信号强度随发射机天顶角偏转角度的变化。由于UWOC采用了高斯光源，该光源围绕其轴线呈中心对称。因此在讨论发射机偏转对接收功率的影响时，可以忽略位角的影响。图5中的4条曲线分别采用与图4相同的信道参数。由图可知，在信道质量较好的链路中，随着发射机出现偏转，接收信号的光功率会显著降低。而在信道质量较差的信道中，接收光功率的衰减相对平缓。图5中，当方位角偏转到18时，4条信道的接受光功率几乎重合。对于散射作用较为强烈的信道，由于抵达接收面的光子中散射光子所占的比例更高，一定程度上放宽了该信道对收发机对齐的要求。
图5 不同水体条件下，接收光强度随发射机天顶角的变化
Fig.5 Received light intensity varies with zenith Angle of transmitter under different water conditions
4 结束语
本文首次提出了纵向UWOC链路的MC仿真方法，通过该方法可获取多层纵向UWOC链路的信道冲击响应曲线以及信道衰减数据。此外本文还考虑了发射机存在角度偏转情况下UWOC链路的信道衰减。本文研究成果对于UWOC链路设计具有重要意义，能够有效推动纵向UWOC系统物理实现。
参考文献
- [1] 殷敬伟,惠俊英,郭龙祥.点对点移动水声通信技术研究[J].物理学报,2008(3):1753-1758.
- [2] YANG M,LI X K,MENG X X.Characteristic analysis of underwater acoustic scattering echoes in the wavelet transform domain[J].船舶与海洋工程学报(英文版)2017(1):93-101.
- [3] ANOUS N,ABDALLAH M,UYSAL M,et al.Performance evaluation of LOS and NLOS vertical inhomogeneous links in underwater visible light communications[J].IEEE Access,2018,6:22408-22420.
- [4] ELAMASSIE M,UYSAL M.Performance characterization of vertical underwater VLC links in the presence of turbulence[C]//2018 11^th International Symposium on Communication Systems,Networks & Digital Signal Processing(CSNDSP).Budapest:IEEE,2018.
- [5] ELAMASSIE M,SAIT S M,UYSAL M.Underwater visible light communications in cascaded gammagamma turbulence[C]//2018 IEEE Globecom Workshops(GC Wkshps).Abu Dhabi:IEEE,2018.
- [6] LOU Y,CHENG J L,NIE D H,et al.Performance of vertical underwater wireless optical communications with cascaded layered modeling[J].IEEE Transactions on Vehicular Technology,2022,71(5):5651-5655.
- [7] YILMAZ A,ELAMASSIE M,UYSAL M.Diversity gain analysis of underwater vertical MIMO VLC links in the presence of turbulence[C]//2019 IEEE International Black Sea Conference on Communications and Networking(BlackSeaCom).Sochi:IEEE,2019.
- [8] CHRISTOPOULOU C,SANDALIDIS H G,VAVOULAS A.Optimal trade-off between depth and coverage for a vertical UWOC link[J].IEEE Photonics Technology Letters,2020,32(4):208-211.
- [9] CHRISTOPOULOU C,SANDALIDIS H G,VAIOPOULOS N.Performance of an underwater optical wireless link with a randomly placed or moving receiver[J].IEEE Journal of Oceanic Engineering,2020,46(3):1068-1079.
- [10] SAIT M,GUO Y,ALKHAZRAGI O,et al.The impact of vertical salinity gradient on non-line-of-sight underwater optical wireless communication[J].IEEE Photonics Journal,2021,13(6):1-9.
- [11] JI X Y,YIN H X,JING L Y,et al.Modeling and performance analysis of oblique underwater optical communication links considering turbulence effects based on seawater depth layering[J].Optics Express,2022,30(11):18874-18888.
- [12] HALTRIN V I.Chlorophyll-based model of seawater optical properties[J].Applied Optics,1999,38(33):6826-6832.
- [13] YAP Y X,JASMAN F,MARCUS T C E.Impact of Chlorophyll concentration on underwater optical wireless communications[C]//2018 7^th International Conference on Computer and Communication Engineering(ICCCE).Kuala Lumper:IEEE,2018.
- [14] JOHNSON L J,GREEN R J,LEESON M S.Underwater optical wireless communications:depth dependent variations in attenuation[J].Applied Optics,2013,52(33):7867-7873.
- [15] HALTRIN V I.Theoretical and empirical phase functions for Monte Carlo calculations of light scattering in seawater[C]//Fourth International Conference on Remote Sensing for Marine and Coastal Environments.Orlando:Environmental Research Institute of Michigan,1997.
- [16] COX W C.Simulation,modeling,and design of underwater optical communication systems[M].Los Angeles:North Carolina State University,2012.

基本信息

中图分类号: TN929.1
文献标识码: A
DOI: 10.19838/j.issn.2096-5753.2022.04.007
文章编号: 2096-5753(2022)04-0322-07

基金信息

国家自然科学基金青年基金资助项目“面向跨海面信息传输的两跳纵向水下无线光通信关键技术研究”(62101152)；

引用信息

孙若凡，刘学蕊，娄毅，等. 分层纵向水下无线光信道建模研究[J]. 数字海洋与水下攻防，2022，5（4）：322-328.

稿件历史

收稿日期: 2022-06-14

参考文献

[1] 殷敬伟,惠俊英,郭龙祥.点对点移动水声通信技术研究[J].物理学报,2008(3):1753-1758.
[2] YANG M,LI X K,MENG X X.Characteristic analysis of underwater acoustic scattering echoes in the wavelet transform domain[J].船舶与海洋工程学报(英文版)2017(1):93-101.
[3] ANOUS N,ABDALLAH M,UYSAL M,et al.Performance evaluation of LOS and NLOS vertical inhomogeneous links in underwater visible light communications[J].IEEE Access,2018,6:22408-22420.
[4] ELAMASSIE M,UYSAL M.Performance characterization of vertical underwater VLC links in the presence of turbulence[C]//2018 11^th International Symposium on Communication Systems,Networks & Digital Signal Processing(CSNDSP).Budapest:IEEE,2018.
[5] ELAMASSIE M,SAIT S M,UYSAL M.Underwater visible light communications in cascaded gammagamma turbulence[C]//2018 IEEE Globecom Workshops(GC Wkshps).Abu Dhabi:IEEE,2018.
[6] LOU Y,CHENG J L,NIE D H,et al.Performance of vertical underwater wireless optical communications with cascaded layered modeling[J].IEEE Transactions on Vehicular Technology,2022,71(5):5651-5655.
[7] YILMAZ A,ELAMASSIE M,UYSAL M.Diversity gain analysis of underwater vertical MIMO VLC links in the presence of turbulence[C]//2019 IEEE International Black Sea Conference on Communications and Networking(BlackSeaCom).Sochi:IEEE,2019.
[8] CHRISTOPOULOU C,SANDALIDIS H G,VAVOULAS A.Optimal trade-off between depth and coverage for a vertical UWOC link[J].IEEE Photonics Technology Letters,2020,32(4):208-211.
[9] CHRISTOPOULOU C,SANDALIDIS H G,VAIOPOULOS N.Performance of an underwater optical wireless link with a randomly placed or moving receiver[J].IEEE Journal of Oceanic Engineering,2020,46(3):1068-1079.
[10] SAIT M,GUO Y,ALKHAZRAGI O,et al.The impact of vertical salinity gradient on non-line-of-sight underwater optical wireless communication[J].IEEE Photonics Journal,2021,13(6):1-9.
[11] JI X Y,YIN H X,JING L Y,et al.Modeling and performance analysis of oblique underwater optical communication links considering turbulence effects based on seawater depth layering[J].Optics Express,2022,30(11):18874-18888.
[12] HALTRIN V I.Chlorophyll-based model of seawater optical properties[J].Applied Optics,1999,38(33):6826-6832.
[13] YAP Y X,JASMAN F,MARCUS T C E.Impact of Chlorophyll concentration on underwater optical wireless communications[C]//2018 7^th International Conference on Computer and Communication Engineering(ICCCE).Kuala Lumper:IEEE,2018.
[14] JOHNSON L J,GREEN R J,LEESON M S.Underwater optical wireless communications:depth dependent variations in attenuation[J].Applied Optics,2013,52(33):7867-7873.
[15] HALTRIN V I.Theoretical and empirical phase functions for Monte Carlo calculations of light scattering in seawater[C]//Fourth International Conference on Remote Sensing for Marine and Coastal Environments.Orlando:Environmental Research Institute of Michigan,1997.
[16] COX W C.Simulation,modeling,and design of underwater optical communication systems[M].Los Angeles:North Carolina State University,2012.

首页

本刊简介

编委会

收录情况

期刊订阅

征稿简则

联系我们

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

分层纵向水下无线光信道建模研究

Modeling of Hierarchical Vertical Underwater Wireless Optical Channel

摘要

Abstract

关键词

Keywords

0 引言

1 系统和信道模型

1.1 系统模型

1.2 吸收与散射模型

1.2.1 吸收系数

(1)

(2)

(3)

1.2.2 散射系数

(4)

(5)

(6)

1.2.3 衰减系数

(7)

1.3 散射相位函数

(8)

(9)

2 蒙特卡洛仿真

2.1 光子的初始状态

(10)

(11)

2.2 光子的传播

(12)

(13)

(14)

(15)

2.3 光子的接收

(16)

(17)

3 仿真结果与讨论

4 结束语

参考文献

基本信息

基金信息

引用信息

稿件历史

参考文献