0 引　言

随着海洋在经济和军事应用中重要性的不断凸显，水声通信作为开发和利用海洋科学勘探、海洋环境监测和水下机器人通信等相关领域的有效工具，受到高度重视^[1-3]。在水声通信网络中，一般均采用正交多址接入(orthogonal multiple access, OMA)技术提供数据服务^[4]，虽然OMA技术通过用户独占时域/频域资源块的方式，极大地避免了数据冲突和用户间干扰，但这种资源接入方式也极大地限制了用户接入数量和资源利用率的提升。而随着水声通信用户不断增加，水声通信网络需要在有限的系统带宽和功率范围内，为越来越多的用户提供服务^[5]。因此，探索新的多址接入方式对进一步提高水声通信网络的通信质量和资源利用率至关重要。

与OMA技术相比，非正交多址接入技术(non-orthogonal multiple access，NOMA)通过复用功率域，增加非正交的资源控制，同时发送多个用户信号，达到了大规模提高频谱利用率和用户接入数量的目的^[6-7]。近年来，国内外学者从不同的角度，对NOMA技术在水声通信网络中的应用做了初步的探索研究。Goutham V等^[8]在非精确串行干扰消除(successive interference cancellation，SIC)技术假设下，研究了基于NOMA技术的水声协作通信网络的可达容量和中断概率，揭示了NOMA技术在水声协作系统资源利用率上的优势。Cheon J等^[9]在基于NOMA技术的下行水声通信网络中，提出一种等传输时间的功率分配方案，以避免NOMA节点之间的信息传输时间不等而造成资源浪费。Goutham V等^[10]将全双工技术和NOMA技术引入水声协作传感器网络中，分析推导自干扰和非精确SIC假设下系统的遍历容量和能量效率，揭示了中继节点位置、风速和航运程度对系统性能的影响。

虽然NOMA技术应用于水声通信网络可以明显提高系统资源利用率和可达性能，但这些性能改善是在NOMA节点使用相同频段的前提下取得的。实际中，通信距离越远，水声节点适用的最优频率越小，3 dB带宽越窄^[11]。值得注意的是，虽然通信距离不同的节点适用的频段不会完全重叠，但分布在同一目的节点周围的节点通信距离不会相差特别大。因此，水声通信场景中，节点的适用通信频段虽然不会完全重叠，但会部分重叠。

Kim B等^[12]通过调整NOMA用户间的频段重叠率，提出的部分NOMA (partial NOMA, P-NOMA)技术以达到减少NOMA用户间干扰和进一步提升NOMA系统可达性能的目的。考虑到水声通信场景的独特特点，在文献[11]和文献[12]基础上，本文在上行水声通信中引入P-NOMA技术，并在功率均匀分布和非精确SIC假设下，分析推导各P-NOMA节点在通信频段内的可达容量和能量效率表达式，以期通过P-NOMA技术的引入，进一步有效利用有限的水声频段资源，提高系统资源利用率，为水声通信网络在下一代移动通信网络中服务更多的水声节点/用户奠定基础。

1 系统模型

如图1所示，本文所提的面向非精确SIC的P-NOMA水声上行通信网络由一个目的节点和M $\left( {M \geqslant 2} \right)$ 个水声监测节点组成，监测节点根据任务需要分布在目的节点周围，各监测节点的通信频率和3 dB带宽由其与目的节点之间的通信距离决定。在监测节点和目的节点均配置单天线的假设下，本文假设随机选择2个通信距离差异较大的节点i和j组成一个P-NOMA组，该组水声节点利用P-NOMA技术同时传输信号，如图2所示。由于通信距离不同节点的频段会部分重合，该P-NOMA组节点在频段非重叠范围内使用OMA技术发送各自监测信号；在频段重叠范围内( $ B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right) $ )，利用最大或受限发送功率同时同频发送监测信息给目的节点。目的节点再将接收到的信号通过电磁波转发到控制中心/船载通信基站进行解码，为了简化分析，本文假设目的节点具有解码能力。

利用P-NOMA技术，目标节点处接收到的信号为：

$ y\left( {{f_m}} \right) = \sum\limits_{m = i,j}^{} {\frac{{\sqrt {{P_m}} }}{{\sqrt {A\left( {{l_m},{f_m}} \right)} }}} {X_m}\left( {{f_m}} \right) + N\left( {{f_m}} \right)。$

