海浪背景下舰船尾流感应磁场特性

魏宝君 陈涛 刘健 李静

摘要 ：基于麥克斯韦方程组以及海浪和尾流的速度场表达式得到无限深海水中海浪背景下舰船尾流感应电磁场的表达式，根据分界面处感应电磁场的连续性条件确定表达式中的待定系数。利用所得到的表达式分析无限深海水中尾流感应磁场、海浪感应磁场和海浪背景下舰船尾流感应磁场的特性。结果表明：舰船尾流感应磁场各分量的幅度在靠近海面的海水中均存在极大值点，且这些极大值点的数值和对应的海水深度均随船速的增大而增加；在海水中某一深度处舰船尾流感应磁场主要集中分布在一个“V”形平面区域内，且船速越大，“V”形区域内感应磁场周期性振荡的幅度越大，沿尾流运动方向振荡的空间间隔越大；海浪对舰船尾流感应磁场 z分量的影响大于对x 分量的影响。风向对尾流感应磁场的影响存在周期为π的周期性；当风速小于船速时，海浪对舰船尾流感应磁场的影响较弱，在“V”形区域内尾流感应磁场占主导作用；当风速继续增加时，海水中海浪感应磁场的强度显著增大，对尾流感应磁场的影响明显增强，在“V”形区域内尾流感应磁场逐渐被海浪感应磁场所覆盖；与风速变化所产生的影响差异比较明显相比，风向变化对尾流感应磁场影响的差别相对要小。

关键词 ：舰船尾流； 海浪； 感应磁场； 三维海水运动

中图分类号 ：P 631.9    文献标志码 ：A

引用格式 ：魏宝君，陈涛，刘健，等.海浪背景下舰船尾流感应磁场特性［J］.中国石油大学学报（自然科学版），2024，48（1）：104-114.

WEI Baojun， CHEN Tao， LIU Jian， et al. Characteristics of  ship' s wake induction magnetic field under  background of

waves［J］. Journal of China University of Petroleum（Edition of Natural Science）， 2024，48（1）：104-114.

Characteristics of  ship' s wake induction magnetic

field under  background of waves

WEI Baojun， CHEN Tao， LIU Jian， LI Jing

（College of Science in China University of Petroleum （East China）， Qingdao 266580， China）

Abstract ：The expression detailing the ship' s wake-induced electromagnetic field in the presence of waves within infinite deep sea water was obtained based on Maxwell' s equations and the velocity field expressions of waves and wakes. The coefficients within the derived expression were determined by ensuring the continuity condition of induced electromagnetic fields at the interface. The obtained expressions were then employed to analyze the distinct characteristics of the ship' s wake-induced magnetic field， the wave' s induction magnetic field and the ship' s wake-induced magnetic field in the presence of waves. The findings revealed that the amplitude of each component of the ship' s wake-inducted magnetic field reaches the peak value near the sea surface， and both the peak value and the corresponding sea depth increase with the ship' s speed. The ship' s wake-induced magnetic field is mainly distributed over a "V" shaped plane area at a certain depth in the seawater. Higher ship speeds correlate with the increased amplitude of the periodic vibration of the induced magnetic field， and larger spatial interval of the vibration along the direction of the ship' s wake movement within the "V" shaped area. The influence of waves on the  z -component of the ship' s wake-induced magnetic field is greater than that on the  x -component. The influence of wind direction on the ship' s wake-induced magnetic field is periodic with a period of π. When the wind speed is smaller than the ship speed， the influence of waves on the ship' s wake-induced magnetic field is weak， and the wake-induced magnetic field plays a dominant role in the "V" shaped area. When the wind speed increases continuously， the strength of the wave' s induction magnetic field in the seawater increases significantly and its influence on the ship' s wake induced magnetic field is significantly enhanced， and the ship' s wake-induced magnetic field is gradually covered by the wave' s induction magnetic field in the "V" shaped area. Compared to the more pronounced influence caused by changes in wind speed， the impact on the ship' s wake-induced magnetic field produced by changes in wind directions is relatively small.

Keywords ：ship' s wake; waves; induced magnetic fields; three-dimensional sea water movements

