[发明专利]瞬变电磁信号叠前去噪方法

摘要

本发明公开了一种瞬变电磁勘探的时间序列信号的叠前去噪方法，包括以下步骤：S1、存贮瞬变电磁时间序列；S2、对瞬变电磁测量的场分量时间序列的每个周期衰减曲线进行倒序时域递归滤波；S3、对滤波后每个周期的场分量衰减曲线进行加窗处理，获取对数等间隔的场值衰减曲线；S4、采取设定阈值的统计叠加处理方法，对对数等间隔的场值衰减曲线中的多个观测周期的同一时窗场值进行去飞点叠加。本发明使得LOATEM方法资料处理快速有效，信噪比明显改善，资料畸变明显减小，资料品质明显提高，为视电阻率曲线计算，高分辨的精细反演奠定了坚实基础。

主权项

1.一种瞬变电磁勘探的时间序列信号的叠前去噪方法，其特征在于，包括以下步骤：S1、存贮瞬变电磁时间序列；S2、对瞬变电磁测量的场分量时间序列的每个周期衰减曲线进行倒序时域递归滤波；S3、对滤波后每个周期的场分量衰减曲线进行加窗处理，获取对数等间隔的场值衰减曲线；S4、采取设定阈值的统计叠加处理方法，对对数等间隔的场值衰减曲线中的多个观测周期的同一时窗场值进行去飞点叠加；步骤S1中，仅提取关断后的信号，且将每个周期内的两次关断后的二次场衰减曲线相加后再平均，形成一条衰减曲线，在一个测点进行K个周期观测时，则为K个周期衰减曲线的同一延迟时间的场值相加后平均；步骤S3中进行加窗处理时，将关断后的时间分成N个对数等间隔时窗，对数域里的间隔宽度为Δ，每个时窗的中心时间点记为t₁,t₂......t_N，其所对应的场值分别计为EMF(t₁),EMF(t₂),.......EMF(t_N)，N为自然数；将每个时窗分成4份，依次记为t_i+Δ/4,t_i+Δ/2,t_i+3Δ/4,t_i+Δ；为了计算中心时间点处的场分量值EMF(t_i)，将时窗宽度扩大到3Δ/2，以保证每个时窗有Δ/2的重叠；将落入时窗中的K个采样点进行最小二乘拟合；设y是采样时间ts及参量的函数，(y_i,ts_i)是K对观测值，其中y_i是观测场值，ts_i是对应的采样时间，i＝1,......,K，K、M均为自然数；求参量使目标函数达到最小，其中σ_i为观测场值y_i误差的标准差；则时窗中心点所对应的场值为其中i＝1,......,N，N为自然数；步骤S4中设瞬变电磁场分量衰减曲线的时窗数为N，观测周期数为M，加窗后的第j个周期中的场值记为EMF^j(t_i)，其中i＝1,......,N；j＝1,......,M，M、N均为自然数；并设阈值为ε，则M个周期叠加处理如下：①计算不同周期同一时窗中心点场值的平均值：②设定阈值ε，剔除的场值，其中j＝1,......,M；③对余下的场值按①中公式重新计算时窗中心点场值。