时间域航空电磁数据约束多项式拟合调平方法技术

本发明专利技术涉及一种时间域航空电磁数据约束多项式拟合调平方法，主要针对传统拟合方法忽略各时间道间关系，导致调平结果各道之间扭曲问题。首先去除电磁数据的高度影响；其次对于各时间道数据选择基准测线，依次拟合相邻测线调平误差，得到一次调平结果；选择基准道，通过各测点衰减关系建立各时间道间的不等式关系，以此不等式为约束，拟合相邻测线调平误差，得到调平结果。本发明专利技术的多项式拟合调平方法，可快速、有效的去除测区数据的条带状误差，在调平过程中加入了测点各时间道衰减关系不等式约束，调平结果的剖面曲线得到很大改善，避免调平结果各道数据之间的扭曲，保证了调平结果的准确性；不需要测线切割线辅助校正，节省了大笔费用。

【技术实现步骤摘要】

本专利技术涉及一种时间域航空电磁资料处理方法，尤其是带约束的多项式拟合航空电磁数据调平方法。
技术介绍
时间域航空电磁法(ATEM)是一种以飞机为载体的地球物理探测方法，广泛用于矿产资源、地下水及油气探测等方面，具有速度快、效率高等优点。但是航空电磁系统在飞行探测过程中，由于受到诸多的不可避免的干扰因素，如飞行方向、飞行高度、温度等，致使测区的各测线数据基准处于不同的水平上，在平面等值线图上呈现条带状误差，区域与局部异常被扭曲、畸变，产生虚假异常，数据质量变差。为了去除条带状误差，得到准确的异常分布信息，需要对航空电磁数据进行调平(leveling)。CN105068131A公开了“一种航磁数据调平方法”。该方法沿测线方向和穿越测线方向对航磁数据进行总体平均经验模态分解，得到若干个子数据块，然后在该若干个子数据块中找出干扰子数据块并消除，将其余的子数据块相加得到调平后的航磁数据。该方法通过可视化操作找出调平误差所在，但是并未涉及航空电磁数据各时间道之间的关联约束。航空数据调平方法在磁数据和电磁数据方面均有应用。Huang(2008)基于航空电磁数据各测线间相关性，以各时间道数据为一个整体，选取数据基准，依次拟合相邻测线调平误差，取得较好的效果；Majid等(2010)提出在一维和二维窗中进行多项式拟合，通过移动窗对数据进行逐点调平，该调平方法能更好的保留地质特征，且通用于规则测线和不规则测线的调平。到目前为止，仍未见航空电磁数据各时间道相互约束的调平方法。所以，本专利技术提出一种带约束的航空电磁数据调平方法。
技术实现思路
本文的专利技术目的在于提供一种准确的时间域航空电磁系统带约束的多项式拟合调平方法，该方法主要针对传统拟合方法忽略各时间道间关系，可能导致调平结果各道之间扭曲的问题。首先由于航空电磁数据的高度敏感性，去除电磁数据的高度影响；其次对于各时间道数据，选择基准测线，依次拟合相邻测线调平误差，得到一次调平结果；最后取一次调平结果，选择基准道，通过各测点的衰减关系建立各时间道间的不等式关系，以此不等式为约束，拟合相邻测线调平误差，得到调平结果。本方法不仅考虑了各测线间的连续性，而且加入了各时间道间的相关性，使整体调平过程更加的完整。时间域航空电磁数据约束多项式调平方法，包括以下步骤：a、录入航空电磁数据；b、去除数据高度影响；c、选择基准测线，对各道数据分别进行一次调平；d、输出一次调平结果；e、判断一次调平结果各测线不符合衰减规律测点数量是否达到阈值，否，进行下一步；f、选择基准道，对不符合衰减规律测线进行二次调平；g、判断调平结果，达到阈值要求，否，返回到上一步；h、是，把调平结果转换到原数据等级；i、输出数据调平结果。步骤a所述的录入时间域航空电磁数据，基于叠加偶极子假设，步骤b所述的为去除原始数据dt高度影响，将系统接收的瞬态数据变换为数据dd＝dt×(h/R)3 (1)其中，h为飞机飞行高度，R为发射线圈半径。若电磁探测系统发射线圈非圆形，取等效发射线圈半径其中，S为探测系统发射线圈半径。步骤c所述的基准测线选择，根据数据的平面等值线图，在无明显条带状处，选择基准测线，假设选择测区数据d的第r条测线为基准测线，则基准测线表示为：其中，N表示各测线测点数，M表示各测点的时间道数，各测点表示为x＝(1,2,…,N)，一次调平是对各时间道数据进行独立调平，以第m道为例，基准测线dr和其相邻测线dr+1的第m道数据分别可以表示为和选择测线中包含调平误差f(x)为多项式形式，即，f(x)＝a0+a1x+a2x2+…+akxk， (2)其中，ak为多项式系数，k为多项式拟合阶数，相邻测线的调平结果可获得： dl m r + 1 = d m r + 1 - f ( x ) , - - - ( 3 ) ]]>由于相邻测线间数据的相关性，存在 | | dl m r + 1 - d m r | | 2 = | | d m r + 1 - f ( x ) - d m r | | 2 = min - - - ( 4 ) ]]>联合式(2)和式(4)，解最小二乘方程，得到调平误差系数a＝(a0,a1,…,ak)T，代入式(3)，得到调平结果循环取出基准测线dr和其相邻测线dr+1的第m(m＝1,2,...,M)道数据，代入公式(2)～(4)得到相邻测线所有道调平结果dlr+1，以dlr+1为基准测线，依照此步骤调平其相邻测线dr+2，从而调平整个测区数据，得到一次调平结果dl。步骤e判断一次调平结果各测线不符合衰减规律测点数量是否达到阈值，设定阈值为测点数的p倍(N*p,p＜1)，判断一次调平结果dl各测线不符合衰减规律测点数是否小于阈值。若大于设定阈值，对此测线数据进行二次调平，根据初始数据和一次调平结果衰减曲线，选择较为符合衰减规律的道为基准道，基准道不需要进行二次调平，直接取一次调平结果，以第r+1条测线为例，假设选择各测点的第m-1道作为基准道，则其相邻第m道的二次调平结果数据应该满足： dl m - 1 r + 1 ≥ dls m r + 1 - - - ( 5 ) ]]>同时，该测线和其基准测线(第r条测线)的第m道数据还需要满足： | | dls m r + 1 - dl m r | | 2 = | | dl m r + 1 本文档来自技高网...

