-
液体中的激光致声,即在液体中激光能量达到一定强度激发产生声波的过程。根据激光相互作用区域的能量密度与时空分布,可把液体中光声效应机制主要分为热膨胀、汽化与光击穿3种。与汽化机制和光击穿机制相比,热膨胀机制物理过程简单,实验上容易实现。相对于另外两种致声机制,热膨胀机制致声原理的研究相对比较成熟,本文中关于液体中激光致声空-水跨介质研究将针对热膨胀机制进行。
在热膨胀机制下,假设液体为均匀的理想流体,激光激发水下声波的转换方程[21]:
$ \nabla^2 s-\frac{1}{c^2} \frac{\partial^2 s}{\partial t^2}=-\frac{\beta}{c_p} \frac{\partial \boldsymbol{H}(x, y, z, t)}{\partial t} $
(1) 其中:
$ s=\left(\alpha c^2 \beta E_0 / c_p\right) \cdot \exp (-\alpha z) $
(2) $ \beta=-\frac{1}{\rho_0}\left(\frac{\partial \rho_0}{\partial T}\right) $
(3) 式中:s为声压; c为声速; x、y、z为空间坐标方向距离; t为时间; α为水的吸收系数; E0为表面激光能量密度;cp为单位质量定压比热容;H (x, y, z, t)为单位时间内液体吸收并进行能量转化的电磁能量密度;β为液体的体积热膨胀系数(ρ0为密度,T为温度,假设T在过程中不变,β为常数)。由参考文献[16],将式(1)傅里叶变换取实部,得到其激光垂直入射表达式为:
$ s(r, \omega)=-\frac{A I_0 \beta a^2}{2 c_p} E(\omega) I(\omega) $
(4) 其中,
$ A=1-\left[\frac{n_2 \cos i-n_1 \cos \phi}{n_2 \cos i+n_1 \cos \phi}\right]^2 $
(5) 式中:r为接收点与声源之间的距离; ω为频率; A为给定激光偏振情况下激光对液体的透射率; i为入射角; ϕ为折射角; n1为空气中的折射率; n2为液体中的折射率; I0为光强; a为光束半径; I(ω)为激光频谱; E(ω)为单位时间吸收并转化为热量的电磁波能量密度频谱,可表示为:
$ \left\{\begin{array}{l} E(\omega)=\frac{\exp (\mathrm{i} k r)}{r} \frac{\omega^2 \tau_0}{1+\omega^2 \tau_0^2} \exp \left(\frac{-\omega^2 \tau_1^2}{4}\right) \\ \tau_0=\cos \theta /(\mu c) \\ \tau_1=a \sin \theta / c \end{array}\right. $
(6) 式中:k为波数,与激光波长有关; μ为水对激光的吸收系数; θ为观测角,即r与垂直向下方向的夹角; τ0为垂直特性延迟时间; τ1为水平特性延迟时间。
-
n-MFSK调制是在多进制频移键控(n-frequency shift keying, n-FSK)调制基础上,进一步提高通信速率并改善误比特率。本文中针对激光致声空-水跨介质长脉冲法和重复频率法,提供n-MFSK调制方案,整体激光致声通信原理图如图 1所示。
-
长脉冲法是利用某种调制方法对长脉冲激光进行调制的方法,本文中根据参考文献[16]来选取热膨胀效应下的光声转换仿真参数,具体参数如表 1所示。使用Nd∶glass激光器,长脉冲激光波形近似为方形脉冲,在边缘逐渐下降,该波形可在极端聚焦情况下减少非线性效应波形突变的影响。
parameter value laser wavelength 1053 nm pulse width 2 ms peak power 25 kW coefficient of thermal expansion 2×10-4 ℃-1 laser beam radius 6 mm specific heat capacity of liquids 4.17×103 J·kg-1·K-1 speed of sound in liquids 1500 m/s optical absorption coefficient 15.7 m-1 optimal center frequency 10 kHz 由于水声信道存在多径效应,需选取合适的符号间隔以减轻符号间干扰,故本文中选取n-FSK调制,该调制系统由多台激光器构成激光阵列,通过控制激光器发射频率减少由多径引起的符号间干扰。以2-FSK为例,将其中一个调制频率记为1,另一个调制频率记为0,两个调制频率中心频率都为最佳中心频率。每台激光器,符号持续时间为2 ms,符号间隔为98 ms,通过机械斩波器控制激光频率变化,例如,传输速率为100 bits/s时,需由10台激光器组成阵列,以激光器发射频率并行发射。该方案可通过增加调制频率数量,以减少激光器数量并提升传输速率,但所需带宽随传输速率呈非线性增加。理论上可采用增加符号持续时间,减少所占带宽,为此需减少激光器发射功率,以增加脉冲持续时间,这会导致声波在水中的声压级(sound pressure level,SPL)和传输距离减少。此外,随着调制频率数量增加,会降低对多径效应的抑制,可通过信道编码减少多径效应,但会降低传输速率。
为此,本文作者提出了另一种调制方式,即n-MFSK调制,在一个符号间隔期间发射多个频率,以增加调制频率数量。以2-MFSK为例,将两台相同激光器基频设为f1和f2,其中,当两台激光器均不发射激光时,记为00;当激光器1发射频率为f1、激光器2不发射时,记为01;当激光器1不发射、激光器2发射频率为f2时,记为10;当激光器1发射频率为f1、激光器2发射频率为f2时,记为11。当传输速率达到100 bits/s时,需要5台激光器组成阵列。2-FSK和2-MFSK调制的时域和频域结果分别如图 2a和图 2b所示。
图 2 长脉冲法时域和频域仿真结果
Figure 2. Long-pulse-duration laser time domain and frequency domain simulation results
该方案的优点为通过激光器频率组合,增加每个符号传输比特数来提高传输速率,并且在使用开关键控时,可提高带宽利用率。
根据表 1中的参数,计算不同调制方法下长脉冲法水下通信距离为[17]:
$ \left\{\begin{array}{l} X_{\mathrm{s}}-X_{\mathrm{t}}-X_{\mathrm{n}}=R_{\mathrm{SNR}} \\ X_{\mathrm{t}}=20 \lg r+\alpha r \times 10^{-3} \end{array}\right. $
(7) 式中:Xs为声源级; Xt为声传播损失; Xn为背景噪声; RSNR为接收器检测信噪比; α为海水吸收系数。
表 2为长脉冲法n-FSK和n-MFSK调制方法通信速率、带宽、声压级和水中传输距离的比较。其中背景噪声取3级海况,RSNR=10 dB,α=0.02 dB/km,声压级作用范围以137 dB为基准(依据表 1计算)。随着调制频率数量增加,n-MFSK调制可提高通信速率,并改善频带利用率,但声压级和水中传输距离会随符号持续时间内调制频率数量增加而减小。比较4-FSK和2-MFSK调制,两者可达到相同的通信速率,且2-MFSK只需要1 kHz的带宽,然而在最坏情况下,当一个符号持续时间传输两个频率时,激光能量必须除以2,因此声压级减小到131 dB。这导致水中传输距离从4-FSK的708 m减小到355 m,随着一个符号时间内传输频率数量增加,水中传输距离将大幅度减少。因此,当采用n-MFSK调制时,需权衡通信速率与水中传输距离的关系,根据实际需求选择合适的调制方式。
表 2 长脉冲法n-FSK和n-MFSK调制的比较
Table 2. Comparison of n-FSK and n-MFSK modulation by long-pulse-duration laser method
modulation
methoddata rate/
(bit·s-1)bandwidth/
kHzSPL/
dBin-water
range/m2-FSK 10 1 137 708 4-FSK 20 2 137 708 8-FSK 30 4 137 708 16-FSK 40 8 137 708 2-MFSK 20 1 131 355 3-MFSK 30 1.5 128 251 4-MFSK 40 2 125 178 -
重复频率法是基于激光致声的热膨胀效应通过控制激光器的重复频率产生调制声信号的方法。由于单激光脉冲产生的声信号频谱所占带宽大,频带利用率低,其信号能量会在水下发射严重衰减,易产生畸变,不利于传输编码。因此, 需采用产生窄带且中心频率可控的声信号,即高重复率的激光脉冲信号。
当激光垂直射向水面时,激光脉冲重复率为fr=1/T,T为激光脉冲周期时,激光脉冲的频谱可表示为[20]:
$ \left\{\begin{array}{l} I(\omega)=I_0(\omega) G(\omega) \\ G(\omega)=\frac{\sin (N \omega T / 2)}{\sin (\omega T / 2)} \end{array}\right. $
(8) 式中:I0(ω)为单脉冲的频谱; G(ω)频谱由于谐波关系呈现梳状结构; N为阶数; fr为相邻谱峰的间距,且其周期可控。将式(8)代入式(4)可得[21]:
$ s(r, \omega)=-\frac{A I_0 \beta a^2}{2 c_p} E(\omega) I_0(\omega) G(\omega) $
(9) 式中:E(ω)与μ、a、θ等参数有关,其中激光器参数μ和a是固定的,故只需控制fr使s(r, ω)窄带特性最好即可。选择激光器发射重复率fr时,要使声信号高次谐波被抑制且中心频率处的能量强,因此,激光器的重复频率需满足[20]:
$ \left\{\begin{array}{l} 20 \lg \left|\frac{s\left(f_{\mathrm{r}}\right)}{s\left(f_{\max }\right)}\right| \leqslant 3 \mathrm{~dB} \\ 20 \lg \frac{\left|s\left(f_{\mathrm{r}}\right)\right|}{\max \left|s\left(n f_{\mathrm{r}}\right)\right|} \geqslant 3 \mathrm{~dB}, (n=1, 2, \cdots) \end{array}\right. $
(10) 式中:s(fr)为中心频率的幅度; s(fmax)为频谱幅度的最大值; max|s(nfr)|为高次谐波频谱幅度的最大值, 由此,可确定激光器发射重复率的范围。根据表 3中的参数,采用调Q Nd∶YAG高重频激光器,并采用灵敏度为-170 dB的光纤水听器对声波进行接收,联合式(8)~ 式(10)计算fr的范围可得:
parameter value laser wavelength 1060 nm pulse width 8 ns single pulse energy 100 mJ~800 mJ optical absorption coefficient 13.7 m-1 angle of observation from the vertical 1.31 rad $ 12 \mathrm{kHz} \leqslant f_{\mathrm{r}} \leqslant 31 \mathrm{kHz} $
(11) 目前现有的激光器难以从实验的角度进行高速通信研究,由于相邻激光脉冲产生的声信号具有很好的重复性,故可采用单个或多个实验信号进行重频组合,获得所需的高重复率激光。图 3a中将单脉冲激光致声信号以fr=21 kHz重频组合,图 3b是依据式(9)计算后得到结果和单脉冲重频组合信号频谱进行对比。两者频谱特征基本吻合,验证了该方案的可行性。
图 3 a—单脉冲激光信号时域波形b—重频组合信号时域波形c—I0(ω)G(ω)与单脉冲重频组合频域波形
Figure 3. a—single pulse laser signal time domain waveform b—re-frequency combination signal time domain waveform c—spectrogram with single laser signal re-frequency combination and I0(ω)G(ω)
通过控制激光器的重复频率使其有规律地变化,从而完成n-FSK和n-MFSK调制以搭载信息,n-FSK调制中一个码元信号以重复频率fr的N个激光脉冲激发声信号表示,在符号间隔(0, T0)中,码元信号可表示为:
$ s_i(t)=\operatorname{Re}\left[\operatorname{IFFT}\left(s\left(f_{\mathrm{r}, i}\right)\right)\right], (i=1, 2, \cdots, n) $
(12) 故n-FSK调制信号可表示为:
$ \begin{aligned} s(t)= & \sum\limits_{m=0}^{\infty} s_i(t) g_\tau\left(t-m T_0\right), \\ & (i=1, 2, \cdots, n) \end{aligned} $
(13) 式中:IFFT(s(fr, i))表示对s(fr, i)进行的傅里叶逆变换(inverse fast Fourier transform, IFFT);gτ(t- mT0)是脉宽为τ、周期为T0的门函数。
n-MFSK调制中每个码元含有零个、一个或多个调制频率,单个码元信号可表示为:
$ \begin{gathered} s_{\mathrm{MFSK}, i}(t)=\sum\limits_{i=1}^n \cos \left(\frac{n_i \pi}{2}\right) s_i(t-(i-1) \tau), \\ (i=1, 2, \cdots, n) \end{gathered} $
(14) 式中:ni取值为0或1,取0代表该码元含有调制频率fr, i,取1代表不含调制频率fr,i。
n-MFSK调制信号可表示为:
$ \begin{aligned} s_{\mathrm{MFSK}}(t)= & \sum\limits_{m=0}^{\infty} s_{\mathrm{MFSK}, i}(t) g_\tau\left(t-m T_0\right), \\ & (i=1, 2, \cdots, n) \end{aligned} $
(15) 通过FSK和MFSK调制信号搭载信息,实现激光信号的编码传输。图 4a和图 4b分别对应4-FSK和2-MFSK调制时域和频域波形,符号间隔为50 ms。
4-FSK调制频率fr,i根据式(12)重复频率范围等间隔取值,分别取为13 kHz、17 kHz、21 kHz和25 kHz,代表 4个码元,每个码元用二进制依次表示为{00,01,10,11},码元带宽约为500 Hz。2-MFSK调制频率fr,i分别取为17 kHz和25 kHz,第1个码元不包含频率信息为fr,0,第2个码元含有1个17 kHz频率信息为fr,1,第3个码元含有1个25 kHz频率信息为fr,2,第4个码元含有两个17 kHz和25 kHz频率信息为fr,1,2,每个码元用二进制依次表示为{00,01,10,11}。如图 3所示,2-MFSK调制可减少基频数量, 同时提高通信速率并提高带宽利用率。
表 4为重复频率法n-FSK和n-MFSK调制的比较。参数计算方式与长脉冲法相同,声压级作用范围以139 dB为基准(依据表 2计算)。对比n-FSK和n-MFSK调制,在相同调制频率数量时,n-MFSK调制可提高通信速率,并提高频带利用率,与此同时代价是水中通信距离的损失。
表 4 重复频率法n-FSK和n-MFSK调制的比较
Table 4. Comparison of high repetitive rate method n-FSK and n-MFSK modulation
modulation
methoddata rate/
(bit·s-1)bandwidth/
kHzSPL/
dBin-water
range/m2-FSK 20 1 139 897 4-FSK 40 2 139 897 8-FSK 60 4 139 897 16-FSK 80 8 139 897 2-MFSK 40 1 133 459 3-MFSK 60 1.5 130 316 4-MFSK 80 2 127 224
基于n-MFSK调制的激光致声空-水跨介质通信方法
Method of laser-generated sound with n-MFSK modulation for air-water trans-media communication
-
摘要: 为了改善空中平台和水下目标之间的激光致声通信技术, 采用了一种多进制多频移键控(n-MFSK)调制方式来提升激光致声空-水跨介质通信速率的方法。在借助激光致声热膨胀效应实现空-水界面处的光声转换基础上, 分别采用长脉冲法和重复频率法进行了调制及仿真验证, 得到了在调制频率数量为2的情况下, 2-MFSK调制可在2-FSK调制基础上将通信速率提高1倍的结果。结果表明, 长脉冲法主要通过激光阵列结合频率叠加的方式实现n-MFSK调制, 重复频率法则通过控制激光器频率变化进而在时域上以分配时间段的方式实现n-MFSK调制; 随着调制频率数量增加, 相对于n-FSK调制, n-MFSK调制下通信速率更高, 并可改善频带利用率, 但声压级和水中传输距离会随符号码元持续时间内调制频率数量增加而减小。该研究为未来激光致声空-水跨介质通信实际应用提供了参考。Abstract: In order to improve the laser acoustic communication technology between air platforms and underwater objects, the n-multifrequency shift keying (n-MFSK) modulation method was proposed to further enhance the laser acoustic air-water cross-media communication rate. Based on the thermal expansion effect of laser acoustic, the modulation methods of the long-pulse-duration laser method and the high repetitive rate method were defined and simulated. The result was obtained that 2-MFSK modulation doubles the communication rate on top of 2-FSK modulation for a modulation frequency number of 2. The n-MFSK modulation was achieved by the long-pulse-duration laser method, mainly through the laser array by means of frequency superposition, and the high repetitive rate method by controlling the allocation of time periods on the time domain of the laser frequency change. The results show that as the number of modulation frequencies increases, n-MFSK modulation can increase the communication rate and improve the band utilization relative to n-frequency shift keying (n-FSK) modulation at the cost of loss of the sound pressure level (SPL) and in-water-range. This study provides a reference for the practical application of laser acoustic air-water cross-media communication in the future.
-
parameter value laser wavelength 1053 nm pulse width 2 ms peak power 25 kW coefficient of thermal expansion 2×10-4 ℃-1 laser beam radius 6 mm specific heat capacity of liquids 4.17×103 J·kg-1·K-1 speed of sound in liquids 1500 m/s optical absorption coefficient 15.7 m-1 optimal center frequency 10 kHz 表 2 长脉冲法n-FSK和n-MFSK调制的比较
Table 2. Comparison of n-FSK and n-MFSK modulation by long-pulse-duration laser method
modulation
methoddata rate/
(bit·s-1)bandwidth/
kHzSPL/
dBin-water
range/m2-FSK 10 1 137 708 4-FSK 20 2 137 708 8-FSK 30 4 137 708 16-FSK 40 8 137 708 2-MFSK 20 1 131 355 3-MFSK 30 1.5 128 251 4-MFSK 40 2 125 178 parameter value laser wavelength 1060 nm pulse width 8 ns single pulse energy 100 mJ~800 mJ optical absorption coefficient 13.7 m-1 angle of observation from the vertical 1.31 rad 表 4 重复频率法n-FSK和n-MFSK调制的比较
Table 4. Comparison of high repetitive rate method n-FSK and n-MFSK modulation
modulation
methoddata rate/
(bit·s-1)bandwidth/
kHzSPL/
dBin-water
range/m2-FSK 20 1 139 897 4-FSK 40 2 139 897 8-FSK 60 4 139 897 16-FSK 80 8 139 897 2-MFSK 40 1 133 459 3-MFSK 60 1.5 130 316 4-MFSK 80 2 127 224 -
[1] SPEARS A, HOWARD A M, WEST M, et al. Determining underwater vehicle movement from sonar data in relatively featureless seafloor tracking missions[C]//IEEE Winter Conference on Applications of Computer Vision. New York, USA: IEEE, 2014: 909-916. [2] 方尔正, 李宗儒, 桂晨阳. 穿海牵天提升对潜通信保障能力——跨介质通信技术现状及展望[J]. 国防科技工业, 2022, 262(2): 59-62. FANG E Zh, LI Z R, GUI Ch Y. Crossing the sea and holding the sky to enhance the submarine communication security capability with the current situation and outlook of cross-media communication technology[J]. Defence Science & Technology Industry, 2022, 262(2): 59-62(in Chinese). [3] NASH P, STRUDLEY A, CRICKMORE R, et al. High efficiency TDM/WDM architectures for seismic reservoir monitoring[J]. Proceedings of the SPIE, 2009, 7503: 1125-1128. [4] WESTERVELT P J, LARSON R S. Laser-excited broadside array[J]. The Journal of the Acoustical Society of America, 1973, 54(1): 121-122. doi: 10.1121/1.1913551 [5] BUNKIN F V, KARLOV N V, KOMISSAROV V M, et al. Excitation of sound when a surface layer of a liquid absorbs a laser pulse[J]. Soviet Journal of Experimental and Theoretical Physics Letters, 1971, 13(9): 341-343. [6] ZHANG Y D, ZHAO Zh Y, YAO Ch B, et al. The nonlinear absorption and optical limiting in phenoxy-phthalocyanines liquid in nano-and femto-second regime: Experimental studies[J]. Optics & Laser Technology, 2014, 58(11): 207-214. [7] CARVER C J, TIAN Z, ZHANG H, et al. Amphilight: Direct air-water communication with laser light[J]. GetMobile: Mobile Computing & Communications review, 2021, 24(3): 26-29. [8] LEI L, JU Z, LEI Z, et al. Research on laser induced acoustic source based underwater communication system[J]. Proceedings of the SPIE, 2016, 10157: 477-483. [9] PARTON S, DALGLEISH F, BEAUJEAN P P, et al. Acquisition and registration of bathymetric acoustic data and MOFSLI (multiple overlapping field of view serial laser imager)[C]// 2013 MTS/IEEE OCEANS-Bergen. New York, USA : IEEE, 2013: 1-8. [10] JIANG H, QIU H, HE N, et al. Research on the optoacoustic communication system for speech transmission by variable laser-pulse repetition rates[J]. Results in Physics, 2018, 9(4): 1291-1296. [11] 任亚辉, 林菊平, 童勇, 等. 高能脉冲氙灯驱动电路特性分析[J]. 激光技术, 2022, 46(3): 397-401. REN Y H, LIN J P, TONG Y, et al. Analysis of driving circuit characteristics of high-power pulsed xenon lamp[J]. Laser Techno-logy, 2022, 46(3): 397-401(in Chinese) [12] 陈清明. 脉冲激光在液体中激发的声波特性研究[J]. 光学与光电技术, 2006, 4(3): 28-31. CHEN Q M. Characteristics of photoacoustics induced by pulse laser in liquid[J]. Optics & Optoelectronic Technology, 2006, 4(3): 28-31(in Chinese). [13] 王晓宇, 王江安, 宗思光, 等. 光击穿机制下的光声能量转换效率[J]. 强激光与粒子束, 2013, 25(3): 579-583. WANG X Y, WANG J A, ZONG S G, et al. Laser acoustic energy conversion efficiency in optical breakdown mechanism[J]. High Power Laser and Particle Beams, 2013, 25(3): 579-583(in Chinese). [14] WHITE R M. Elastic wave generation by electron bombardment or electromagnetic wave absorption[J]. Journal of Applied Physics, 1963, 34(7): 2123-2124. doi: 10.1063/1.1729762 [15] WHITE R M. Generation of elastic waves by transient surface heating[J]. Journal of Applied Physics, 1963, 34(12): 3559-3567. doi: 10.1063/1.1729258 [16] BLACKMON F, ESTES L, FAIN G. Linear optoacoustic underwater communication[J]. Applied Optics, 2005, 44(18): 3833-3845. doi: 10.1364/AO.44.003833 [17] BLACKMON F A, ANTONELLI L. Experimental detection and reception performance for uplink underwater acoustic communication using a remote, in-air, acousto-optic sensor[J]. IEEE Journal of Oceanic Engineering, 2006, 31(1): 179-187. [18] BLACKMON F, ANTONELLI L. Remote, aerial, trans-layer, linear and nonlinear downlink underwater acoustic communication[C]// OCEANS 2006. New York, USA: IEEE, 2006: 1-7. [19] ANTONELLI L, BLACKMON F. Laser-based method for docking unmanned underwater vehicle to a submarine: US, 8047149B1[P]. 2011-11-01. [20] 彭水, 张明敏, 王江安. 一种采用高重复频率激光进行水声通信的方法[J]. 中国激光, 2012, 39(7): 0705005. PENG Sh, ZHANG M M, WANG J A. Method of laser-generated sound with high repetition rate for underwater acoustic communication[J]. Chinese Journal of Lasers, 2012, 39(7): 0705005 (in Chinese). [21] LI P, ZHAO Y, ZHOU Zh Q, et al. Research on laser induced acoustic detection of trans-media aerial-underwater[J]. Infrared and Laser Engineering, 2021, 50(5): 20200310 (in Chinese).