海水作为良导体,对高频电磁信号具有强烈屏蔽作用,而低频信号不易被吸收,信号强度一定的情况下,频率越低,衰减越慢,传播距离越远,因此可以利用极低频电磁信号(本文为0.1~300 Hz范围内)进行深水通信[1]。
卓贤军等[2]对海水中极低频信号传播过程中所受影响因素进行了研究,认为接收端得到的电、磁信号与海水本身电阻率、深度以及海底的海床地质条件有关。本文从极化模式的角度出发,推导出TE和TM极化模式下的极低频电磁波在海水中的传播公式,利用有限元算法,研究2种极化模式下极低频电磁波在不同深度的衰减特性。
1 两种极化模式公式推导海水介质的电阻率一般为0.25
$ \nabla \times \vec E = i\omega \mu \vec H, $ | (1) |
$ \nabla \times \vec H = \sigma \vec E, $ | (2) |
$ \nabla ·\overrightarrow{E}=0, $ | (3) |
$ \nabla ·\overrightarrow{H}=0。$ | (4) |
在笛卡尔坐标系中,令z轴垂直向下,x和y轴在地表水平面内,可以把式(1)写为:
$ \frac{{\partial {E_z}}}{{\partial y}}\; - \frac{{\partial {E_y}}}{{\partial z}} = i\omega \mu {H_x}, $ | (5) |
$ \frac{{\partial {E_x}}}{{\partial z}}\; - \frac{{\partial {E_z}}}{{\partial x}} = i\omega \mu {H_y}, $ | (6) |
$ \frac{{\partial {E_y}}}{{\partial x}}\; - \frac{{\partial {E_x}}}{{\partial y}} = i\omega \mu {H_z}。$ | (7) |
把式(2)写为:
$ \frac{{\partial {H_z}}}{{\partial y}}\; - \frac{{\partial {H_y}}}{{\partial z}} = \sigma {E_x}, $ | (8) |
$ \frac{{\partial {H_x}}}{{\partial z}}\; - \frac{{\partial {H_z}}}{{\partial x}} = \sigma {E_y}, $ | (9) |
$ \frac{{\partial {H_y}}}{{\partial x}}\; - \frac{{\partial {H_x}}}{{\partial y}} = \sigma {E_z}。$ | (10) |
当极低频电磁波通过“地-电离层”空腔中进行传导[4],抵达近海海面时可以将其近似看作平面电磁波垂直入射到海面,将海水简单近似为均匀各向同性介质,此时其电磁场在水平方向上是均匀的,可得:
$ \frac{{\partial {E_z}}}{{\partial x}} = \frac{{\partial {E_z}}}{{\partial y}} = 0\text{,}\frac{{\partial {H_z}}}{{\partial x}} = \frac{{\partial {H_z}}}{{\partial y}} = 0, $ | (11) |
此时将式(11)代入式(5)~式(10),可以得到:
$ - \frac{{\partial {E_y}}}{{\partial z}} = i\omega \mu {H_x}, $ | (12) |
$ \frac{{\partial {E_x}}}{{\partial z}}\; = i\omega \mu {H_y}, $ | (13) |
$ {H_z} = 0, $ | (14) |
$ - \frac{{\partial {H_y}}}{{\partial z}} = \sigma {E_x}, $ | (15) |
$ \frac{{\partial {H_x}}}{{\partial z}}\; = \sigma {E_y}, $ | (16) |
$ {E_z} = 0。$ | (17) |
从式(12)~式(17)可以看出,
TE极化模式(
$ \frac{{\partial {E_x}}}{{\partial z}} = i\omega \mu {H_y}, $ | (18) |
$ - \frac{{\partial {H_y}}}{{\partial z}} = \sigma {E_x}, $ | (19) |
$ \frac{{\partial {E_x}^2}}{{\partial {z^2}}} - {k^2}{E_x} = 0\;\;\;\;\;\;{\rm{or}}\;\;\;\;\;\;\frac{{\partial {H_y}^2}}{{\partial {z^2}}} - {k^2}{H_y} = 0, $ | (20) |
TM极化模式(
$ \frac{{\partial {H_x}}}{{\partial z}} = \sigma {E_y}, $ | (21) |
$ - \frac{{\partial {E_y}}}{{\partial z}} = i\omega \mu {H_x}, $ | (22) |
$ \frac{{\partial {H_x}^2}}{{\partial {z^2}}} - {k^2}{H_x} = 0\;\;\;\;\;\;{\rm{or}}\;\;\;\;\;\;\frac{{\partial {E_y}^2}}{{\partial {z^2}}} - {k^2}{E_y} = 0。$ | (23) |
式中:
首先以TE模式为例,给定的电场
$ \frac{{\partial {E_x}^2}}{{\partial {z^2}}} - i\omega \mu \sigma (z){E_x} = 0。$ | (24) |
电场所满足的边界条件为:
1)在上边界AB处假设
2)下边界CD以下为均质岩石,电磁波在CD以下将按照负指数衰减,即
该边值问题所对应的变分问题为:
$ F({E_x}) = \int_{z1}^{zN} {\left[ {{{\left( {\frac{{\partial {E_x}}}{{\partial z}}} \right)}^2} - i\omega \mu \sigma {{\left( {{E_x}} \right)}^2}} \right]} {\rm{d}}z + a{E_x}^2|{z_N}, $ | (25) |
$ \delta F({E}_{x})=0。$ | (26) |
采用有限元方法进行计算,首先应将区域剖分为若干网格,在网格内电导率必须是连续的。在保证计算精度的前提下,采用三次插值。最终,计算出节点处的电场值[5-6]。
2.2 TM模式磁场
$ \frac{\partial }{\partial z}\left(\frac{1}{\sigma (z)}\frac{\partial {H}_{x}}{\partial z}\right)+i\omega \mu {H}_{x}=0。$ | (27) |
不管海水介质电性如何分布,空气中的
本文的主要研究内容为海水中极低频电磁波的衰减,因此构建为3层模型,第1层为海水层,水深400 m,电阻率0.25
根据前文的公式及边界条件,经过有限元算法代码实现,得到TE和TM模式下,4 Hz,17 Hz,88 Hz和178 Hz的电、磁场随深度的衰减情况,其中Ex0代表海面处的水平电场值,Ex代表某一深度处的水平电场值,Hy0代表海面处的水平磁场值,Hy代表某一深度处的水平电场值,如图2和图3所示。
1)从TE和TM模式均可以看出,极低频水平电场和磁场在海水中的衰减特征均表现为频率越高,衰减越快,同时海水中磁场衰减快,电场相对衰减慢。
2)水平电场分量
极低频电磁波在海水中传播,海底的海床地质条件也会影响到电磁波在海水中的传播特性[7-9]。为进一步研究TE和TM极化模式下海底介质电阻率值变化对海水中极低频电磁波衰减造成的影响,为此以0.5 Hz为例,改变海底岩石电阻率值和厚度值,如表1所示。再据此参数表进行仿真计算,得到0.5 Hz时水平电场分量和水平磁场分量随深度的的衰减特征,如图4所示。
由图4和图5可以看出:1)当模型参数发生变化时,TE和TM模式的极低频电、磁场衰减特征几乎相同;2)从TE和TM模式的水平电场和磁场分量计算结果均看出,模型1和模型3的计算结果一样,在图上出现曲线重合的现象。说明海水下低阻岩石的厚度增加对电场、磁场的衰减几乎无影响;3)模型2计算的结果相比较模型1和模型3发现,当低阻岩层电阻率增大时,电场衰减速度变慢,而磁场衰减速度变快;4)频率降低,变为0.5 Hz时,电场衰减速度变慢的现象更加明显,而磁场这一现象并不明显。
通过对TE和TM模式的极低频电磁波传播过程对比分析发现:
1)在海水中2种模式下的水平电场分量和磁场分量衰减特性基本一致,说明极化模式对极低频电磁波的传播几乎无影响;
2)仿真结果表明,当信号频率越低时,电场和磁场分量衰减速度越慢,这和电磁波传播特征完全一致;
3)当频率在0.5 Hz时,磁场衰减速度明显快于电场,同时当海底岩石电阻率值越大,电场衰减速度变慢,而磁场衰减速度变快。因此在实际布设中可以尝试考虑利用磁力梯度仪开展水平磁场分量在垂向上的信号采集工作。
[1] |
卓贤军, 陆建勋. “极低频探地工程”在资源探测和地震预测的应用与展望[J]. 舰船科学技术, 2010, 32(6): 3-7+30. |
[2] |
卓贤军, 张佳炜. 极低频电磁波在海水中的传播特性[J]. 舰船科学技术, 2011, 33(6): 117-120. |
[3] |
柳建新. 大地电磁测深法勘探: 资料处理, 反演与解释 [M]. 北京: 科学出版社, 2012.
|
[4] |
底青云. “地-电离层”模式电磁波传播特征研究 [M]. 北京: 科学出版社, 2013.
|
[5] |
徐世浙. 地球物理中的有限单元法[M]. 北京: 科学出版社, 1994.
|
[6] |
陈润滋, 安志国, 杨良勇. 基于Matlab语言的二维大地电磁OCCAM快速反演[J]. 地球物理学进展, 2018, 33(4): 1461-1468. |
[7] |
刘胜道, 周骏, 龚沈光. 海水中电磁波的弯曲[J]. 海军工程大学学报, 2001, 13(4): 75-78. DOI:10.3969/j.issn.1009-3486.2001.04.017 |
[8] |
陈聪, 周骏, 龚沈光. 海水中电磁波传播特性的研究[J]. 海军工程大学学报, 2004, 18(2): 61-64. DOI:10.3969/j.issn.1009-3486.2004.02.016 |
[9] |
许建荣, 柳建新, 李爱勇, 等. 水域大地电磁测深中水体对电磁场的影响[J]. 中南大学学报(自然科学版), 2007, 38(3): 567-573. |