海洋电磁探测作为一种重要的海洋地球物理探测技术，在海洋油气资源勘探和海底地质结构研究等领域具有广阔的应用前景。在海洋电磁探测过程中，常常以舰船为载体释放与回收海底电磁采集站。舰船在海水中行驶会产生尾流，海面也经常产生海浪，它们都会因为切割地磁场而产生感应电磁场，此类感应电磁场是海洋电磁探测的主要噪声来源  ［1-2］ 。此外，随着科技的发展，舰船的隐身性能逐渐受到重视，从而使舰船的特征信号很难捕捉。当隐身舰船的航速低于一定数值时，很难利用声呐技术对其进行探测和跟踪，因此发展非声波探测手段就变得尤为重要。在此情况下舰船尾流因为难以消除成为新的舰船信号源，而舰船在海水表面行驶时由于是以海浪作为背景，因此海浪信号是尾流信号的主要噪声之一。舰船尾流探测技术主要包括声特性探测、电磁特性探测、光特性探测和热特性探测等。海水水域环境、海水表面风速、舰船航速和舰船几何尺寸等因素均可对舰船尾流产生影响  ［3］ ，且不同条件下舰船尾流的电磁学、光学、声学和热学等特性截然不同  ［4］ 。在舰船尾流探测和舰船定位中，电磁探测具有其独特的优势，如识别能力强、效率高、定位精度高、命中率高和成本低等。综上所述，研究尾流感应电磁场和海浪感应电磁场的特性以及海浪背景下尾流感应電磁场的特性，对于海洋电磁探测和舰船尾流电磁探测均具有重要理论意义和实际应用价值。海浪感应电磁场的研究始于二十世纪五、六十年代。Weaver  ［5］ 基于表面波理论建立了二维无限深海水中海浪感应磁场模型，研究结果表明海浪感应磁场大小与海浪波高和波长有关。Groskaya等  ［6］ 推导出了浅海表面重力波感应磁场的解析解，张自力等  ［7］ 分析了二维无限深海水中海浪感应磁场的频谱特征。朱晓建等  ［8］ 将不同频率的海浪感应磁场叠加，得到了三维无限深多频海浪感应磁场的响应。目前关于三维无限深海水中海浪感应电磁场的研究尚不完善，多数研究针对的是单频海浪感应电磁场。事实上海浪是一种多频三维海水运动，针对多频海浪感应电磁场计算方法的研究及特性分析需要继续深入。在舰船尾流感应电磁场研究方面，Dan  ［9］ 建立了无限深水面船舶尾流感应电磁场模型，结果表明尾流感应电磁场是可测的。张伽伟等  ［10］ 模拟了浅海中水面船舶产生的感应电磁场，Zou等  ［11］ 研究了航空磁测尾流感应磁场的可行性，Yaakobi等  ［12］ 推导出了有限深水面船舶感应磁场在空气层中的表达式。王向坤等  ［1］ 考虑了有限深海浪以及海底介质电导率的影响，丰富了浅海中船舶尾流感应电磁场模型。笔者在前人研究的基础上推导出了无限深海水中海浪和尾流感应电磁场新的解析表达式，进一步完善了无限深海水中海浪感应电磁场和尾流感应电磁场的计算方法。 在海浪感应电磁场模拟中考虑到海浪的多频性和随机性，将各种频率随机初相位的单频海浪感应电磁场线性叠加，得到多频海浪感应电磁场的响应，并给出空气层和海水层海浪感应电磁场的表达式。在推导舰船尾流的感应电磁场时借鉴海浪感应电磁场的推导思路。利用上述推导结果针对无限深海水中尾流和海浪感应磁场各自的传播特性以及无限深海水背景下海浪对尾流感应磁场的影响进行计算和分析，并研究风速和船速变化对舰船尾流感应磁场的影响规律。

1  无限深海水中海浪和尾流感应电磁场理论

为了研究的方便，针对运动中的水面舰船建立如图1所示的直角坐标系。 设z方向为垂直方向， z= 0代表海面，z>0的区域为空气层，z<0的区域为海水层，并假定舰船以恒定速度U在海面沿x轴负方向运动，则尾流的运动方向为x轴正方向。对小尺度的海浪和舰船尾流而言，可假定研究区域内的地磁场为固定值，则由图1的几何关系可以得到总地磁场的各分量为

B   E =  B   E  （ cos  I cos  α e  x+ cos  I sin  α e  y- sin  I e  z）.  （1）

式中， B   E 为地磁场的磁感应强度；I为地磁倾角；α为地磁北极与x轴正方向的夹角。