(1)

式中： $ {X_m}\left( {{f_m}} \right) $ 为监测节点m $\left( {m = i,j} \right)$ 在发射功率 $ {P_m} $ 条件下的传输信号； $ A\left( {{l_m},{f_m}} \right) $ 为依赖于传播距离 $ {l_m} $ 和频率 $ {f_m} $ 的声传播路径损耗函数，可写为^[13-14]：

$ 10\lg \left( {A\left( {{l_m},{f_m}} \right)} \right) = {k_m} \times 10\lg {l_m} + \left( {{{{l_m}} \mathord{\left/ {\vphantom {{{l_m}} {{{10}^3}}}} \right. } {{{10}^3}}}} \right) \times 10\lg a\left( {{f_m}} \right)。$

(2)

其中： $ {k_m} $ 为能量扩散因子， $ {k_m} $ =1，1.5和2分别表示浅水通信时的柱面扩散，实际通信扩散和深水通信时的球面扩散， $ a\left( {{f_m}} \right) $ 为声信号吸收系数，可写为^{[13 ,14]}：

$ 10\lg a\left( {{f_m}} \right) = \frac{{0.11f_m^2}}{{1 + f_m^2}} + \frac{{44f_m^2}}{{4100 + f_m^2}} + 2.75 \times {10^{ - 4}}f_m^2 + 0.003 。$

(3)

由式(2)~式(3)可知，无论是通信距离还是频率的增加，都会增大水声通信的路径损耗，这就限制了不同通信距离下，水声节点适用的3 dB带宽和最优频率。

$ N\left( {{f_m}} \right) $ 表示海洋环境总噪声谱密度(dB re μ Pa²/Hz)，可写为：

$ N\left( {{f_m}} \right) = {N_{tb}}\left( {{f_m}} \right) + {N_s}\left( {{f_m}} \right) + {N_w}\left( {{f_m}} \right) + {N_{th}}\left( {{f_m}} \right) 。$

(4)

式中： $ \lg {N_{tb}}\left( {{f_m}} \right) = 1.7 - 3\lg \left( {{f_m}} \right) $ 为海洋湍流噪声； $\lg {N_s} ( {{f_m}} ) = 4 + 2\left( {s - 0.5} \right) + 2.6\lg \left( {{f_m}} \right) - 6\lg \left( {0.03 + {f_m}} \right)$ 为航运噪声， $ s = 0,0.5,1 $ 分别表示轻微，中等和较强航运密度； $ \lg {N_{th}}\left( {{f_m}} \right) = - 1.5 + 2\lg \left( {{f_m}} \right) $ 表示热噪声； $\lg {N_w}\left( {{f_m}} \right) = 5 + 0.75\sqrt w + 2\lg \left( {{f_m}} \right) - 4\lg \left( {0.4 + {f_m}} \right)$ 为风浪噪声， $w \geqslant 0$ m/s为海面风速。

2 P-NOMA和OMA技术下的可达容量和能量效率 2.1 P-NOMA技术非精确SIC假设下的可达容量和能量效率

假设 $ {l_i} < {l_j} $ ，由于目标节点处节点 $ i $ 的信号较强，在频段 $ B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right) $ 内，目标节点将节点j的信号当作干扰，直接解码节点i的信号。因此，节点 $ i $ 的信干噪比（signal-to-interference-plus-noise ratio，SINR）为：

