CN102638308A - 一种水下光传输特性的分析方法 - Google Patents
一种水下光传输特性的分析方法 Download PDFInfo
- Publication number
- CN102638308A CN102638308A CN2012100415091A CN201210041509A CN102638308A CN 102638308 A CN102638308 A CN 102638308A CN 2012100415091 A CN2012100415091 A CN 2012100415091A CN 201210041509 A CN201210041509 A CN 201210041509A CN 102638308 A CN102638308 A CN 102638308A
- Authority
- CN
- China
- Prior art keywords
- light
- layer
- water body
- attenuation
- prime
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种水下光传输特性的分析方法,其特征在于包括:1:对水体沿深度方向进行分层,获取不同层相应的吸收系数值和散射系数值,根据每一层的吸收系数值和散射系数值计算得到相应层的漫射衰减系数;2:根据每一层的漫射衰减系数,计算每一层的透过率,将每一层的透过率相乘,得到水体通信信道的整体透过率,之后采用截断误差对该整体透过率进行修正,得到修正后的整体透过率;3:根据水体的光衰减情况,输出水体通信信道传输特性。本发明的水下光传输特性的分析方法利用了绿叶素浓度的垂直分布模型实现对光在水中通信信道的传输特性的分析,相对于现有的衰减系数模型,由于无需通过实验测出不同深度下的衰减系数,方法简便,易于应用。
Description
技术领域
本发明属于激光通信领域,特别是一种水下光传输特性的分析方法。
背景技术
地球上百分之七十的面积都被海水覆盖着,在海上活跃着大量的民用船只和进行海战的舰船、潜艇等。随着现代高新技术的发展,水下捕捞、探测、控制等的需要,逐渐形成了水下通信这一个特殊的应用领域。传统的通信手段,如无线电波通信、水声通信等已经无法满足水下通信的需求,成了限制水下通信发展的技术瓶颈。与传统的通信技术相比,光学通信技术可以克服传统通信的带宽窄、受环境影响大、可适用的载波频率低、传输的时延大等缺陷。例如,空中飞行器可以利用激光作为载波工具而向水下一定深度处的潜艇发送通信数据。
然而由于光在水中的衰减很严重,因而有必要预先对光在水中传输的通信信号进行分析。目前,关于光在海水介质中传输的研究一般有三种方法:分析方法、实验方法和数值模拟方法,其中的分析方法又包括两种模型,一种是基于生物光学的衰减模型,一种是衰减系数模型。基于生物光学的衰减模型主要是通过分析海水中黄色物质、叶绿素、悬浮颗粒等物质对光吸收和散射的影响来研究光在水下的衰减特性,其只能用于光在水平方向上的传输,而当光垂直传输时这种方法就不再适用了;衰减系数模型主要是通过海水的衰减系数来研究光在水下的衰减特性,其虽然是基于实测得来数据进行分析,比较有说服力,但是当光垂直传输时就必须通过实验测出不同深度下的衰减系数,如果要求误差越小,那么所测的值就会越多,分析方法也就越复杂,应用不方便。
发明内容
针对现有水下光传输特性的两种分析方法分别存在的上述问题,本发明提出了一种水下光传输特性的分析方法。本发明采用的技术手段如下:
一种水下光传输特性的分析方法,其特征在于包括:
步骤1:对水体沿深度方向进行分层,获取不同层相应的吸收系数值和散射系数值,根据每一层的吸收系数值和散射系数值计算得到相应层的漫射衰减系数;
步骤2:根据每一层的漫射衰减系数,计算每一层的透过率,将每一层的透过率相乘,得到水体通信信道的整体透过率,之后采用截断误差对该整体透过率进行修正,得到修正后的整体透过率,该修正后的整体透过率表征了水体的光衰减情况;其中的截断误差表示为:
步骤3:根据水体的光衰减情况,输出水体通信信道的传输特性。
本发明的水下光传输特性的分析方法利用了绿叶素浓度的垂直分布模型实现对光在水中通信信道的传输特性的分析,相对于现有的衰减系数模型,由于无需通过实验测出不同深度下的衰减系数,方法简便,易于应用。
附图说明
图1为本发明的截断误差E与层数j之间的关系图。
图2为叶绿素浓度在南海几个站点的实测图。
图3为本发明采用的叶绿素浓度垂直分布模型图。
图4为采用叶绿素浓度模型时海水深度与功率衰减的关系图。
图5为采用衰减系数模型时海水深度与功率衰减的关系图。
图6为叶绿素浓度模型与衰减系数模型的比较图。
图7为利用本发明的分析方法得到的理想天气条件下功率衰减与传输距离的关系曲线图。
图8为利用本发明的分析方法得到的理想天气条件下误码率与传输距离的关系曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。
本发明的水下光传输特性的分析方法包括:
步骤1:对水体沿深度方向进行分层,获取不同层相应的吸收系数值和散射系数值,并根据每一层的吸收系数值和散射系数值计算得到相应层的漫射衰减系数。
由于光信号在水体中的衰减主要是由光吸收和光散射引起的,光吸收主要包括纯水体对光的吸收、溶解有机物对光的吸收、浮游植物对光的吸收、悬浮颗粒及碎屑对光的吸收;光散射主要包括纯水体对光的散射、浮游植物对光的散射、小颗粒对光的散射、大颗粒对光的散射。假设光吸收系数为A(λ),纯水体对光的吸收系数为Aw(λ),溶解有机物对光的吸收系数为ACDOM(λ),浮游植物对光的吸收系数为Ac(λ),悬浮颗粒及碎屑对光的吸收系数为Ad(λ),则有:
A(λ)=Aw(λ)+ACDOM(λ)+Ac(λ)+Ad(λ) (1)
又假设光散射系数为B(λ),纯水体对光的散射系数为Bw(λ)、浮游植物对光的散射系数为Bc(λ)、小颗粒对光的散射系数为Bs(λ)、大颗粒对光的散射系数为Bl(λ),则有:
B(λ)=Bw(λ)+Bc(λ)+CsBs(λ)+ClBl(λ) (2)
又假设光信号在水体中的衰减系数是K(λ),则有:
K(λ)=A(λ)+B(λ)
由于水体散射作用的存在,描述准直光束所形成的水下光场的功率衰减不能简单使用光束衰减系数K(λ)表示,为此,本发明采用光信号在水体中的漫射衰减系数k,假设μ=A(λ)+B(λ),则漫射衰减系数k可表示为:
步骤2:根据每一层的漫射衰减系数,计算每一层的透过率,将每一层的透过率相乘,得到水体通信信道的整体透过率,之后采用截断误差对该整体透过率进行修正,得到修正后的整体透过率,该修正后的整体透过率表征了水体的光衰减情况。
假设将水体平均分为j层,每一层的叶绿素浓度不变,z为光信号的水下垂直传输距离,则每一层光透过率Lw表示为:
Lw=exp(-kz) (4)
而在实际情况下,水体中叶绿素浓度的分布式连续的,将其分层后会产生误差,且所分层数越多,误差越小,但所分层数越多,计算量越大,为了获取最佳分层数,下面通过截断误差来分析所分层数与误差之间的关系:
复合梯形公式截断误差E满足:
其中,a、b分别为水体深度的边界值,f″(η)为叶绿素高斯分布的表达式的二阶导数,如图1示出了误差E与层数j之间的关系。由图1可以看出,在所分层数j达到800层之前,误差E随着层数j的增加而急剧的减小,其变化范围跨度了3个数量级,而当所分层数在800至2000层之间变化时,误差E的减小逐渐趋于平缓,稳定在一个数量级以内。所以,在误差允许的范围内考虑计算量,那么,水深200米时,所分的最佳层数应该为800层左右。
步骤3:根据水体的光衰减情况,输出水体通信信道的传输特性,该传输特性包括误码率、传输速率、传输距离等。
当上述方法应用于对机载激光探测系统的水通信链路传输特性分析时,在步骤3之后,还可以包括:获取机载激光发射机发出的激光束在到达水面时的激光发射功率;根据激光发射功率及整体透过率,得出光在水下传输时功率衰减的变化情况;根据光功率的衰减、接收机的灵敏度,计算光在水下传输的距离。
其中,根据激光发射功率及整体透过率,得出光在水下传输时功率衰减的变化情况的步骤可表示为:
其中,根据光功率的衰减、接收机的灵敏度,计算光在水下传输的距离的步骤可表示为:
式中,n表示光功率刚好能被接收机接收时所达到的层数,z0表示光所传输的距离。
利用式(7)和式(8),还可以在已知激光接收功率Pr的前提下,逆向求得激光发射功率Pt,以为合理的设置机载激光发射机发射功率提供参考。
为了验证上述采用叶绿素浓度模型的分析方法的可行性,下面采取相同海域的参数对衰减系数模型与叶绿素浓度模型进行了比较分析:
对于叶绿素浓度模型,假设叶绿素浓度的实测值如图2所示,某一深度下的叶绿素浓度值可表示为:
其中,c0是背景生物量(mg/m3);zmax是最大叶绿素浓度对应的深度(m);σ是标准偏差(m);h为背景生物量以上的总生物量(mg/m2),表示基线以上的峰值高度。Cc(z)是深度z处的叶绿素浓度值。由式(9)和图2中06年的叶绿素浓度的实测值,可得到如图3所示的叶绿素浓度垂直分布图。若选取灵敏度为-59dBm的接收机,发射功率Pt=106W,仿真结果如图4所示。
对于衰减系数模型,假设系统的有效衰减系数为Γ,系统参数为Pmax,最小可检测到的光场功率为Psmin,水体最大探测深度Lmax,则有:
Lmax=ln(Pmax/Psmin)/(2Γ) (10)
如下表选择的不同深度的海水漫射衰减系数,其与上面选取的不同深度海水的叶绿素浓度来自同一监测站的实测数据:
深度/m | 5 | 10 | 20 | 30 | 50 | 100 | 200 | 240 |
衰减系数/m-1 | 0.111 | 0.113 | 0.116 | 0.114 | 0.111 | 0.078 | 0.052 | 0.051 |
根据不同深度海水的衰减系数和公式(10),可以得到不同深度海水的功率衰减,若选取灵敏度为-59dBm的接收机,发射功率Pt=106W,仿真结果如图5所示。
图6为两种模型对比的仿真图,可以看出,随着激光在海水中传输距离的增加,接收功率呈指数形式衰减。当选用灵敏度为-59dBm的接收机,发射功率Pt=106W时,可以看出,水下的传输距离为200m左右,两种模型得到的功率衰减特性基本相同,但是用叶绿素浓度来分析明显要更加的方便,适用性也更加的广泛。
下面通过一仿真实例进一步说明上述方法,该仿真实例用以实现功率衰减及误码率仿真分析:
仿真条件:定义一波长为λ的光波,其传输通过距离为z的过程中所产生的空间损耗Lp满足:
考虑纯大气情况,卫星轨道高610km,此时空间链路的功率损耗可由公式(11)计算得出,采用波长为488nm的激光进行通信,那么通过大气空间链路的功率损耗为263.7744dB。设采用的激光器的发射功率为36dBm,发射天线增益为109dB,接受器的灵敏度为-59dBm,接收机天线增益为135dB,大气背景噪声取经验值2.17dBm,大气海水界面的透过率取0.83dB,通过这些参数,我们就可以得到理想天气条件下的功率衰减,误码率的变化规律。
通过上面的参数,我们能够计算出激光到达海面时的功率大约为13.17dBm,此时采用的叶绿素浓度模型如图3所示,根据上述方法进行仿真分析后的结果如图7所示。
同样,通过叶绿素浓度高斯分布模型我们还可以分析当叶绿素浓度的分布如图3所示时,不同调制的误码率的变化规律。我们分别仿真分析了OOK、2DPSK、4PPM等几种常用调制的误码率,仿真结果如图8所示。其误码率的计算公式如下:
OOK调制的误码率为:
2DPSK调制的误码率为:
4PPM调制的误码率为:
式中,S/N表示信噪比,k表示阈值门限,一般取0.5。
本发明的水下光传输特性的分析方法利用了绿叶素浓度的垂直分布模型实现对光在水中通信信道的传输特性的分析,相对于现有的衰减系数模型,由于无需通过实验测出不同深度下的衰减系数,方法简便,易于应用。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (6)
2.根据权利要求1所述的方法,其特征在于计算得到相应层的漫射衰减系数的步骤表示为:
其中,k为漫射衰减系数,Aw(λ)为光吸收系数,B(λ)为光散射系数,μ=A(λ)+B(λ)。
3.根据权利要求2所述的方法,其特征在于计算每一层的透过率的步骤表示为:
Lw=exp(-kz)
其中,z为光信号的水下垂直传输距离,Lw为每一层光透过率。
5.根据权利要求4所述的方法,其特征在于步骤3之后还包括:
获取机载激光发射机发出的激光束在到达水面时的激光发射功率;
根据激光发射功率以及整体透过率,得出光在水下传输时功率衰减的变化情况,其步骤表示为:
其中,Pt为激光发射功率,Pr为激光接收功率,为光在水下的透射率。
6.据权利要求5所述的方法,其特征在于根据光功率的衰减、接收机的灵敏度,计算光在水下传输的距离,其步骤表示为:
式中,n表示光功率刚好能被接收机接收时所达到的层数,z0表示光所传输的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100415091A CN102638308A (zh) | 2012-02-22 | 2012-02-22 | 一种水下光传输特性的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100415091A CN102638308A (zh) | 2012-02-22 | 2012-02-22 | 一种水下光传输特性的分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102638308A true CN102638308A (zh) | 2012-08-15 |
Family
ID=46622565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100415091A Pending CN102638308A (zh) | 2012-02-22 | 2012-02-22 | 一种水下光传输特性的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102638308A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108429591A (zh) * | 2017-11-13 | 2018-08-21 | 西北工业大学 | 一种适用于深海信道的多载波水声通信方法 |
CN109238993A (zh) * | 2018-11-28 | 2019-01-18 | 南昌航空大学 | 海水中叶绿素含量对于水下光传播特性影响的探测方法 |
CN110095784A (zh) * | 2019-05-09 | 2019-08-06 | 北京航空航天大学 | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 |
CN110166120A (zh) * | 2019-04-08 | 2019-08-23 | 中山大学 | 适用于多种水质的智能双向可见光通信系统 |
CN110650054A (zh) * | 2019-09-29 | 2020-01-03 | 华为技术有限公司 | 配置调制器光工作点的方法、装置和系统 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN113746544A (zh) * | 2021-09-09 | 2021-12-03 | 桂林电子科技大学 | 一种非均匀海水的垂直信道链路分层方法 |
CN113810137A (zh) * | 2021-08-23 | 2021-12-17 | 华北电力大学(保定) | 一种基于垂直叶绿素浓度分布的水下光子传输仿真方法 |
-
2012
- 2012-02-22 CN CN2012100415091A patent/CN102638308A/zh active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108429591A (zh) * | 2017-11-13 | 2018-08-21 | 西北工业大学 | 一种适用于深海信道的多载波水声通信方法 |
CN108429591B (zh) * | 2017-11-13 | 2021-05-18 | 西北工业大学 | 一种适用于深海信道的多载波水声通信方法 |
CN109238993A (zh) * | 2018-11-28 | 2019-01-18 | 南昌航空大学 | 海水中叶绿素含量对于水下光传播特性影响的探测方法 |
CN110166120A (zh) * | 2019-04-08 | 2019-08-23 | 中山大学 | 适用于多种水质的智能双向可见光通信系统 |
CN110095784A (zh) * | 2019-05-09 | 2019-08-06 | 北京航空航天大学 | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 |
CN110650054A (zh) * | 2019-09-29 | 2020-01-03 | 华为技术有限公司 | 配置调制器光工作点的方法、装置和系统 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN111965608B (zh) * | 2020-07-16 | 2024-01-12 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN113810137A (zh) * | 2021-08-23 | 2021-12-17 | 华北电力大学(保定) | 一种基于垂直叶绿素浓度分布的水下光子传输仿真方法 |
CN113810137B (zh) * | 2021-08-23 | 2023-10-17 | 华北电力大学(保定) | 一种基于垂直叶绿素浓度分布的水下光子传输仿真方法 |
CN113746544A (zh) * | 2021-09-09 | 2021-12-03 | 桂林电子科技大学 | 一种非均匀海水的垂直信道链路分层方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102638308A (zh) | 一种水下光传输特性的分析方法 | |
Salam et al. | Di-Sense: In situ real-time permittivity estimation and soil moisture sensing using wireless underground communications | |
CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
Guenther et al. | Meeting the accuracy challenge in airborne lidar bathymetry | |
Greinert et al. | Hydroacoustic experiments to establish a method for the determination of methane bubble fluxes at cold seeps | |
CN106021864B (zh) | 星载散射计后向散射系数的检验方法和装置 | |
CN104360329B (zh) | 一种全数字阵列相控阵天气雷达的强度定标方法 | |
CN102012524B (zh) | 一种海上三维地震观测系统羽状漂移定量评估方法 | |
Suberg et al. | Assessing the potential of autonomous submarine gliders for ecosystem monitoring across multiple trophic levels (plankton to cetaceans) and pollutants in shallow shelf seas | |
CN103076290A (zh) | 确定云的液态水含量的方法 | |
Matolak et al. | Initial results for air-ground channel measurements & modeling for unmanned aircraft systems: Over-sea | |
CN106526611A (zh) | 一种利用激光雷达系统测量海洋固有光学参数的方法 | |
CN107976686A (zh) | 一种多视场角海洋激光雷达及其视场角择优方法 | |
KR20130044551A (ko) | 해상풍 탐지시스템 및 이를 이용한 해상풍 탐지방법 | |
Chaichitehrani et al. | Colored dissolved organic matter dynamics in the northern Gulf of Mexico from ocean color and numerical model results | |
CN106442233A (zh) | 利用卫星数据估测海水悬浮颗粒物粒径的混合型新方法 | |
CN103278804A (zh) | 一种波导超视距雷达 | |
Jou-Claus et al. | Applicability of Landsat 8 thermal infrared sensor for identifying submarine groundwater discharge springs in the Mediterranean Sea basin | |
CN106556559A (zh) | 一种云/气溶胶的监测方法、装置和系统 | |
Mason | An integrated assessment of oil and gas release into the marine environment at the former Taylor Energy MC20 site | |
Belov et al. | Atmospheric bistatic communication channels with scattering. Part 2. Field experiments in 2013 | |
CN114720426A (zh) | 新的星载gnss反射信号的溢油检测方法 | |
Chen et al. | Review of airborne oceanic lidar remote sensing | |
Morison et al. | ICESat-2 ATL12 ocean surface height and ATL19 gridded dynamic ocean topography | |
Saputra et al. | Effect of Turbidity, Temperature and Salinity of Waters on Depth Data from Airborne LiDAR Bathymetry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C05 | Deemed withdrawal (patent law before 1993) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120815 |