假设海浪和尾流运动的速度分别为 V   w 和 V   s ，忽略海浪和尾流之间的相互作用，则海浪和尾流切割磁场所产生的感应电场可表示为 E   i =（ V   w + V   s ）×（ B + B   E ），此处 B 为感应磁场的磁感应强度。由于地磁场的磁感应强度 B   E 远大于感应磁场的磁感应强度 B ，所以海浪和尾流的感应电场可近似为 E   i =（ V   w + V   s ）× B   E 。尾流和海浪的感应电磁场满足以下麦克斯韦方程组：

× E =-μ   H   t  ， （2）

× H =σ［ E +（ V   w + V   s ）× B   E ］+ε   E   t  . （3）

式中， E 为电场强度； H 为磁场强度；μ、σ、ε分别为介质的磁导率、电导率和介电常数。

由式（2）、（3）可知，若已知任一时刻海洋中某点的海浪速度 V   w （x，y，z，t）和尾流速度 V   s （x，y，z，t）以及 B、E的边界条件，就可以通过计算得到空气和海洋中某点的B、E。因此，求解B、E的前提是对速度函数V   w （x，y，z，t）和 V   s （x，y，z，t）的求解。 假设海水是不可压缩无旋流体，根据不可压缩流动的定义，海水密度不随海水流动发生变化，根据流体的连續性方程有 ·（ V   w + V   s ）=0。利用该条件，由式（2）和（3）可得到尾流与海浪感应磁场所满足的微分方程为

2 H -μσ   H   t -με   2 H   t 2 = 0， z>0;

-σ（  B   E ·  ）（ V   w + V   s ）， z≤0.  （4）

Higgins等  ［13］ 将三维海浪视为由 M×N个振幅为a  mn 、角频率为ω m、波数为k m、随机初始相位为φ  mn 、传播方向和x轴夹角为θ n的正弦波线性叠加而得到，其中M为频率数目、N为传播方向数目，则无限深海水中多频三维海浪速度场可以表示为

V   w （x，y，z，t）=  Re  ∑ M  m=1  ∑ N  n=1

q   w （z） exp ［ i （ω mt-k mx cos  θ n- k my sin  θ n+φ  mn ）］ . （5）

其中

q   w （z）=a  mn ω m exp （k mz）（ i  cos  θ n e  x+ isin  θ n e  y- e  z），

a  mn = 2S（ω m，θ n）ΔωΔθ  ，

S（ω m，θ n）= 8.1×10  -3  4k 3 m  exp   -0.74g 2 k 2 mU 4  19.5   ·

1  π  ［1+ a cos  2（θ n-θ）+b cos  4（θ n-θ）］，

a=0.5+0.82 exp ［-ω 4 m/（2ω 4  p ）］，

b=0.32 exp ［-ω 4 m/（2ω 4  p ）］，

ω  p =0.877g/U  19.5 .

式中，a  mn 为海浪振幅函数，与海浪谱S（ω m，θ n）有关，本文中采用的海浪谱由 PM谱  ［12］ 与Cote  ［14］ 提出的方向谱组成; ω m为角频率，满足频散关系式ω 2 m=gk m;g为重力加速度;φ  mn 为初始相位，随机分布在（0，2 π ）范围内;U  19.5 为海平面上方19.5  m处的风速;

Δ ω和 Δ θ分别为角频率和传播方向的分割间隔;θ为风向与x轴的夹角，称为风向角;θ n分布在海平面风向两侧 π/2 的范围内。

舰船的尾流速度场可以表示为  ［10］

V   s （x，y，z，t）=

Re  ∫   π/2    - π/2   q   s （z） exp ［ i （ω 0t- k 0x cos  θ-k 0y sin  θ）］ d θ . （6）

其中

q   s （z）=S（θ） exp （k 0z）（ i  cos  θ e  x+ isin  θ e  y- e  z），

S（θ）=-  2i k 2 0U π    Y （ x′，z′ ）   x′  exp （k 0z′+  i k 0x′ cos  θ） d x′ d z′，

ω 0=k 0U cos  θ， k 0=g  sec   2θ/U 2.

式中，ω 0为尾流速度场角频率；k 0为尾流波数；Y（x′，z′）为船形函数，与船舶的几何特性有关。根据船体理论，细长型船的船形函数用如下公式  ［9］ 表示：

Y（x′，z′）= B 2  1+ z′ D   1+b 1 z′ D +b 2 z′ 2 D 2   1- x′ 2 l 2  ×  1+c 1 x′ l +c 2 x′ 2 l 2  ，

x′ ≤l，-D≤z′≤0. （7）

式中，B为船宽；l为半船长；D为舰船吃水深度；b 1、b 2、c 1、c 2为权重系数，与船体主尺度（包括船长、船宽和吃水深度等参数）有关，一般都是小于1的常数。为了简化计算，本文中仅考虑舰船主尺度对尾流特性的影响。对于细长型船体，假设舰船是抛物线型，其船型函数可表示为  ［8］

Y（x′）= B 2  1- x′ 2 l 2  ，  x′ ≤l，-D≤z′≤0. （8）

下面分别求解尾流和三维海浪感应电磁场所满足的微分方程（4）。尾流感应电磁场是由尾流速度场引起的，由于尾流速度场是简谐波，根据线性波动理论，其感应电磁场也一定具有简谐波函数的形式  ［5］ ，故尾流感应电磁场可表示为

H   s （x，y，z，t）= Re  ∫   π/2    - π/2   h  θ（z） exp ［ i （ω 0t- k 0x cos  θ-k 0y sin  θ）］ d θ ，  （9）

E   s （x，y，z，t）= Re  ∫   π/2    - π/2   e  θ（z） exp ［ i （ω 0t-k 0x cos  θ-k 0y sin  θ）］ d θ . （10）

类似地，海浪感应电磁场可以表示为

H   w （x，y，z，t）= Re  ∑ M  m=1  ∑ N  n=1   h   θ ， mn （z） exp ［ i （ω mt- k mx cos  θ n-k my sin  θ n+φ  mn ）］ ， （11）

E   w （x，y，z，t）= Re  ∑ M  m=1  ∑ N  n=1   e   θ ， mn （z） exp ［ i （ω mt- k mx cos  θ n-k my sin  θ n+φ  mn ）］ . （12）

式中， h  θ（z）和 e  θ（z）为尾流感应磁场与感应电场表达式中的待求解部分，

h   θ，mn （z）和 e   θ，mn （z）为海浪感应磁场与感应电场表达式中的待求解部分，均需要基于速度场通过求解麦克斯韦方程组得到。

由于尾流速度场和海浪速度场与它们激发的感应电磁场具有相同的简谐波函数形式，且感应电磁场指数项中的随机初始相位φ  mn 不參与运算，所以可以统一令

q = q （z） exp ［ i （ωt-kx cos  θ-ky sin  θ）］， （13）

h = h （z） exp ［ i （ωt-kx cos  θ-ky sin  θ）］， （14）

e = e （z） exp ［ i （ωt-kx cos  θ-ky sin  θ）］. （15）

式中， q （z）为 q   s （z）或 q   w （z）； h （z）为尾流和海浪感应磁场表达式中待求解部分 h  θ（z）或 h   θ，mn （z）； e （z）为尾流和海浪感应电场表达式中待求解部分 e  θ（z）或 e   θ，mn （z）；角频率ω为尾流速度场的角频率ω 0或海浪的角频率ω m；波数k为尾流的波数k 0或海浪的波数k m。

将式（14）代入式（4），可得到深海条件下空气和海水中感应磁场微分方程的通解为

h （z）=  a  1 exp （-δ 0z）， z>0;

a  2 exp （δ 1z）-  P  k 2-δ 2 1 ， z≤0.  （16）

式中，  a  1、  a  2为待定系数；δ 2 i=k 2+ i ωμσ′ i，σ′ i=σ i+ i ωε i（i=0，1）为复电导率，下标i=0，1分别对应空气层和海水层；

P 为包含速度场和地磁场信息的源项，其表达式为

P =［σ 1（ B   E · ） q ］ exp ［- i （ωt-kx cos  θ-ky sin  θ）］.

将式（15）和（16）代入式（3），可得到感应电场的表达式为

e （z）=  （ × h ） exp ［- i （ωt-kx cos  θ-ky sin  θ）］ σ′  0   ， z>0;

（ × h ） exp ［- i （ωt-kx cos  θ-ky sin  θ）］- Q  σ′  1   ，  z≤0.  （17）

其中

Q =σ 1（ q × B   E ） exp ［- i （ωt-kx cos  θ-ky sin  θ）］.

式中， Q 为另一个包含速度场和地磁场信息的源项。

利用空气和海水分界面处电磁场的连续性条件可求得待定系数

a 1、 a 2， 进而由式（16）和式（17）得到海浪以及尾流感应电磁场的表达式。考虑到空气层和海水层均可视为非磁性介质，其磁导率相同，则分界面处感应磁场的各分量连续，由式（16）可得

