屈文璋，许 诚，孙 溥，李丁山，闫建峰

(中国舰船研究院，北京 100192)

0 引 言

海水作为良导体，对高频电磁信号具有强烈屏蔽作用，而低频信号不易被吸收，信号强度一定的情况下，频率越低，衰减越慢，传播距离越远，因此可以利用极低频电磁信号（本文为0.1~300 Hz 范围内）进行深水通信[1]。

卓贤军等[2]对海水中极低频信号传播过程中所受影响因素进行了研究，认为接收端得到的电、磁信号与海水本身电阻率、深度以及海底的海床地质条件有关。本文从极化模式的角度出发，推导出TE 和TM 极化模式下的极低频电磁波在海水中的传播公式，利用有限元算法，研究2 种极化模式下极低频电磁波在不同深度的衰减特性。

1 两种极化模式公式推导

海水介质的电阻率一般为0.25 Ω ·m，根据极低频所处频段范围，估算出位移电流与传导电流的最大比值ωε/σ ≈ρ/2T×10-10=4×10-9，因此在海水介质中可以忽略位移电流对场分布的影响，即此时研究的是似稳电磁场问题[3]，于是，谐变场的Maxwell 方程组表示为：

在笛卡尔坐标系中，令z轴垂直向下，x和y轴在地表水平面内，可以把式（1）写为：

把式（2）写为：

当极低频电磁波通过“地-电离层”空腔中进行传导[4]，抵达近海海面时可以将其近似看作平面电磁波垂直入射到海面，将海水简单近似为均匀各向同性介质，此时其电磁场在水平方向上是均匀的，可得:

此时将式（11）代入式（5）~式（10），可以得到：

从式（12）~式（17）可以看出，Ex只和Hy有关，Ey只 和Hx有关，且两者均沿着z轴传播。假设真空中电磁波波前与x轴平行，在y和z坐标平面内考虑问题，这时的电磁波可以解耦为电场仅有水平分量的极化模式横电（TE）波型和磁场仅有水平分量的极化模式横磁（TM）波型，同时，可以看出2 种波型中电场和磁场没有垂向分量。

TE 极化模式（Ex-Hy）：

TM 极化模式（Ey-Hx）：

图1 TE 和TM 模式的研究区域Fig. 1 The research area of TE and TM polarization modes

2 有限元数值算法

2.1 TE 模式

首先以TE 模式为例，给定的电场Ex所满足的微分方程为：

电场所满足的边界条件为：

1）在上边界AB 处假设Ex|AB=1。

2）下边界CD 以下为均质岩石，电磁波在CD 以下将按照负指数衰减，即Ex=Ex0e-az。式中：Ex0为常数；，对u求导，因为在C D 处 有， 所以C D 处的边界条件为：

该边值问题所对应的变分问题为：

采用有限元方法进行计算，首先应将区域剖分为若干网格，在网格内电导率必须是连续的。在保证计算精度的前提下，采用三次插值。最终，计算出节点处的电场值[5-6]。

2.2 TM 模式

磁场Hx满足的微分方程如下式：

不管海水介质电性如何分布，空气中的Hx近似为常量，其上边界条件可取海面AB 的Hx|AB=1，下边界CD 的边界条件同TE 极化模式。

3 模型仿真计算

3.1 TE 和TM 模式仿真计算

本文的主要研究内容为海水中极低频电磁波的衰减，因此构建为3 层模型，第1 层为海水层，水深400 m，电阻率0.25 Ω·m第2 层为低阻岩石层，厚度200 m，电阻率为50 Ω ·m，第3 层之下为半空间的均质高阻岩层，电阻率为2 000 Ω ·m。

根据前文的公式及边界条件，经过有限元算法代码实现，得到TE 和TM 模式下，4 Hz，17 Hz，88 Hz 和178 Hz 的电、磁场随深度的衰减情况，其中Ex0代表海面处的水平电场值，Ex代表某一深度处的水平电场值，Hy0代表海面处的水平磁场值，Hy代表某一深度处的水平电场值，如图2 和图3 所示。

图2 TE 模式下海水中Ex 和HY 随深度衰减Fig. 2 The attenuation curves of Ex and HY in the seawater with depth in TE polarization mode

图3 TM 模式下海水中Ey 和Hx 随深度衰减Fig. 3 The attenuation curves of Ey and Hx in the seawater with depth in TM polarization mode

1）从TE 和TM 模式均可以看出，极低频水平电场和磁场在海水中的衰减特征均表现为频率越高，衰减越快，同时海水中磁场衰减快，电场相对衰减慢。

2）水平电场分量Ex和Ey、水平磁场分量Hx和Hy，两者在TE 模式或TM 模式下均有相同的衰减特征：衰减曲线走势一致、幅度一致。

3.2 电阻率和厚度对衰减特征的影响

极低频电磁波在海水中传播，海底的海床地质条件也会影响到电磁波在海水中的传播特性[7-9]。为进一步研究TE 和TM 极化模式下海底介质电阻率值变化对海水中极低频电磁波衰减造成的影响，为此以0.5 Hz为例，改变海底岩石电阻率值和厚度值，如表1 所示。再据此参数表进行仿真计算，得到0.5 Hz 时水平电场分量和水平磁场分量随深度的的衰减特征，如图4所示。

表1 理论模型参数表Tab. 1 The table of theoretical model parameter

由图4 和图5 可以看出：1）当模型参数发生变化时，TE 和TM 模式的极低频电、磁场衰减特征几乎相同；2）从TE 和TM 模式的水平电场和磁场分量计算结果均看出，模型1 和模型3 的计算结果一样，在图上出现曲线重合的现象。说明海水下低阻岩石的厚度增加对电场、磁场的衰减几乎无影响；3）模型2 计算的结果相比较模型1 和模型3 发现，当低阻岩层电阻率增大时，电场衰减速度变慢，而磁场衰减速度变快；4）频率降低，变为0.5 Hz 时，电场衰减速度变慢的现象更加明显，而磁场这一现象并不明显。

图4 频率为0.5 Hz 时水平电场分量随深度的衰减曲线Fig. 4 The attenuation curve of horizontal electric field components with depth at frequency 0.5 Hz

图5 频率为0.5 Hz 时水平磁场分量随深度的衰减曲线Fig. 5 The attenuation curve of horizontal magnetic field components with depth at frequency 0.5 Hz

4 结 语

通过对TE 和TM 模式的极低频电磁波传播过程对比分析发现：

1）在海水中2 种模式下的水平电场分量和磁场分量衰减特性基本一致，说明极化模式对极低频电磁波的传播几乎无影响；

2）仿真结果表明，当信号频率越低时，电场和磁场分量衰减速度越慢，这和电磁波传播特征完全一致；

3）当频率在0.5 Hz 时，磁场衰减速度明显快于电场，同时当海底岩石电阻率值越大，电场衰减速度变慢，而磁场衰减速度变快。因此在实际布设中可以尝试考虑利用磁力梯度仪开展水平磁场分量在垂向上的信号采集工作。