基于同步提取變換的地空頻率域電磁信號幅度提取方法

摘要：地空頻率域電磁法探測信號為多頻非平穩信號，為了解決應用傳統傅里葉變換方法提取其幅度時分辨率較差的問題，本文提出了基于同步提取變換（SET）的地空頻域電磁信號幅度提取方法。該方法對電磁數據進行SET，得到高分辨率時頻圖，并利用能量算子使電磁數據時頻譜能量更為集中；采用貪心算法提取脊線，得到時頻圖的高能量帶；通過自回歸模型自適應地補充了脊線中的0值，解決了由窗函數引起的端點效應問題。根據脊線位置的時頻圖復數值，得到各頻率分量幅度隨時間的變化，研究了不同信噪比情況下基于SET的地空頻率域電磁信號幅度提取結果的準確性。結果表明：當信噪比≥10 dB時，幅度提取結果的平均相對均方根誤差均小于5%；當信噪比＜10 dB時，幅度提取結果的平均相對均方根誤差在10%以內。提取效果良好。將此方法應用于新疆霍拉山隧道工程地空頻率域電磁法探測中，成功提取了多頻電磁信號各頻率分量的幅度，與采用傅里葉變換提取非平穩信號幅度的方法相比，該方法有效地提高了幅度提取結果的分辨率。

關鍵詞：地空頻率域電磁法；時頻分析；非平穩信號；參數估計；同步提取變換

doi：10.13278/j.cnki.jjuese.20230084

中圖分類號：P319.3

文獻標志碼：A

0引言

地空頻率域電磁法（ground-airborne frequency-domain electromagnetic method， GAFDEM）采用地面發射、空中接收的形式。地面大功率發射系統以大地為負載，采用電偶/磁偶極矩超大功率電性/磁偶源向大地發射多個不同頻率的偽隨機波，空中磁感應傳感器采集磁場變化。通過分析和解釋不同頻率的空中磁場信號，獲得不同深度地下結構的電導率信息2]。由于GAFDEM具有地面發射功率強、勘探深度大、工作效率高等優勢，目前被廣泛用于地熱調查、火山結構調查、地下水鹽漬化研究、礦產資源勘探等領域。

空中接收系統受到飛機顛簸等因素的影響，導致數據存在大量運動噪聲和隨機噪聲，嚴重影響數據質量；因此，如何準確、有效地從中提取、分析地空頻率域電磁法信號各頻率分量的幅度變化，是研究者研究的重點之一。例如：康利利提出了基于小波的低頻運動基線去除算法，抑制了低頻運動噪聲在高頻段內的泄露，進而提高數據的信噪比和測量精度；Gao等提出了基于Robust-M 估計的時窗疊加濾波方法，通過抑制GAFDEM數據的環境噪聲提高數據的信噪比。地空頻率域電磁信號為隨時間變化的非平穩信號，傳統的頻域分析方法快速傅里葉變換（fast Fourier transform， FFT）無法準確有效地提取其各頻率分量的幅度。在非平穩信號分析領域，時頻分析（time-frequency analysis， TFA）起著重要的作用，它可以將一維信號映射到二維，采用時頻表示（time-frequency representation， TFR），用以表征信號的時變特性。近年來，研究人員提出了如重分配法（reassignment method， RM）、同步壓縮變換（Synchrosqueezing transform， SST）、高階SST、同步提取變換（synchroextracting transform， SET）、高階SET等時頻分析方法。SET方法利用同步提取算子，使信號能量更為集中，并且，SET信號恢復需要較少的參數，可以更方便地實現信號重建。因此，本文提出基于SET的地空頻率域電磁信號幅度提取方法，以解決應用傳統傅里葉變換方法提取地空頻率域電磁法探測信號幅度時分辨率較差的問題。

1方法原理

1.1地空頻率域電磁法數據特征

地空頻率域電磁系統主要由空中磁場收錄系統、地面大功率電性源發射系統和地面參考站三部分組成。其中：空中磁場收錄系統主要由機艙內接收機、雷達高度計、GPS等輔助設備、機艙外磁感應傳感器等構成，用于收集系統在測區三維坐標下磁場變化信息；地面大功率電性源發射系統由地面發電機、大功率多頻發射機、電偶極矩和大地構成，該系統通過向大地發射2n（n為正整數）個不同頻率的偽隨機波，得到不同深度的地下結構信息；地面參考站主要由接收機、三分量磁感應傳感器構成，負責收錄周圍空間磁場的變化。與傳統地面電磁法相比，地空頻率域電磁法的靈活性高、環境適應性更強18]?？罩写鸥袘獋鞲衅鞑杉讲煌l率對應的地空頻率域電磁信號，包含多個頻率分量，各分量的幅度隨時間變化。通過數據預處理后，該信號可表示為多個不同幅度和頻率的疊加：

由于采用了貪心算法，提取到的脊線為局部最優解。若在設定閾值范圍內沒有找到最優解，則可能會導致某個點無法提取到脊線，從而使得脊線矩陣出現預期外的0值，影響幅度提取的過程。由于信號幅度變化為線性的，通常情況下，選取合適的搜索閾值可以避免脊線0值點的出現，即使脊線中出現0值點，一般為稀疏的0值，難以出現連續的0值。為了消除脊線0值對信號幅度提取造成的影響，將脊線中的0值替換為空值“NaN”，采用自回歸建模方法，對NaN點進行填補。此方法從剩余樣本的正向和反向自回歸擬合推斷出估計值，替換信號中存在的所有NaN，消除了脊線0值的影響。優化效果如圖2所示。其中：圖2a為采用貪心算法實現的提取方法，可以看出存在若干斷點；圖2b為采用此方法處理過后的信號，可見斷點均已得到補充，曲線平滑。

自回歸建模方法同樣可以用以解決SET、STFT由于加窗函數而引起信號兩端出現的端點效應。設信號窗寬為h，在信號兩端各補充長度為h/2的NaN，對補充后的信號做SET、STFT，輸出長度取h2]即可。脊線提取效果如圖3所示。其中，圖3a為處理前的提取效果，圖3b為處理后的效果，可以明顯地看出端點效應得到了解決，提取信號兩端完整，符合預期。

圖4為采用貪心算法對圖1b中仿真信號進行脊線提取到的結果。在提取到的數據矩陣中，能量點位置值為1，其余位置值為0，有效地表示了混疊信號中各頻率分量的信息。從圖4可以看出，貪心算法可以準確提取GAFDEM數據經SET后的脊線。

2仿真分析

仿真信號由幅度緩慢變化的倍頻信號、50 Hz工頻以及噪聲信號混疊而成，仿真時間取10 s，采樣頻率1 kHz。其中，倍頻信號由16、32、64、128、256 Hz等頻率分量組成；幅度變換采用一次函數、正弦函數、對數函數、e指數函數等實現，混疊信號的表達式如下：

經過上述處理后，我們提取到了SET的脊線，根據脊線所在的位置，可以找到對應的能量點復矩陣，對復矩陣中的元素求模即可得到其幅度。采用SET提取到的各信號分量振幅如圖6所示，可以看出基于SET方法得到的不同頻率的振幅（實線）與仿真信號幅度（虛線）基本重合，驗證了基于SET的地空頻率域電磁信號幅度提取的有效性。

為了進一步評估本方法的性能，本文研究了不同信噪比條件下SET方法對地空頻率域電磁信號幅度提取結果的影響。仿真產生1組地空頻率域電磁信號，分別加入4組不同的高斯噪聲，使其信噪比分別為20、15、10和5 dB。采用式（11）分別計算了不同頻率幅度提取結果的相對均方根誤差，結果見表1。

由表1可以看出：當信噪比大于等于10 dB時，提取結果的平均相對均方根誤差均小于5%，且各頻率分量最大相對均方根誤差不大于10%；當信噪比小于10 dB時，提取結果的平均相對均方根誤差在10%以內。

3實測數據驗證

實測信號為新疆霍拉山隧道工程勘察項目的地空頻率域電磁系統測量數據，信號采樣頻率為12 kHz，測量時間為30 s，發射信號頻率分量包含32、128、208、256、272、336、448、688、752、1 280 Hz等10個頻率分量，接收到的信號包含對應頻率信號及其高次諧波，同時，信號中也含有50 Hz工頻噪聲、隨機高斯噪聲和尖峰噪聲等，其時域波形如圖7所示。

應用上述方法，對該信號進行SET，提取結果如圖8a所示。為了進一步提取各頻率分量的幅度信息，提取其脊線，結果如圖8b所示。從圖8中可以看出，除發射信號的10個頻率分量外，在0～1 400 Hz頻帶內還含有其他諧波分量。

從圖8b的脊線中選取欲提取的10個頻率，先計算其振幅，得到SET方法的振幅提取結果，如圖9a所示。圖9b為圖9a在0.1 s內的局部放大圖。然后將幅度轉換為頻率域電磁法的磁感應強度，與FFT方法進行對比，結果如圖9c所示?？梢钥闯?，相較于FFT方法，SET方法在保證提取精度的同時，顯著提高了提取結果的分辨率。

4結論

本文針對地空頻率域電磁數據特征，提出了基于同步提取變換的地空頻率域電磁信號幅度提取方法，得到如下結論。

1） 當信噪比大于等于10 dB時，幅度提取結果的平均相對均方根誤差均小于5%，且各分量的最大相對均方根誤差不大于10%；

當信噪比小于10 dB時，幅度提取結果的平均相對均方根誤差在10%以內。提取效果較好。

2）將同步提取變換方法應用于新疆霍拉山隧道工程勘查項目的實測數據處理中，成功提取了信號的幅度，且同步提取變換方法比快速傅里葉變換方法具有更好的空間分辨率。

參考文獻（References）：

1.周海根. 多場源地空頻率域電磁探測方法研究 . 長春：吉林大學， 2017.