a  2-  P    z=0  k 2-δ 2 1 = a  1. （18）

根据均匀介质中感应磁场的散度为零的条件，由式（14）可得

d  d z （ h · e  z）= i k h · e  ρ. （19）

其中

e  ρ= cos  θ e  x+ sin  θ e  y.

将式（16）中的第二个表达式代入式（19），并考虑到源项 P 的特点，利用正交关系和坐标变换可得

a  2= a  2z δ 1  i k （ cos  θ e  x+ sin  θ e  y）+a  2z  e  z. （20）

在海水分界面处，感应磁场的垂直分量以及该分量对z的偏导数连续，由式（16）可得

a  1z =a  2z - P z   z=0  k 2-δ 2 1  ， （21）

-δ 0a  1z =δ 1a  2z - （ P z/ z）   z=0  k 2-δ 2 1  . （22）

由式（21）、（22）解得

a  2z = δ 0P z   z=0 +（ P z/ z）   z=0  （k 2-δ 2 1）（δ 0+δ 1）  . （23）

将式（23）代入式（20）得到  a  2，再由式（18）可得到 a  1。

对海浪感应电磁场，需要将上述表达式中的k和ω替换为海浪的波数k m和角频率ω m，将源项 P和Q替换为由海浪速度场q   w （z）求取的 P   w 和 Q   w 后进行叠加运算，从而求得海浪感应电磁场的响应。对尾流感应电磁场，需要将表达式中的k和ω替换为尾流波数k 0和角频率ω 0，将源项 P和Q替换为由尾流速度场q   s （z）求取的 P   s 和 Q   s 后对表达式进行积分运算，从而求得尾流感应电磁场的响应。

2 海浪和尾流感应磁场的数值模拟

在进行数值模拟时假设空气和海水的电导率分別为σ 0=0  S·m  -1  和σ 1=5  S·m  -1  ，空气的相对介电常数为1，海水的相对介电常数为81。取地磁场磁感应强度为  B   E =5×10 4 nT，地磁倾角 I= π /6，地磁北极与x轴正方向的夹角α= π /10。

2.1 尾流感应磁场

假定半船长 l=20  m ，船宽B=5  m ，舰船吃水深度D=1.5  m 。当船速U取9、10和11  m/s 三种不同的数值时，海平面坐标x=1  m 、y=1  m 处尾流感应磁场水平分量的大小 B x

和垂直分量的大小 B z 随海水深度的变化关系如图2所示。由图2 （a） 可以看出，舰船尾流感应磁场的水平分量 B z 有两个极值点，分别出现在海面z=0处和海水表面下方的海水中。尾流感应磁场水平分量的符号在稍低于海水表面的海水中发生改变，其大小先随海水深度的增加从海面z=0处的最大值快速减小至零，然后随海水深度的进一步增加而增大到反向极大值点，在达到极大值点后再随海水深度的增加而逐渐衰减至零。由图2 （b） 可以看出，海水中尾流感应磁场的垂直分量 B z 先随海水深度的增加从海面处的某一数值快速达到极大值，然后随海水深度的进一步增加而呈现较为缓慢的指数衰减。对比 图2（a）和（b），尾流感应磁场水平分量在海水中的符号改变点（即零值点）与垂直分量在海水中的极大值点对应的深度相同，该深度和尾流感应磁场水平分量在海水中的极大值点对应的深度均随船速的增大而增加，

只不过后一个深度随船速的改变更加明显。另由图2可以看出，海面上方空气中尾流感应磁场的水平分量和垂直分量均随高度的增加呈指数衰减，且在空气中比在海水中衰减得更快。舰船尾流感应磁场 B x与B z这两个分量的幅度差别不大，其数值均在同一数量级。在其他条件相同的情况下，舰船速度越大，其尾流感应磁场的幅度越大。

取船速 U仍分别为9、10和11  m/s，海面下方1 m处舰船中心线上（ y=0、z=-1  m ）尾流感应磁场的水平分量B x和垂直分量B z随舰船后方坐标x的变化关系如图3所示。由图3可以看出，不同船速情况下B x和B z均随x的变化呈周期性振荡，B z的振荡幅度均稍大于B x的振荡幅度，且距离舰船越远它们的振荡幅度越小。另由图3可知，船速U越大，B x和B z的振幅越大。

