CN105116448A - 一种转换波方位各向异性校正方法及装置 - Google Patents
一种转换波方位各向异性校正方法及装置 Download PDFInfo
- Publication number
- CN105116448A CN105116448A CN201510490184.9A CN201510490184A CN105116448A CN 105116448 A CN105116448 A CN 105116448A CN 201510490184 A CN201510490184 A CN 201510490184A CN 105116448 A CN105116448 A CN 105116448A
- Authority
- CN
- China
- Prior art keywords
- wave
- slow
- component data
- shear
- slow shear
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种转换波方位各向异性校正方法及装置,包括利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;利用所述快横波道集和慢横波道集以及预设的滑动时窗,计算每个样点所对应的快横波和慢横波之间随时间变化的时延;利用计算出的时延,对所述慢横波道集进行校正,得到慢横波校正道集;对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到转换波径向分量数据和横向分量校正数据。通过本申请实施例所公开的技术方案,可以实现正交介质中转换波方位各向异性的校正。
Description
技术领域
本申请涉及转换波地震勘探资料处理技术领域,特别涉及一种基于横波分裂分析的转换波方位各向异性校正方法及装置。
背景技术
转换横波分裂分析和校正是转换波方位各向异性处理的一项关键技术。转换横波分裂分析和校正可以是指在横波分裂分析的基础上,对转换波方位各向异性进行校正。利用该方法不仅能改善转换波成像质量,还可以获得地层相关的裂隙发育方向和密度,这可以用于储层预测。
转换波方位各向异性校正可以是指根据横波分裂分析所确定的裂缝发育方位,对转换波进行快、慢波分离,并利用分析的快、慢时差参数对慢波做时差校正,并且在目的层为多层时,可以采用层剥离法来实现转换波方位各向异性校正。
目前,主要是利用互相关法、能量比值法、横向能量最小的双扫描法或径向能量叠加能量最大的双扫描法等角度扫描法来对所采集的转换波地震数据进行横波分裂分析,实现快慢、横波的分离并确定出裂缝发育方位;然后利用所确定的裂缝发育方位,获得裂缝发育方向上的快、慢横波并计算快、慢横波的时延;再根据所得到的固定时延,对慢横波进行各向异性校正,将快横波分量和校正后的慢横波分量反转回原来的径向和横向分量,实现转换波方位各向异性校正。
在实现本申请过程中,发明人发现现有技术中至少存在如下问题:
现有技术中基于横波分裂分析的转换波方位各向异性校正方法,是基于单层或多层水平层状HTI(HorizontalTransverseIsotropy水平横向各向同性)介质模型,但在广泛存在垂直裂隙或定向分布的裂隙的薄互层分布的地区,地下介质的裂缝各向异性变化复杂,这种情况下可以将介质模型等效为正交各向异性介质(可以简称为正交介质)模型,对于这种介质模型在实际应用中很难再用单层或多层HTI介质来描述。因此,需要提供一种新的技术方案,以实现正交各向异性介质中转换波方位各向异性的校正。
发明内容
本申请实施例的目的是提供一种转换波方位各向异性校正方法及装置,以实现正交介质中转换波方位各向异性的校正。
为解决上述技术问题,本申请实施例提供一种转换波方位各向异性校正方法及装置是这样实现的:
本申请实施例提供了一种转换波方位各向异性校正方法,包括:
利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;
将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;
利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延;
利用计算出的所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;
对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
本申请实施例还提供了一种转换波方位各向异性校正装置,该装置包括:
确定单元,用于利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;
变换单元,用于将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;
计算单元,用于利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延;
校正单元,用于利用计算出的所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;
逆变换单元,用于对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
与现有技术相比,本申请实施例通过计算预设滑动时窗内每个样点所对应的快横波和慢横波之间随时间变化的时延;利用所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据,这使得所得到的转换波径向分量数据和横向分量校正数据更加精确,可以反映出裂缝各向异性变化复杂的介质(即正交介质)的特性,因此,这实现了提供适用于正交介质的转换波方位各向异性校正方法的目的。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例所提供的一种转换波方位各向异性校正方法的流程图;
图2是步骤S110所包括的子步骤流程图。
图3是正交各向异性介质中转换横波分裂分析和校正前的径向分量。
图4是正交各向异性介质中转换横波分裂分析和校正后的径向分量。
图5是本申请实施例所提供的一种转换波方位各向异性校正装置的模块图。
具体实施方式
本申请实施例提供一种转换波方位各向异性校正方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图对本申请所述的一种转换波方位各向异性校正方法进行详细的说明。虽然本申请提供了如下述实施例或流程图所述的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中,这些步骤的执行顺序不限于本申请实施例提供的执行顺序。图1是本申请所述转换波方位各向异性校正方法的一种实施例的方法流程图。该方法包括如下步骤:
S110:利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向。
在采集到目标区域中的转换波地震资料后,可以获取目标区域中第一目的层所对应的含有转换波径向分量数据和横向分量数据的地震道集,然后利用所述转换波径向分量数据和横向分量数据,确定第一目的层的裂缝发育方向。
所述转换波地震资料可以包括转换波在野外采集坐标系中各方向(X、Y和Z方向)上的野外地震数据以及观测系统定义文件等。所述观测系统定义文件可以包括观测系统参数,例如炮线距、接收线距、最大偏移距以及覆盖次数等。
所述目标区域可以使整个勘探区域,也可以是勘探区域中的部分区域。所述第一目的层可以是所述目标区域中的上套地层,即目的层埋深比较浅的地层。在所述目标区域中只有一层目的层时,所述第一目的层就是所述目标区域中唯一的目的层。每一目的层对应一套正交各向异性介质。
如图2所示,该步骤具体的可以包括如下子步骤:
S111:获取第一目的层的转换波地震资料。
在所述目标区域中只有一层目的层时,所述获取第一目的层的转换波地震资料可以是直接利用所采集的目标区域中的转换波地震资料;也可以是对所采集的目标区域中的转换波地震资料进行正演来得到第一目的层的转换波地震资料。
在所述目标区域中含有多层目的层时,所述获取第一目的层的转换波地震资料可以是从所采集的目标区域中的转换波地震资料中选取第一目的层的转换波地震资料;也可以是先对所采集的目标区域中的转换波地震资料进行正演,然后在正演得到的转换波地震资料中选取第一目的层的转换波地震资料。
可以利用现有技术中的方法对转换波地震资料进行正演,在此不再赘叙。
S112:对所述转换波地震资料进行预处理,得到含有转换波径向分量数据和横向分量数据的地震道集。
在得到第一目的层的转换波地震资料后,可以对所得到的转换波地震资料进行预处理,得到含有转换波径向分量数据和横向分量数据的地震道集。
所述预处理可以包括坐标旋转、加载观测系统、静校正、波场分离、叠前去噪、振幅处理和/或叠前时间偏移等。其中,
所述坐标旋转可以是指利用炮检方向与检波器X轴方向的夹角,将野外采集的转换波X分量和Y分量旋转到处理坐标,即径向分量和横向分量。
所述加载观测系统可以是指根据观测系统定义文件,确定炮点、检波点、共中心点、转换点之间的几何关系,并将所确定的几何关系存储在地震数据的道头中。
所述静校正可以包括采用基于层位控制的转换波静校正方法计算转换波的长波长静校正量以及根据有效信号的连续性和光滑性,采用互相关法计算转换波的短波长静校正量。
所述波场分离可以是指将所述转换地震波资料中的纵波和横波进行分离。
所述叠前去噪可以是指压制所述转换地震波资料中的各种干扰。
所述振幅处理可以包括球面扩散补偿、地表一致性振幅补偿以及地层吸收衰减等能量补偿。
叠前时间偏移可以是指在保证覆盖次数的条件下,将转换波的横向分量和径向分量分为不同的方位扇区,对每个方位扇区内的横向分量和径向分量进行叠前时间偏移处理,以提高转换波的成像精度。
所述地震道集可以是共转换点道集、共检波点道集,或是完成分方位角的叠前偏移获得分方位角的共成像点道集。所述共转换点道集的共转换点号以及所述共检波点道集的共检波点号均预置在地震数据的道头中。
S113:在所选择的横波分裂分析时窗内,利用所述转换波径向分量数据和横向分量数据,确定所述第一目的层的裂缝发育方向。
在得到含有转换波径向分量数据和横向分量数据的地震道集后,可以在所选择的横波分裂分析时窗内,利用所述地震道集中的转换波径向分量数据和横向分量数据,分析横波分裂方位角,确定所述第一目的层的裂缝发育方向。
可以利用现有技术中的径向能量最小法来确定第一目的层的裂缝发育方位角,但并不限于该方法,具体分析过程不再赘叙。
S120:将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集。
在确定出第一目的层的裂缝发育方向后,可以将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集。
在一具体实现方式中,可以对所述转换波径向分量数据和横向分量数据进行三角旋转变换,将所述地震道集中的所有方位上的转换波径向分量数据和横向分量数据旋转到裂缝发育方向上,得到对应的快横波道集和慢横波道集。
对所述转换波径向分量数据和横向分量数据进行三角旋转变换可以用下述公式表示:
上式中,θ为炮检方向与裂缝发育方向之间的夹角;S1和S2分别为单个地震道上的快横波和慢横波;R和T分别为转换波的径向分量和横向分量。
所得到的所有地震道所有方位上的快横波分量数据构成了快横波道集,所得到的所有地震道所有方位上的慢横波分量数据构成了慢横波道集。
需要说明的是,所述预设变换并不限于三角旋转变换,也可以是其他可行的坐标变换。
S130:利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延。
在得到裂缝发育方向上的快横波道集和慢横波道集后,可以将快横波道集中所有地震道的快横波进行叠加,得到快波地震道;将慢横波道集中所有地震道的慢横波进行叠加,得到慢波地震道;然后利用预设的滑动时窗,逐点计算所述快波地震道和所述慢波地震道之间的时延。
所述滑动时窗可以是整个横波分裂分析时窗,也可以是横波分裂分析时窗内的子时窗。
所述时延可以是指随着时间变化而变化的时间延迟。所述时间延迟是指接收快横波和接收慢横波之间的时间差。
所述快波地震道和所述慢波地震道的计算公式可以表示如下:
上式中,为快波地震道;为慢波地震道。
逐点计算所述快波地震道和所述慢波地震道之间的时延可以是采用互相关法,利用滑动时窗,计算第一目的层内每个样点所对应的快波地震道和所述慢波地震道之间的时延。具体的,
可以利用下述公式(4)计算每个样点所对应的快波地震道和所述慢波地震道的互相关函数的极大值;然后利用互相关函数极大值所对应的k值,计算每个样点所对应的快波地震道和所述慢波地震道之间的时延,即快横波和慢横波之间的时延。
可以用下述公式计算每个样点所对应的快波地震道和所述慢波地震道之间的时延:
V(t)=k'*dt
(5)
t=i*dt
上式(4)中,fi(k)为第i个样点的互相关函数,L和M为横波分裂分析时窗的样点范围,i和j为正整数,i的取值范围为1~NT,NT为第一目的层内的样点个数,k的取值范围为0~K。
上式(5)中,V(t)为时延,k'为互相关函数极大值所对应的k值,dt为采样间隔,t为时间。
L和M的计算公式如下:
上面两式中,N为滑动时窗内的样点个数。
时窗内样点个数为N通常可以采用下式计算:
N=2K+1(8)
K为预设参数,其值可以由需要预先定义的快横波和慢横波之间的最大时延τmax来确定,即
K=τmax/dt(9)
根据上面式(5)可以看出,每个样点所对应的快波地震道和所述慢波地震道之间的时延是随时间变化的,而不是固定的,这可以提高后续计算结果的精确度。
S140:利用计算出的所述时延,对所述慢横波道集进行校正,得到慢横波校正道集。
在计算出每个样点所对应的快横波和慢横波之间随时间变化的时延之后,可以利用所得到的时延对慢横波道集中各方位上的慢横波分量数据进行补偿,将各方位上的慢横波分量数据都校正到快横波分量数据所在时间上,得到慢横波分量校正数据。这样慢横波分量校正数据的“接收时间”与快横波分量数据的接收时间相同。需要说明的是,此处“接收时间”并不是慢横波分量数据的真实接收时间,而是指成像时所显示的时间。
每个方位的慢横波分量校正数据可以表示如下:
所有方位上的慢横波分量数据构成了慢横波校正道集。
S150:对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
在得到慢横波校正道集后,可以对所述快横波道集中各方位上的快横波分量数据以及所述慢横波校正道集中各方位上的慢横波分量校正数据进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
对所述快横波分量数据以及所述慢横波分量校正数据进行所述预设变换的逆变换可以用公式表示如下:
上式中,R'和T'分别为校正后的转换波径向分量数据和横向分量校正数据。
由于快横波和慢横波之间的时延是随时间变化的,而不是固定时延,这可以使计算出的转换波径向分量数据和横向分量校正数据更加精确。此外,由于快横波分量数据和慢横波分量校正数据之间不存在时间差,因此校正后转换波径向分量数据和横向分量校正数据可以在同一时刻成像,这可以提高成像品质。
图3和图4分别示出了正交各向异性介质中转换波方位各向异性校正前和校正后的径向分量。对比这两图可以看出,利用本申请实施例所提供的转换波方位各向异性校正方法对径向分量校正后,各方位上的转换波径向分量之间的时间差大幅减小,提高了成像品质,并且转换波径向分量更加精确。
通过上述步骤可以看出,本申请实施例通过利用滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延;利用所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据,这使得所得到的转换波径向分量数据和横向分量校正数据更加精确,可以反映出裂缝各向异性变化复杂的介质(即正交介质)的特性,因此,这实现了提供适用于正交介质的转换波方位各向异性校正方法的目的。
在所述目标区域中存在多套各向异性介质(即多层目的层)时,例如H套,所述方法还包括以下步骤:
S160:将所述第一目的层的转换波径向分量数据和横向分量校正数据作为第二目的层的转换波径向分量数据和横向分量数据,参照步骤S110-S150,获取第二目的层的转换波径向分量数据和横向分量校正数据。
S170:按照获取第二目的层的转换波径向分量数据和横向分量校正数据的方法,依次获取目标区域中剩余H-2个目的层的转换波径向分量数据和横向分量校正数据。
需要说明的是,第I层目的层的埋深要大于第I-1层目的层的埋深。H和I均为正整数。此外,本申请中大写字母和小写字母的取值可以不同。
上述步骤的具体执行过程可以通过采用现有技术中的层剥离法来实现,在此不再赘叙。
本申请实施例中还提供了一种转换波方位各向异性校正装置,如图5所示,该装置包括确定单元510、变换单元520、计算单元530、校正单元540以及逆变换单元550。其中,确定单元510用于利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;变换单元520用于将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;计算单元530用于利用所获得的快横波道集和慢横波道集,计算预设滑动时窗内每个样点所对应的快横波和慢横波之间随时间变化的时延;校正单元540用于利用计算出的快横波和慢横波之间随时间变化的时延,对所述慢横波道集进行校正,得到慢横波校正道集;逆变换单元550用于对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
在一实施例中,确定单元510可以包括(图中未示出):
获取子单元,用于获取第一目的层的转换波地震资料;
预处理子单元,用于对所述转换波地震资料进行预处理,得到含有转换波径向分量数据和横向分量数据的地震道集;
确定子单元,用于在所选择的横波分裂分析时窗内,利用所述转换波径向分量数据和横向分量数据,确定所述第一目的层的裂缝发育方向。
在一实施例中,计算单元530可以包括(图中未示出):
第一叠加子单元,用于将快横波道集中所有地震道的快横波进行叠加,得到快波地震道;
第二叠加子单元,用于将慢横波道集中所有地震道的慢横波进行叠加,得到慢波地震道;
计算子单元,用于利用预设的滑动时窗,逐点计算所述快波地震道和所述慢波地震道之间的时延。
在一实施例中,所述计算子单元还用于计算第一目的层内每个样点所对应的快波地震道和所述慢波地震道的互相关函数的极大值;然后利用互相关函数极大值所对应的k,计算每个样点所对应的快波地震道和慢波地震道之间的时延。
在一实施例中,所述装置还可以包括(图中未示出):
第一获取单元,用于将所述第一目的层的转换波径向分量数据和横向分量校正数据作为第二目的层的转换波径向分量数据和横向分量数据,参照步骤S110-S150,获取第二目的层的转换波径向分量数据和横向分量校正数据。
第二获取单元,用于按照获取第二目的层的转换波径向分量数据和横向分量校正数据的方法,依次获取目标区域中剩余H-2个目的层的转换波径向分量数据和横向分量校正数据。
上述实施例阐明的装置或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM,ROM,EEPROM,CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种转换波方位各向异性校正方法,其特征在于,包括:
利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;
将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;
利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延;
利用计算出的所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;
对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
2.根据权利要求1所述的方法,其特征在于,所述利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向包括:
获取第一目的层的转换波地震资料;
对所述转换波地震资料进行预处理,得到含有转换波径向分量数据和横向分量数据的地震道集;
在所选择的横波分裂分析时窗内,利用所述转换波径向分量数据和横向分量数据,确定所述第一目的层的裂缝发育方向。
3.根据权利要求1所述的方法,其特征在于,所述利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算每个样点所对应的快横波和慢横波之间随时间变化的时延包括:
将快横波道集中所有地震道的快横波进行叠加,得到快波地震道;
将慢横波道集中所有地震道的慢横波进行叠加,得到慢波地震道;
利用预设的滑动时窗,逐点计算所述快波地震道和所述慢波地震道之间的时延。
4.根据权利要求3所述的方法,其特征在于,所述逐点计算所述快波地震道和所述慢波地震道之间的时延包括:
利用下述公式计算每个样点所对应的快波地震道和所述慢波地震道的互相关函数的极大值:
利用互相关函数极大值所对应的k值以及下述公式,计算每个样点所对应的快波地震道
和慢波地震道之间的时延:
V(t)=k'*dt,
t=i*dt
其中,fi(k)为第i个样点的互相关函数,L和M为横波分裂分析时窗的样点范围,
L和M的计算公式如下:
i和j为正整数,i的取值范围为1~NT,N为滑动时窗内的样点个数,NT为第一目的层内的样点个数,k的取值范围为0~K,K为预设参数,V(t)为时延,k'为互相关函数极大值所对应的k值,dt为采样间隔。
5.根据权利要求1所述的方法,其特征在于,所述利用计算出的所述时延,对所述慢横波道集进行校正包括利用下述公式来对所述慢横波道集中各方位上的慢横波分量数据进行校正,得到慢横波分量校正数据:
上式中,S2(t)为慢横波分量数据;为慢横波分量校正数据,t为时间。
6.一种转换波方位各向异性校正装置,其特征在于,该装置包括:
确定单元,用于利用所获取的转换波径向分量数据和横向分量数据,确定目标区域中第一目的层的裂缝发育方向;
变换单元,用于将所述转换波径向分量数据和横向分量数据进行预设变换,获得所述裂缝发育方向上的快横波道集和慢横波道集;
计算单元,用于利用所获得的快横波道集和慢横波道集以及预设的滑动时窗,计算第一目的层内每个样点所对应的快横波和慢横波之间随时间变化的时延;
校正单元,用于利用计算出的所述时延,对所述慢横波道集进行校正,得到慢横波校正道集;
逆变换单元,用于对所述快横波道集和所述慢横波校正道集进行所述预设变换的逆变换,得到横波分裂校正后所述第一目的层所对应的转换波径向分量数据和横向分量校正数据。
7.根据权利要求6所述的装置,其特征在于,所述确定单元包括:
获取子单元,用于获取第一目的层的转换波地震资料;
预处理子单元,用于对所述转换波地震资料进行预处理,得到含有转换波径向分量数据和横向分量数据的地震道集;
确定子单元,用于在所选择的横波分裂分析时窗内,利用所述转换波径向分量数据和横向分量数据,确定所述第一目的层的裂缝发育方向。
8.根据权利要求6所述的装置,其特征在于,所述计算单元包括:
第一叠加子单元,用于将快横波道集中所有地震道的快横波进行叠加,得到快波地震道;
第二叠加子单元,用于将慢横波道集中所有地震道的慢横波进行叠加,得到慢波地震道;
计算子单元,用于利用预设的滑动时窗,逐点计算所述快波地震道和所述慢波地震道之间的时延。
9.根据权利要求8所述的装置,其特征在于,所述计算子单元还用于:
利用下述公式计算每个样点所对应的快波地震道和所述慢波地震道的互相关函数的极大值:
利用互相关函数极大值所对应的k值以及下述公式,计算每个样点所对应的快波地震道和慢波地震道之间的时延:
V(t)=k'*dt,
t=i*dt
其中,fi(k)为第i个样点的互相关函数,L和M为横波分裂分析时窗的样点范围,
L和M的计算公式如下:
i和j为正整数,i的取值范围为1~NT,N为滑动时窗内的样点个数,NT为第一目的层内的样点个数,k的取值范围为0~K,K为预设参数,V(t)为时延,k'为互相关函数极大值所对应的k值,dt为采样间隔。
10.根据权利要求6所述的装置,其特征在于,所述校正单元还用于利用下述公式来对所述慢横波道集中各方位上的慢横波分量数据进行校正,得到慢横波分量校正数据:
上式中,S2(t)为慢横波分量数据;为慢横波分量校正数据,t为时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490184.9A CN105116448B (zh) | 2015-08-11 | 2015-08-11 | 一种转换波方位各向异性校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510490184.9A CN105116448B (zh) | 2015-08-11 | 2015-08-11 | 一种转换波方位各向异性校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105116448A true CN105116448A (zh) | 2015-12-02 |
CN105116448B CN105116448B (zh) | 2018-03-13 |
Family
ID=54664480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510490184.9A Active CN105116448B (zh) | 2015-08-11 | 2015-08-11 | 一种转换波方位各向异性校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105116448B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168678A (zh) * | 2016-07-15 | 2016-11-30 | 中国矿业大学(北京) | 一种裂隙中传播横波的分离方法及系统 |
CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
CN107894616A (zh) * | 2017-11-14 | 2018-04-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 多分量转换波裂缝预测方法 |
CN108957544A (zh) * | 2018-08-16 | 2018-12-07 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
CN110737020A (zh) * | 2018-07-20 | 2020-01-31 | 中国石油化工股份有限公司 | 页岩气储层裂缝参数计算方法及系统 |
CN111830563A (zh) * | 2019-04-17 | 2020-10-27 | 中国石油天然气集团有限公司 | 纯横波多层裂缝裂缝方向及快慢波时差确定方法及装置 |
CN112558154A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种正交各向异性介质速度建模方法及系统 |
WO2024093855A1 (zh) * | 2022-11-04 | 2024-05-10 | 中国石油天然气集团有限公司 | 一种横波分裂校正方法、系统、存储介质及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4881209A (en) * | 1988-10-12 | 1989-11-14 | Mobil Oil Corporation | Method of processing converted mode seismic data |
CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
CN102053262A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 获得地震转换波的方位速度的方法及处理地震数据的方法 |
CN103576197A (zh) * | 2012-08-02 | 2014-02-12 | 中国石油天然气集团公司 | 一种转换波角道集抽取方法 |
-
2015
- 2015-08-11 CN CN201510490184.9A patent/CN105116448B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4881209A (en) * | 1988-10-12 | 1989-11-14 | Mobil Oil Corporation | Method of processing converted mode seismic data |
CN102053262A (zh) * | 2009-10-29 | 2011-05-11 | 中国石油化工股份有限公司 | 获得地震转换波的方位速度的方法及处理地震数据的方法 |
CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
CN103576197A (zh) * | 2012-08-02 | 2014-02-12 | 中国石油天然气集团公司 | 一种转换波角道集抽取方法 |
Non-Patent Citations (1)
Title |
---|
刘红爱: "转换横波分裂分析及校正技术研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106168678B (zh) * | 2016-07-15 | 2018-08-17 | 中国矿业大学(北京) | 一种裂隙中传播横波的分离方法及系统 |
CN106168678A (zh) * | 2016-07-15 | 2016-11-30 | 中国矿业大学(北京) | 一种裂隙中传播横波的分离方法及系统 |
CN107678063B (zh) * | 2017-09-25 | 2019-11-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
CN107678063A (zh) * | 2017-09-25 | 2018-02-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于等级相关分析的多分量转换波裂缝预测方法 |
CN107894616A (zh) * | 2017-11-14 | 2018-04-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 多分量转换波裂缝预测方法 |
CN107894616B (zh) * | 2017-11-14 | 2020-01-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 多分量转换波裂缝预测方法 |
CN110737020A (zh) * | 2018-07-20 | 2020-01-31 | 中国石油化工股份有限公司 | 页岩气储层裂缝参数计算方法及系统 |
CN108957544A (zh) * | 2018-08-16 | 2018-12-07 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
CN108957544B (zh) * | 2018-08-16 | 2019-07-26 | 中国石油大学(北京) | 近地表各向异性参数的测量方法、装置、地震计及介质 |
CN111830563A (zh) * | 2019-04-17 | 2020-10-27 | 中国石油天然气集团有限公司 | 纯横波多层裂缝裂缝方向及快慢波时差确定方法及装置 |
CN111830563B (zh) * | 2019-04-17 | 2023-06-30 | 中国石油天然气集团有限公司 | 纯横波多层裂缝裂缝方向及快慢波时差确定方法及装置 |
CN112558154A (zh) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | 一种正交各向异性介质速度建模方法及系统 |
CN112558154B (zh) * | 2019-09-25 | 2023-10-13 | 中国石油化工股份有限公司 | 一种正交各向异性介质速度建模方法及系统 |
WO2024093855A1 (zh) * | 2022-11-04 | 2024-05-10 | 中国石油天然气集团有限公司 | 一种横波分裂校正方法、系统、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105116448B (zh) | 2018-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105116448A (zh) | 一种转换波方位各向异性校正方法及装置 | |
Tang | Imaging near-borehole structure using directional acoustic-wave measurement | |
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN102928873B (zh) | 基于四维能量聚焦的地面微地震定位方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105301641A (zh) | 一种方位各向异性速度反演的方法及装置 | |
CN105093281A (zh) | 一种反演框架下的地震多波建模方法 | |
CN102692648A (zh) | 一种基于大炮初至的折射波剩余静校正方法 | |
CN102073067A (zh) | 一种提高地震数据初至波自动拾取效率的方法 | |
US7397427B1 (en) | Phase event detection and direction of arrival estimation | |
CN104820242A (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN103852785A (zh) | 地层各向异性的评价方法 | |
CN103308941A (zh) | 一种基于任意广角波动方程的成像方法及装置 | |
CN109696705A (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
CN113238280A (zh) | 一种基于格林函数的地震监测方法 | |
CN104216018A (zh) | 一种基于小道距层析反演的表层结构调查方法 | |
Le Stunff et al. | Depth-domain velocity analysis in VTI media using surface P-wave data: Is it feasible? | |
CN106950600A (zh) | 一种近地表散射面波的去除方法 | |
CN102778691B (zh) | 一种计算检波器组内静校正时差的方法 | |
CN111999770A (zh) | 一种tti介质转换ps波精确束偏移成像方法及系统 | |
CN1308703C (zh) | 用于生成相似性图和计算反射体倾角的地震探测方法及设备 | |
Wu et al. | Offset VSP P wave and converted S wave salt proximity migration study | |
CN111474580B (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 |