CN106226811A - 一种基于属性的纵波数据和转换波数据匹配的方法及装置 - Google Patents

一种基于属性的纵波数据和转换波数据匹配的方法及装置 Download PDF

Info

Publication number
CN106226811A
CN106226811A CN201610524422.8A CN201610524422A CN106226811A CN 106226811 A CN106226811 A CN 106226811A CN 201610524422 A CN201610524422 A CN 201610524422A CN 106226811 A CN106226811 A CN 106226811A
Authority
CN
China
Prior art keywords
wave
converted
data
attribute
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.)
Pending
Application number
CN201610524422.8A
Other languages
English (en)
Inventor
张铁强
孙鹏远
钱忠平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201610524422.8A priority Critical patent/CN106226811A/zh
Publication of CN106226811A publication Critical patent/CN106226811A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing 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

本申请涉及地震数据处理领域,尤其涉及一种基于属性的纵波数据和转换波数据匹配的方法及装置,该方法包括:获取目标区域的振幅随偏移距变化(AVO)的纵波道集和转换波道集;将所述纵波道集进行AVO反演,获取所述纵波的拟横波阻抗属性;将所述转换波道集进行叠加处理,获得转换波叠后数据;根据所述拟横波阻抗属性和所述转换波叠后的相关性获取满足预设条件的速度比;根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。本申请实施例的拟横波阻抗属性与转换波数据具有更高的相关性,从而可以提高纵波和转换波时间匹配的精度,为多波地震资料联合反演与解释提供准确的基础。

Description

