[发明专利]时间域航空电磁数据约束多项式拟合调平方法

说明书

技术领域

本发明涉及一种时间域航空电磁资料处理方法，尤其是带约束的多项式拟合航空电磁数据调平方法。

背景技术

时间域航空电磁法(ATEM)是一种以飞机为载体的地球物理探测方法，广泛用于矿产资源、地下水及油气探测等方面，具有速度快、效率高等优点。但是航空电磁系统在飞行探测过程中，由于受到诸多的不可避免的干扰因素，如飞行方向、飞行高度、温度等，致使测区的各测线数据基准处于不同的水平上，在平面等值线图上呈现条带状误差，区域与局部异常被扭曲、畸变，产生虚假异常，数据质量变差。为了去除条带状误差，得到准确的异常分布信息，需要对航空电磁数据进行调平(leveling)。

CN105068131A公开了“一种航磁数据调平方法”。该方法沿测线方向和穿越测线方向对航磁数据进行总体平均经验模态分解，得到若干个子数据块，然后在该若干个子数据块中找出干扰子数据块并消除，将其余的子数据块相加得到调平后的航磁数据。该方法通过可视化操作找出调平误差所在，但是并未涉及航空电磁数据各时间道之间的关联约束。

航空数据调平方法在磁数据和电磁数据方面均有应用。Huang(2008)基于航空电磁数据各测线间相关性，以各时间道数据为一个整体，选取数据基准，依次拟合相邻测线调平误差，取得较好的效果；Majid等(2010)提出在一维和二维窗中进行多项式拟合，通过移动窗对数据进行逐点调平，该调平方法能更好的保留地质特征，且通用于规则测线和不规则测线的调平。到目前为止，仍未见航空电磁数据各时间道相互约束的调平方法。所以，本发明提出一种带约束的航空电磁数据调平方法。

发明内容

本文的发明目的在于提供一种准确的时间域航空电磁系统带约束的多项式拟合调平方法，该方法主要针对传统拟合方法忽略各时间道间关系，可能导致调平结果各道之间扭曲的问题。首先由于航空电磁数据的高度敏感性，去除电磁数据的高度影响；其次对于各时间道数据，选择基准测线，依次拟合相邻测线调平误差，得到一次调平结果；最后取一次调平结果，选择基准道，通过各测点的衰减关系建立各时间道间的不等式关系，以此不等式为约束，拟合相邻测线调平误差，得到调平结果。

本方法不仅考虑了各测线间的连续性，而且加入了各时间道间的相关性，使整体调平过程更加的完整。

时间域航空电磁数据约束多项式调平方法，包括以下步骤：

a、录入航空电磁数据；

b、去除数据高度影响；

c、选择基准测线，对各道数据分别进行一次调平；

d、输出一次调平结果；

e、判断一次调平结果各测线不符合衰减规律测点数量是否达到阈值，否，进行下一步；

f、选择基准道，对不符合衰减规律测线进行二次调平；

g、判断调平结果，达到阈值要求，否，返回到上一步；

h、是，把调平结果转换到原数据等级；

i、输出数据调平结果。

步骤a所述的录入时间域航空电磁数据，基于叠加偶极子假设，步骤b所述的为去除原始数据d_t高度影响，将系统接收的瞬态数据变换为数据d

d＝d_t×(h/R)³ (1)

其中，h为飞机飞行高度，R为发射线圈半径。若电磁探测系统发射线圈非圆形，取等效发射线圈半径其中，S为探测系统发射线圈半径。

步骤c所述的基准测线选择，根据数据的平面等值线图，在无明显条带状处，选择基准测线，假设选择测区数据d的第r条测线为基准测线，则基准测线表示为：其中，N表示各测线测点数，M表示各测点的时间道数，各测点表示为x＝(1,2,…,N)，一次调平是对各时间道数据进行独立调平，以第m道为例，基准测线d^r和其相邻测线d^r+1的第m道数据分别可以表示为和

选择测线中包含调平误差f(x)为多项式形式，即，

f(x)＝a₀+a₁x+a₂x²+…+a_kx^k，(2)

其中，a_k为多项式系数，k为多项式拟合阶数，相邻测线的调平结果可获得：

由于相邻测线间数据的相关性，存在