[发明专利]一种基于稀疏反演的滑动扫描谐波压制方法

说明书

技术领域

本发明涉及石油勘探中可控震源滑动扫描采集技术领域，是用于去除可控震源滑动扫描采集数据中存在的谐波污染的一种基于稀疏反演的滑动扫描谐波压制方法。 

背景技术

在陆上采用可控震源进行采集时，由于可控震源与地表的非线性耦合、可控震源本身的非线性因素以及不充分的回馈系统等原因，往往可观测到力信号和接收记录均存在谐波畸变。如果采用降频扫描方式，力信号的基波分量与接收记录互相关，谐波干扰存在于记录的续至位置，那么强初至的谐波干扰往往会影响反射波的识别；而采用升频扫描方式，力信号的基波分量与接收记录互相关，谐波干扰存在于记录的前端，这样强初至的谐波干扰就不会影响反射波。因此在野外施工时，往往采用升频线性扫描方式。 

1970年A.J.Seriff等人提出滑动扫描记录中谐波的性质与特点以及各种压制谐波的尝试，其中大部分针对理论数据进行处理。滑动扫描采集就是在前一组可控震源扫描的数据记录结束之前，后一组准备就绪的可控震源在仪器设定的时间间隔(滑动时间)之后开始扫描，其中滑动时间不能小于记录长度，仪器对接收到的原始母记录进行处理，首先根据TB时间按照每炮的记录长度切分出每炮的原始记录，然后用参考信号与各炮的原始记录相关，输出相关后的单炮记录。因此，它是一种利用多组可控震源进行 连续或重叠扫描的采集技术，因其能够缩短单一振次的采集时间而提高采集效率。 

然而，可控震源常规采集方法通过升频扫描在很大程度上可以减弱谐波污染的影响，但在可控震源滑动扫描作业时则必须正视后续扫描对前一扫描产生的谐波污染影响，被迫在采集效率与数据品质之间寻求平衡——一方面延长滑动时间，减少谐波干扰的影响，但降低了采集效率；另一方面则容忍一定程度的谐波污染，可以显著提高生产效率，但谐波污染使数据品质受到影响。因而，成为影响其广泛应用的瓶颈问题。 

发明内容

本发明目的在于提过一种解决可控震源滑动扫描记录中存在的谐波污染问题，提高数据品质的基于稀疏反演的滑动扫描谐波压制方法。 

本发明通过以下具体步骤实现： 

1)可控震源采集地震数据，滑动扫描相关前数据和力信号，首先从最后一炮开始处理，对最后一炮的力信号利用短时傅立叶变换在时频域进行分离，分离出基波、二次谐波、三次谐波成分，将分离出的基波分量与相关前记录作互相关，得到相关后记录x₀(t)； 

2)以一道为例，从相关后记录中取绝对值最大x_max，从序列中搜寻|x₀(t)|＞αx_max，形成以x_max为中心的一段区域，确定可控震源的强初至； 

其中：x_max为绝对值最大，α为＜1的常数； 

3)利用从力信号中分离出的基波分量进行自相关得到零相位子波，利用该子波在相关后记录中进行匹配追踪，得到初至的稀疏化反射系数，得 到相关前记录强初至的基波分量； 

4)从相关前记录中减去强初至的基波分量得到剩余记录 将其与力信号的1阶谐波s₁(t)进行互相关，得到相关后记录x₁(t)； 

5)对相关后记录x₁(t)，用1阶谐波自相关得到的子波w₁(t)在强初至区域内采用匹配追踪方法，得到一阶谐波的初至系数，从而得到强初至在相关前记录中的1阶谐波分量； 

6)从步骤4)中得到的剩余记录 中减去步骤5)得到的强初至一阶谐波分量，得到另一个剩余记录，按照从步骤4)到步骤5)提取一阶谐波分量的方法，提取出强初至在相关前记录中的2阶、3阶等更高阶谐波分量； 

7)采用步骤3)得到的相关前记录强初至的基波分量和步骤5)和步骤6)得到的1阶、2阶等高阶谐波分量，按照采集的滑动时间，从前炮中在相应位置减去这些谐波分量，得到没有谐波污染的相关前记录； 

8)把步骤7)得到的没有谐波污染的相关前记录和从前炮力信号中提取出的基波分量进行互相关，就得到前炮压制谐波后的相关后正常单炮记录。 

本发明通过谐波压制达到了和模拟的无谐波理论数据相差无几的程度，能很好的把主要的谐波污染成分提取出来，并很好的压制掉，从时频分析结果看，1次、2次一直到5次以内的高频的谐波成分都被明显压制(图2所示)。经过压制处理后的单炮记录(图4所示)与原始单炮比质量品质显著提高，消除了中深层的强的高频谐波污染，这些谐波污染都集中在近偏移距范围(如图5所示)。 

附图说明