一种基于属性的纵波数据和转换波数据匹配的方法及装置
技术领域
本发明涉及地震数据处理领域,尤其涉及一种基于属性的纵波数据和转换波数据匹配的方法及装置。
背景技术
在多分量地震勘探中,反射纵波(PP波)是由下行P波和上行P波组成,反射转换波(PS波)是由下行P波和上行S波组成。由于纵波(P波)和横波(S波)的传播机理不同,造成接收到的纵波数据和转换波数据在旅行时和波形等方面都存在差异,其主要表现是:对于相同的反射界面,纵波(PP波)和转换波(PS波)在反射时间上存在差异,并且在各自的地震反射剖面上的同相轴个数不相同;与纵波(PP波)相比,转换波(PS波)能量弱、主频低。
纵波(PP波)和转换波(PS波)联合分析可以提高构造成像、储层描述等的精度和可信度。在联合分析时,需要将来自地下相同反射界面的纵波(PP波)和转换波(PS波)时间进行匹配。由于纵波(PP波)和转换波(PS波)的动力学特征差异较大(振幅、频率、相位等),使得直接利用纵波(PP波)和转换波(PS波)数据的波形相似性进行匹配时精度低,效果较差。
发明内容
本申请实施例提供一种基于属性的纵波数据和转换波数据匹配的方法及装置,以提高纵波数据和转换波数据时间匹配的精度,为多波地震资料联合反演与解释提供准确的基础。
为达到上述目的,一方面,本申请实施例提供了一种基于属性的纵波数据和转换波数据匹配的方法,该方法包括:
获取目标区域的振幅随偏移距变化(AVO)的纵波道集和转换波道集;
对所述转换波道集进行叠加处理,获得转换波叠后数据;
将所述纵波道集进行AVO反演,获取拟横波阻抗属性;
根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比;
根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
进一步地,所述方法还包括:
根据所述速度比将所述转换波道集从转换波时间域转换为纵波时间域,获得匹配后的转换波道集。
进一步地,根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比,包括:
对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,获取波组特性一致性最好的层位;
拾取所述波组特性一致性最好的层位相对应的纵波层位时间和转换波层位时间;
根据所述纵波层位时间和转换波层位时间计算速度比。
进一步地,根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比,包括:
采用一系列的速度比计算所述拟横波阻抗属性和所述转换波叠后数据之间的相关系数;
比较所述相关系数,获取所述相关系数的最大值对应的速度比,所述相关系数的最大值对应的速度比为满足预设条件的速度比。
进一步地,采用以下公式对所述纵波道集进行AVO反演:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - [ 1 2 tan 2 θ - 2 γ 2 sin 2 θ ] Δ ρ ρ ‾
式中,RPP为纵波反射系数;θ为入射角;ΔIP为纵波阻抗差;为纵波阻抗的平均值;ΔIS为横波阻抗差;为横波阻抗的平均值;γ为速度比,表达式为VS为横波速度,VP为纵波速度。
进一步地,采用以下公式对所述转换波叠后数据从转换波时间域转换为纵波时间域:
TPP=2*TPS/(1+γ)
其中,TPP为纵波时间,TPS为转换波时间,γ为速度比。
另一方面,本申请实施例还提供了一种基于属性的纵波数据和转换波数据匹配的装置,该装置包括:
道集获取单元,用于获取目标区域的振幅随偏移距变化AVO的纵波道集和转换波道集;
叠加处理单元,用于对所述转换波道集进行叠加处理,获得转换波叠后数据;
纵波反演单元,用于将所述纵波道集进行AVO反演,获取所述纵波的拟横波阻抗属性;
速度比获取单元,用于根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比;
匹配单元,用于根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
进一步地,所述匹配单元用于根据所述速度比将所述转换波道集从转换波时间域转换为纵波时间域,获得匹配后的转换波道集。
进一步地,所述速度比获取单元包括:
层位波组特性比较子单元,用于对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,获取波组特性一致性最好的层位;
层位时间拾取子单元,用于拾取所述波组特性一致性最好的层位相对应的纵波层位时间和转换波层位时间;
速度比计算子单元,用于根据所述纵波层位时间和转换波层位时间计算速度比。
进一步地,所述速度比获取单元包括:
相关系数计算子单元,用于采用一系列的速度比计算所述拟横波阻抗属性和所述转换波叠后数据之间的相关系数;
比较获取子单元,用于比较所述相关系数,获取所述相关系数的最大值对应的速度比,所述相关系数的最大值对应的速度比为满足预设条件的速度比。
进一步地,采用以下公式对所述纵波道集进行AVO反演:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - [ 1 2 tan 2 θ - 2 γ 2 sin 2 θ ] Δ ρ ρ ‾
式中,RPP为纵波反射系数;θ为入射角;ΔIP为纵波阻抗差;为纵波阻抗的平均值;ΔIS为横波阻抗差;为横波阻抗的平均值;γ为速度比,表达式为VS为横波速度,VP为纵波速度。
进一步地,采用以下公式将所述转换波叠后数据从转换波时间域转换为纵波时间域:
TPP=2*TPS/(1+γ)
其中,TPP为纵波时间,TPS为转换波时间,γ为速度比。
本申请实施例通过将纵波道集进行反演获得拟横波阻抗属性,再利用拟横波阻抗属性和转换波叠后数据之间的相关性获得速度比,最后根据速度比将转换波叠后数据从转换波时间域转换为纵波时间域,从而获得匹配后的转换波数据。本申请实施例的拟横波阻抗属性与转换波叠后数据具有更高的相关性,从而用这种属性和转换波叠后数据进行时间匹配,可以提高纵波和转换波时间匹配的精度,为多波地震资料联合反映与解释提供准确的基础。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例的纵波和转换波的传播路径示意图;
图2是本申请实施例的基于属性的纵波数据和转换波数据匹配的方法示意图;
图3是本申请实施例中的基于属性的纵波数据和转换波数据匹配的装置结构示意图;
图4(a)~4(b)是本申请一实施例的纵波AVO道集和转换波AVO道集示意图;
图5(a)~5(e)是纵波和转换波匹配前后数据对比示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面结合附图,对本申请实施例的具体实施方式作进一步的详细说明。
如图1所示为纵波和转换波的传播路径示意图,图中O为震源点,R’为转换波接收点,R为纵波接收点。从图中可以看出,纵波由于入射为纵波,反射也为纵波,其传播速度是一样的。而对于转换波地震资料,入射为纵波,反射为横波,其传播速度和路径都不一样。因此,反射纵波地震资料与转换波地震资料对于同一反射界面的同向轴在时间域内不同。因此,在进行纵波与转换波联合分析时,需要将来自地下相同反射界面的纵波(PP波)和转换波(PS波)时间进行匹配。本申请实施例提出了一种基于属性驱动的方法,从而提高纵波与转换波匹配效果。属性驱动方法的关键是找到纵波和转换波之间最大相关性的属性,相关程度是提高匹配有效性的基础。
参考图2,本申请实施例提供了一种基于属性的纵波数据和转换波数据匹配的方法,所述方法包括:
S1、获取目标区域的振幅随偏移距变化(AVO)的纵波道集和转换波道集。
本申请实施例中,首先在野外利用纵波震源激发地震波并利用检波器记录地震波,然后按照常规地震资料处理流程对采集的数据进行相对振幅保持的高保真处理,从而获取正常时差校正(NMO)后的振幅随偏移距变化(AVO)的纵波道集和转换波道集。具体的,正常时窗时差校正是对地震资料的一个处理步骤,进行正常时差校正后,一个道集内相同时间的数据都来自地下相同的深度,以便进行后续的反演处理。
S2、对所述转换波道集进行叠加处理,获得转换波叠后数据。
本申请实施例中,对转换波道集进行叠加处理,从而可以获得转换波叠后数据。同时,对纵波道集进行叠加处理,从而可以获得纵波叠后数据。
S3、将所述纵波道集进行AVO反演,获取拟横波阻抗属性。
本申请实施例中,可以采用以下反射系数近似公式对所述纵波道集进行AVO反演:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - [ 1 2 tan 2 θ - 2 γ 2 sin 2 θ ] Δ ρ ρ ‾ - - - ( 1 )
式中,RPP为纵波反射系数;θ为入射角;ΔIP为纵波阻抗差;为纵波阻抗的平均值;ΔIS为横波阻抗差;为横波阻抗的平均值;γ表示速度比,即界面两侧地层的平均横波速度与平均纵波速度的比值表达式为VS为横波速度,VP为纵波速度。
公式(1)中的第三项可以忽略,从而获得以下公式:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - - - ( 2 )
本申请实施例中,通过公式(2),从纵波道集中可以直接反演获得拟横波阻抗属性,即本申请实施例中获取拟横波阻抗属性的方法在反演过程中不需要速度比,而且反射系数近似公式的精度高。
本申请实施例中,转换波反射系数具有如下近似关系:
R P S ( θ ) ≈ 2 γ ΔI S I ‾ S sin θ - - - ( 3 )
R P S ( θ ) ≈ 2 γ ΔI S I ‾ S sinθcos 2 θ - - - ( 4 )
其中,RPS为转换波反射系数,θ,ΔIS和γ的意义同上。
由(3)式和(4)式可以看出,反演得到的拟横波阻抗属性与转换波反射系数之间只存在比例关系,因而与转换波数据的相似性更高,进而有利于层位对比,可以提高层位对比精度。
S4、根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比。
本申请中速度比表示界面两侧地层的平均横波速度与平均纵波速度的比值,本申请实施例中横波阻抗属性与转换波叠后数据的相似性高,通过两者的相关性可以获取满足预设条件的速度比,从而可以再根据满足预设条件的速度比对转换波叠后数据进行转换处理,从而可以得到匹配后的转换波数据。
在本申请实施例中,根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比,包括:
对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,获取波组特性一致性最好的层位;
拾取所述波组特性一致性最好的层位相对应的纵波层位时间和转换波层位时间;
根据所述纵波层位时间和转换波层位时间计算速度比。
本申请实施例中,利用拟横波阻抗属性和转换波叠后数据的相似性,对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,在两个数据集上寻找波组特性一致性最好的层位,分别进行层位时间拾取,取得相对应的纵波和转换波层位时间,再通过层位时间和速度比的转换公式,获取速度比。本申请实施例中,层位时间和速度比的转换公式如下:
TPP=2*TPS/(1+γ) (5)
其中,TPP为纵波时间,TPS为转换波时间,γ为速度比。
在本申请另一实施例中,所述根据拟横波阻抗属性和所述转换波叠后数据获取速度比,包括:
采用一系列的速度比计算所述拟横波阻抗属性和所述转换波叠后数据之间的相关系数;
比较所述相关系数,获取所述相关系数的最大值对应的速度比,所述相关系数的最大值对应的速度比为满足预设条件的速度比。
本申请实施例中,采用一系列的平均横波速度与平均纵波速度的比值计算拟横波阻抗属性和转换波叠后数据之间的相关系数,寻找相关系数最大值对应的速度比。
S5、根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
本申请实施例中,根据步骤S4获取的速度比对转换波叠后数据进行压缩,将所述转换波叠后数据从转换波时间域转换为纵波时间域,得到纵波时间域的转换波数据,即在纵波时间匹配后的转换波数据,从而可以进行后续的纵波数据和转换波数据联合分析。
在本申请另一实施例中,根据对数据处理的需要,可以根据步骤S4获取的速度比对转换波道集从转换波时间域转换为纵波时间域,从而可以获得在纵波时间匹配后的转换波道集,可以进行后续的纵波道集和转换波道集联合分析。
本申请实施例中,原始的转换波数据是在转换波时间域,公式(5)同时也表达了转换波时间和纵波时间的关系。在获得速度比γ后,可以通过公式(5)将转换波数据进行压缩,即将转换波数据中的转换波时间根据公式(5)转换为纵波时间,从而得到纵波时间域的转换波数据。
本申请实施例通过将纵波道集进行反演获得拟横波阻抗属性,再利用拟横波阻抗属性和转换波叠后数据的相关性获得速度比,最后根据速度比对转换波叠后数据进行压缩,从而获得在纵波时间匹配后的转换波数据。本申请实施例的拟横波阻抗属性与转换波叠后数据具有更高的相关性,从而用这种属性和转换波数据进行层位对比和时间匹配,可以提高纵波和转换波时间匹配的精度,为多波地震资料联合反映与解释提供准确的基础。
为了清楚的说明本申请实施例的有益效果,下面结合附图进行说明:
如图4(a)、4(b)分别为纵波AVO道集和转换波AVO道集,从图中可以看出,反射纵波地震资料与转换波地震资料对于同一反射界面的同向轴在时间域内不同。图5(a)为转换波叠后数据,图5(b)为拟横波阻抗属性,图5(c)为速度比,图5(d)为纵波时间域的转换波数据,即匹配后的转换波数据,图5(e)为纵波叠后数据。经过对数据的相关性分析,纵波叠后数据与匹配后的转换波数据之间的相关系数为0.742,拟横波阻抗属性与匹配后的转换波数据之间的相关系数为0.996。
如图3所示,本申请实施例还提供了一种基于属性的纵波数据和转换波数据匹配的装置,所述装置包括:
道集获取单元21,用于获取目标区域的振幅随偏移距变化AVO的纵波道集和转换波道集;
叠加处理单元22,用于对所述转换波道集进行叠加处理,获得转换波叠后数据;
纵波反演单元23,用于将所述纵波道集进行AVO反演,获取所述纵波的拟横波阻抗属性;
速度比获取单元24,用于根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比;
匹配单元25,用于根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
本实施例的装置的各组成部分分别用于实现前述实施例的方法的各步骤,由于在方法实施例中,已经对各步骤进行了详细说明,在此不再赘述。
本申请实施例通过将纵波道集进行反演获得拟横波阻抗属性,再利用拟横波阻抗属性和转换波叠后数据的相关性获得速度比,最后根据速度比对转换波叠后数据进行压缩,从而获得在纵波时间匹配后的转换波数据。本申请实施例的拟横波阻抗属性与转换波数据具有更高的相关性,从而用这种属性和转换波数据进行层位对比和时间匹配,可以提高纵波和转换波时间匹配的精度,为多波地震资料联合反映与解释提供准确的基础。
在一个或多个示例性的设计中,本申请实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。
以上所述的具体实施例,对本申请的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本申请实施例的具体实施例而已,并不用于限定本申请的保护范围,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (12)

