基于包絡(luò)信號的軸頻電場深度換算方法

靳 雄, 姜潤翔, 程錦房, 陳新剛

基于包絡(luò)信號的軸頻電場深度換算方法

靳 雄1, 姜潤翔2, 程錦房1, 陳新剛3

(1. 海軍工程大學 兵器工程學院, 湖北 武漢, 430033; 2. 海軍工程大學 電氣工程學院, 湖北 武漢, 430033; 3. 中國人民解放軍92942部隊, 北京, 100161)

在僅利用水平時諧電偶極子對艦船軸頻電場信號進行換算時, 存在近場深度誤差較大的問題。文中首先在仿真數(shù)據(jù)的基礎(chǔ)上, 明確了船體表面不同位置電流密度與軸地電阻值波動時的變化規(guī)律; 其次, 提出了基于包絡(luò)信號的軸頻電場深度換算方法, 該方法在希爾伯特(Hilbert)變換計算軸頻電場信號包絡(luò)的基礎(chǔ)上, 將軸頻電場近場的深度換算問題轉(zhuǎn)化為包絡(luò)信號等效靜電場的近場換算問題, 并利用點電荷模型建立了軸頻電場信號包絡(luò)在“空氣—海水—海床”3層介質(zhì)條件下的正演及反演模型。最后, 分別利用4種陰極防腐狀態(tài)下的船模試驗數(shù)據(jù)對所提方法的有效性進行了檢驗, 結(jié)果表明, 所提方法能夠較好地實現(xiàn)對軸頻電場信號包絡(luò)值的準確換算, 以相對均方根誤差作為評價準則, 水深為1倍船模寬度的換算誤差小于15%。文中方法可為艦船近場的軸頻電場反演提供新的途徑。

艦船; 軸頻電場; 點電荷; 反演

0 引言

艦船軸頻電場因其頻率低、傳播距離遠、線譜特征明顯, 可被用做水下遠距離探測的被動信號源[1-7]。為了能夠?qū)悍脚灤S頻電場信號的量級及特性進行評估, 同時能夠有效評價軸頻電場抑制裝置的有效性, 均需對軸頻電場的源強度及不同深度的軸頻電場信號進行反演計算。

在遠場(1.0倍船長)附近, 軸頻電場可近似視為水平時諧電偶極子產(chǎn)生的電場, 程銳等[8]、Daya等[9]利用水平時諧電偶極子對軸頻電場信號進行了反演計算; 熊露等[10]利用單個水平時諧電偶極子對淺海中的船舶軸頻電場進行了建模, 以上都是針對等效電偶極矩水平方向的分量進行的研究。由于艦船等效電偶極矩還存在垂直方向的分量, 該分量雖然遠小于水平方向的等效電偶極矩, 但其在近場附近仍然可以產(chǎn)生較大量值的信號, 因此, 在僅利用水平時諧電偶極子對軸頻電場換算時, 存在近場誤差較大的問題, 為此, 需要研究新的軸頻電場反演算法。

1 軸頻電場信號

1.1 產(chǎn)生機理

1.2 軸地電阻波動的隨機性

典型艦船尾部電路結(jié)構(gòu)如圖1所示。圖中, ICCP (impressed current cathodic protection)指外加電流陰極保護。

圖1 艦船尾部電路結(jié)構(gòu)示意圖

由圖1可知, 船體內(nèi)部回路電阻并聯(lián)可得

文獻[12]研究結(jié)果表明, 軸接地裝置的電阻(碳刷和滑環(huán)的接觸電阻)明顯小于其他接地電阻, 即有。與軸的載荷、偏心度、碳刷彈簧的彈性等因素密切相關(guān), 艦船同種轉(zhuǎn)速條件下的軸地等效電阻的波動情況如圖2所示, 其中由軸地電壓值除以軸電流(非接觸電流傳感器測量結(jié)果)得到。由圖2可知, 軸地電阻的波動值具有明顯的隨機性, 且左軸與右軸的差異較大。

軸地電阻波動的隨機性必將導致水下軸頻電場信號的隨機性, 因此, 在對軸頻電場信號建模時, 一種合適的方法是對某一時間段內(nèi)軸頻電場信號的包絡(luò)進行建模, 以得到軸頻電場信號的最大幅值。

1.3 軸地電阻對船體表面電流分布的影響

利用邊界元法對典型艦船(如圖3所示)在不同條件下的靜電場信號進行仿真計算, 可得船體在自然腐蝕狀態(tài)及陰極保護狀態(tài)(4對輔助陽極恒電流輸出、減搖鰭、美人架及舵板局部犧牲陽極防腐)條件下軸地電阻變化時船體表面電位、電流密度及水下電場的特性。

圖4 自然腐蝕狀態(tài)條件下船體表面電流密度曲線

圖5 陰極防腐條件下船體表面電流密度曲線

圖4和圖5對應(yīng)的1倍船寬深度, 正橫距為5 m處的水下電場信號分別如圖7和圖8所示。由圖可知, 軸地等效電阻的變化也將導致水下靜電場信號的變化, 且不同位置的電場值變化不同, 即無法利用某一固定的波動系數(shù)在靜電場的基礎(chǔ)上完成對軸頻電場的建模。

圖8 陰極防腐狀態(tài)下不同軸地等效電阻靜電場信號變化曲線

2 軸頻電場包絡(luò)反演算法

2.1 基本原理

2.2 反演步驟

根據(jù)基本原理, 確定反演步驟如下。

1) 利用Hilbert變換計算已知深度不同測線軸頻電場信號的峰值包絡(luò), 并對峰值包絡(luò)進行0.5 Hz的低通濾波, 在峰值包絡(luò)過零點處對信號進行反相, 得到軸頻電場信號包絡(luò)。

的求解。

對式(5)增加控制方程

同時, 根據(jù)電中性條件有

利用最小值條件

3) 在步驟2)的基礎(chǔ)上, 利用式(3)即可實現(xiàn)對目標深度的軸頻電場信號進行正演, 得到不同位置處軸頻電場信號的包絡(luò)。

4) 對比目標深度軸頻電場信號包絡(luò)與反演出的軸頻電場信號的包絡(luò)以檢驗反演效果, 并計算其相對均方根誤差

需要說明的是, 在利用Hilbert變換計算得到的軸頻電場包絡(luò)均是正值, 這與實際艦船等效電荷電中性的前提條件相矛盾, 為此, 在包絡(luò)信號過零點處, 將軸頻電場信號的包絡(luò)反相, 以保證包絡(luò)信號具有正、負電荷的特性。

3 算例驗證

鑒于縮比模型法已成為應(yīng)用于艦船外加電流陰極保護系統(tǒng)設(shè)計、腐蝕相關(guān)電場特性、涂層破損和流速對腐蝕和防腐的影響等研究領(lǐng)域的標準方法[14-16], 文中利用縮比模型船模(如圖9所示, 模型長度=2.72 m、船寬=34 cm、吃水深度=9 cm, 船體材料為涂層鋼, 螺旋槳為鎳鋁青銅)試驗對所提方法的有效性進行了檢驗。

圖9 縮比船模試驗

分別對船模在自然腐蝕狀態(tài)、局部犧牲陽極防腐(美人架、舵板及減搖鰭附近安裝有犧牲陽極)、局部犧牲陽極與外加電流陰極(2對輔助陽極)保護、全船犧牲陽極防腐4種狀態(tài)下的軸頻電場信號進行了測量, 并利用參比電極和非接觸電流傳感器對船體中心的電位及軸電流進行了監(jiān)測。

值得注意的是, 由于模型等比例縮小, 受測量電極尺寸的影響, 難以準確獲知軸頻電場的E、E和E分量。因此, 試驗中選擇對軸頻電場信號的電位進行測量, Ag/AgCl電極分別置于水深24 cm、34 cm和51 cm 等3個不同深度, 每個深度平面3個Ag/AgCl電極的正橫距分別為0 cm、17 cm和34 cm, 基準Ag/AgCl電極置于水深1.2 m, 基準電極距離測量電極的最小距離為4 m。電機拖動裝置控制船模以10.08 cm/s的速度勻速通過測量電極上方時, 利用微弱電場測量裝置實時記錄軸頻電場的電位信號, 系統(tǒng)帶寬為0~20 Hz, 采樣頻率=100 Hz, 螺旋槳轉(zhuǎn)動頻率為280r/min左右。

圖10為全船犧牲陽極防腐狀態(tài)下24 cm水深的軸頻電場信號及其包絡(luò)值。圖中: 軸頻電場信號在74.5 s過零點; 在74.5 s之前, 軸頻電場信號包絡(luò)與Hilbert包絡(luò)重合; 在74.5 s之后, 兩者波形反相。

圖10 水深24 cm時軸頻電場信號及其包絡(luò)值

將軸頻電場包絡(luò)的等效36個點電荷分別置于船體中軸線、左弦0.5、右弦0.5這3條線上, 每條線的點電荷數(shù)為12個, 點電荷的垂直位置為0.5。圖11和圖12為采用文獻[8]所提靜電場反演的算法對等效源強度進行計算, 并對34 cm和51 cm深度的軸頻電場信號包絡(luò)進行反演結(jié)果。

圖11 水深34 cm時軸頻電場信號及其包絡(luò)值

由圖11和圖12可知, 反演所得軸頻電場信號的包絡(luò)與實測軸頻電場信號的包絡(luò)具有良好的一致性, 證明了算法的有效性。按式(8)計算所得不同狀態(tài)下軸頻電場的反演誤差如表1所示。由表1可知, 在不同陰極保護狀態(tài)下, 所提方法均能對軸頻電場信號的包絡(luò)進行準確反演, 最大反演誤差小于25%, 且隨著測量水深的增加, 軸電流的增大, 誤差有減小趨勢, 在1.0深度平面上, 反演誤差小于15%, 這與靜電場反演所得結(jié)論一致, 即反演算法由淺向深反演具有較高的精度, 隨著軸電流的增大, 信噪比提高, 反演誤差也減小。

表1 不同船體狀態(tài)條件下軸頻電場信號包絡(luò)反演誤差

4 結(jié)束語

為了解決近場軸頻電場反演精度不高的問題, 文中提出一種基于軸頻電場包絡(luò)的反演算法, 該方法將軸頻電場的反演問題轉(zhuǎn)化為與靜電場反演相似的問題, 避免了求解軸地等效電阻波動時引起船體表面各部分電流變化率的難題, 通過船模實測數(shù)據(jù)檢驗發(fā)現(xiàn), 該算法換算效果較好, 精度較高, 為近場的軸頻電場反演提供了一種新的途徑。下一步研究將進行實船實驗, 以驗證文中所提方法的有效性。

[1] Bostick F, Smith H, Boehl J. The Detection of ULF–ELF Emissions from Moving Ships[R]．New York: State Aca- demic Educational Institutions, 1977: 13-24．．

[2] Dymarkowski K, Uczciwek J. The Extremely Low Frequency Electromagnetic Signature of the Electric Field of the Ship[C]//Underwater Defence Technology Conference, London: IEEE, 2001, 1-6．

[3] Zolotarevskii Y, Bulygin F, Ponomarev A, et al．Methods of Measuring the Low-frequency Electric and Magnetic Fields of Ships[J]．Measurement Techniques, 2005, 48(11): 1140-1144．

[4] Lü X C, Chen H Q, Ma S J, et al．The Electromagnetic Fi- eld Generated by Time-harmonic Electric Current Element in Shallow Sea[C]//Hardy T, Undersea Defence Techno- logy, Hamburg: IEEE, 2006．

[5] Jia Y Z, Jiang R , Gong S G．Research on Wavelet Modulus Maximum-based Detection Algorithm of Ship’s Sh- aft-rate Electric Field[J]．Acta Armamentarii, 2013, 34(5): 579-584．

[6] Gheorghe S, Serghei R, Camelia C. Mesurements of Elect- rical and Magnetic Fields on Board Container Ships[C]//Pro- ceedings of the Scientific Conference AFASES. Brasov.Rumania:Air Force Academy: 2013.

[7] 賈亦卓, 姜潤翔, 龔沈光．基于小波模極大值的船舶軸頻電場檢測算法研究[J]．兵工學報, 2013, 34(5): 579-584．Jia Yi-zhuo, Jiang Run-xiang, Gong Shen-guang．Research on Wavelet Modulus Maximum-based Detection Algori- thm of Ship’s Shaft-rate Electric Field[J]. Acta Electro- nica Sinica, 2013: 34(5): 549-584.

[8] 程銳, 姜潤翔, 龔沈光．船舶軸頻電場等效源強度計算[J]．國防科技大學學報, 2016, 38(2): 138-143．Cheng Rui, Jiang Run-xiang, Gong Shen-guang. Calculation Method of Vessel’s Shaft Rate Electric Field Equivalent Source Magnitude[J]. Journal of National University of Defense Technology, 2016, 38(2): 138-143.

