一种多震源地震资料处理中的谐波去除方法技术

本发明专利技术提供了一种多震源地震资料处理中的谐波去除方法，属于石油勘探领域。本方法包括：(1)对有谐波干扰的多震源单炮资料，以每个排列为单位，按排序增量进行道重排，重排后谐波在资料上的显示为随机异常状分布；(2)对重排过的炮集资料进行区域异常衰减测试及处理，将呈随机异常状分布的谐波衰减掉；(3)按炮序和道序将步骤(2)得到的炮资料反分选回去，得到的炮集资料即为将谐波去除干净的炮集资料。

【技术实现步骤摘要】

本专利技术属于石油勘探领域，具体涉及，用于油气勘探中多震源同步激发(可控震源同步扫描)地震资料的处理。
技术介绍
多震源同步激发地震技术是利用常规采集技术采集一炮资料的时间采集多炮地震记录，为了达到减少相邻震次间的时间间隔，通过采用多组震元交替扫描来减少搬迁时间即交替扫描方式，其对采集效率的提高是非常明显的。但在一个很短的时间内有多个震源同步激发，同步震源信号之间的相互干涉是难免的。其首要面对的谐波畸变干扰问题。作为一种重要的可控震源高效采集方式，当采用滑动扫描方式时，前一次扫描还没结束，后一次扫描就已开始，由于等待时间较短，也就引发了可控震源采集的新问题，其中之一就是谐波畸变问题。畸变的谐波窜入信号会造成干扰，让正常的信号识别困难甚至出现差错。特别对勘探目标资料信噪比较低或者目标层较弱时，谐波的干扰尤其不能忽视。因此，找到一种高效去除谐波干扰的方法非常重要。目前提的较多的有线性预测法(请参考:张贤达，现代信号处理。北京:清华大学出版社。1994.68-142)，谐波归一化频率法(请参考:孙燕等，一种消除强谐波干扰的有效方法数据采集与处理vol.13.N0.31998)，涉及谐波算子相关法。在处理反面，主要有Martin等利用类似陷波的方法及其他主要通过设计算子预测相关等方法，但各种方法有的计算繁琐，有的谐波去除不干净，而且都有一定的适用条件。另外，在常规处理中，有一种区域异常衰减技术，但这个技术对成片的异常信号是束手无策的，如何把成片出现的谐波信号分散在单炮资料的整个剖面上，利用谐波与周围信号的差异，结合区域异常衰减的方法来消除谐波对资料的影响是现在还没有解决的问题。
技术实现思路
本专利技术的目的在于解决上述现有技术中存在的难题，提供，解决多震源地震资料中饱受谐波干扰的低信噪比资料的处理的问题。本专利技术是通过以下技术方案实现的:—种多震源地震资料处理中的谐波去除方法，包括:(1)对有谐波干扰的多震源单炮资料，以每个排列为单位，按排序增量进行道重排，重排后谐波在资料上的显示为随机异常状分布；(2)对重排过的炮集资料进行区域异常衰减测试及处理，将呈随机异常状分布的谐波衰减掉；(3)按炮序和道序将步骤(2)得到的炮资料反分选回去，得到的炮集资料即为将谐波去除干净的炮集资料。 所述步骤(3)进一步包括:如果得到的炮集资料的谐波去除还不能达到要求，则返回步骤(2)。所述步骤(3)进一步包括:如果得到的炮集资料的谐波去除还不能达到要求，则分析所述排序增量是否合适，设定新的排序增量，然后返回步骤(1)。所述步骤(1)是这样实现的:对有谐波干扰的多震源单炮资料，将道头字defined_trC_p0S加1，按排序增量为20进行求余，即MOD(defined_trc_pos+l, 20)，然后将求取的余统一减个一成为新的道序trace_num.c,即 trace_num.c = MOD(defined_trc_pos+l, 20)-1,重新按炮和新的道序选排，这样实现对每个排列的道序重排。与现有技术相比，本专利技术的有益效果是:利用本方法就可以用常规的区域异常衰减来解决扰人的谐波干扰，省却很多繁琐的运算，而且方便可行。【附图说明】图1-1是多震源单炮资料的一个排列。图1-2是将多震源单炮资料的一个排列按一定方式重排后的显示。图1-3是图1-2得到的资料，通过区域异常处理后的效果。图1-4是去掉谐波后，重新按炮道分选回来的单炮。图2本专利技术的步骤框图。【具体实施方式】下面结合附图对本专利技术作进一步详细描述:作为一种重要的可控震源高效采集方式，滑动扫描技术已较为成熟，并广泛应用。但是米用滑动扫描方式时，由于前一次扫描还没有结束，后一次扫描就已经开始，这样后续扫描所激发的谐波就会污染前面几次扫描的记录，而形成对相邻炮的谐波干扰。为此探求高效可行的谐波干扰压制技术是高效可控震源滑动扫描地震资料数据处理的关键技术之一。在本专利技术中主要是利用谐波干扰的振幅特征及分布规律，进行道转换及异常衰减。当然，道转换的增量及衰减的参数测试，都需要以面临的实际资料进行不同的参数组合。由于谐波在单炮上出现形式是大面积的相对有规律的分布区域，这为按区域异常方式衰减提出了难题。因为区域异常衰减是基于能量统计的地表一致性的去噪手段，可以在共炮点、共检波点，共偏移距和共深度点多个方面进行信号统计。通过均方根振幅、平均振幅、最大绝对值震幅，涉及分析时窗、门槛参数，拾取振幅能量对能量分析计算并分解，通过统计分析能量的差异来压制异常振幅及相关噪音(请参考:苑益华，牛滨华，王焕第等。浅谈地震资料处理中的提高信噪比处理技术。石油地球物理勘探，2005，40(2);168-171)。如何把具有谐波干扰的多震源弹炮资料上的谐波干扰转换成符合区域异常衰减去除干扰的资料分布特点，是本方法的重点。本专利技术将单炮资料在每个排列内抽到重新排序，将谐波干扰分散化、重新分选，区域异常处理；反分选；通过将成片状成块分布的谐波干扰分散化排列，使之适用于区域异常衰减的条件范围，再进行区域异常处理，将处理完成区域异常后的资料反分选回去。通过这样一个过程来完成谐波干扰的去除，运算简单，方便易行。如图2所示，本专利技术的，包括:对需要进行谐波去除处理的单炮，以每个排列为单位，按一定的增量(根据干扰的分布情况，目的是把成堆的干扰分散化，一般选取20左右)进行道重排，即进行道头字操作过的数据，重新按炮和新生成的道序选排分选。对重新排序过的炮集资料(这时，谐波在资料上的显示为随机异常状分布)，进行区域异常衰减测试及处理，也就是将呈随机异常状分布的谐波衰减掉；将完成了区域异常衰减的炮资料反分选回去(按炮序和道序分选)，此时，得到的炮集资料已是将谐波去除干净的炮集资料了。如果此时得到的炮集资料谐波去除还不能达到要求，可以重新测试区域异常衰减参数，甚至再回到第一步，分析排序增量是否合适，循环刚才的步骤，直到满意为止。本专利技术的一个实施例如下:1，将有谐波干扰的多震源单炮资料，以道头字“defined_trC_p0S” (每个排列内的道序)进行操作，将defined_trc_p0S加1，按20 (根据谐波分布情况，一般选20足够了)进行求余，即MOD (defined_trc_pos+l, 20),将求取的余统一减个一成为新的道序trace_num.c,即trace_num.c = MOD(defined_trc_pos+l, 20)-1,通过这个对每个排列的道序重排。2，将进行道头字操作过的数据，重新按炮和新生成的道序选排分选。3，经过重新分选后的单炮资料，此时，谐波的干扰分布为散点随机的，适合运用区域异常衰减的方法了。根据区域异常衰减的特点，进行时窗等参数测试(按照常规处理即可)，及常规的区域异常处理。4，把经过区域异常处理过后的单炮资料，重新按照炮，道序(道头traCe_num)进行分选，返回到开始的炮道集。谐波干扰去除完成。图1-1是多震源单炮资料的一个排列，从图中可以看到明显的谐波干扰黑乎乎的呈片状分布。图1-2是将多震源单炮资料的一个排列按一定方式重排后的显示，从图中可以看到重拍后谐波干扰随机分布在排列里，没有规律，在区域上表现为异常震幅。图1-3是图1-2得到的资料，通过区域异常处理后的效果。可本文档来自技高网...

【技术保护点】
一种多震源地震资料处理中的谐波去除方法，其特征在于：所述方法包括：(1)对有谐波干扰的多震源单炮资料，以每个排列为单位，按排序增量进行道重排，重排后谐波在资料上的显示为随机异常状分布；(2)对重排过的炮集资料进行区域异常衰减测试及处理，将呈随机异常状分布的谐波衰减掉；(3)按炮序和道序将步骤(2)得到的炮资料反分选回去，得到的炮集资料即为将谐波去除干净的炮集资料。

【技术特征摘要】

【专利技术属性】
技术研发人员：陈胜红，佘德平，武永山，杨丽，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11