CN106483573B - 基于双极性超光滑宽带激发源的航空电磁方法 - Google Patents
基于双极性超光滑宽带激发源的航空电磁方法 Download PDFInfo
- Publication number
- CN106483573B CN106483573B CN201611054894.8A CN201611054894A CN106483573B CN 106483573 B CN106483573 B CN 106483573B CN 201611054894 A CN201611054894 A CN 201611054894A CN 106483573 B CN106483573 B CN 106483573B
- Authority
- CN
- China
- Prior art keywords
- signal
- ultra
- wave
- bipolarity
- smooth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
- G01V3/165—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat operating with magnetic or electric fields produced or modified by the object or by the detecting device
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
基于双极性超光滑宽带激发源的航空电磁方法。包括以下步骤：选择与记录实际发射波形发射；实际信号的接收与处理；计算反演算法参数。利用设计的发射波形对地下分层情况进行反演。在发射波形为超光滑宽带激发源，其一阶及二阶导数在开关时刻均为零，并利用发射超光滑激发波来对地下电导率参数进行反演。适用于全航空瞬变电磁系统。克服现有双极性方波的下降沿信号由于混叠误差而难以准确探测的问题，同时抑制双极性方波信号正演过程中出现混叠误差现象。以BHW波为例的超光滑宽带激发源为宽带信号，包含大量的频率信息，同时在带外衰减迅速，能量集中在带内，可使发射信号能量更集中，同时提高接收信号频带内的信噪比。
Description
技术领域
本发明属于地球物理勘探技术领域，尤其是涉及一种基于双极性超光滑宽带激发源的航空电磁方法。
背景技术
瞬变电磁法(Transient Electromagnetic Methods)又称时间域电磁方法(TimeDomain Electromagnetic Methods)，简称TDEM或TEM。近年来瞬变电磁发展为航空瞬变电磁法(Airborne Transient Electromagnetics Methods，ATEM)和半航空瞬变电磁法(Semi-airborne Transient Electromagnetic Methods，SATEM)，也被称为地面电源航空瞬变电磁方法(Ground Electrical Source Airborne Transient ElectromagneticsMethods，GREATEM)。其中ATEM和GREATEM，由于接收机在直升机或者无人机中，探测范围大，并且可以探测到一些地形复杂的区域，近年来被广泛应用于矿产普查、地下水探测、石油勘探等领域。
航空电磁法是基于地下岩石的差异性，如电性或磁性差异，利用电磁感应原理，以固定翼或无人机等飞行器为载体，发射时域波形，根据接收信号来分析或反演出地下情况的勘探方法，其具有高效、经济、试用性强等优点。
然而对于ATEM其发射波形为双极性的半正弦波，在实际发射过程中容易产生吉布斯效应，导致实际测量结果与发射波形之间存在误差，从而影响实际反演结果；GREATEM发射波形为双极性方波信号，由于方波信号在实际产生时难以达到理想的状态，导致其上升沿或下降沿的时间太长，准确的上升时间或下降时间要用很高的采样率才能精确记录，但是也很难达到很高的精度，从而对反演结果有很大的影响。
发明内容
本发明的目的在于提供一种基于双极性超光滑宽带激发源的航空电磁方法。
本发明包括以下步骤：
(1)选择与记录实际发射波形发射，具体方法如下：
发射波形选择双极性的BHW波形，在发射机附近设有一个记录发射波形的接收机，用于确定发射信号的频谱信息，将记录得到的发射信号做快速傅里叶变换(FFT)，得到频域的增益系数；
(2)实际信号的接收与处理，具体方法如下：
接收机位于无人机或固定翼飞行器上，对无人机或固定翼飞行器的飞行姿态进项校正，对信号进行降噪处理；
在步骤(2)中，所述对无人机或固定翼飞行器的飞行姿态进项校正，对信号进行降噪处理的具体方法可为：先记录接收机的姿态，然后进行角度校正；采用多个周期累加求平均方法对信号进行降噪处理，当信号处理后，进行时频转换，得到接收信号的频谱信息和发射信号的频谱信息。
(3)计算反演算法参数，具体方法如下：
得到接收信号的频谱信息和发射信号的频谱信息后，用反演方法对地下电导率等信息进行反演，得到地下的电导率分布情况。
本发明提出了用发射超光滑宽带激发源作为一种新型的航空电磁和半航空电磁的发射波形。例如基于BHW波，其波形表达式如下：
其中a0＝0.353222，a1＝-0.488,a2＝0.145,a3＝-0.0102222；信号的特征频率为fc＝1/T，对于零阶BHW波信号的最大频率近似为fmax≈3.351fc。这种波形的一阶及二阶导数在开关时刻均为零，故具有良好的频谱特性，消除了传统半正弦波及方波的混叠误差。
对于上述BHW波其一阶导数表达式为：
二阶的表达式为:
由表达式可知，BHW波的一阶与二阶导数在0时刻为0，这可以避免实际应用中双极性方波在0时刻导数不为零所引起的不连续问题。并且BHW波为余弦信号按照不同的比例因子叠加而成，所以在实际应用中容易实现，且不会出现过载现象，在正演仿真时可以避免因为频率分量过高而出现的混叠现象，导致仿真结果不稳定；在实际应用中要去掉去直流分量，同时为了避免零点偏移现象，所以采用双极性BHW波作为发射波形；BHW波的波谱能量集中，带外衰减快，有利于实际接收信号在带内的信噪比提高；同时可以避免在实际应用中出现的过载现象，有利于实际接收机的应用与提高发射信号的记录精度；其双极性BHW波表达式如下：
对于T和T0可以取不同的数值，这里取T＝0.005，T0＝0.005；
本发明提出一种新型的航空电磁与半航空电磁的发射波形，同时利用设计的发射波形对地下分层情况进行反演。在瞬变电磁系统新型的发射波形和利用新型发射波形准确反演地层电导率的算法。本发明的最大特点是在发射波形为超光滑宽带激发源，其一阶及二阶导数在开关时刻均为零(如BHW波)，并利用发射超光滑激发波来对地下电导率参数进行反演。
本发明使用超光滑宽带激发源(其一阶及二阶导数在开关时刻均为零)，如BHW(Blackman-Harris Window)波是简化汉明窗的一种，具有频谱能量集中，带外衰减快的优点；且是宽带信号包含频率成分大，对于带内频率能量大，可以很好的提高带内频率的信噪比；在实际产生过程中由于其是由正弦信号叠加而成，所以实际电路容易产生，且可以避免其他波形的过载和混叠效应等。
附图说明
图1为BHW波形示意图。
图2为BHW波频谱衰减示意图。
图3为BHW波一阶导数示意图。
图4为BHW波二阶导数示意图。
图5为双极性BHW波的波形图。
图6为双极性BHW波频谱衰减示意图。
图7为利用BHW波反演出的一维分层模型结果。
具体实施方式
以下将结合附图对本发明的技术方案及其突出效果作进一步说明。
(1)实际发射波形发射选择与记录
因为实际接收机受到环境、仪器等多种情况的影响，所以一般的发射波形选择双极性的BHW波形(如图2)。这样可以减少直流分量对接收信号的影响；同时在发射机附近要有一个记录发射波形的接收机，以便与确定信号的频谱信息。将记录得到的发射信号做FFT变换得到其频域的增益系数；
(2)实际信号的接收与处理
实际接收机一般位于无人机或者固定翼飞行器上，所以就需要对其飞行的姿态进项校正和对信号进行降噪处理。这里需要记录接收机的姿态，然后对其进行角度校正；为了对接收信号进行降噪，一般采用多个周期累加求平均方法。当信号处理好后，需要对其进行时频转换，得到其频谱信息。
(3)反演算法参数的计算
得到接收信号的频谱信息和发射信号的频谱信息后，就可用反演方法对地下电导率等信息进行反演，来得到地下的电导率分布情况。
本发明不仅适用于半航空瞬变电磁系统，而且适用于全航空瞬变电磁系统。本发明从瞬变电磁的发射波出发，提出了一种新型的航空电磁发射波形：超光滑宽带激发源(其一阶及二阶导数在开关时刻均为零)，并以BHW波为例指出了光滑激发源的优点。超光滑宽带激发源波形克服现有双极性方波的下降沿信号由于混叠误差而难以准确探测的问题，同时抑制双极性方波信号正演过程中出现混叠误差现象。以BHW波为例的超光滑宽带激发源为宽带信号，包含大量的频率信息，同时在带外衰减迅速，能量集中在带内，可以使发射信号能量更集中，同时提高接收信号频带内的信噪比。最后给出1D的BHW波波形信息反演出来的结果。
本反演算法在理论上具有较高的反演精度：
图1为0阶BHW波示意图，图中的BHW波的特征频率为fc＝200Hz，T＝0.005s。
图2为0阶BHW波频谱示意图，从图中可以看出0阶的BHW波在特征频率内的增益衰小，但在特征频率的3.351倍的最大频率下其衰减可以达到-40dB以上，说明BHW波的能量主要集中在0～fc频带内。
图3为1阶BHW波波形示意图，此波形为图1BHW波的一阶导数，也成为一阶BHW波，从图中可以看出其在0时刻时的值为0。
图4为2阶BHW波波形示意图，此波形为图1BHW波的二阶导数，也成为二阶BHW波，从图中可以看出其在0时刻时的值为0。
图5为双极性BHW波波形示意图，图中的双极性BHW波，其中单个的BHW波的特征周期为200Hz，其中T＝0.005s，中间的间隔时间T0＝0.005s，整个双极性BHW波的周期为50Hz。
图6为双极性BHW波频谱示意图，从图中可以看出双极性BHW波在单极性波的特征频率内的增益衰减由于受周期性的影响在频带内变化比原来大，但是其双极性波没有直流分量，并且在单极性BHW波特征频率的3.0855倍的最大频率下其衰减可以达到-40dB以上，说明BHW波的波谱能量主要集中在其频带内。
图7利用BHW波反演结果图，利用特征频率为200Hz的BHW波所反演出的一维数据结果，可以看出模型为5层的分层模型，反演过程中分10层来对模型进行反演，最后得到的结果和实际的真实值结果很接近，基本上可以认为反演结果和真实结果相吻合。
Claims (2)
1.基于双极性超光滑宽带激发源的航空电磁方法，其特征在于，包括以下步骤：
(1)选择与记录实际发射波形
在发射机附近设置一记录发射波形的接收机，用于确定发射信号的频谱信息，将记录得到的发射信号的频谱信息做快速傅里叶变换，得到频谱信息的增益系数；其中，发射机为双极性超光滑宽带激发源，发射波形为双极性超光滑的BHW波，具体表达式如下：
其中，对于T和T0取不同的数值，取T＝0.005，T0＝0.005；
(2)实际信号的接收与处理
将接收机设置于无人机或固定机翼飞行器上，对无人机或固定翼飞行器的飞行姿态进行校正，对信号进行降噪处理；
(3)计算反演算法参数
将得到的接收信号的频谱信息和发射信号的频谱信息，用反演方法对地下电导率信息进行反演，得到地下的电导率分布情况。
2.如权利要求1所述的基于双极性超光滑宽带激发源的航空电磁方法，其特征在于步骤(2)中，对无人机或固定翼飞行器的飞行姿态进行校正，对信号进行降噪处理的方法为：先记录接收机的姿态，然后进行角度校正；采用多个周期累加求平均方法对信号进行降噪处理；当信号处理完成后，进行时频转换，得到接收信号的频谱信息和发射信号的频谱信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611054894.8A CN106483573B (zh) | 2016-11-25 | 2016-11-25 | 基于双极性超光滑宽带激发源的航空电磁方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611054894.8A CN106483573B (zh) | 2016-11-25 | 2016-11-25 | 基于双极性超光滑宽带激发源的航空电磁方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106483573A CN106483573A (zh) | 2017-03-08 |
| CN106483573B true CN106483573B (zh) | 2019-01-22 |
Family
ID=58274201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611054894.8A Expired - Fee Related CN106483573B (zh) | 2016-11-25 | 2016-11-25 | 基于双极性超光滑宽带激发源的航空电磁方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106483573B (zh) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2006297766B2 (en) * | 2005-09-19 | 2011-04-28 | Exxonmobil Upstream Research Company | Method for performing controlled source electromagnetic surveying with multiple transmitters |
| WO2014035543A1 (en) * | 2012-08-30 | 2014-03-06 | Exxonmobil Upstream Research Company | Processing methods for time division csem data |
| CN105807325B (zh) * | 2014-12-31 | 2018-07-13 | 中国船舶重工集团公司第七研究院 | 一种频率域航空极低频电磁法 |
| CN104953884A (zh) * | 2015-07-14 | 2015-09-30 | 中国科学院电子学研究所 | 全航空瞬变电磁系统双极性半正弦电流产生装置和方法 |
| CN106125145A (zh) * | 2016-06-30 | 2016-11-16 | 北京工业大学 | 基于宽频带lfm激励源的大地电磁脉冲响应辨识方法及系统 |
-
2016
- 2016-11-25 CN CN201611054894.8A patent/CN106483573B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN106483573A (zh) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107576943B (zh) | 基于瑞利熵的自适应时频同步压缩方法 | |
| Chen et al. | Effect of transmitter current waveform on TEM response | |
| AU2011100086A4 (en) | A method for measuring biomass of tree-roots by ground-penetrating radar | |
| Ji et al. | Noise reduction of grounded electrical source airborne transient electromagnetic data using an exponential fitting-adaptive Kalman filter | |
| CN105629317B (zh) | 一种基于站间传递函数的大地电磁噪声压制方法 | |
| YIN et al. | Exploration capability of airborne TEM systems for typical targets in the subsurface | |
| US20180038976A1 (en) | Surface Consistent Statics Solution and Amplification Correction | |
| CN106908786B (zh) | 一种基于极化散射矩阵估计的高精度昆虫体轴朝向提取方法 | |
| CN101937104A (zh) | 用于在海洋电磁勘探中估计并除去空气波响应的方法 | |
| Ma et al. | Joint processing of forward and backward extended Prony and weighted spectral semblance methods for robust extraction of velocity dispersion data | |
| Omelchuk et al. | Multichannel target speed estimation by a colocated Doppler-pulse MIMO radar | |
| CN106483573B (zh) | 基于双极性超光滑宽带激发源的航空电磁方法 | |
| CN106405665B (zh) | 基于dbim的瞬变电磁电导率反演方法 | |
| CN110515044A (zh) | 低慢小雷达空频域四维杂波抑制方法 | |
| Zhu et al. | Simultaneous observations of electric field changes, wideband magnetic field pulses, and VHF emissions associated with K processes in lightning discharges | |
| US10274635B2 (en) | Joint inversion of subsurface resistivity and noise parameters | |
| Hennessy et al. | Source-dependent bias of sferics in magnetotelluric responses | |
| Gao et al. | Specific emitter identification based on instantaneous frequency characteristics | |
| US10613246B2 (en) | Joint estimation of electromagnetic earth responses and ambient noise | |
| ZHANG et al. | Calibration study on the damping characteristics of TEM sensor | |
| Eisen et al. | Identifying isochrones in GPR profiles from DEP-based forward modeling | |
| CN104297733B (zh) | 复合背景下时域信号的处理方法及复合信号处理系统 | |
| CN112713907B (zh) | 一种基于字典学习的海洋csem降噪方法及系统 | |
| Shvets et al. | A technique for lightning location and estimation of the lower ionosphere parameters using tweek-atmospherics | |
| CN107024724A (zh) | 反卷积快速优化直升机航空电磁探测早期数据的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190122 Termination date: 20211125 |