海浪背景下艦船尾流感應磁場特性

魏寶君 陳濤 劉健 李靜

摘要 ：基于麥克斯韋方程組以及海浪和尾流的速度場表達式得到無限深海水中海浪背景下艦船尾流感應電磁場的表達式，根據分界面處感應電磁場的連續性條件確定表達式中的待定系數。利用所得到的表達式分析無限深海水中尾流感應磁場、海浪感應磁場和海浪背景下艦船尾流感應磁場的特性。結果表明：艦船尾流感應磁場各分量的幅度在靠近海面的海水中均存在極大值點，且這些極大值點的數值和對應的海水深度均隨船速的增大而增加；在海水中某一深度處艦船尾流感應磁場主要集中分布在一個“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 差別相對較?。▓D2）， 故對艦船尾流施加海風的影響后 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.

（編輯 修榮榮）