取船速 U为10  m/s，图4给出了 t=0时刻海面下方z=-1  m 处B x和B z在舰船后方xy平面上的分布情况。由图4可以看出，尾流感应磁场的水平分量和垂直分量的空间分布规律相似，且各分量的空间分布和舰船尾流速度的分布具有一致性，都呈现 “V”形分布，这种尾流又称为开尔文尾流  ［15］ ，开尔文尾流的感应磁场主要集中分布在“V”形区域内，由向船后方传播的横波和向四周辐散的扩散波组成。在“V”形区域的边界上（又称为开尔文臂  ［15］ ），扩散波具有局部最大值，在“V”形区域外尾流感应磁场快速衰减。由于尾流感应磁场的幅度在开尔文臂区域存在局域最大值，所以在辨别舰船尾流特征信号时应把开尔文臂作为重点研究区域。

取船速 U分别为 12和14  m/s，图5给出了 t=0时刻海面下方z=-1  m 处感应磁场的水平分量B x在舰船后方xy平面上的分布情况（由图4知，B z分量与B x分量的分布状况相似，故只给出了B x分量的分布图）。 对比图5和图4（a）可以看出，在其他条件不变的情况下，船速越大，舰船尾流感应磁场的幅度越大，“V”形区域内感应磁场周期性振荡沿尾流运动方向的空间间隔越大。

2.2 海浪感应磁场

在模拟海浪产生的感应磁场时取风向角 θ= π/4，海浪的频率范围为0.02～1 Hz，所选择的频率个数和方向分割数目为 M=N=20 。取海平面上方19.5 m处的风速 V（V=U  19.5 ）分别为8、10和12  m/s，图6给出了海平面坐标 x =1 m、 y =1 m处海浪感应磁场不同分量的大小随海水深度的变化关系。对比图6与图2可以发现，尽管海浪产生的感应磁场与舰船尾流产生的感应磁场具有完全不同的来源，但两种感应磁场随海水深度 z的 变化规律相似。对比图6（a）和（b）可以看出，与图2所示舰船尾流感应磁场不同的是，海水中海浪感应磁场垂直分量的极大值远大于其水平分量的极大值。

2.3 海浪背景下舰船尾流感应磁场

下面分析无限深海水中有海浪的背景下舰船尾流产生的感应磁场。首先固定风速 V=10  m/s，风向角 θ =π/4，取船速 U 分别为9、10和11 m/s，海平面坐标 x =1 m、 y =1 m处尾流感应磁场水平分量的大小

B x 和垂直分量的大小 B z 随海水深度的变化關系如图7 （a）和（b）所示。其次固定船速 U =10 m/s，风向角保持不变，取风速 V 分别为8、10和12 m/s，海平面坐标 x =1 m、 y =1 m处尾流感应磁场  B x和 B z 随海水深度的变化关系如图7 （c）和（d）所示。

对比图7（a）与图2（a）可以发现，施加海风的影响后，舰船尾流 感应磁场的水平分量 B x 并没有明显的变化。对比 图7（b）与图2（b）可以看出，施加上述方向海风的影响后，舰船尾 流感应磁场的垂直分量 B z 显著增大，变化明显，这说明风速对舰船尾流感应磁场 B z 分量的影响要远大于对 B x 分量的影响。上述结论可以通过对比 图6（a）和（b）得到解释，由于海风引起的海浪感应磁场的垂直 分量 B z 远大于水平分量 B x ，而舰船尾流感应磁场的 B x 与 B z 差别相对较小（图2）， 故对舰船尾流施加海风的影响后 B z 变化明显。也正是由于上述原因， 图7（c）与图6（a）相比变 化明显且不同风速下 B x 的曲线差别较小（对 B x 的主要贡献来自于舰船速度），而图 7（d）与图6（b）相比则变化不大 （对 B z 的主要贡献来自于风速）。另由图7可以看到，在正常的船速和风速取值范围内，海浪背景下舰船尾流感应磁场的数值一般不会超过10  nT。

取船速 U =10 m/s、风速 V =10 m/s，图8和9分别给出了不同风向角情况下 t=0时刻海面下方 z= -1  m处舰船尾流感应磁场的水平分量 B x和垂直分量B z在舰船后方xy平面上的分布情况。由式（5），由于风向对尾流感应磁场的影响存在周期为 π 的周期性， θ=  π 与θ=0的仿真结果完全相同，故图8和9没有给出θ= π的分布。另外，计算发现，风向角为3π/4时的仿真结果与风向角为π/4时的仿真结果相对于 y=0呈对称分布，即风向角为 π/4时尾流感应磁场在图8和9上半区域的分布与风向角为3π/4时尾流感应磁场在图8和9下半区域的分布相同，反之亦然，故图8和9也没有给出 θ= 3π/4的分布。