$ {\gamma _i} = \left\{ {\begin{aligned} &{\dfrac{{\displaystyle\int_{{B_{ionly}}} {{A^{ - 1}}\left( {{l_i},{f_i}} \right){S_i}\left( {{f_i}} \right){\rm{d}}{f_i}} }}{{\displaystyle\int_{{B_{ionly}}} {N\left( {{f_i}} \right){{{\rm{d}}}}{f_i}} }}，{f_i} \in {B_{ionly}}} ，\\ &{\dfrac{{\displaystyle\int_{{B_{mix}}} {{A^{ - 1}}\left( {{l_i},{f_{mix}}} \right){S_i}\left( {{f_{mix}}} \right){\rm{d}}{f_{mix}}} }}{{\displaystyle\int_{{B_{mix}}} {\left[ {{A^{ - 1}}\left( {{l_j},{f_{mix}}} \right){S_j}\left( {{f_{mix}}} \right) + N\left( {{f_{mix}}} \right)} \right]{\rm{d}}{f_{mix}}} }}，{f_{mix}} \in {B_{mix}}}。\end{aligned}} \right .$

(5)

式中： $ {B_{ionly}}{\text{ = }}B\left( {{l_i}} \right) - \left( {B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right)} \right) $ ， ${B_{mix}} = B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right)$ ， $ {S_m}\left( {{f_m}} \right) = {{{P_m}} \mathord{\left/ {\vphantom {{{P_m}} {{B_{3 dBm}}}}} \right. } {{B_{3dBm}}}}\left( {m = i,j} \right) $ 是在功率均匀分布下，节点m信号传输的功率谱密度， $ {B_{3dBm}} $ 是节点m的3dB带宽。

监测节点j由于距离较远，信号较弱，在频段 $ B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right) $ 内，目标节点只能通过SIC技术先解码并减去节点i的信号后，再从剩余信号中解码节点j的信号。利用SIC技术，节点j的SINR为：