1.一种基于属性的纵波数据和转换波数据匹配的方法,其特征在于,所述方法包括:
获取目标区域的振幅随偏移距变化AVO的纵波道集和转换波道集;
对所述转换波道集进行叠加处理,获得转换波叠后数据;
将所述纵波道集进行AVO反演,获取拟横波阻抗属性;
根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比;
根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
2.如权利要求1所述的基于属性的纵波数据和转换波数据匹配的方法,其特征在于,所述方法还包括:
根据所述速度比将所述转换波道集从转换波时间域转换为纵波时间域,获得匹配后的转换波道集。
3.如权利要求1所述的基于属性的纵波数据和转换波数据匹配的方法,其特征在于,根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比,包括:
对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,获取波组特性一致性最好的层位;
拾取所述波组特性一致性最好的层位相对应的纵波层位时间和转换波层位时间;
根据所述纵波层位时间和转换波层位时间计算速度比。
4.如权利要求1所述的基于属性的纵波数据和转换波数据匹配的方法,其特征在于,根据所述拟横波阻抗属性和所述转换波叠后数据的相关性获取满足预设条件的速度比,包括:
采用一系列的速度比计算所述拟横波阻抗属性和所述转换波叠后数据之间的相关系数;
比较所述相关系数,获取所述相关系数的最大值对应的速度比,所述相关系数的最大值对应的速度比为满足预设条件的速度比。
5.如权利要求1所述的基于属性的纵波数据和转换波数据匹配的方法,其特征在于,采用以下公式对所述纵波道集进行AVO反演:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - [ 1 2 tan 2 θ - 2 γ 2 sin 2 θ ] Δ ρ ρ ‾
式中,RPP为纵波反射系数;θ为入射角;ΔIP为纵波阻抗差;为纵波阻抗的平均值;ΔIS为横波阻抗差;为横波阻抗的平均值;γ为速度比,表达式为VS为横波速度,VP为纵波速度。
6.如权利要求1所述的基于属性的纵波数据和转换波数据匹配的方法,其特征在于,采用以下公式将所述转换波叠后数据从转换波时间域转换为纵波时间域:
TPP=2*TPS/(1+γ)
其中,TPP为纵波时间,TPS为转换波时间,γ为速度比。
7.一种基于属性的纵波数据和转换波数据匹配的装置,其特征在于,所述装置包括:
道集获取单元,用于获取目标区域的振幅随偏移距变化AVO的纵波道集和转换波道集;
叠加处理单元,用于对所述转换波道集进行叠加处理,获得转换波叠后数据;
纵波反演单元,用于将所述纵波道集进行AVO反演,获取所述纵波的拟横波阻抗属性;
速度比获取单元,用于根据所述拟横波阻抗属性集和所述转换波叠后数据的相关性获取满足预设条件的速度比;
匹配单元,用于根据所述速度比将所述转换波叠后数据从转换波时间域转换为纵波时间域,获得匹配后的转换波数据。
8.如权利要求7所述的基于属性的纵波数据和转换波数据匹配的装置,其特征在于,所述匹配单元用于根据所述速度比将所述转换波道集从转换波时间域转换为纵波时间域,获得匹配后的转换波道集。
9.如权利要求7所述的基于属性的纵波数据和转换波数据匹配的装置,其特征在于,所述速度比获取单元包括:
层位波组特性比较子单元,用于对所述拟横波阻抗属性和所述转换波叠后数据不同层位的波组特性一致性进行比较,获取波组特性一致性最好的层位;
层位时间拾取子单元,用于拾取所述波组特性一致性最好的层位相对应的纵波层位时间和转换波层位时间;
速度比计算子单元,用于根据所述纵波层位时间和转换波层位时间计算速度比。
10.如权利要求7所述的基于属性的纵波数据和转换波数据匹配的装置,其特征在于,所述速度比获取单元包括:
相关系数计算子单元,用于采用一系列的速度比计算所述拟横波阻抗属性和所述转换波叠后数据之间的相关系数;
比较获取子单元,用于比较所述相关系数,获取所述相关系数的最大值对应的速度比,所述相关系数的最大值对应的速度比为满足预设条件的速度比。
11.如权利要求7所述的基于属性的纵波数据和转换波数据匹配的装置,其特征在于,采用以下公式对所述纵波道集进行AVO反演:
R P P ( θ ) = 1 2 [ 1 + tan 2 θ ] ΔI P I ‾ P - 4 γ 2 sin 2 θ ΔI S I ‾ S - [ 1 2 tan 2 θ - 2 γ 2 sin 2 θ ] Δ ρ ρ ‾
式中,RPP为纵波反射系数;θ为入射角;ΔIP为纵波阻抗差;为纵波阻抗的平均值;ΔIS为横波阻抗差;为横波阻抗的平均值;γ为速度比,表达式为VS为横波速度,VP为纵波速度。
12.如权利要求7所述的基于属性的纵波数据和转换波数据匹配的装置,其特征在于,采用以下公式将所述转换波叠后数据从转换波时间域转换为纵波时间域:
TPP=2*TPS/(1+γ)
其中,TPP为纵波时间,TPS为转换波时间,γ为速度比。
CN201610524422.8A 2016-07-05 2016-07-05 一种基于属性的纵波数据和转换波数据匹配的方法及装置 Pending CN106226811A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610524422.8A CN106226811A (zh) 2016-07-05 2016-07-05 一种基于属性的纵波数据和转换波数据匹配的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610524422.8A CN106226811A (zh) 2016-07-05 2016-07-05 一种基于属性的纵波数据和转换波数据匹配的方法及装置