【技术保护点】
一种时间域航空电磁数据约束多项式拟合调平方法，其特征在于，包括以下步骤：a、录入航空电磁数据；b、去除数据高度影响；c、选择基准测线，对各道数据分别进行一次调平；d、输出一次调平结果；e、判断一次调平结果各测线不符合衰减规律测点数量是否达到阈值，否，进行下一步；f、选择基准道，对不符合衰减规律测线进行二次调平；g、判断调平结果，达到阈值要求，否，返回到上一步；h、是，把调平结果转换到原数据等级；i、输出数据调平结果。

【技术特征摘要】
1.一种时间域航空电磁数据约束多项式拟合调平方法，其特征在于，包括以下步骤：a、录入航空电磁数据；b、去除数据高度影响；c、选择基准测线，对各道数据分别进行一次调平；d、输出一次调平结果；e、判断一次调平结果各测线不符合衰减规律测点数量是否达到阈值，否，进行下一步；f、选择基准道，对不符合衰减规律测线进行二次调平；g、判断调平结果，达到阈值要求，否，返回到上一步；h、是，把调平结果转换到原数据等级；i、输出数据调平结果。2.按照权利要求1所述的时间域航空电磁数据约束多项式拟合调平方法，其特征在于，步骤a所述的录入时间域航空电磁数据，基于叠加偶极子假设，为去除原始数据dt高度影响，将系统接收的瞬态数据变换为数据dd＝dt×(h/R)3 (1)其中，h为飞机飞行高度，R为发射线圈半径。若电磁探测系统发射线圈非圆形，取等效发射线圈半径其中，S为探测系统发射线圈半径。3.按照权利要求1所述的时间域航空电磁数据约束多项式拟合调平方法，其特征在于，步骤c所述的基准测线选择，根据数据的平面等值线图，在无明显条带状处，选择基准测线，假设选择测区数据d的第r条测线为基准测线，则基准测线表示为：其中，N表示各测线测点数，M表示各测点的时间道数，各测点表示为x＝(1,2,…,N)，一次调平是对各时间道数据进行独立调平，以第m道为例，基准测线dr和其相邻测线dr+1的第m道数据分别可以表示为和选择测线中包含调平误差f(x)为多项式形式，即，f(x)＝a0+a1x+a2x2+…+akxk， (2)其中，ak为多项式系数，k为多项式拟合阶数，相邻测线的调平结果可获得： dl m r + 1 = d m r + 1 - f ( x ) , - - - ( 3 ) ]]>由于相邻测线间数据的相关性，存在 | | dl m r + 1 - d m r | | 2 = | | d m r + 1 - f ( x ) - d m r | | 2 = min - - - ( 4 ) ]]>联合式(2)和式(4)，解最小二乘方程，得到调平误差系数a＝(a0,a1,…,ak)T，代入式(3)，得到调平结果循环取出基准测线dr和其相邻测线dr+1的第m(m＝1,2,...,M)道数据，代入公式(2)～(4)得到相邻测线所有道调平结...

【专利技术属性】
技术研发人员：朱凯光，张琼，彭聪，陆依明，于紫凝，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22