$ \gamma _j = \left\{ {\begin{array}{*{20}{c}} {\dfrac{{\displaystyle\int_{{B_{mix}}} {{A^{ - 1}}\left( {{l_j},{f_{mix}}} \right){S_j}\left( {{f_{mix}}} \right){\rm{{\rm{d}}}}{f_{mix}}} }}{{\displaystyle\int_{{B_{mix}}} {\left[ {N\left( {{f_{mix}}} \right) + \xi {A^{ - 1}}\left( {{l_i},{f_{mix}}} \right){S_i}\left( {{f_{mix}}} \right)} \right]{\rm{{\rm{d}}}}{f_{mix}}} }}，{f_{mix}} \in {B_{mix}}} ，\\ {\dfrac{{\displaystyle\int_{{B_{jonly}}} {{A^{ - 1}}\left( {{l_j},{f_{jonly}}} \right){S_j}\left( {{f_{jonly}}} \right){\rm{d}}{f_{jonly}}} }}{{\displaystyle\int_{{B_{jonly}}} {N\left( {{f_{jonly}}} \right){\rm{d}}{f_{jonly}}} }}，{f_{jonly}} \in {B_{jonly}}} 。\end{array}} \right. $

(6)

式中： $ {B_{jonly}} = B\left( {{l_j}} \right) - \left( {B\left( {{l_i}} \right) \cap B\left( {{l_j}} \right)} \right) $ ， $\xi \in \left[ {0,1} \right]$ 是重合频段内SIC解码的非精确程度。当 $0 < \xi < 1$ 时，节点j部分消除节点i对其造成的干扰；当 $\xi = 0$ 时，节点j可完全消除节点i对其造成的干扰。

如式(1)、式(5)和式(6)所示，在非重合频段内，目标节点接收的信号没有来自其他监测节点信号的干扰；在重合频段内，通信距离较远节点的信号虽然受到通信距离较近节点信号的干扰，但利用SIC技术可大大降低其干扰程度。

因此，在功率均匀分布假设下，各P-NOMA节点的可达容量为：

$ \begin{split}& R_i^{P - NOMA} = \int_{{B_{ionly}}} {{{\log }_2}} \left( {1 + \frac{{{{{P_i}} \mathord{\left/ {\vphantom {{{P_i}} {B\left( {{l_i}} \right)}}} \right. } {B\left( {{l_i}} \right)}}}}{{A\left( {{l_i},{f_i}} \right)N\left( {{f_i}} \right)}}} \right){\rm{d}}{f_i} +\\ &\int_{{B_{mix}}} {{{\log }_2}} \left( {1 + \frac{{{{{A^{ - 1}}\left( {{l_i},{f_{mix}}} \right){P_i}} \mathord{\left/ {\vphantom {{{A^{ - 1}}\left( {{l_i},{f_{mix}}} \right){P_i}} {B\left( {{l_i}} \right)}}} \right. } {B\left( {{l_i}} \right)}}}}{{{A^{ - 1}}\left( {{l_j},{f_{mix}}} \right){{{P_j}} \mathord{\left/ {\vphantom {{{P_j}} {B\left( {{l_j}} \right)}}} \right. } {B\left( {{l_j}} \right)}} + N\left( {{f_{mix}}} \right)}}} \right){\rm{d}}{f_{mix}} ，\\[-20pt]\end{split}$

(7)

$ \begin{split} & R_j^{P - NOMA} = \int_{{B_{jonly}}} {{{\log }_2}} \left( {1 + \frac{{{{{P_j}} \mathord{\left/ {\vphantom {{{P_j}} {B\left( {{l_j}} \right)}}} \right. } {B\left( {{l_j}} \right)}}}}{{A\left( {{l_j},{f_j}} \right)N\left( {{f_j}} \right)}}} \right){\rm{d}}{f_j} + \\ &\int_{{B_{mix}}} {{{\log }_2}} \left( {1 + \frac{{{{{A^{ - 1}}\left( {{l_j},{f_{mix}}} \right){P_j}} \mathord{\left/ {\vphantom {{{A^{ - 1}}\left( {{l_j},{f_{mix}}} \right){P_j}} {B\left( {{l_j}} \right)}}} \right. } {B\left( {{l_j}} \right)}}}}{{\xi {A^{ - 1}}\left( {{l_i},{f_{mix}}} \right){{{P_i}} \mathord{\left/ {\vphantom {{{P_i}} {B\left( {{l_i}} \right)}}} \right. } {B\left( {{l_i}} \right)}} + N\left( {{f_{mix}}} \right)}}} \right){\rm{d}}{f_{mix}}.\\[-20pt]\end{split}\;\;$

(8)

水声通信网络中，水声节点一般依靠海洋热能或电池提供能量，尽管海洋热能能量无穷，但电池充放电次数有限。因此，分析水声通信系统的能量效率比单纯分析系统的可达容量更有意义。基于此，对基于P-NOMA技术的水声上行通信系统的能量效率进行分析。在式(7)和式(8)的可达容量条件下，各P-NOMA节点的能量效率可以写为：

$ \eta _m^{P - NOMA} = \dfrac{{R_m^{P - NOMA}}}{{\varsigma {P_m} + {P_c}}} 。$

(9)

式中： $\varsigma $ 为与频率有关的换能器效率， ${P_c}$ 为包含电路功率放大器、信号处理块和声学调制解调器等电路消耗的功率^[14]。

2.2 OMA技术下的可达容量和能量效率

在采用OMA技术的水声系统模型中，例如采用水声通信系统中常用的时分多址接入技术（time division multiple access，TDMA），假定节点利用功率在不同的时隙内发送信号给目标节点，目标节点处的接收信号为：

$ y\left( {{f_m}} \right) = \frac{{\sqrt {{P_m}} }}{{\sqrt {A\left( {{l_m},{f_m}} \right)} }}{X_m}\left( {{f_m}} \right) + N\left( {{f_m}} \right) ，$

(10)

因此，基于TDMA技术的水声节点m的SINR为：

$ \gamma _m^{OMA} = \frac{{\int_{{B_m}} {{A^{ - 1}}\left( {{l_m},{f_m}} \right){S_m}\left( {{f_m}} \right)d{f_m}} }}{{\int_{{B_m}} {N\left( {{f_m}} \right)d{f_m}} }}{\text{ }}{f_m} \in {B_m} ，$

(11)

可达容量为：

$ R_m^{OMA} = 0.5\int_{{B_m}} {{{\log }_2}} \left( {1 + \frac{{{{{P_m}} \mathord{\left/ {\vphantom {{{P_m}} {B\left( {{l_m}} \right)}}} \right. } {B\left( {{l_m}} \right)}}}}{{A\left( {{l_m},{f_m}} \right)N\left( {{f_m}} \right)}}} \right){\rm{d}}{f_m}。$

(12)

式中：系数0.5表示水声节点i和j各占用半个时隙资源。另外，根据式(9)，节点m的能量效率为： $\eta _m^{OMA} = {{R_m^{OMA}} \mathord{\left/ {\vphantom {{R_m^{OMA}} {\left( {\varsigma {P_m} + {P_c}} \right)}}} \right. } {\left( {\varsigma {P_m} + {P_c}} \right)}}$ 。

3 仿真分析

在无风 $ w = 0 $ ，中等航运密度 $ s=0.5 $ ，实际通信能量扩散因子 $ {k}_{m}=1.5 $ ，电路消耗功率 ${P_c} = 0.5$ mW，换能器效率 $\varsigma = 80\text{% }$ 环境下，仿真OMA方案和本文所提方案在不同通信距离、传输功率和非精确SIC程度下，各水声节点和系统的可达容量和能量效率的变化情况，同时分析关键参数对系统性能的影响。

1）在非精确SIC $\xi = 0.01$ 的假定下，图3比较了节点i、节点j和系统分别采用P-NOMA和OMA多址技术下的可达容量。可知，随着发送功率的增加，无论是哪一种接入技术下，各节点和系统的可达容量均随着传输功率的增加而增加。从图3(a)可以发现，节点i/j的通信距离越近/远，节点i的可达容量越高，这个现象可以用式(7)给出的 $R_i^{P - NOMA}$ 来解释：目标节点在解码节点i信号的过程中，将节点j的信号视为干扰，因此节点j信号强度越弱，其对节点i的可达容量影响越小。从图3 (b)可以发现，节点i/j的通信距离越远/近，节点j的可达容量越高，这是源于节点i的通信距离越远，即式(2)表示的 $ A\left({l}_{i}, {f}_{i}\right) $ 越大时，式(8)表示的 $ {R}_{j}^{P-NOMA} $ 在重合频段内来自节点i的干扰噪声功率谱密度越小。特别地，由于SIC非精确解码的存在，节点j的可达容量最终将趋于饱和，这一现象表明应用P-NOMA技术时，必须为通信距离较远的节点j设计合理的传输功率，以同时保障节点j的可达性能和资源利用率。此外，与OMA技术相比，在某些距离组合中，如 ${l_i} = 4$ km, ${l_j} = 6$ km距离组合中，P-NOMA技术的性能优势随着传输功率的增大而变差，结果表明需要根据各节点的功率配置和通信距离，设计合理的OMA或P-NOMA策略切换标准，以及P-NOMA节点组选择标准，以保障水声通信系统的服务质量和资源利用率。

2）根据P-NOMA和OMA下节点i和节点j的可达容量，图4给出了对应节点的能量效率。可以看到，在不同的距离组合中，2种多址接入技术对应的能量效率曲线均先上升后在某个传输功率值处开始下降。可以注意到，增加节点i或节点j的通信距离，会加大能量效率对应的传输功率转折值。同时，与图3相对应，节点i的能量效率随着节点i/j通信距离的增加/减小而减小，节点j的能量

效率随着节点i/节点j通信距离的减小/增加而减小。此外，可以清楚看到节点i和节点j在P-NOMA技术下的能量效率均明显优于其在OMA技术下的能量效率。根据图4的结论，与OMA技术相比，同等条件下，引入P-NOMA技术的水声通信网络可以获得更高的电池使用效率、更长的电池使用时间和更小的电池容量。

3）为了揭示非精确SIC对可达容量的影响，图5给出了不同通信距离组合下，P-NOMA技术下节点j和系统可达容量随SIC解码非精确程度的变化情况，其中 ${P_m} = 145$ dB re μPa。可以看到，随着 $\xi $ 的增加，节点j和系统的可达容量最开始均快速下降，随后下降速度越来越慢。这一现象源于 $\xi $ 越大，重合频段内非精确SIC引入的用户间干扰越大，节点j的可达容量越依赖于非重合频段内所取得的性能。关于如何综合考虑各P-NOMA水声节点的性能需求、功率配置和能量效率要求，以及目标节点SIC精确程度，从系统角度设计一个最优的资源控制方案将作为后续工作。

4 结　语

为了满足日益增长的水声节点接入需求，以及提高水声通信网络的资源利用率，本文在功率均匀分布假设下，研究了面向非精确SIC的P-NOMA上行水声通信网络的容量分析和能量效率分析。仿真结果表明，与OMA相比，选择合适的P-NOMA节点组，所提方案可以大幅提升水声通信网络的可达容量和能量效率。此外，本文研究了SIC解码非精确程度和位置分布对系统可达容量的影响。