Publications (1)

Publication Number Publication Date
CN106226811A true CN106226811A (zh) 2016-12-14

Family

ID=57519906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610524422.8A Pending CN106226811A (zh) 2016-07-05 2016-07-05 一种基于属性的纵波数据和转换波数据匹配的方法及装置

Country Status (1)

Country Link
CN (1) CN106226811A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107589452A (zh) * 2017-07-21 2018-01-16 中国石油天然气集团公司 纵波和转换波的数据匹配方法和装置
CN108508481A (zh) * 2018-05-17 2018-09-07 中国石油天然气集团有限公司 一种纵波转换波地震数据时间匹配的方法、装置及系统
CN110618448A (zh) * 2018-06-19 2019-12-27 中国石油化工股份有限公司 多波动态匹配方法及系统
CN110850504A (zh) * 2019-11-29 2020-02-28 中国地质调查局油气资源调查中心 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法
US11402530B2 (en) 2018-09-30 2022-08-02 Petrochina Company Limited Method for acquiring converted wave, electronic device and readable storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070081422A1 (en) * 2005-10-11 2007-04-12 Keshan Zou PP/PS event matching (registration)
CN102692645A (zh) * 2012-06-01 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用纵波、转换波数据联合反演储层纵横波速度比的方法
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN104570110A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种基于纵横波匹配的多分量资料联合速度分析方法
CN104570079A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种纵波、转换横波地震资料的时间匹配方法
US20150338536A1 (en) * 2014-05-20 2015-11-26 Cgg Services Sa P-wave and s-wave separation of seismic data in the presence of statics and irregular geometry
WO2016027156A1 (en) * 2014-08-19 2016-02-25 Cgg Services Sa Joint inversion of compressional and shear seismic data in native time domains

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070081422A1 (en) * 2005-10-11 2007-04-12 Keshan Zou PP/PS event matching (registration)
CN102692645A (zh) * 2012-06-01 2012-09-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用纵波、转换波数据联合反演储层纵横波速度比的方法
CN103439739A (zh) * 2013-04-08 2013-12-11 中国石油集团东方地球物理勘探有限责任公司 地球物理勘探用纵横波匹配方法及匹配装置
CN104570110A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种基于纵横波匹配的多分量资料联合速度分析方法
CN104570079A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 一种纵波、转换横波地震资料的时间匹配方法
US20150338536A1 (en) * 2014-05-20 2015-11-26 Cgg Services Sa P-wave and s-wave separation of seismic data in the presence of statics and irregular geometry
WO2016027156A1 (en) * 2014-08-19 2016-02-25 Cgg Services Sa Joint inversion of compressional and shear seismic data in native time domains

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐天吉: "纵横波高精度匹配方法研究与应用", 《中国博士学位论文全文数据库•基础科学辑》 *
陈双全 等: "基于叠前反演方法的纵波与转换波时间匹配研究", 《中国地球物理2011》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107589452A (zh) * 2017-07-21 2018-01-16 中国石油天然气集团公司 纵波和转换波的数据匹配方法和装置
CN107589452B (zh) * 2017-07-21 2019-06-11 中国石油天然气集团公司 纵波和转换波的数据匹配方法和装置
CN108508481A (zh) * 2018-05-17 2018-09-07 中国石油天然气集团有限公司 一种纵波转换波地震数据时间匹配的方法、装置及系统
CN108508481B (zh) * 2018-05-17 2019-11-08 中国石油天然气集团有限公司 一种纵波转换波地震数据时间匹配的方法、装置及系统
CN110618448A (zh) * 2018-06-19 2019-12-27 中国石油化工股份有限公司 多波动态匹配方法及系统
US11402530B2 (en) 2018-09-30 2022-08-02 Petrochina Company Limited Method for acquiring converted wave, electronic device and readable storage medium
CN110850504A (zh) * 2019-11-29 2020-02-28 中国地质调查局油气资源调查中心 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法