[9] Daya Z A, Hutt D L, Richards Troy C. Maritime Electromagnetism and DRDC Management Research[R]．Canada: Defence R&D Canada, 2005: 23-25．

[10] 熊露, 王斌, 畢曉文. 艦船軸頻電場場源建模和實驗研究[J]. 武漢理工大學學報, 2015, 37(7): 52-56.Xiong Lu, Wang Bin, Bi Xiao-wen. Modeling and Experimental Study of Ship’s Axial Frequency Electric Field Source[J]. Journal of Wuhan University of Technology, 2015, 37(7): 52-56.

[11] 林春生, 龔沈光．艦船物理場[M]．北京: 兵器工業(yè)出版社, 2007．

[12] Howard L C, Ballston L N. Naval Electro-chemical Corrosion Reducer: United States: 5052962[P]. 1990-10-1．

[13] 姜潤翔, 林春生, 龔沈光．基于點電荷模型的艦船靜電場反演算法研究[J]．兵工學報, 2015, 36(3): 545-551． Jiang Run-xiang, Lin Chun-sheng, Gong Shen-guang. Ele- ctrostatic Electric Field Inversion Method for Ship Based on Point Charge Source Model[J]. Acta Electronica Sinica, 2015: 36(3): 545-551.

[14] Parks A R, Thomas E D, Lucas K E. Physical Scale Modeling Verification with Shipboard Trials[J]. Material Perf- ormance, 1991, 30(5): 26-29.

[15] 邢少華, 張搏, 閆永貴, 等．涂層破損對船體陰極保護電位分布的影響[J]. 材料開發(fā)與應(yīng)用, 2016, 31(1): 69-73．Xing Shao-hua, Zhang Bo, Yan Yong-gui, et al．The Influ- ence of Coating Damage on Cathodic Protection Potential Distribution of Ship[J]. Development and Application of Materials, 2016, 31(1): 69-73.

[16] DeGiorgi V G, Thomas E D II, Lucas K E, et al. A Combined Design Methodology for Impressed Current Catho- dic Protection Systems[J]. Boundary Element Technology, 1996, 12(14): 335-344.

A Depth Conversion Method of Shaft-Rate Electric Field Based on Envelope Signal

JIN Xiong1, JIANG Run-xiang2, CHENG Jin-fang1, CHEN Xin-gang3

(1. Department of Weaponry Engineering, Naval University of Engineering, Wuhan 430033, China; 2. College of Electrical Engineering, Naval University of Engineering, Wuhan 430033, China; 3. 92942thUnit, The People’s Liberation Army of China, Beijing 100161, China)

Only application of the horizontal time-harmonic dipole to the conversion of the shaft-rate electric field signal of ship will result in large near-field depth error. In this study, the variation laws of the current density in different position of hull surface and the shaft-ground resistance are clarified on the basis of simulation data. A depth conversion method of shaft-rate electric field is proposed based on its envelope signal. In this method, the depth conversion problem in the near field of the shaft-rate electric field signal envelope is transformed into the near-field conversion problem of equivalent static electric field after the calculation of the shaft-rate electric field signal envelope via Hilbert transform, and the point charge model is used to establish a forward and inversion model of the shaft-rate electric field signal envelope in the ‘a(chǎn)ir-sea-seabed’ three-layer media condition. The effectiveness of the proposed method is verified by four kinds of cathodic corrosion protection test data of the ship model. The results show that the proposed method can accurately calculate the shaft-rate electric field signal envelope. Taking root square error as evaluation criteria, the conversion error is less than 15% in the depth equaling the ship model width.

ship; shaft-rate electric field; point charge; inversion

U674.70; TP274

A

2096-3920(2020)04-0403-07

10.11993/j.issn.2096-3920.2020.04.008

2019-10-04;

2019-12-09.

青島海洋科學與技術(shù)國家實驗室“問海計劃”項目(2017WHZZB0101).

靳 雄(1991-), 男, 在讀碩士, 主要研究方向為目標特性及信息感知技術(shù).

靳雄, 姜潤翔, 程錦房, 等. 基于包絡(luò)信號的軸頻電場深度換算方法[J]. 水下無人系統(tǒng)學報, 2020, 28(4): 403-409.

(責任編輯: 楊力軍)