与图4相比，由于海风引起的海浪的存在，在图8和9的“V”形区域尤其是开尔文臂区域内，尾流感应磁场因受到海浪感应磁场的影响而产生更加剧烈的震荡，而在“V”形区域之外主要是海浪产生的感应磁场，该磁场在 xy平面上基本呈等幅度变化。另外，对比图9与图8可以看出，海浪感应磁场垂直分量B z的振荡较水平分量B x的振荡剧烈，其对 “V”形区域内尾流感应磁场的影响也更大。对比图8和9不同风向角对应的图像可知，风向角为0（此时尾流运动方向与风向相同，舰船运动方向与风向相反，处于逆风行驶状态）或π（此时尾流运动方向与风向相反，舰船运动方向与风向相同，处于顺风行驶状态）即舰船运动方向与风向平行时，海浪对尾流感应磁场的影响最大，此时海浪感应磁场对尾流感应磁场所在的全部“V”形区域都有比较明显的影响。当风向角为π/2即舰船运动方向与风向垂直时，海浪感应磁场对尾流感应磁场的影响最小。当风向角为π/4时，海浪感应磁场主要分布在图8（b）和图9（b）的上半区域，海浪感应磁场对尾流感应磁场的影响主要表现在“V”形区域的上半部分，由式（5）知这是由于海浪感应磁场的大小沿着风向两侧逐渐衰减造成的。如前所述，风向角为3π/4时的结论正好相反。另由图8和9中不同风向角对应的图像可以看出，在上述船速和风速参数下，尽管不同风向的海浪对舰船尾流感应磁场的影响有差别，但总体而言差别不大。

取船速 U=10  m/s、风向角 θ=0，图10和图11分别给出了不同风速情况下t=0时刻海面下方z=-1  m处舰船尾流感应磁场的水平分量 B x和垂直分量B z在舰船后方xy平面上的分布情况。分别对比图 10（a）与图8（a）、图11（a）与图9（a）可以看出，对于风速为8 m/s（小于10 m/s的船速）的情况，海水中海浪感应磁场的强度更小，影响更弱，此种情况下在“V”形区域内尾流感应磁场占主导作用，尾流感应磁场的整体变化规律与无风速情况下基本一致，因此图10（a）与没有风速影响的图4（a）接近、

图11（a）与没有风速影响的图4（b）接近。由图10（b）和图11（b），当风速增加到12 m/s（超过10 m/s的船速）时，海水中海浪感应磁场的强度显著增大，对尾流感应磁场的影响明显增强。由图10（c）和图11（c），当风速增加到14 m/s时，海水中海浪产生的感应磁场占据主导地位，“V”形区域内尾流感应磁场基本被海浪感应磁场覆盖。

3 结 论

（1） 舰船尾流感应磁场的幅度在靠近海面的海水中均存在极大值点，且这些极大值点的数值和对应的海水深度均随船速的增大而增加，舰船尾流感应磁场在达到极大值点后随海水深度的继续增加呈指数衰减。舰船尾流感应磁场沿尾流运动方向呈周期性振荡，且距离舰船越远振荡幅度越小。在海水中某一深度处舰船尾流感应磁场主要集中分布在一个“V”形平面区域内，且船速越大，“V”形区域内感应磁场周期性振荡的幅度越大，沿尾流运动方向振荡的空间间隔越大。

（2）海浪感应磁场随海水深度的变化关系与尾流感应磁场相似，且风速越大海浪感应磁场的幅度越大。与舰船尾流感应磁场不同的是，海浪感应磁场 z分量的幅度远大于x分量。

（3）海浪对舰船尾流感应磁场z分量的影响大于对x 分量的影响。风向对尾流感应磁场的影响存在周期为π的周期性，当舰船运动方向与风向平行时海浪对尾流感应磁场的影响最大，当舰船运动方向与风向垂直时海浪对尾流感应磁场的影响最小。当风速小于船速时，海浪对舰船尾流感应磁场的影响较弱，在“V”形区域内尾流感应磁场占主导作用；当风速大于船速时，随着风速的继续增加，海水中海浪感应磁场的强度显著增大，对尾流感应磁场的影响明显增强，在“V”形区域内尾流感应磁场逐渐被海浪感应磁场所覆盖。与风速变化所产生的影响差异比较明显相比，风向变化对尾流感应磁场影响的差别相对要小。

