CN106646611A - 一种纵波和转换波地震数据时间匹配的方法及装置 - Google Patents
一种纵波和转换波地震数据时间匹配的方法及装置 Download PDFInfo
- Publication number
- CN106646611A CN106646611A CN201611122196.7A CN201611122196A CN106646611A CN 106646611 A CN106646611 A CN 106646611A CN 201611122196 A CN201611122196 A CN 201611122196A CN 106646611 A CN106646611 A CN 106646611A
- Authority
- CN
- China
- Prior art keywords
- wave
- data
- converted
- converted wave
- compressional
- 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
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000007906 compression Methods 0.000 claims abstract description 31
- 230000006835 compression Effects 0.000 claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 241000761557 Lamina Species 0.000 claims description 3
- 239000003550 marker Substances 0.000 abstract 2
- 238000003860 storage Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000005096 rolling process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000010979 ruby Substances 0.000 description 1
- 229910001750 ruby Inorganic materials 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种纵波和转换波地震数据时间匹配的方法及装置。所述方法包括:计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比;利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据;根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据;基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。本申请实施例提供的一种纵波和转换波地震数据时间匹配的方法及装置,可以提高目标工区中小层进行纵波和转换波地震数据时间匹配的精度。
Description
技术领域
本申请涉及石油地球物理勘探技术领域,特别涉及一种纵波和转换波地震数据时间匹配的方法及装置。
背景技术
转换波地震勘探开始于二十世纪八十年代,直到二十世纪九十年代中后期,随着三分量数字检波器、万道地震仪器和海量数据存储技术的发展和采集技术的进步,转换波地震数据的采集成本不断下降,转换波地震勘探技术得到越来越广泛的应用和发展。转换波是下行的地震纵波在地下层界面发生发射、由入射的下行地震纵波转换成的出射的上行地震横波。基于纵波和转换波地震数据联合解释的纵横波地震数据解释在岩性识别、地质体空间成像、流体检测和储层裂缝描述中也起到越来越重要的作用。然而,进行纵波和转换波地震数据联合解释的前提是完成同一地质层位在纵波和转换波地震数据上的时间匹配工作。
纵波和转换波地震数据时间匹配是将转换波双程反射传播时间转换为纵波双程反射传播时间刻度,或者将纵波双程反射传播时间转换为转换波双程反射传播时间刻度,以便于纵波与转换波地震剖面的联合对比分析。一旦将目的层位的反射同相轴分别在纵波和转换波剖面上拾取完成后,即可以压缩转换波剖面,使拾取的转换波同相轴位于相应的纵波同相轴的时间刻度上,时间匹配的精度取决于层位拾取的可靠性。时间匹配精度的高低将决定着联合解释和联合反演的成败。
现有技术中主要是采用属性驱动的方法进行纵波和转换波地震数据时间匹配,其主要过程是确定纵波地震数据和转换波地震数据之间最大相关性的属性,利用该属性分别拾取纵波地震数据和转换波地震数据的标志层的层位信息,得到目的工区内各个标志层上的纵波地震数据和转换波地震数据的初始速度比,以初始速度比为参考中心,对最大相关性的属性和转换波地震数据进行滑动时窗速度比扫描,确定各个标志层上的纵波地震数据和转换波地震数据的目标速度比,基于目标速度比,将转换波地震数据压缩到纵波地震数据的时间刻度上,完成在目的工区内各个标志层上纵波地震数据和转换波地震数据的时间匹配。
发明人发现现有技术中至少存在如下问题:现有技术主要是针对标志层进行纵波和转换波地震数据时间匹配。而对于标志层之间的小层,由于其反射信号较弱以及噪音的干扰,现有技术所采用的方法很难分辨小层的反射特征,从而很难确定在小层上纵波和转换波之间的相关属性,以至于很难对小层进行精确地纵波和转换波地震数据时间匹配。
发明内容
本申请实施例的目的是提供一种纵波和转换波地震数据时间匹配的方法及装置,以提高对目标工区中小层进行纵波和转换波地震数据时间匹配的精度。
为解决上述技术问题,本申请实施例提供一种纵波和转换波地震数据时间匹配的方法及装置是这样实现的:
一种纵波和转换波地震数据时间匹配的方法,包括:
计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比;
利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据;
根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据;
基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
优选方案中,所述根据纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据,包括:
分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部;
将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部;
根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值;
根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
优选方案中,所述目标工区包括:至少两个标志层和各标志层之间的小层。
优选方案中,所述基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据,包括:
针对所述至少一个小层中的一个小层,对所述纵波属性数据和所述转换波属性数据进行时间匹配处理,确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量;
根据所述时移量,计算所述纵波地震数据和第一匹配的转换波地震数据在所述小层上的速度比;
利用所述小层上的速度比对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
优选方案中,所述时间匹配处理包括:动态时间调整。
优选方案中,所述对所述纵波属性数据和所述转换波属性数据进行动态时间调整包括:
计算所述目标纵波属性数据和所述转换波属性数据之间的对齐误差;
根据所述对齐误差计算累积距离;
基于所述累积距离,确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量。
优选方案中,所述计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比,包括:
在所述纵波地震数据和所述转换波地震数据上分别确定所述目标工区中标志层对应的纵波层位时间和转换波层位时间;
根据所述纵波层位时间和所述转换波层位时间,计算所述目标工区标志层上的纵波地震数据和转换波地震数据的速度比。
优选方案中,第一压缩处理包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述转换波地震数据进行压缩;具体地,所述压缩后的转换波地震数据的时间等于所述纵波地震数据的时间。
优选方案中,第二压缩处理包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述第一匹配的转换波地震数据进行压缩;具体地,压缩后的第一匹配的转换波地震数据的时间等于所述纵波地震数据的时间。
一种纵波和转换波地震数据时间匹配的装置,所述装置包括:标志层速度比计算模块、第一压缩处理模块、属性数据确定模块和第二压缩处理模块;其中,
所述标志层速度比计算模块,用于确定目标工区中标志层上的纵波地震数据和转换波地震数据的速度比;
所述第一压缩处理模块,用于利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据;
所述属性数据确定模块,用于根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据;
所述第二压缩处理模块,用于基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
优选方案中,所述属性数据确定模块,包括:复信号虚部确定模块、复信号实部确定模块、幅值确定模块和波属性数据确定模块;其中,
所述复信号虚部确定模块,用于分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部;
所述复信号实部确定模块,用于将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部;
所述幅值确定模块,用于根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值;
所述波属性数据确定模块,用于根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
本申请实施例提供了一种纵波和转换波地震数据时间匹配的方法及装置,分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特(Hilbert)变换,得到变换后的纵波属性数据和转换波属性数据,这些数据可以突出小层的反射特征,使小层对应弱同相轴信号变强,提高了小层的对比精度,从而可以提高对目标工区中小层上的纵波和转换波地震数据时间匹配的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种纵波和转换波地震数据时间匹配的方法实施例的流程图;
图2是本申请实施例中纵波地震数据和转换波地震数据的地震剖面示意图;
图3是本申请实施例中第一匹配的转换波地震数据的地震剖面示意图;
图4是本申请实施例中纵波属性数据和第一匹配的转换波属性数据的地震剖面示意图;
图5是本申请实施例中目标转换波地震数据的地震剖面示意图;
图6是本申请纵波和转换波地震数据时间匹配的装置实施例的组成结构图;
图7是本申请纵波和转换波地震数据时间匹配的装置实施例中属性数据确定模块的组成结构图。
具体实施方式
本申请实施例提供一种纵波和转换波地震数据时间匹配的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请一种纵波和转换波地震数据时间匹配的方法实施例的流程图。如图1所示,所述纵波和转换波地震数据时间匹配的方法,包括以下步骤。
步骤S101:计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比。
所述目标工区可以包括:至少两个标志层和各标志层之间的小层。
具体地,可以分别在所述纵波地震数据和所述转换波地震数据上确定所述目标工区中标志层对应的纵波层位时间和转换波层位时间。根据所述纵波层位时间和所述转换波层位时间,可以计算所述目标工区中标志层上的纵波地震数据和转换波地震数据速度比。进一步地,可以采用下述公式计算所述目标工区的纵波地震数据和转换波地震数据在所述目标工区中标志层上的速度比:
公式(1)中,TPP表示所述纵波层位时间,TPS表示所述转换波层位时间,γ表示所述目标工区中标志层上的纵波地震数据和转换波地震数据的速度比。
进一步地,所述分别在所述纵波地震数据和所述转换波地震数据上确定所述目标工区中标志层对应的纵波层位时间和转换波层位时间,可以包括以下步骤。
1)利用纵波震源激发地震波并利用三分量检波器采集所述纵波地震数据和横波地震数据,对采集的纵波地震数据进行相对振幅保持的高保真处理。
2)对步骤1)中的纵波地震数据中的共反射点道集数据进行不同的偏移距谱均衡处理。
3)对谱均衡处理后的纵波地震数据进行叠前振幅随偏移距变化(AmplitudeVersus Offset,AVO)处理,提取截距和梯度的值,对梯度的值和截距的值进行求差运算得到差异值。
4)对步骤1)中的转换波地震数据和步骤3中得到的差异值进行子波一致性处理。
5)基于子波一致性处理后的差异值和转换波地震数据,进行所述标志层拾取,分别确定所述目标工区中标志层对应的纵波层位时间和转换波层位时间。
例如,图2是本申请实施例中纵波地震数据和转换波地震数据的地震剖面示意图。图2中(a)图和(b)图分别为纵波地震数据和转换波地震数据的地震剖面示意图。图2中(a)图和(b)图的黑色虚线所处位置分别表示所述标志层分别在纵波地震数据和转换波地震数据中的层位时间。如图2所示,所述标志层分别在纵波地震数据和转换波地震数据的层位时间不一致。
步骤S102:利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据。
具体地,可以利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据。进一步地,第一压缩处理可以包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述转换波地震数据进行压缩。具体地,所述压缩后的转换波地震数据的时间等于所述纵波地震数据的时间。例如,图3是本申请实施例中第一匹配的转换波地震数据的地震剖面示意图。如图3所示,与图2中(a)图对比,在所述目标工区中标志层上第一匹配的转换波地震数据与纵波地震数据的层位时间匹配较好。
步骤S103:根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据。
具体地,可以分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特(Hilbert)变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部。可以将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部。根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值。根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
例如,采用下述公式对所述纵波地震数据进行Hilbert变换,得到所述纵波复信号的虚部:
h(t)=H(s(t)) (2)
公式(2)中,s(t)表示所述纵波地震数据,h(t)表示所述纵波复信号的虚部,H(t)表示Hilbert变换。可以将s(t)作为所述纵波复信号的实部。根据所述纵波复信号的实部和虚部,可以采用下述公式得到所述纵波复信号的幅值:
公式(3)中,A(t)表示所述纵波复信号的幅值。根据所述纵波复信号的实部和幅值,可以采用下述公式确定所述纵波属性数据:
cosφ(t)=s(t)/A(t) (4)
公式中cosφ(t)表示所述纵波属性数据。
对所述转换波地震数据的各项操作与上述步骤相同,在此不再一一赘述。
图4是本申请实施例中纵波属性数据和第一匹配的转换波属性数据的地震剖面示意图。图4中(a)图和(b)图分别为纵波属性数据和第一匹配的转换波属性数据的地震剖面示意图。从图4中可以看出,在所述纵波属性数据和所述转换波属性数据中所述小层的同相轴可以清楚地识别,表明通过Hilbert变换所得到的纵波属性数据和转换波属性数据可以增强所述目标工区中标志层之间的小层的反射特征,有利于对所述小层进行识别。
步骤S104:基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
具体地,针对所述至少一个小层中的一个小层,对所述纵波属性数据和所述转换波属性数据进行时间匹配处理,可以确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量。根据所述时移量,可以计算所述纵波地震数据和第一匹配的转换波地震数据在所述小层上的速度比。利用所述小层上的速度比对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。其中,所述时间匹配处理具体可以为动态时间调整。进一步地,第二压缩处理可以包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述第一匹配的转换波地震数据进行压缩。具体地,压缩后的第一匹配的转换波地震数据的时间等于所述纵波地震数据的时间。
进一步地,所述对所述纵波属性数据和所述转换波属性数据进行动态时间调整可以包括:可以计算所述目标纵波属性数据和所述转换波属性数据之间的对齐误差;根据所述对齐误差可以计算得到累积距离;基于所述累积距离,可以确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量。例如,可以采用下述公式计算所述目标纵波属性数据和所述转换波属性数据之间的对齐误差:
e[i,l]=(f[i]-g[i+l])2 (5)
公式(5)中,e[i,l]表示所述对齐误差,f[i]表示所述纵波属性数据中第i个采样点对应的采样时间,g[i+l]表示第一匹配的转换波地震数据中第i+l采样点对应的采样时间,l表示采样点的时延量。
所述根据对齐误差计算得到累积距离,可以利用下述公式对所述对齐误差e[i,l]进行迭代处理来实现,得到累积距离:
d[0,l]=e[0,l]
公式(6)中,d[i,l]表示所述累积距离,i=1,2,...,N-1,N表示所述迭代处理过程中的采样点总个数。
可以采用下述公式在所述累积距离中进行反向追踪,确定最小路径,即时移量序列u[0:N-1]:
公式(7)中,u[N-1]为追踪的第一个时移量,u[0]为追踪的最后一个时移量,j=1,2,...,N,N表示所述迭代处理过程中的采样点总个数。
图5是本申请实施例中目标转换波地震数据的地震剖面示意图。如图5所示,所述目的工区中各标志层以及各标志层之间的小层上的目标转换波地震数据与纵波地震数据的时间匹配的精度较高。
所述纵波和转换波地震数据时间匹配的方法实施例,分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特(Hilbert)变换,得到变换后的纵波属性数据和转换波属性数据,这些数据可以突出小层的反射特征,使小层对应弱同相轴信号变强,提高了小层的对比精度,从而可以提高对目标工区中小层上的纵波和转换波地震数据时间匹配的精度。
图6是本申请纵波和转换波地震数据时间匹配的装置实施例的组成结构图。所述纵波和转换波地震数据时间匹配的装置可以包括:标志层速度比计算模块100、第一压缩处理模块200、属性数据确定模块300和第二压缩处理模块400。
所述标志层速度比计算模块100,可以用于计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比。
所述第一压缩处理模块200,可以用于利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据。
所述属性数据确定模块300,可以用于根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据。
所述第二压缩处理模块400,可以用于基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
图7是本申请纵波和转换波地震数据时间匹配的装置实施例中属性数据确定模块的组成结构图。如图7所示,图6中的属性数据确定模块300,可以包括:复信号虚部确定模块310、复信号实部确定模块320、幅值确定模块330和波属性数据确定模块340。
所述复信号虚部确定模块310,可以用于分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部。
所述复信号实部确定模块320,可以用于将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部。
所述幅值确定模块330,可以用于根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值。
所述波属性数据确定模块340,可以用于根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
所述纵波和转换波地震数据时间匹配的装置实施例与所述纵波和转换波地震数据时间匹配的方法实施例相对应,可以突出小层的反射特征,使小层对应弱同相轴信号变强,提高了小层的对比精度,从而可以提高对目标工区中小层上的纵波和转换波地震数据时间匹配的精度。
在20世纪90年代,对于一个技术的改进可以很明显地区分是硬件上的改进(例如,对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而,随着技术的发展,当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此,不能说一个方法流程的改进就不能用硬件实体模块来实现。例如,可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray,FPGA))就是这样一种集成电路,其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上,而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且,如今,取代手工地制作集成电路芯片,这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现,它与程序开发撰写时所用的软件编译器相类似,而要编译之前的原始代码也得用特定的编程语言来撰写,此称之为硬件描述语言(Hardware Description Language,HDL),而HDL也并非仅有一种,而是有许多种,如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等,目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚,只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中,就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中,内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括短暂电脑可读媒体(transitory media),如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (11)
1.一种纵波和转换波地震数据时间匹配的方法,其特征在于,包括:
计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比;
利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据;
根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据;
基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
2.根据权利要求1所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述根据纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据,包括:
分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部;
将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部;
根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值;
根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
3.根据权利要求1所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述目标工区包括:至少两个标志层和各标志层之间的小层。
4.根据权利要求3所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据,包括:
针对所述至少一个小层中的一个小层,对所述纵波属性数据和所述转换波属性数据进行时间匹配处理,确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量;
根据所述时移量,计算所述纵波地震数据和第一匹配的转换波地震数据在所述小层上的速度比;
利用所述小层上的速度比对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
5.根据权利要求4所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述时间匹配处理包括:动态时间调整。
6.根据权利要求5所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述对所述纵波属性数据和所述转换波属性数据进行动态时间调整包括:
计算所述目标纵波属性数据和所述转换波属性数据之间的对齐误差;
根据所述对齐误差计算累积距离;
基于所述累积距离,确定在所述小层上所述纵波属性数据和所述转换波属性数据之间的时移量。
7.根据权利要求1所述的纵波和转换波地震数据时间匹配的方法,其特征在于,所述计算目标工区中标志层上的纵波地震数据和转换波地震数据的速度比,包括:
在所述纵波地震数据和所述转换波地震数据上分别确定所述目标工区中标志层对应的纵波层位时间和转换波层位时间;
根据所述纵波层位时间和所述转换波层位时间,计算所述目标工区标志层上的纵波地震数据和转换波地震数据的速度比。
8.根据权利要求1所述的纵波和转换波地震数据时间匹配的方法,其特征在于,第一压缩处理包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述转换波地震数据进行压缩;具体地,所述压缩后的转换波地震数据的时间等于所述纵波地震数据的时间。
9.根据权利要求1所述的纵波和转换波地震数据时间匹配的方法,其特征在于,第二压缩处理包括:根据所述转换波地震数据的时间和所述纵波地震数据的时间,对所述第一匹配的转换波地震数据进行压缩;具体地,压缩后的第一匹配的转换波地震数据的时间等于所述纵波地震数据的时间。
10.一种纵波和转换波地震数据时间匹配的装置,其特征在于,所述装置包括:标志层速度比计算模块、第一压缩处理模块、属性数据确定模块和第二压缩处理模块;其中,
所述标志层速度比计算模块,用于确定目标工区中标志层上的纵波地震数据和转换波地震数据的速度比;
所述第一压缩处理模块,用于利用所述标志层上的速度比对所述转换波地震数据进行第一压缩处理,得到第一匹配的转换波地震数据;
所述属性数据确定模块,用于根据所述纵波地震数据和第一匹配的转换波地震数据,分别确定纵波属性数据和转换波属性数据;
所述第二压缩处理模块,用于基于所述纵波属性数据和所述转换波属性数据,对第一匹配的转换波地震数据进行第二压缩处理,得到目标转换波地震数据。
11.根据权利要求10所述的一种纵波和转换波地震数据时间匹配的装置,其特征在于,所述属性数据确定模块,包括:复信号虚部确定模块、复信号实部确定模块、幅值确定模块和波属性数据确定模块;其中,
所述复信号虚部确定模块,用于分别对所述纵波地震数据和第一匹配的转换波地震数据进行希尔伯特变换,分别得到所述纵波地震数据对应的纵波复信号的虚部和第一匹配的转换波地震数据对应的转换波复信号的虚部;
所述复信号实部确定模块,用于将所述纵波地震数据和第一匹配的转换波地震数据分别作为所述纵波复信号的实部和所述转换波复信号的实部;
所述幅值确定模块,用于根据所述纵波复信号的实部和虚部可以得到所述纵波复信号的幅值,以及根据所述转换波复信号的实部和虚部可以得到所述转换波复信号的幅值;
所述波属性数据确定模块,用于根据所述纵波复信号的幅值和实部确定所述纵波属性数据,以及根据所述转换波复信号的幅值和实部确定所述转换波属性数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611122196.7A CN106646611B (zh) | 2016-12-08 | 2016-12-08 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611122196.7A CN106646611B (zh) | 2016-12-08 | 2016-12-08 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646611A true CN106646611A (zh) | 2017-05-10 |
CN106646611B CN106646611B (zh) | 2018-11-16 |
Family
ID=58818782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611122196.7A Active CN106646611B (zh) | 2016-12-08 | 2016-12-08 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646611B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107167843A (zh) * | 2017-05-25 | 2017-09-15 | 中国石油天然气集团公司 | 多波时间域匹配方法及装置 |
CN107589452A (zh) * | 2017-07-21 | 2018-01-16 | 中国石油天然气集团公司 | 纵波和转换波的数据匹配方法和装置 |
CN108761533A (zh) * | 2018-05-17 | 2018-11-06 | 中国石油天然气集团有限公司 | 一种确定纵横波速度比的方法、装置及系统 |
CN108897039A (zh) * | 2018-08-27 | 2018-11-27 | 电子科技大学 | 基于层位控制的dtw多波匹配方法 |
CN110618448A (zh) * | 2018-06-19 | 2019-12-27 | 中国石油化工股份有限公司 | 多波动态匹配方法及系统 |
CN111665558A (zh) * | 2019-03-07 | 2020-09-15 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
CN112379439A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 地震数据中纵波与横波匹配的方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021184A1 (en) * | 2001-06-08 | 2003-01-30 | Yaohui Zhang | Extraction of P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion |
CN101359056A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气集团公司 | 一种生成纵波时间域高精度转换波剖面的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
CN104237938A (zh) * | 2014-09-04 | 2014-12-24 | 中国石油天然气集团公司 | 利用动态时间调整进行纵波和转换波数据匹配的方法 |
CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
CN105445785A (zh) * | 2014-06-12 | 2016-03-30 | 中国石油化工股份有限公司 | 一种横波偏移速度建模方法 |
-
2016
- 2016-12-08 CN CN201611122196.7A patent/CN106646611B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030021184A1 (en) * | 2001-06-08 | 2003-01-30 | Yaohui Zhang | Extraction of P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion |
CN101359056A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气集团公司 | 一种生成纵波时间域高精度转换波剖面的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
CN105445785A (zh) * | 2014-06-12 | 2016-03-30 | 中国石油化工股份有限公司 | 一种横波偏移速度建模方法 |
CN104237938A (zh) * | 2014-09-04 | 2014-12-24 | 中国石油天然气集团公司 | 利用动态时间调整进行纵波和转换波数据匹配的方法 |
Non-Patent Citations (1)
Title |
---|
刘兰锋 等: ""动态图像变形的纵横波自动时间匹配"", 《石油地球物理勘探》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107167843A (zh) * | 2017-05-25 | 2017-09-15 | 中国石油天然气集团公司 | 多波时间域匹配方法及装置 |
CN107167843B (zh) * | 2017-05-25 | 2019-05-07 | 中国石油天然气集团公司 | 多波时间域匹配方法及装置 |
CN107589452A (zh) * | 2017-07-21 | 2018-01-16 | 中国石油天然气集团公司 | 纵波和转换波的数据匹配方法和装置 |
CN107589452B (zh) * | 2017-07-21 | 2019-06-11 | 中国石油天然气集团公司 | 纵波和转换波的数据匹配方法和装置 |
CN108761533A (zh) * | 2018-05-17 | 2018-11-06 | 中国石油天然气集团有限公司 | 一种确定纵横波速度比的方法、装置及系统 |
CN108761533B (zh) * | 2018-05-17 | 2019-10-11 | 中国石油天然气集团有限公司 | 一种确定纵横波速度比的方法、装置及系统 |
CN110618448A (zh) * | 2018-06-19 | 2019-12-27 | 中国石油化工股份有限公司 | 多波动态匹配方法及系统 |
CN108897039A (zh) * | 2018-08-27 | 2018-11-27 | 电子科技大学 | 基于层位控制的dtw多波匹配方法 |
CN108897039B (zh) * | 2018-08-27 | 2019-11-12 | 电子科技大学 | 基于层位控制的dtw多波匹配方法 |
CN111665558A (zh) * | 2019-03-07 | 2020-09-15 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
CN111665558B (zh) * | 2019-03-07 | 2023-02-28 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
CN112379439A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 地震数据中纵波与横波匹配的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106646611B (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646611A (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
CN106054251B (zh) | 一种初至波拾取方法及装置 | |
CN106772589B (zh) | 一种叠前地震反演方法及装置 | |
CN107092031B (zh) | 一种声波测井数据校正方法及装置 | |
CN106353818B (zh) | 一种地震数据时差校正的方法及系统 | |
CN107193045A (zh) | 一种地震数据处理方法及装置 | |
CN106842304A (zh) | 一种叠前深度偏移方法及装置 | |
CN109061734A (zh) | 一种地震数据的处理方法及装置 | |
CN110794458B (zh) | 基于时频分析的含气性检测方法、装置及存储介质 | |
CN107179546B (zh) | 一种确定速度场的方法及装置 | |
CN105353406B (zh) | 一种生成角道集的方法和装置 | |
CN106291691A (zh) | 一种地震偏移成像方法及装置 | |
CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
CN105044773B (zh) | 转换波井震标定方法及装置 | |
CN109100797A (zh) | 一种确定致密储层含气性的方法及装置 | |
CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
Decker et al. | Quantifying and correcting residual azimuthal anisotropic moveout in image gathers using dynamic time warping | |
CN106707345A (zh) | 一种角度弹性参数岩性识别方法和装置 | |
CN106772596B (zh) | 一种确定叠前时间偏移速度场的方法及装置 | |
CN111665558B (zh) | 一种油气分布检测方法和系统 | |
CN107329169A (zh) | 一种角度道集的提取方法及装置 | |
CN106125137B (zh) | 一种转换波地震数据处理方法及装置 | |
CN109932748A (zh) | 一种地表一致性振幅补偿处理方法、装置及存储介质 | |
CN107290783A (zh) | 一种油气检测方法及装置 | |
US20230184982A1 (en) | Full-wavefield angle gather for high-contrast inter thin-bed models |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |