[发明专利]一种拖曳式瞬变电磁数据预处理方法

权利要求书

1.一种拖曳式瞬变电磁数据预处理方法，其特征在于，该方法包括以下步骤：

S1、获取瞬变电磁仪器连续拖曳移动发射时采集的全时拖曳瞬变电磁数据，根据发射频率与接收机采样率提取每个发射周期下接收的瞬变电磁数据，每个发射周期采集的数据记为一个测点数据；

S2、设置相关度阈值，从第一个测点开始，依次计算该测点数据与相邻第一至第一百的其它测点的相关度值并记录，随后滑至下一个测点数据进行相关度计算，直至滑动完成所有测点数据的相关度计算；所述步骤S2中相关度阈值是经验阈值系数，只有当计算的其它相邻测点相关度高于这个阈值时，才进行相关加权叠加，通过降低阈值增加相关叠加数目，提高视电阻率成像纵向分辨率，通过提高阈值增加成像横向异常分辨率，阈值选取的范围为0.5到0.9之间；

S3、从第一个测点开始，将与该测点相关度值大于设置阈值的其他测点数据进行相关加权叠加，权值为对应的相关度值，将叠加后的数据记为新的测点数据，随后滑至下一个测点数据进行相关加权叠加计算，直至完成所有测点数据的相关加权叠加计算；

S4、所有相关叠加后的新测点数据为双极性梯形波发射的数据，将其正极性周期发射关断后对应的数据减去负极性周期发射关断后对应的数据后除以二，得到的数据记为拖曳式瞬变电磁预处理后的数据；

S5、所有测点完成上述计算后，完成拖曳式瞬变电磁数据预处理。

2.按照权利要求1所述的方法，其特征在于，所述步骤S1中的瞬变电磁仪器是由人工或者全地形车牵引的拖曳式瞬变电磁系统，探测线圈放置于移动平台上，采用中心回线配置，发射系统控制发射线圈在拖曳过程中连续发射，接收系统通过接收线圈同步采集全波形信号数据。

3.按照权利要求1所述的方法，其特征在于，所述步骤S1中的测点数据，其数据长度L_m由发射频率f_v与接收机采样率f_s定义，将拖曳式数据中每L_m个数据记为一个测点数据，表示为：

4.按照权利要求1所述的方法，其特征在于，

所述步骤S2相关度计算公式满足：

其中A与B为两个不同测点数据序列，A_m、B_m分别为对应测点数据序列里的第m个点，为A测点数据序列的均值，为B测点数据的均值，r为A与B测点的相关度计算值。

5.按照权利要求1所述的方法，其特征在于，所述步骤S3中相关加权叠加具体包括：通过计算目标测点与全部测点的相关度值，选取相关度值高于设置阈值的测点，将其作为权系数分别乘以对应的测点数据，叠加到目标测点，最后除以叠加的权系数之和进行数据归一化，得到目标测点叠加值U_m满足：

其中U为测点数据，U_k为对应的第k个测点数据序列，N为总测点数目，r_k为第k测点数据与第m测点数据的相关度值，r_min为相关度阈值，U_m为加权叠加后的第m测点的数据。

6.按照权利要求1所述的方法，其特征在于，所述步骤S4中双极性梯形波发射为：拖曳式瞬变电磁发射系统的连续发射波形，双极性梯形波发射周期为工频周期的偶数倍以抑制工频噪声的干扰，正负极性周期发射关断后对应的数据相减满足：

其中为相关加权叠加处理后测点数据的正极性周期发射关断后对应的接收数据，为负极性周期发射关断后对应的接收数据，U_m'为预处理后的第m测点拖曳式瞬变电磁数据。