Similar Documents

Publication Publication Date Title
CN106226811A (zh) 一种基于属性的纵波数据和转换波数据匹配的方法及装置
CN103149586B (zh) 一种倾斜层状粘弹性介质中波场正演模拟方法
CN104730579B (zh) 一种基于表层横波速度反演的纵横波联合静校正方法
CN106855640B (zh) 一种海底电缆地震数据初至计算方法及装置
CN103645497B (zh) 一种基于出射角的多分量波场分离方法
CN106970417B (zh) 椭圆展开转换波速度分析方法与系统
CN105093281A (zh) 一种反演框架下的地震多波建模方法
CN103576197B (zh) 一种转换波角道集抽取方法
CN106154319A (zh) 一种成像道集的分选方法
CN105629299A (zh) 角度域叠前深度偏移的走时、角度表获取方法及成像方法
CN102401908A (zh) 一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法
Yan et al. Angle gather extraction for acoustic and isotropic elastic RTM
CN106199704A (zh) 一种三维三分量海底电缆地震资料速度建模方法
CN102798888B (zh) 一种利用非零井源距数据计算纵横波速度比的方法
CN104182651B (zh) 用于三分量检波器接收的微地震事件方位角自动质控方法
CN106353798A (zh) 多分量联合高斯束叠前逆时偏移成像方法
CN106338766A (zh) 基于分步傅里叶算法的叠前时间偏移方法
CN103135134B (zh) 三维地震弹性偏移检波波场中标量横波的确定方法及装置
CN107589452A (zh) 纵波和转换波的数据匹配方法和装置
CN106990434B (zh) 椭圆展开转换波成像方法及系统
CN105510958A (zh) 一种适用复杂介质的三维vsp观测系统设计方法
CN103217707B (zh) 一种直接抽取纵波时间域转换波角道集的方法
CN101963674A (zh) 共转换点抽道集方法
US11686870B2 (en) Interpretive-guided velocity modeling seismic imaging method and system, medium and device
CN103472482A (zh) 基于基因排序体系的多波地震资料时间域匹配方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20161214