CN104898166B - 一种处理地震数据的方法及装置 - Google Patents
一种处理地震数据的方法及装置 Download PDFInfo
- Publication number
- CN104898166B CN104898166B CN201510254138.9A CN201510254138A CN104898166B CN 104898166 B CN104898166 B CN 104898166B CN 201510254138 A CN201510254138 A CN 201510254138A CN 104898166 B CN104898166 B CN 104898166B
- Authority
- CN
- China
- Prior art keywords
- geophone station
- water
- data
- energy
- correction data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种处理地震数据的方法及装置。该方法包括对所采集的地震数据进行校正,获得检波点水深;利用所述检波点水深,获取所述地震数据的水面反射系数;利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。通过利用本申请实施例所公开的一种处理地震数据的方法能构造出确定性反褶积因子。
Description
技术领域
本申请涉及地震探勘处理技术领域,特别涉及一种处理地震数据的方法及装置。
背景技术
利用拖缆进行海上地震探勘时,一般需要将震源和内置检波器的电缆放置于海水面以下的一定深度处。由于海水与空气的接触面(即海平面)是一个良好的反射面,因此在震源激发出地震波后,部分地震波自激发点(即炮点)传播到检波器所在的检波点之后继续上行到海平面,在海平面处发生反射,反射之后地震波又被检波器接收。这种被检波器接收的反射地震波通常称为虚反射波(也可以称为鬼波或伴随波)。
虚反射波的存在会造成地震记录的低频响应,降低有效波成像的分辨率,有时甚至还会产生假的同相轴,给地震解释造成困难。现有技术中,通常采用反褶积技术来压制虚反射波。反褶积技术通过采用一系列假设条件,例如反射系数序列谱白噪化假设、反射系数时间序列随机假设、子波最小相位假设以及子波时不变假设等,来构造反褶积因子,然后利用所构造的反褶积因子对采集的地震数据进行褶积运算,从而实现对虚反射波的压制。但在实际勘探中通过拖缆采集到的地震数据并不一定能满足这些假设条件,导致所构造的反褶积因子具有不确定性,进而导致不能有效的用于压制虚反射波。
发明内容
本申请实施例的目的是提供一种处理地震数据的方法及装置,以实现构造确定性反褶积因子的目的。
为实现上述目的,本申请实施例通过以下技术方案来实现:
本申请提供了一种处理地震数据的方法,包括:
对所采集的地震数据进行校正,获得检波点水深;
利用所述检波点水深,获取对应所述地震数据的水面反射系数;
利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
优选的,所述对所采集的地震数据进行校正,获得检波点水深包括:
利用预先获取的炮点水深对所述地震数据进行校正,得到炮点校正数据;
获取所述炮点校正数据所对应的一组检波点校正能量;
利用所述一组检波点校正能量确定出检波点水深。
优选的,所述获取所述炮点校正数据所对应的一组检波点校正能量包括:
利用获取的所有检波点因子水深对所述炮点校正数据进行校正,得到一组检波点校正数据;
对所述一组检波点校正数据进行能量计算,得到一组检波点校正能量。
优选的,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的所有检波点因子水深中的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
按照获取所述第一检波点校正能量的方法,依次获取所述所有检波点因子水深中剩余检波点因子水深各自所对应的检波点校正能量。
优选的,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
利用预设的第一比例因子,依次计算第一预设参数范围内剩余比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正数据的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正数据;
对所得到的所有检波点校正数据进行能量计算,得到一组检波点校正能量。
优选的,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
利用预设的第一比例因子,依次计算第一预设参数范围内剩余比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正能量的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正能量。
优选的,所述利用所述检波点水深,获取对应所述地震数据的水面反射系数包括:
获得与所述检波点水深对应的检波点校正数据;
获取与所述检波点校正数据对应的一组修正能量;
根据所述一组修正能量,确定水面反射系数。
优选的,所述获取与所述检波点校正数据对应的一组修正能量包括:
获取第二预设参数范围内的N个水面反射记录系数,N为正整数;
依次利用每个水面反射记录系数对所述检波点校正数据进行修正,得到一组修正数据;
对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。
优选的,所述获取与所述检波点校正数据对应的一组修正能量还包括:
获取第二预设参数范围内的N个水面反射记录系数;
利用所述第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
对所述第一修正数据中每个第一样点修正数据进行平方运算,得到第一修正能量;
按照获得所述第一修正能量的方法,依次获得剩余N-1个水面反射记录系数各自所对应的修正能量。
优选的,所述获取与所述检波点校正数据对应的一组修正能量还包括:
利用预设的第一水面反射记录系数对所述检波点校正数据进行修正,得到第一修正数据;
利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正数据的方法,依次获得所述N-1个水面反射记录系数各自对应的N-1个修正数据;
对所获得的每个修正数据进行能量计算,得到一组修正能量。
优选的,所述获取与所述检波点校正数据对应的一组修正能量还包括:
利用预设的第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
对所述第一修正数据中的所有样点修正数据进行平方运算,得到第一修正能量;
利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正能量的方法,依次获得所述N-1个水面反射记录系数各自对应的N-1个修正能量。
本申请实施例还提供了一种处理地震数据的装置,该装置包括:
校正单元,用于对所采集的地震数据进行校正,获得检波点水深;
获取单元,用于利用所述检波点水深,获取所述地震数据的水面反射系数;
构造单元,用于利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
与现有技术相比,本申请实施例通过对实际采集的地震数据进行处理来构造反褶积因子,而不是采用假设条件来构造反褶积因子,这能实现构造确定性反褶积因子的目的。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例所提供的一种处理地震数据的方法的步骤流程图;
图2是一种获得检波点水深的方法的步骤流程图。
图3是一种获得一组检波点校正能量的方法的步骤流程图。
图4是另一种获得一组检波点校正能量的方法的步骤流程图。
图5是另一种获得一组检波点校正能量的方法的步骤流程图。
图6是另一种获得一组检波点校正能量的方法的步骤流程图。
图7是一种获得水面反射记录系数的方法的步骤流程图。
图8是拖缆采集地震数据的原始剖面。
图9是图2中的地震数据经过处理之后的剖面。
图10是本申请实施例所提供的一种处理地震数据的装置的模块图。
具体实施方式
本申请实施例提供一种处理地震数据的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图对本申请实施例所提供的一种处理地震数据的方法进行详细说明。虽然本申请提供了如下述实施例或流程图所述的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤,操作步骤之间的执行顺序并没有限制。如图1所示,本申请实施例所提供的一种处理地震数据的方法包括如下步骤:
S100:对所采集的地震数据进行校正,获得检波点水深;
所采集的地震数据可以是整个研究区域中所有样点的地震数据;也可以是某个预设时窗内所有样点的地震数据。所述地震数据中可以含有虚反射波和有效波。
所述检波点水深可以表示检波点离水平面的真实距离,即检波器的真实沉放深度。
所述对所采集的地震数据进行校正,获得检波点水深可以包括以下步骤,如图2所示:
S101:利用预先获取的炮点水深对所述地震数据进行校正,得到炮点校正数据。
所述炮点水深可以是指地震激发源离水平面的距离,例如气枪的沉放深度。
所述利用预先获取的炮点水深对所采集的地震数据进行校正,得到炮点校正数据可以包括对所述炮点水深进行单程旅行时的Z变换;利用变换后的炮点水深对每个样点的地震数据进行校正,得到炮点校正数据。所述单程旅行时可以是指地震波从炮点到水平面的时间。
所述炮点水深的Z变换可以表示为:
其中,hs为炮点水深,v为预先获取的水速,w为角频率。
对于第j个样点,利用变换后的炮点水深对所采集的地震数据进行校正,得到的炮点样点校正数据可以表示为:
X(tj)s=ZsX(tj) (2)
其中,X(tj)s为第j个样点的炮点样点校正数据,X(tj)为所采集的第j个样点的地震数据,j为正整数。
相应的,所有样点的炮点样点校正数据构成了所述炮点校正数据,其可以表示为
X(t)s=[X(t1)s,X(t2)s...X(tm)s] (3)
其中,X(t)s为炮点校正数据;m为样点数。
S102:获取所述炮点校正数据所对应的一组检波点校正能量。
所述检波点校正能量可以是指对所述炮点校正数据进行校正后所得到数据的平方值。
在一实施例中,如图3所示,所述获取所述炮点校正数据所对应的一组检波点校正能量可以包括以下步骤。:
S1021A:利用获取的所有检波点因子水深对所述炮点校正数据进行校正,得到一组检波点校正数据。
所述检波点因子水深可以表示第一预设参数范围内所有比例因子与预先获取的检波点测量水深进行乘积运算所得到的数值。第i检波点因子水深的计算公式如下:
hci=hrchi (4)
其中,hci为第i检波点因子水深,hr为检波点测量水深,chi为第i比例因子。
所述所有检波点因子水深可以即时获取,也可以预先获取。获取所有检波点因子水深可以是首先获取第一预设参数范围内所有比例因子,然后获取每个比例因子所对应的检波点因子水深;还可以是获取一个比例因子后获取该比例因子的检波点因子水深;然后再依次获取第一预设参数范围内剩余比例因子以及对应的检波点因子水深。
所述获取第一预设参数范围内所有比例因子可以是利用第一预设步长(例如0.1)以及预设的第一比例因子,计算出第一预设参数范围内所有的比例因子。计算第i比例因子的公式可以如下式所示:
chi=ch1+(i-1)Δ (5)
上式中,chi为第i比例因子;ch1为第一比例因子;Δ为第一预设步长;i=1~n,n为比例因子个数。
所述第一预设参数范围可以是根据海水涨落变化深度来设定的,例如0.5~1.5。需要说明的是,所有比例因子的取值均在第一预设参数范围之内。
所述利用获取的所有检波点因子水深对所述炮点校正数据进行校正,得到一组检波点校正数据可以是利用所获取的所有检波点因子水深对所述炮点校正数据中每个样点的炮点样点校正数据进行校正,得到一组检波点校正数据。具体的,可以包括对每个检波点因子水深进行Z变换;利用变换后的每个检波点因子水深对每个样点的炮点样点校正数据分别进行正向校正和反向校正,得到相应的检波点样点正向校正数据和检波点样点反向校正数据。
一个样点的检波点样点正向校正数据和检波点样点反向校正数据构成了一个检波点样点校正数据;一个检波点因子水深所对应的所有样点的检波点样点校正数据构成了一个检波点校正数据;n个检波点因子水深所对应的n个检波点校正数据构成了所述一组检波点校正数据。
所述正向校正可以是指增大所述炮点校正数据所在时刻,即对于同一计算水深处的炮点校正数据,其校正后所在时刻相比校正前所在时刻有所增大。所述反向校正可以是指减小所述炮点校正数据所在时刻,即对于同一计算水深处的炮点校正数据,其校正后所在时刻相比校正前所在时刻有所减小。
所述对每个检波点因子水深进行Z变换可以表示为:
其中,hci为第i检波点因子水深,v为水速。
利用变换后的每个检波点因子水深对每个样点的炮点样点校正数据进行正向校正,得到该样点的检波点样点正向校正数据可以表示为:
其中,为第i检波点因子水深所对应的第j个样点的检波点样点正向校正数据,X(tj)s为第j个样点的炮点样点校正数据。
利用变换后的每个检波点因子水深对每个样点的炮点样点校正数据进行反向校正,得到该样点的检波点样点反向校正数据可以表示为:
其中,为第i检波点因子水深所对应的第j个样点的检波点样点反向校正数据。
第i检波点因子水深所对应的所有样点的检波点样点正向校正数据可以表示为:
所有检波点因子水深所对应的检波点正向校正数据X(t)rhd可以表示为:
第i检波点因子水深所对应的所有样点的检波点样点反向校正数据可以表示为:
所有检波点因子水深所对应的检波点反向校正数据X(t)rhu可以表示为:
第i检波点因子水深所对应的所有样点的检波点样点校正数据可以表示为:
所有检波点因子水深所对应的所有样点的检波点校正数据X(t)rh可以表示为:
或
S1022A:对所述一组检波点校正数据进行能量计算,得到一组检波点校正能量。
在得到所有检波点因子水深所对应的一组检波点校正数据后,对所述一组检波点校正数据进行能量计算,得到对应的一组检波点校正能量。具体的,可以在第一预设时窗内依次将n个检波点校正数据中的所有检波点样点正向校正数据和所有检波点样点反向校正数据进行对应相加,然后进行平方运算,得到每个检波点校正数据所对应的检波点校正能量。所得到的n个检波点校正能量构成了所述一组检波点校正能量。例如,计算第i个检波点因子水深所对应的第i检波点校正能量的公式如下式所示:
相应的,所述一组检波点校正能量可以表示如下:
在获得所述一组检波点校正数据和所述一组检波点校正能量后,还可以对检波点因子水深、检波点校正数据和检波点校正能量进行保存,优选的将检波点因子水深、检波点校正数据与检波点校正能量进行关联存储,例如可以将三者保存在同一位置。此外,还可以对检波点因子水深、检波点校正数据与检波点校正能量进行标记,以便于区分。
在另一实施例中,如图4所示,所述获取所述炮点校正数据所对应的一组检波点校正能量还可以包括以下步骤:
S1021B:利用所获取的n个检波点因子水深中的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据。
所述n个检波点因子水深与第一预设参数范围内的n个比例因子相对应。所述n个检波点因子水深可以即时获取,也可以预先获取。获取所述n个检波点因子水深的过程可以参考S1021A,在此不再赘叙。
在得到所有n个检波点因子水深后,可以对这n个检波点因子水深中的第一检波点因子水深进行Z变换,然后利用变换后的第一检波点因子水深对所有样点的炮点校正数据进行校正,得到相应的第一检波点校正数据。
该步骤的具体执行过程,可参考S1021A,在此不再赘述。
S1022B:对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量。
在得到所述第一检波点校正数据后,可以利用公式(17)来得到第一检波点校正能量。
该步骤的具体执行过程,可参考S1022A,在此不再赘述。
S1023B:参照步骤S1021B-S1022B,按照获取所述第一检波点校正能量的方法,依次获取剩余n-1个检波点因子水深各自所对应的检波点校正能量。
在得到第一检波点校正能量后,按照获取所述第一检波点校正能量的方法,重复步骤S1021B-S1022B,依次获取剩余n-1个检波点因子水深所对应的n-1个检波点校正能量,所得到的n个检波点校正能量构成了所述一组检波点校正能量。
在另一实施例中,如图5所示,所述获取所述炮点校正数据所对应的一组检波点校正能量也可以包括以下步骤:
S1021C:利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据。
所述第一检波点因子水深可以预先获取,也可以即时获取。
获取第一检波点因子水深可以是将第一预设参数范围内预设的第一比例因子与预先获取的检波点水深进行乘积运算,得到第一检波点因子水深。计算所述第一因子水深的公式可参考公式(4)。
该步骤的具体执行过程可参考S1022B,在此不再赘述。
S1022C:利用所述第一比例因子,依次计算第一预设参数范围内剩余的n-1个比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正数据的方法,依次获得所述n-1个检波点因子水深各自对应的检波点校正数据。
在得到第一检波点校正数据后,利用第一预设步长和第一比例因子,计算第一预设参数范围内的第二比例因子,然后根据所得到的第二比例因子,计算对应的第二检波点因子水深,再参照步骤S1021C,得到对应的第二检波点校正数据。以此类推,重复该步骤,依次计算出第一预设参数范围内剩余n-2个比例因子以及对应的n-2个检波点校正数据。所得到的第一检波点校正数据、第二检波点校正数据以及后续得到的n-2个检波点校正数据构成了所述一组检波点校正数据。
计算第一预设参数范围内剩余所有比例因子的过程可参考S1021A,在此不再赘叙。
S1023C:对所得到的所有检波点校正数据进行能量计算,得到一组检波点校正能量。
该步骤与S1022A相同,在此不再赘叙。
在另一实施例中,如图6所示,所述获取所述炮点校正数据所对应的一组检波点校正能量也还可以包括以下步骤:
S1021D:利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据。
上述步骤S1021D与步骤S1021C相同,在此不再赘叙。
S1022D:对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量。
可以利用公式(17)来计算所述第一检波点校正能量,即:
S1023D:利用所述第一比例因子,依次计算第一预设参数范围内剩余的n-1个比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正能量的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正能量。
在得到第一检波点校正能量后,利用第一预设步长和第一比例因子,计算第一预设参数范围内的第二比例因子,然后根据所得到的第二比例因子,计算对应的第二检波点因子水深,再参照步骤S1021D,得到对应的第二检波点校正数据;最后参照步骤S1022D,得到对应的第二检波点校正能量。以此类推,重复该步骤,依次计算出第一预设参数范围内的其他n-2个比例因子、n-2个检波点校正数据以及对应的n-2个检波点校正能量。所得到的第一检波点校正数据、第二检波点校正数据以及后续得到的n-2个检波点校正数据构成了所述一组检波点校正数据;所得到的第一检波点校正能量、第二检波点校正能量以及后续得到的n-2个检波点校正能量构成了所述一组检波点校正能量。
计算第一预设参数范围内所有剩余比例因子的过程可参考S1021A,在此不再赘叙。
S103:利用所述一组检波点校正能量获得检波点水深。
在得到所述一组检波点校正数据和所述一组检波点校正能量后,可以利用所述一组检波点校正能量来获得检波点水深。具体的,可以将所述一组检波点校正能量中的每个检波点校正能量进行对比,找出最小检波点校正能量;然后根据所得到的最小检波点校正能量,查找出对应的检波点因子水深。所查找到的检波点因子水深即为检波器所在的真实水深,即检波点水深。
S110:利用所述检波点水深,获取对应所述地震数据的水面反射系数。
水面反射系数可以表示地震数据中虚反射波与一次波之间的振幅比值。由于受天气、气温等外界环境变化的影响,测量得到的水面反射记录系数可能会在某一范围内产生变化,所以实际勘探中测量得到的水面反射记录系数并不能代表真实的水面反射系数,因此需要获取真实的水面反射系数。
在获得检波点水深后,可以利用所述检波点水深,获得对应的真实水面反射系数。如图7所示,具体的,可以包括以下步骤:
S111:获取与所述检波点水深对应的检波点校正数据。
所述获取与所述检波点水深对应的检波点校正数据可以是利用所述检波点水深,在存储结果中查找到对应的检波点校正数据,也还可以是利用所述检波点水深,再次计算出对应的检波点校正数据。所述检波点校正数据包括对应所述检波点水深的所有样点的检波点样点校正数据,所述检波点样点校正数据包括检波点样点正向校正数据和检波点样点反向校正数据。
具体计算检波点校正数据的过程可以参考S1021A,在此不再赘叙。
S112:获取与所述检波点校正数据对应的一组修正能量。
在获得对应所述检波点水深的检波点校正数据后,可以通过利用第二预设参数范围内所有的水面反射记录系数对所述检波点校正数据进行修正来得到一组修正数据,通过对所述一组修正数据中每个修正数据进行能量计算来得到对应的一组修正能量。具体的,
在一实施例中,S112可以包括如下子步骤:
S1121A:获取第二预设参数范围内的N个水面反射记录系数。
所述第二预设参数范围可以是考虑到外界环境的影响,根据实际采集的地震数据来进行设定的,其可以与所述第一预设参数范围相同,也可以不同,例如其为0.5~1。
所述水面反射记录系数包括水面反射系数。
所述获取第二预设参数范围内N个水面反射记录系数可以是利用第二预设步长和预设的第一水面反射记录系数,计算出第二预设参数范围内所有的水面反射记录系数。所述第二预设步长可以与所述第一预设步长相同,也可以不同。
计算第k水面反射记录系数的公式如下:
rk=r1+(k-1)Δ' (18)
其中,rk为第k水面反射记录系数,r1为第一水面反射记录系数,Δ'为第二预设步长,k为1~N,N为正整数。所有N个水面反射记录系数的取值均在所述第二预设参数范围内。
可以利用上述公式(18)来计算出所有的水面反射记录系数。
S1122A:利用获取的所有水面反射记录系数对所述检波点校正数据进行修正计算,得到一组修正能量。
在一实施例中,所述利用所获取的所有水面反射记录系数对所述检波点校正数据进行修正计算,得到一组修正能量可以是在得到所有的水面反射记录系数后,依次利用每个水面反射记录系数对所述检波点校正数据进行修正,得到一组修正数据,然后可以在第二预设时窗内对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。所述第二预设时窗可以与所述第一预设时窗相同,也可以不同。
所述利用每个水面反射记录系数对所述检波点校正数据进行修正,得到一组修正数据可以是依次将每个水面反射记录系数与所述检波点校正数据中每个样点的检波点样点正向校正数据进行乘积运算,对应每个水面反射记录系数均得到一组样点正向修正数据;然后依次将所得到的每个样点正向修正数据与所述检波点校正数据中对应的检波点样点反向校正数据相加,对应每个水面反射记录系数均得到一组样点修正数据,即一个修正数据。所有水面反射记录系数所对应的所有修正数据构成了所述一组修正数据。
第j个样点的第k样点正向修正数据的计算公式可以表示如下:
其中,为第j个样点的第k样点正向修正数据;为第j个样点的检波点样点正向校正数据。
相应的,第j个样点的第k样点修正数据的计算公式为:
其中,为第j个样点的第k样点修正数据;为第j个样点的检波点样点反向校正数据。
第k修正数据是所有样点的第k样点修正数据的组合,其可以表示为:
所述一组修正数据X(t)rdru可以表示为:
上面两式中,M为第二预设时窗内的样点数,N为水面反射记录系数的个数。需要说明的是M和m、N和n的数值各自可以分别不同。
所述在第二预设时窗内对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量可以是针对每个修正数据,将该修正数据在第二预设时窗内所含有的所有样点的样点修正数据先各自进行平方运算,然后将每个样点修正数据的平方运算结果相加,得到该修正数据所对应的修正能量。所有修正数据所对应的所有修正能量,构成了所述一组修正能量。
可以利用下述公式来得到所述一组修正能量中的每个修正能量:
其中,为第k修正能量。
在另一实施例中,所述利用所获取的所有水面反射记录系数对所述校正数据进行修正计算,得到一组修正能量也可以是在得到所有的水面反射记录系数后,利用所述第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据,然后在第二预设时窗内对所述第一修正数据中每个第一样点修正数据进行平方运算,得到第一修正能量,再按照获得第一修正数据和第一修正能量的方法,依次获得其他水面反射记录系数所对应的修正数据和修正能量。所获得的所有修正数据构成了所述一组修正数据,所获得的所有修正能量构成了所述一组修正能量。
在本实施例中,同样可以利用公式(19)-(23)来得到所述第一修正数据、所述第一修正能量以及其他水面反射记录系数所对应的修正数据和修正能量。具体计算过程可参考该步骤中的上一实施例,在此不再赘叙。
在另一实施例中,S121还可以包括如下子步骤:
S1121B:利用预设的第一水面反射记录系数对所述检波点校正数据进行修正,得到第一修正数据。
可以利用公式(19)-(21)来得到第一修正数据,所得到的第一修正数据表示如下
S1122B:利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正数据的方法,依次获得对应的N-1个修正数据。
在得到第一修正数据后,利用第而预设步长和第一水面反射记录系数,计算第而预设参数范围内的第二水面反射记录系数,然后根据所得到的第二水面反射记录系数,参照S1211’,得到对应的第二修正数据。以此类推,重复该步骤,依次计算出第二预设参数范围内的其他N-2个比例因子以及对应的N-2个修正数据。所得到的第一修正数据、第二修正数据以及后续得到的N-2个修正数据构成了所述一组修正数据。
可以利用公式(18)来计算出第二预设参数范围内剩余N-1个水面反射记录系数,具体过程可参考S1121A,在此不再赘叙。
S1123B:对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。
该步骤的具体过程可参考S1121A,在此不再赘叙。
在另一实施例中,S112还可以包括如下子步骤:
S1121C:利用预设的第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据。
该步骤与S1121B相同,在此不再赘叙。
S1122C:对所述第一修正数据中的所有样点修正数据进行平方运算,得到第一修正能量。
可以根据公式(23)来计算第一修正能量,所得到的第一修正能量可以表示如下:
该步骤的具体过程可参考S1122A,在此不再赘叙。
S1123C:利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正能量的方法,依次获得对应的N-1个修正能量。
在得到第一修正数据后,利用第而预设步长和第一水面反射记录系数,计算第而预设参数范围内的第二水面反射记录系数,然后根据所得到的第二水面反射记录系数,参照S1121C,得到对应的第二修正数据,再参照S1122C,得到对应的第二修正能量。以此类推,重复该步骤,依次计算出第二预设参数范围内的其他N-2个比例因子以及对应的N-2个修正能量。所得到的第一修正能量、第二修正能量以及后续得到的N-2个修正能量构成了所述一组修正能量。
可以利用公式(18)来计算出第二预设参数范围内剩余N-1个水面反射记录系数,具体过程可参考S1121,在此不再赘叙。
S113:根据所述一组修正能量,获得对应所述地震数据的水面反射系数。
在得到所述一组修正能量后,将所述一组修正能量中的每个修正能量进行对比,找出这组修正能量中的最小修正能量;然后根据所得到的最小修正能量,查找到对应该最小修正能量的水面反射记录系数。所查找到的水面反射记录系数即为对应所述地震数据的真实水面反射系数。
S120:根据所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
在得到水面反射系数后,可以利用所述水面反射系数来构造鬼波因子,然后利用所构造的鬼波因子在频率域中计算出反褶积因子。所得到的反褶积因子可以用于对校正后的地震数据进行褶积运算,压制所述地震数据中的虚反射波。
所构造的鬼波因子Wg以及计算出的反褶积因子Dcong可以分别表示如下:
Wg=1-Z2rg (26)
其中,hc为检波点水深;rg为对应所述地震数据的水面反射系数。
将所得到的反褶积因子对校正后的地震数据进行褶积运算,压制所述地震数据中的虚反射波可以是将所述反褶积因子与在步骤S111中所获得的检波点校正数据进行褶积运算,优选的是将所述反褶积因子与所述检波点校正数据中所有样点的检波点样点正向校正数据进行褶积运算,得到运算结果。进行褶积运算的公式如下:
其中,X(t)out为运算结果;为检波点样点正向校正数据。
通过上述步骤可以看出,本申请实施例通过对采集的地震数据进行校正来获取检波点水深,然后通过对所述检波点水深所对应的检波点校正数据进行修正来获取水面反射系数,最后通过利用所获取的水面反射系数来构造用于压制虚反射波的反褶积因子。由于所构造的反褶积因子是通过对实际采集的地震数据进行处理来得到的,而不是采用假设条件来得到的,因而具有确定性。此外,由于所构造的反褶积因子是通过对实际采集到的地震数据进行处理得到的,因此其能反映出虚反射波的真实存在情形,所以利用该反褶积因子来对地震数据进行褶积运算,能更好的实现有效压制虚反射波的目的。
下面以本申请实施例所提供的处理地震数据的方法在实际探勘区域中的应用来说明该方法的实际应用效果。
图8和图9分别为拖缆采集地震数据的原始剖面以及利用本申请实施例所提供的反褶积因子对地震数据进行褶积运算后所得到的剖面。通过对比这两幅图可以看出,利用本申请实施例所提供的数据处理方法处理所采集的地震数据之后,所述地震数据中的虚反射(即鬼波)得到了有效的压制。
本申请实施例还提供了一种处理地震数据的装置,如图10所示。该装置包括:校正单元210、获取单元220以及构造单元230。其中,校正单元210用于对所采集的地震数据进行校正,获得检波点水深;获取单元220用于利用所述检波点水深,获取对应所述地震数据的水面反射系数;构造单元230用于利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
在一实施例中,校正单元210可以包括(图中未示出):
第一校正子单元,用于利用预先获取的炮点水深对所述地震数据进行校正,得到炮点校正数据;
能量获取子单元,用于获取所述炮点校正数据所对应的一组检波点校正能量;
真实水深获取子单元,用于利用所述一组检波点校正能量获得检波点水深。
在一实施例中,所述能量获取单元可以包括:
第二校正子单元,用于利用获取的所有检波点因子水深对所述炮点校正数据进行校正,得到一组检波点校正数据;
一组能量计算子单元,用于对所述一组检波点校正数据进行能量计算,得到一组检波点校正能量。
在另一实施例中,所述能量获取子单元还可以包括:
第三校正子单元,用于利用所获取的n个检波点因子水深中的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
第一能量计算子单元,用于对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
第二能量计算子单元,用于按照获取所述第一检波点校正能量的方法,依次获取剩余的n-1个检波点因子水深各自所对应的检波点校正能量。
在另一实施例中,所述能量获取子单元还可以包括:
第三校正子单元,用于利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
第四校正子单元,用于利用预设的第一比例因子,依次计算第一预设参数范围内剩余的n-1个比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正数据的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正数据;
一组能量计算子单元,用于对所得到的所有检波点校正数据进行能量计算,得到一组检波点校正能量。
在另一实施例中,所述能量获取子单元还可以包括:
第三校正子单元,用于利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
第一能量计算子单元,用于对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
第二能量计算子单元,用于利用预设的第一比例因子,依次计算第一预设参数范围内剩余的n-1个比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正能量的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正能量。
在一实施例中,真实水深获取子单元可以包括:
第一获取子单元,用于获取与所述检波点水深对应的检波点校正数据;
第二获取子单元,用于获取与所述检波点校正数据对应的一组修正能量;
第三获取子单元,用于根据所述一组修正能量,获得水面反射系数。
在一实施例中,所述第二获取子单元可以包括:
因子获取子单元,用于获取第二预设参数范围内的N个水面反射记录系数,N为正整数;
修正数据获取子单元,用于依次利用每个水面反射记录系数对所述检波点校正数据进行修正,得到一组修正数据;
修正能量获取子单元,用于对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。
在另一实施例中,所述第二获取子单元还可以包括:
因子获取子单元,用于获取第二预设参数范围内的N个水面反射记录系数;
第一修正数据获取子单元,用于利用所述第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
第一修正能量获取子单元,用于对所述第一修正数据中每个第一样点修正数据进行平方运算,得到第一修正能量;
第二修正能量获取子单元,用于按照获得所述第一修正能量的方法,依次获得剩余N-1水面反射记录系数各自所对应的修正能量。
需要说明的是,第一修正能量获取子单元和第二修正能量获取子单元可以构成上述修正能量获取子单元。
在另一实施例中,所述第二获取子单元还可以包括:
第一修正数据获取子单元,用于利用预设的第一水面反射记录系数对所述检波点校正数据进行修正,得到第一修正数据;
第二修正数据获取子单元,用于利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正数据的方法,依次获得对应的N-1个修正数据;
修正能量获取子单元,用于对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。
在另一实施例中,所述第二获取子单元还可以包括:
第一修正数据获取子单元,用于利用预设的第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
第一修正能量获取子单元,用于对所述第一修正数据中的所有样点修正数据进行平方运算,得到第一修正能量;
第三修正能量获取子单元,用于利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正能量的方法,依次获得对应的N-1个修正能量。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (12)
1.一种处理地震数据的方法,其特征在于,包括:
对所采集的地震数据进行校正,获得检波点水深;
利用所述检波点水深,获取对应所述地震数据的水面反射系数;
利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
2.根据权利要求1所述的方法,其特征在于,所述对所采集的地震数据进行校正,获得检波点水深包括:
利用预先获取的炮点水深对所述地震数据进行校正,得到炮点校正数据;
获取所述炮点校正数据所对应的一组检波点校正能量;
利用所述一组检波点校正能量确定出检波点水深。
3.根据权利要求2所述的方法,其特征在于,所述获取所述炮点校正数据所对应的一组检波点校正能量包括:
利用获取的所有检波点因子水深对所述炮点校正数据进行校正,得到一组检波点校正数据;
对所述一组检波点校正数据进行能量计算,得到一组检波点校正能量。
4.根据权利要求2所述的方法,其特征在于,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的n个检波点因子水深中的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据,n为正整数;
对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
按照获取所述第一检波点校正能量的方法,依次获取剩余的n-1个检波点因子水深各自所对应的检波点校正能量。
5.根据权利要求2所述的方法,其特征在于,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
利用预设的第一比例因子,依次计算第一预设参数范围内剩余比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正数据的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正数据;
对所得到的所有检波点校正数据进行能量计算,得到一组检波点校正能量。
6.根据权利要求2所述的方法,其特征在于,所述获取所述炮点校正数据所对应的一组检波点校正能量还包括:
利用所获取的第一检波点因子水深对所述炮点校正数据进行校正,得到第一检波点校正数据;
对所述第一检波点校正数据进行能量计算,得到第一检波点校正能量;
利用预设的第一比例因子,依次计算第一预设参数范围内剩余比例因子所对应的n-1个检波点因子水深,并按照获得所述第一检波点校正能量的方法,依次获得所述n-1个检波点因子水深各自所对应的检波点校正能量。
7.根据权利要求1所述的方法,其特征在于,所述利用所述检波点水深,获取对应所述地震数据的水面反射系数包括:
获得与所述检波点水深对应的检波点校正数据;
获取与所述检波点校正数据对应的一组修正能量;
根据所述一组修正能量,确定水面反射系数。
8.根据权利要求7所述的方法,其特征在于,所述获取与所述检波点校正数据对应的一组修正能量包括:
获取第二预设参数范围内的N个水面反射记录系数,N为正整数;
依次利用每个水面反射记录系数对所述检波点校正数据进行修正,得到一组修正数据;
对所述一组修正数据中每个修正数据进行能量计算,得到一组修正能量。
9.根据权利要求7所述的方法,其特征在于,所述获取与所述检波点校正数据对应的一组修正能量还包括:
获取第二预设参数范围内的N个水面反射记录系数;
利用第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
对所述第一修正数据中每个第一样点修正数据进行平方运算,得到第一修正能量;
按照获得所述第一修正能量的方法,依次获得剩余的N-1个水面反射记录系数各自所对应的修正能量。
10.根据权利要求7所述的方法,其特征在于,所述获取与所述检波点校正数据对应的一组修正能量还包括:
利用预设的第一水面反射记录系数对所述检波点校正数据进行修正,得到第一修正数据;
利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正数据的方法,依次获得所述N-1个水面反射记录系数各自所对应的N-1个修正数据;
对所获得的每个修正数据进行能量计算,得到一组修正能量。
11.根据权利要求7所述的方法,其特征在于,所述获取与所述检波点校正数据对应的一组修正能量还包括:
利用预设的第一水面反射记录系数对所述校正数据进行修正,得到第一修正数据;
对所述第一修正数据中的所有样点修正数据进行平方运算,得到第一修正能量;
利用所述第一水面反射记录系数,依次计算第二预设参数范围内剩余的N-1个水面反射记录系数,并按照获得所述第一修正能量的方法,依次获得所述N-1个水面反射记录系数各自所对应的N-1个修正能量。
12.一种处理地震数据的装置,其特征在于,包括:
校正单元,用于对所采集的地震数据进行校正,获得检波点水深;
获取单元,用于利用所述检波点水深,获取所述地震数据的水面反射系数;
构造单元,用于利用所述水面反射系数,构造用于对校正后的地震数据进行褶积运算的反褶积因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510254138.9A CN104898166B (zh) | 2015-05-19 | 2015-05-19 | 一种处理地震数据的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510254138.9A CN104898166B (zh) | 2015-05-19 | 2015-05-19 | 一种处理地震数据的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898166A CN104898166A (zh) | 2015-09-09 |
CN104898166B true CN104898166B (zh) | 2017-09-15 |
Family
ID=54030926
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510254138.9A Active CN104898166B (zh) | 2015-05-19 | 2015-05-19 | 一种处理地震数据的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104898166B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425297B (zh) * | 2015-10-29 | 2018-02-02 | 中国石油天然气集团公司 | 一种压制虚反射信号的方法及装置 |
CN105301656B (zh) * | 2015-10-29 | 2018-02-02 | 中国石油天然气集团公司 | 一种压制虚反射信号的方法及装置 |
CN107193045B (zh) * | 2017-07-07 | 2019-02-15 | 中国石油天然气集团公司 | 一种地震数据处理方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2961316A1 (fr) * | 2010-06-10 | 2011-12-16 | Cggveritas Services Sa | Procede de traitement de donnees sismiques marines |
CN102937720B (zh) * | 2011-08-15 | 2016-09-28 | 中国石油化工股份有限公司 | 井控提高地震资料分辨率的方法 |
US9671511B2 (en) * | 2012-08-31 | 2017-06-06 | Cgg Services Sas | Horizontal streamer broadband marine seismic acquisition configuration and processing |
-
2015
- 2015-05-19 CN CN201510254138.9A patent/CN104898166B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104898166A (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US4380059A (en) | F-K Filtering of multiple reflections from a seismic section | |
US9575196B2 (en) | Coherent noise attenuation | |
WO2016076917A1 (en) | Creating a high-resolution earth model using seismic tomography and impedance inversion | |
KR20110057124A (ko) | 지진 표면파들의 파형들을 사용하는 토양 특성들의 추정 | |
CN104898166B (zh) | 一种处理地震数据的方法及装置 | |
CN103210323B (zh) | 处理地震数据的系统和方法 | |
WO2011103297A2 (en) | Estimating internal multiples in seismic data | |
Hlebnikov et al. | Noise types and their attenuation in towed marine seismic: A tutorial | |
US20140269186A1 (en) | System and method for isolating signal in seismic data | |
CN112285775B (zh) | 一种海底浅部沉积层剪切波速度探测装置及方法 | |
Oren et al. | An overview of reproducible 3D seismic data processing and imaging using Madagascar | |
CN114185095B (zh) | 一种三维平面波域地震数据多次波压制的方法 | |
CN111257935A (zh) | 一种加快层析反演速度收敛的速度融合方法及处理终端 | |
CN113687417B (zh) | 一种三维叠前地震数据层间多次波预测和压制方法 | |
WO2014176116A1 (en) | Surface multiple prediction | |
CN112666610B (zh) | 地层有机碳含量计算方法及系统 | |
CN111665574B (zh) | 岩石热学参数测井解释方法及系统 | |
Raeisdana et al. | Oriented NMO correction of VTI data in the absence of near-offset traces | |
Fei et al. | 3D common-reflection-point-based seismic migration velocity analysis | |
Zhang et al. | Using FWI and deep learning to characterize velocity anomalies in crosswell seismic data | |
CN112925020B (zh) | 一种参量阵型浅剖智能能量补偿方法 | |
US20160238727A1 (en) | Removal of an Estimated Acquisition Effect from a Marine Survey Measurement | |
CN113568040B (zh) | 一种时移地震采集数据的可重复性分析方法及系统 | |
CN111624660B (zh) | 海底节点地震数据上下行波场分离的方法及装置 | |
WO2024016572A1 (zh) | 噪声识别方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |