陳 醒，岳瑞永，姜楷娜，趙 哲

(大連測控技術(shù)研究所，遼寧 大連 116013)

0 引 言

海洋環(huán)境電場作為艦船測試的主要背景場有著重要的研究價值。海洋中存在多種電場效應(yīng)，根據(jù)場源的形式可分為天然電場和人為因素形成的電場。天然電場主要包括海洋大地電場、海水在地場中運動產(chǎn)生的感應(yīng)電場兩部分；人為因素形成的電場主要來源于沿海工業(yè)設(shè)施產(chǎn)生的工頻干擾。為了充分掌握海洋環(huán)境電場特性，支撐艦船電場信號檢測，本文基于歷史觀測的海洋環(huán)境電場數(shù)據(jù)，對海洋環(huán)境電場時頻典型特性及空間相關(guān)性進行分析，確定了典型頻段環(huán)境電場主要特征，可為艦船信號提取提供依據(jù)[1–4]。

國內(nèi)外學者對海洋環(huán)境電場做了相關(guān)研究。蔡岳豐等[5]根據(jù)點電極、電偶極子等不同種類模型對海洋環(huán)境電場分類進行論述。崔培[6]利用三分量水下電場測量裝置對大連附近海域進行觀測，得到了海洋環(huán)境電場幅值、峰值的主要分布特征。Hermand[7]對海洋中產(chǎn)生電場的各種場源從頻率和振幅特征進行了定性分析。Lil-ley[8]研究了澳大利亞南部海域的海浪電場功率譜情況。瑞典國防研究局人員利用三分量水下電場傳感器對瑞典西海岸海域的海洋環(huán)境水下電場進行了長達10個月的觀測，測量覆蓋了數(shù)MHz～20 Hz頻率范圍，觀測到了地磁脈動和雷電產(chǎn)生的舒曼諧振現(xiàn)象[9]。從上述研究成果來看，目前國內(nèi)外環(huán)境電場研究主要限于時頻域的分析，對不同空間位置的電場相關(guān)性研究相對較少。本文將從時頻域、極化方向、空間相關(guān)性多角度研究海洋環(huán)境電場，初步得到淺海海域環(huán)境電場空間相關(guān)性的基本特征。

1 海洋環(huán)境電場數(shù)據(jù)預處理

利用水下電場測量系統(tǒng)對北黃海海域進行長時間觀測。該系統(tǒng)由2個水下測量節(jié)點組成，通過采用分布式采集和有纜式傳輸方式，可同步獲取海底多節(jié)點的水下電場矢量測量數(shù)據(jù)，并可獲取測量體方位、姿態(tài)等信息參數(shù)。環(huán)境電場數(shù)據(jù)取自2019年7月30日下午4點，測量環(huán)境為海況2級，海面5～6級風，電導率為4.55 S/m，水溫為23.3 ℃，試驗海域水深10 m，兩測量體間距為20 m，選取測量數(shù)據(jù)中典型時間段的樣本進行分析。

為了便于不同測量節(jié)點數(shù)據(jù)的統(tǒng)一分析，通過坐標系轉(zhuǎn)換將傳感器測量數(shù)據(jù)轉(zhuǎn)換到北東坐標系，方向為地理北方向，方向為地理東方向。另外考慮到海洋環(huán)境電場多源特征，基于相關(guān)理論分析結(jié)果，根據(jù)場源不同大致可分為海流信號、海面波浪信號、地磁擾動、雷電干擾等，將電場信號劃分為DC～0.1 Hz，0.1 Hz～1 Hz，1 Hz～5 Hz，5 Hz～10 Hz這 4 個頻帶,對這 4 個頻段采用不同頻率分辨率進行精細化頻譜分析。

2 環(huán)境電場信號的時頻特性、矢量極化方向特征分析

經(jīng)過預處理得到的環(huán)境信號單獨從時域或頻域觀察信號不能充分體現(xiàn)環(huán)境電場信號特征，只有從時間和頻率以及極化方向的角度同時表征信號，才能更全面、更細致地反映信號的特性，因此需要對信號進行時域、頻域、極化特性分析。

2.1 信號的時域特征分析

環(huán)境電場幅值是測量得到的最基本物理量，其變化規(guī)律反映了海洋環(huán)境電場最普遍的特征。圖1為海洋環(huán)境電場不同測量體不同頻段不同分量的典型樣本數(shù)據(jù)時域分析圖。

圖1 不同測量體、不同頻段、不同分量樣本數(shù)據(jù)時域分析曲線Fig. 1 Time domain analysis curves of sample data of different measuring bodies, different frequency bands,and different components

由圖可知，雖然時域低頻電場信號隨機性都很強，但是其幅值振動的情況是較為穩(wěn)定，其中DC～0.1 Hz，0.1 ～1 Hz頻段中2個測量體的環(huán)境電場信號有一定相位差，但振幅趨勢一致，DC–0.1 Hz相關(guān)性較差，0.1 ～1 Hz相關(guān)性較好，1 ～5 Hz兩環(huán)境電場區(qū)別較大，5 ～10 Hz兩環(huán)境電場數(shù)據(jù)重合度較高，相關(guān)性較好。為了研究以上現(xiàn)象產(chǎn)生的原因，需要進行環(huán)境電場信號的頻域特征分析。

2.2 信號的頻域特征分析