Zhou Haigen.Research on Ground-Airborne Frequeney-Domain Electromagnetie Methods with Multiple Sources . Changchun： Jilin University， 2017.

Ward S H， Hohmann G W. Electromagnetic Theory for Geophysical Applications. Electromagnetic Methods in Applied Geophysics， 1988， 1： 131311.

2.康利利. 地空頻率域電磁探測系統運動噪聲抑制技術研究 . 長春：吉林大學， 2019.

Kang Lili. Research on the Motion Noise Suppression Technology of the Electromagnetic Detection System in the Ground-Air Frequency Domain. Changchun： Jilin University， 2019.

Ito H， Kaieda H， Mogi T， et， al. Grounded Electrical-Source Airborne Transient Electromagnetics （GREATEM） Survey of Aso Volcano， Japan. Exploration Geophysics， 2014， 45（1）： 4348.

Gao L， Yu S， Zhou H， et， al. Depth-Focused Waveform Based on SHEPWM Method for Ground-Airborne Frequency-Domain Electromagnetic Survey. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing， 2019， 12（6）： 19811990.

3.劉長勝， 朱文杰， 馬金發， 等. 地空頻率域電磁法探測范圍與深度分析. 中國礦業大學學報， 2020， 49（5）： 10061012.

Liu Changsheng， Zhu Wenjie， Ma Jinfa， et al. Analysis of the Detection Range and Depth of the Electromagnetic Method in the Ground-Air Frequency Domain. Journal of China University of Mining and Technology， 2020， 49（5）： 10061012.

Pons L J， Antonino D J A， Riera G M， et al. Advanced Induction Motor Rotor Fault Diagnosis via Continuous and Discrete Time-Frequency Tools. IEEE Transactions on Industrial Electronics， 2015， 62（3）： 17911802.

Yang Y， Peng Z K， Dong X J， et al. General Parameterized Time-Frequency Transform. IEEE Transactions on Signal Processing， 2014， 62（11）： 27512764.

Auger F， Flandrin P. Improving the Readability of Time-Frequency and Time-Scale Representations by the Reassignment Method. IEEE Transactions on Signal Processing， 1995， 43（5）： 1068–1089.

Daubechies I， Lu J， Wu H T. Synchrosqueezed Wavelet Transforms： An Empirical Mode Decomposition-Like Tool. Applied and Computational Harmonic Analysis， 2011， 30（2）： 243261.

Oberlin T， Meignen S， Perrier V. Second-Order Synchrosqueezing Transform or Invertible Reassignment？ Towards Ideal Time-Frequency Representations. IEEE Transactions on Signal Processing， 2015， 63（5）： 13351344.

Yu G， Yu M， Xu C. Synchroextracting Transform. IEEE Transactions on Industrial Electronics， 2017， 64（10）： 80428054.

Bao W， Li F， Tu X， et al. Second-Order Synchroextracting Transform with Application to Fault Diagnosis. IEEE Transactions on Instrumentation and Measurement， 2020， 70. doi： 10.1109/TIM.2020.3045841.

Huang Z L， Zhang J H， Zhao T H， et al.Synchrosqueezing S-Transform and Its Application in Seismic Spectral Decomposition. IEEE Transactions on Geoscience and Remote Sensing， 2016， 54（2）： 817825.

4.馮敏，陳輝，鄧居智，等.地空頻率傾子測深法及一維正演特征分析.石油地球物理勘探，2022，57（3）：728737， 497498.

Feng Min， Chen Hui， Deng Juzhi， et al. Ground-Air Frequency Dipper Sounding Method and Analysis of One-Dimensional Forward Modeling Characteristics. Petroleum Geophysical Exploration， 2022， 57（3）：728737， 497498.

5.齊彥福，李貅，孫乃泉，等. 電性源短偏移距瞬變電磁地形影響特征分析. 吉林大學學報（地球科學版），2022，52（1）：247260.

Qi Yanfu， Li Xiu， Sun Naiquan， et al. Analysis of Influence Characteristics of Topography on Grounded-Source Short-Offset Transient Electromagnetic Responses. Journal of Jilin University （Earth Science Edition）， 2022， 52 （1）： 247260.

6.李建華，林品榮，張強，等. 可控源電磁法中關鍵技術研究與應用.吉林大學學報（地球科學版）， 2022， 52（3）： 713724.

Li Jianhua， Lin Pinrong， Zhang Qiang， et al.Research and Application of Key Technology of Controlled Source Electromagnetic Method. Journal of Jilin University （Earth Science Edition）， 2022， 52（3）： 713724.

7.王新宇，嚴良俊，毛玉蓉，等. 起伏地形條件下長偏移距瞬變電磁三維正演. 吉林大學學報（地球科學版），2022，52（3）：754765.

Wang Xinyu， Yan Liangjun， Mao Yurong， et al. Three-Dimensional Forward Modeling of Long-Offset Transient Electromagnetic Method over Topography. Journal of Jilin University （Earth Science Edition）， 2022， 52 （3）： 754765.