参考文献 ：

［1］  王向坤，李予国，李建凯，等.有限水深海浪背景下船舶尾流感应电磁场数值模拟［J］.中国海洋大学学报（自然科学版），2020，50（12）：107-114.

WANG Xiangkun， LI Yuguo， LI Jiankai， et al. Numerical simulation of electromagnetic fields induced by ship wakes in ocean waves of finite depth water［J］. Periodical of Ocean University of China（Edition of Natural Science），2020，50（12）：107-114.

［2］  刘健，王江涛，王子鸣，等.洋流和潮波产生的感应电磁场对比［J］.中国石油大学学报（自然科学版），2022，46（4）：52-62.

LIU Jian， WANG Jiangtao， WANG Ziming， et al. Comparison on induced electromagnetic fields generated by ocean current and tidal wave［J］.Journal of China University of Petroleum （Edition of Natural Science）， 2022，46（4）：52-62.

［3］  冀邦杰，周德善，张建生.基于舰船尾流光效应的制导鱼雷［J］.鱼雷技术，2000，8（3）：44-47.

JI Bangjie， ZHOU Deshan， ZHANG Jiansheng. Guided torpedo based on wake light effect［J］.Torpedo Technology，2000，8（3）：44-47.

［4］  程學峰，杨宝民，冀邦杰，等.水下声成像仿真分析及图像处理［J］.舰船科学技术，2008，30（1）：153-156.

CHENG Xuefeng， YANG Baomin， JI Bangjie， et al. Simulation analysis and image processing of underwater acoustic imaging［J］.Ship Science and Technology，2008，30（1）：153-156.

［5］  WEAVER J T. Magnetic variations associated with ocean waves and swell ［J］. Journal of Geophysical Research Atmospheres， 1965，70（8）：1921-1929.

［6］  GROSKAYA Y M， SKRYNNIKOV R G， SOKOLOV G V. Magnetic field variations induced by the motion of sea waves in shallow water［J］.Geomagnetism and Aeronomy， 1972，12（9）：131.

［7］  张自力，魏文博，李庚伟，等.海浪感应磁场的频谱特性［J］.渤海大学学报（自然科学版），2006，27（4）：354-357.

ZHANG Zili， WEI Wenbo， LI Gengwei， et al. Spectral character of magnetic field generated by ocean waves［J］. Journal of Bohai University （Natural Science Edition）， 2006，27（4）：354-357.

［8］  ZHU X， XIA M. Magnetic field induced by wake of moving body in wind waves［J］. Progress in Electromagnetics Research， 2014，149（1）：109-118.

［9］  MADURASINGHE D. Induced electromagnetic fields associated with large ship wakes［J］. Wave Motion， 1994，20（3）：283-292.

［10］  张伽伟，姜润翔，龚沈光.浅海中船舶尾流产生的感应电磁场［J］.哈尔滨工程大学学报，2014，35（8）：931-935.

ZHANG Jiawei， JIANG Runxiang， GONG Shenguang. Study of the electric field induced by the wake of a moving ship ［J］. Journal of Harbin Engineering University， 2014，35（8）：931-935.

［11］  ZOUN，NEHORAI A. Detection of ship wakes using an airborne magnetic transducer［J］.IEEE Transactions onGeoscience and Remote Sensing，2000，38（1）：532-539.

［12］  YAAKOBI O， ZILMAN G， MILOH T.Detection of the electromagnetic field induced by the wake of a ship moving in a moderate sea state of finite depth［J］. Journal of Engineering Mathematics， 2011，70（1）：17-27.

［13］  LONGUET-HIGGINS M S， STERN M E， STOMMEL H M. The electrical field induced by ocean currents and waves， with applications to the method of towed electrodes［J］. Papers in Physical Oceanography and Meteorology， 1954，13（1）：1-37.

［14］  張宝强.基于小波分析的海洋Mt资料海浪感应电磁噪声压制方法研究［D］.青岛：中国海洋大学，2018.

ZHANG Baoqiang. Removal of induced electromagnetic noise due to ocean waves in marine MT data by using wavelet analysis method ［D］. Qingdao： Ocean University of China， 2018.

［15］  孙荣庆.海面舰船尾迹电磁散射研究［D］.西安：西安电子科技大学，2013.

SUN Rongqing. Research on electromagnetic scattering from ship wakes［D］. Xi' an： Xidian University， 2013.

（编辑 修荣荣）