與環(huán)境時域特征分析相對應(yīng)，信號的頻域特征需要通過計算信號功率譜進行分析與研究。由于不同測量體電場分量、電場分量的頻譜特征基本一致，僅對比不同頻段下環(huán)境電場分量的頻譜圖，得出環(huán)境電場頻域特征。選取測量數(shù)據(jù)中長400 s的樣本進行頻域分析，頻域分辨率接近0.001 Hz，分析曲線如圖2所示。

圖2 不同頻段樣本數(shù)據(jù)頻域分析曲線Fig. 2 Frequency domain analysis curves of sample data of different measuring bodies, different frequency bands, and different components

從圖 2 可知，DC～0.1 Hz、0.1 ～1 Hz、1 ～5 Hz這3個頻段均存在明顯的譜峰特征，其中DC～0.1 Hz,0.1 –1 Hz頻段的電場峰值頻率分別為0.003 Hz，0.15 Hz左右。經(jīng)分析是海洋涌浪和表面波浪運動感應(yīng)產(chǎn)生的電場[10]。在1 ～5 Hz的環(huán)境電場頻率峰值出現(xiàn)在1.5 Hz左右，經(jīng)分析是地磁擾動引起的電場變化以及一些環(huán)境隨機信號[11]。5～10 Hz的環(huán)境電場頻率峰值出現(xiàn)在左右，符合舒曼諧振特征。舒曼諧振基頻以7.5 Hz為中心，在地球表層與電離層之間形成諧振腔，對雷電某些特定頻點的電磁波能量具有放大作用[9]。從圖2結(jié)合圖1綜合分析，在DC～0.1 Hz，0.1 ～1 Hz這2個頻段易受海洋波浪和涌浪影響，兩環(huán)境電場信號有一定相位差。0.1 ～1 Hz相關(guān)性較好；1 ～5 Hz兩環(huán)境電場相關(guān)性較差，主要是因為該頻段地磁擾動、環(huán)境信號帶來的隨機性。5 ～10 Hz環(huán)境電場相關(guān)性較好，經(jīng)分析是因為“舒曼諧振”成分為該頻段主要能量，該信號具有一定的空間均勻性。

2.3 信號的極化方向特征分析

環(huán)境電場信號具有一定的方向性，本文對不同頻段信號水平分量的極化特性進行研究，如圖3所示。

從圖3可知，DC～0.1 Hz測量體2的環(huán)境電場極化方向發(fā)散，與測量體1電場極化方向有一定偏差，相關(guān)性較差；0.1～1 Hz，5 ～10 Hz頻段，兩電場信號存在明顯的極化方向且極化方向一致，相關(guān)性較好；1 ～5 Hz頻段內(nèi)，兩電場信號極化方向相反，且測量體2的環(huán)境電場極化方向發(fā)散，相關(guān)性較差。

3 環(huán)境電場信號歸一化相關(guān)性分析

為了準確給出4個頻段的空間相關(guān)性，研究電場各測量體各頻段各分量之間的內(nèi)在關(guān)聯(lián)，對環(huán)境電場信號進行歸一化相關(guān)性分析是十分重要的，以下給出相關(guān)性定義。

定義信號x(t)和y(t)的互相關(guān)函數(shù)為：

通過對環(huán)境電場互相關(guān)分析可以得到圖4曲線，橫坐標代表的是偏移時間，縱坐標代表2個測量體的歸一化互相關(guān)度值。當互相關(guān)度峰值對應(yīng)的偏移時間越接近零，表示兩環(huán)境電場信號相關(guān)性越好。

圖3 不同測量體、不同頻段、不同分量樣本數(shù)據(jù)極化圖Fig. 3 Polarization diagram of sample data of different measuring bodies, different frequency bands, and different components

4 環(huán)境電場信號的相位補償

相位補償是讓信號不再超前或延時，實現(xiàn)在電場同步下的細微調(diào)節(jié)。從環(huán)境電場相關(guān)性分析可知，兩環(huán)境電場信號相位有所偏差，通過選取環(huán)境電場相關(guān)度峰值點對應(yīng)的偏移時間對環(huán)境電場信號進行相位補償，分析結(jié)果如圖5所示。

圖4 不同分量環(huán)境電場歸一化相關(guān)性曲線Fig. 4 Normalized correlation curves of environmental electric fields with different components

從圖5對比兩測量體、不同分量補償前后的環(huán)境電場信號可以發(fā)現(xiàn)：

圖5 三個頻段的相位補償Fig. 5 Phase compensation for three frequency bands

5 結(jié) 語

本文通過對中國北黃海某淺水海域不同測量體測得的環(huán)境電場數(shù)據(jù)進行了分析。通過坐標系轉(zhuǎn)換和分頻濾波對進行預處理，分析其各頻段的時頻域、極化方向的信號特征。對不同測量體、不同頻段、不同分量的信號特征進行描述，對產(chǎn)生的峰值信號分析其原因。對海洋環(huán)境電場信號進行空間相關(guān)性分析，結(jié)果表明，頻段的環(huán)境電場相關(guān)性較好，頻段環(huán)境電場相關(guān)性較差，驗證了環(huán)境電場相關(guān)性易受海浪、“舒曼諧振”影響。對環(huán)境電場信號按照環(huán)境電場相關(guān)性偏移時間進行相位補償，可以發(fā)現(xiàn)其相關(guān)性得到提高，為后期基于陣列處理的環(huán)境電場抑制提供依據(jù)。