CN107589452B - 纵波和转换波的数据匹配方法和装置 - Google Patents
纵波和转换波的数据匹配方法和装置 Download PDFInfo
- Publication number
- CN107589452B CN107589452B CN201710598851.4A CN201710598851A CN107589452B CN 107589452 B CN107589452 B CN 107589452B CN 201710598851 A CN201710598851 A CN 201710598851A CN 107589452 B CN107589452 B CN 107589452B
- Authority
- CN
- China
- Prior art keywords
- wave
- converted
- time
- indicate
- converted wave
- 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
技术领域
本发明涉及地质勘探技术领域,特别涉及一种纵波和转换波的数据匹配方法和装置。
背景技术
纵波(PP波)和转换波(PS波)的联合分析处理可以提高储层构造成像、岩性分析、各向异性分析、流体预测以及储层描述的精度。在PP波和PS波的联合分析处理中,进行PP波和PS波的同相轴时间匹配是至关重要的。
现有技术中可以直接利用PP波和PS波中波形资料的波形相似性,或者基于属性驱动的方法进行PP波和PS波的数据匹配。这两种方法都可以是利用互相关的方法进行相似性运算,从而完成数据匹配的。然而,互相关方法本身会受到很多条件的限制,诸如:相关时窗、门槛值、时移量、信噪比、带限有限长度信号、分辨率等条件,因此,在实际应用中,采用上述两种方式进行时间匹配时精度不高,进一步地,这极大地影响着后续继续对储层进行联合反演、联合解释等分析工作。
发明内容
本发明提供了一种纵波和转换波的数据匹配方法和装置,以解决现有技术中采用互相关方法对纵波和转换波进行数据匹配时精度较低的问题。
本发明实施例提供了一种纵波和转换波的数据匹配方法,提供有:匹配数据集,所述匹配数据集可以包括但不限于以下至少之一:纵波和转换波之间的时移量、与所述纵波相对应的纵波时间、与所述转换波相对应的转换波时间,所述方法可以包括:根据所述时移量、所述纵波时间计算所述纵波和转换波之间的平均速度比;利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
在一个实施例中,在得到匹配后的纵波时间域的第一转换波数据之后,所述方法还可以包括:采用插值的方式对所述第一转换波匹配时间进行重采样,得到重采样处理后的第二转换波匹配时间。
在一个实施例中,可以按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
在一个实施例中,所述匹配数据集还可以包括:转换波属性,所述转换波属性是对转换波道集进行叠前反演得到的,所述转换波道集是对获取的地震资料进行正常时差校正后得到的。
在一个实施例中,根据所述时移量、所述纵波时间计算所述纵波属性和转换波属性之间的平均速度比,可以包括:
确定所述转换波属性是否进行过常量压缩处理;
当所述转换波属性进行过常量压缩处理时,按照以下公式计算所述平均速度比:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间,f表示第三参量,g表示第四参量,h表示第五参量;
当所述转换波属性没有进行过常量压缩处理时,按照以下公式计算所述平均速度比:
上式中,k表示第六参量。
在一个实施例中,所述压缩系数可以为大于1的常量。
本发明实施例还提供了一种纵波和转换波的数据匹配装置,提供有:匹配数据集,所述匹配数据集包括以下至少之一:纵波和转换波之间的时移量、与所述纵波相对应的纵波时间、与所述转换波相对应的转换波时间,所述装置可以包括:速度比计算模块,用于根据所述时移量、所述纵波时间计算所述纵波和转换波之间的平均速度比;数据转换模块,用于利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
在一个实施例中,所述装置还可以包括:重采样模块,用于在得到匹配后的纵波时间域的第一转换波数据之后,采用插值的方式对所述第一转换波数据进行重采样,得到重采样处理后的第二转换波数据。
在一个实施例中,所述数据转换模块具体可以用于按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
在一个实施例中,所述匹配数据集还可以包括:转换波属性,所述转换波属性是对转换波道集进行叠前反演得到的,所述转换波道集是对获取的地震资料进行正常时差校正后得到的。
在本发明实施例中,基于所述纵波和转换波之间的时移量计算得到所述平均速度比,利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。将所述时移量转换为所述平均速度比之后,再利用所述平均速度比对所述纵波时间以及所述转换波时间进行数据匹配,解决了现有技术中将所述时移量直接应用于转换波波形数据时,由于时移量变化比较剧烈,导致波形产生不同程度畸变的缺陷,提高了纵波和转换波的数据匹配精度,进一步地,对所得到的第一转换波匹配时间进行重采样,保证了重采样处理后的第二转换波匹配时间的高保真性,提高了联合利用所述纵波时间、所述第一转换波匹配时间进行储层反演、层位解释等分析工作的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的一种纵波和转换波的数据匹配方法流程图;
图2是本申请提供的纵波道集示意图;
图3是本申请提供的转换波道集示意图;
图4是本申请提供的纵波和转换波之间的时移量示意图;
图5是本申请提供的平均速度比示意图;
图6是本申请提供的第一转换波匹配时间示意图;
图7是本申请提供的纵波和转换波的数据匹配装置的一种结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
纵波可以指的是入射波为纵波、反射波也为纵波的波形,简称为PP波;转换波可以指的是入射波为纵波,反射波为横波的波形,简称为PS波。由于纵波和横波在传播和反射中的机理不同,因而,地下反射界面在纵波剖面和转换波剖面上的旅行时间以及波形都存在差异,这些差异主要表现在以下几个方面:1)纵波和转换波对于同一地层而言具有不同的反射特征,从而导致纵波和转换波在地震反射剖面上同相轴的个数不相同;2)横波的能量衰减很快,一般纵波和转换波能量相差几个数量级,并且转换波的相位畸变严重,同时,可能会存在相位反转的问题;3)转换波主频较低,并且转换波的频带一般只有纵波的一半,分辨率也比纵波低得多。考虑到现有技术中匹配纵波、转换波数据时,直接将时移量应用于波形数据会引起波形不同程度畸变的缺陷,本申请中,在得到所述时移量之后,先根据所述时移量计算得到平均速度比,再根据所述平均速度比对所述转换波时间进行压缩,得到匹配后的纵波时间域的第一转换波匹配时间。基于此,提出了一种纵波和转换波的数据匹配方法,如图1所示,可以包括以下步骤:
步骤101:根据所述时移量、所述纵波时间计算所述纵波和转换波之间的平均速度比。
在本实施例中,可以根据所述时移量、所述纵波时间,计算得到所述纵波和转换波之间的平均速度比。纵波可以指的是入射波为纵波、反射波也为纵波的波形,简称为PP波;转换波可以指的是入射波为纵波,反射波为横波的波形,简称为PS波。
在本实施例中,在计算所述平均速度比之前,可以先确定所述匹配数据集,其中,所述匹配数据集可以包括但不限于以下至少之一:纵波和转换波之间的时移量、与所述纵波相对应的纵波时间、与所述转换波相对应的转换波时间,以及所述转换波属性。
具体的,可以按照以下方式确定所述纵波和转换波之间的时移量以及所述转换波属性:
S1-1-1:获取振幅随偏移距变化的纵波道集和转换波道集。
可以将采集的地震资料进行正常时差校正等处理,并生成振幅随偏移距变化的纵波道集和转换波道集。在野外利用纵波震源激发地震波并利用检波器记录地震波,按照常规地震资料处理流程对采集的地震数据进行相对振幅保持的高保真处理,形成振幅随偏移距变化的正常时差校正后的纵波道集和转换波道集。
S1-1-2:分别对所述纵波道集和所述转换波道集进行叠前反演,得到所述纵波道集的多个待测纵波属性和所述转换波道集的多个待测转换波属性。
其中,所述待测纵波属性和所述待测转换波属性均可以包括截距、梯度、密度、纵波速度、横波速度、波阻抗、拉梅常数或弹性参数中的至少一种。所述反演方法可以是线性的或非线性的。可以利用Fatti公式对所述纵波道集和所述转换波道集进行叠前反演。
S1-1-3:计算所述多个待测纵波属性中第ii个待测纵波属性与所述多个待测转换波属性中第jj个待测转换波属性的互相关值,并将所得到的多个互相关值中的最大值对应的待测纵波属性和待测转换波属性,作为所述纵波属性和所述转换波属性。
其中,ii=1,2,…,n,jj=1,2,…,m,n表示待测纵波属性的个数,m表示待测转换波属性的个数。
S1-1-4:计算所述纵波和所述转换波之间的时移量。
基于动态时间调整的方法分别对所述纵波属性和所述转换波属性进行计算,得到所述纵波与所述转换波之间的时移量。具体的,基于动态时间调整的计算获得纵波与转换波之间的时移量的过程可以包括以下步骤:
步骤1:计算获得所述纵波属性与所述转换波属性之间的对齐误差;
步骤2:对所述对齐误差进行迭代计算,获得累积距离;
步骤3:将所述累积距离通过反向追踪获得的最小路径确定为时移量序列。
具体的,即为:首先,对于属性f[i]和g[i],首先计算f[i]和g[i]之间的对齐误差e[i,a]。
e[i,l]=(f[i]-g[i+a])2
上式中,a是整数的样点时延量。
然后,对对齐误差e[i,a]做迭代计算,得到累积距离d[i,a],计算公式如下:
d[0,a]=e[0,a]
上式中,i=1,2,…,N-1。
最后,在累积距离d[i,a]中反向追踪寻找最小路径,即时移量u[0:N-1]。开始追踪的第一个时延是u[N-1],最后一个时延是u[0]。
u[N-1]=argmind[N-1,a]
上式中,i=N-1,N-2,…,1。
在得到所述时移量之后,可以根据所述时移量计算所述纵波属性和转换波属性之间的平均速度比。具体的,根据所述时移量、所述纵波时间计算所述纵波属性和转换波属性之间的平均速度比,可以包括以下步骤:
S1-2-1:确定所述转换波属性是否进行过常量压缩处理;
S1-2-2:当所述转换波属性进行过常量压缩处理时,可以按照以下公式计算所述平均速度比:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间,f表示第三参量,g表示第四参量,h表示第五参量;
在本申请的一个实施例中,f=2,g=1,h=1,即,该式可以表示为:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间;
S1-2-3:当所述转换波属性没有进行过常量压缩处理时,可以按照以下公式计算所述平均速度比:
上式中,k表示第六参量。
在本申请的一个实施例中,k=1,即,该式可以表示为:
上式中,所述压缩系数可以为大于1的常量。
由于时移量变化比较剧烈,直接应用时容易对波形进行拉伸和压缩,从而引起波形产生不同程度的畸变。因此,在本申请中先将所述时移量转换为所述平均速度比,再应用所述平均速度比对纵波、转换波进行数据匹配,从而实现纵波和转换波数据的无畸变匹配。
步骤102:利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
在本申请的一个实施例中,在得到所述平均速度比之后,可以根据所述平均速度比对所述转换波时间进行压缩,得到纵波时间域的第一转换波匹配时间,即匹配后的转换波数据,从而可以进行后续的纵波和转换波联合分析。
可以按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
在本申请的一个实施例中,b=2,c=1,即:该表达式可以表示为:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比。
进一步的,在得到匹配后的纵波时间域的第一转换波数据之后,还可以采用插值的方式对所述第一转换波匹配时间进行重采样,得到重采样处理后的第二转换波匹配时间。例如:转换波时间域中的第10个采样点进行纵波和转换波的数据匹配处理之后,转换至纵波时间域的第一转换波匹配时间中的第7.5个采样点,因此,可以通过插值的方式使得转换后的采样点在整数样点上。采样点的插值可以利用sinc函数采用高保真插值的方法确定。采用sinc函数所设计的滤波器可以采用等波纹逼近的方法,从而可以保证纵波时间的转换波数据的高保真、波形不畸变。
下面结合一个具体的实施例对上述纵波和转换波的数据匹配方法进行具体说明,然而值得注意的是,该具体实施例仅是为了更好地说明本发明,并不构成对本发明的不当限定。
步骤1:在野外利用纵波震源激发地震波并利用检波器记录地震波,按照常规地震资料处理流程对采集的数据进行相对振幅保持的高保真处理,生成振幅随偏移距变化的正常时差校正后的纵波道集和转换波道集。图2是纵波道集示意图,图3是转换波道集示意图。
步骤2:对步骤1所形成的纵波道集和转换波道集进行反演,得到纵波属性和转换波属性。
步骤3:利用动态时间调整对步骤2中的纵波属性和转换波属性进行计算,得到如图4所示的纵波和转换波之间的时移量示意图。
步骤4:利用步骤3中的时移量计算得到如图5所示的平均速度比示意图。
步骤5:利用平均速度比对转换波数据进行压缩,得到纵波时间域的第一转换波匹配时间,即如图6所示的与纵波匹配后的第一转换波匹配时间示意图,从而可以进行后续的纵波和转换波联合分析。
基于同一发明构思,本发明实施例中还提供了一种纵波和转换波的数据匹配装置,如下面的实施例所述。由于纵波和转换波的数据匹配装置解决问题的原理与纵波和转换波的数据匹配方法相似,因此纵波和转换波的数据匹配装置的实施可以参见纵波和转换波的数据匹配方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图7是本发明实施例的纵波和转换波的数据匹配装置的一种结构框图,如图7所示,可以包括:速度比计算模块701、数据转换模块702,下面对该结构进行说明。
速度比计算模块701,可以用于根据所述时移量、所述纵波时间计算所述纵波和转换波之间的平均速度比;
数据转换模块702,可以用于利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
在一个实施例中,所述纵波和转换波的数据匹配装置还可以包括:重采样模块,用于在得到匹配后的纵波时间域的第一转换波数据之后,采用插值的方式对所述第一转换波数据进行重采样,得到重采样处理后的第二转换波数据。
在一个实施例中,所述数据转换模块具体可以用于按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
在一个实施例中,所述匹配数据集还可以包括:转换波属性,所述转换波属性可以是对转换波道集进行叠前反演得到的,所述转换波道集可以是对获取的地震资料进行正常时差校正后得到的。
在一个实施例中,所述速度比计算模块可以包括:
判断单元,可以用于确定所述转换波属性是否进行过常量压缩处理;
第一确定单元,可以用于在确定所述转换波属性进行过常量压缩处理的情况下,按照以下公式计算所述平均速度比:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间,f表示第三参量,g表示第四参量,h表示第五参量;
第二确定单元,用于在确定所述转换波属性没有进行过常量压缩处理的情况下,按照以下公式计算所述平均速度比:
上式中,k表示第六参量。
在一个实施例中,所述压缩系数可以为大于1的常量。
尽管本申请内容中提到第一转换波匹配时间的计算方式、匹配数据集包含的内容、平均速度比的计算方式、压缩系数的取值范围等描述,但是,本申请并不局限于必须是符合本申请实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据计算/选取等获取的实施例,仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的单元、装置或模块等,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (7)
1.一种纵波和转换波的数据匹配方法,其特征在于,提供有:匹配数据集,所述匹配数据集包括以下至少之一:纵波和转换波之间的时移量、与所述纵波相对应的纵波时间、与所述转换波相对应的转换波时间、以及转换波属性,所述转换波属性是对转换波道集进行叠前反演得到的,所述转换波道集是对获取的地震资料进行正常时差校正后得到的,所述方法包括:
根据所述时移量、所述纵波时间计算所述纵波属性和转换波属性之间的平均速度比,包括:
确定所述转换波属性是否进行过常量压缩处理;
当所述转换波属性进行过常量压缩处理时,按照以下公式计算所述平均速度比:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间,f表示第三参量,g表示第四参量,h表示第五参量;
当所述转换波属性没有进行过常量压缩处理时,按照以下公式计算所述平均速度比:
上式中,k表示第六参量;
利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
2.如权利要求1所述的方法,其特征在于,在得到匹配后的纵波时间域的第一转换波数据之后,所述方法还包括:
采用插值的方式对所述第一转换波匹配时间进行重采样,得到重采样处理后的第二转换波匹配时间。
3.如权利要求1所述的方法,其特征在于,按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
4.如权利要求1所述的方法,其特征在于,所述压缩系数为大于1的常量。
5.一种纵波和转换波的数据匹配装置,其特征在于,提供有:匹配数据集,所述匹配数据集包括以下至少之一:纵波和转换波之间的时移量、与所述纵波相对应的纵波时间、与所述转换波相对应的转换波时间、以及转换波属性,所述转换波属性是对转换波道集进行叠前反演得到的,所述转换波道集是对获取的地震资料进行正常时差校正后得到的,所述装置包括:
速度比计算模块,用于根据所述时移量、所述纵波时间计算所述纵波属性和转换波属性之间的平均速度比,包括:
判断单元,用于确定所述转换波属性是否进行过常量压缩处理;
第一确定单元,用于在确定所述转换波属性进行过常量压缩处理的情况下,按照以下公式计算所述平均速度比:
上式中,γ0表示所述平均速度比,C表示压缩系数,Δt表示所述时移量,tpp表示所述纵波时间,f表示第三参量,g表示第四参量,h表示第五参量;
第二确定单元,用于在确定所述转换波属性没有进行过常量压缩处理的情况下,按照以下公式计算所述平均速度比:
上式中,k表示第六参量;
数据转换模块,用于利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间。
6.如权利要求5所述的装置,其特征在于,还包括:
重采样模块,用于在得到匹配后的纵波时间域的第一转换波数据之后,采用插值的方式对所述第一转换波数据进行重采样,得到重采样处理后的第二转换波数据。
7.如权利要求5所述的装置,其特征在于,所述数据转换模块具体用于按照以下公式利用所述平均速度比,将转换波时间域的所述转换波时间转换为纵波时间域的第一转换波匹配时间:
上式中,tpp表示所述第一转换波匹配时间,tps表示所述转换波时间,γ0表示所述平均速度比,b表示第一参量,c表示第二参量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710598851.4A CN107589452B (zh) | 2017-07-21 | 2017-07-21 | 纵波和转换波的数据匹配方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710598851.4A CN107589452B (zh) | 2017-07-21 | 2017-07-21 | 纵波和转换波的数据匹配方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107589452A CN107589452A (zh) | 2018-01-16 |
CN107589452B true CN107589452B (zh) | 2019-06-11 |
Family
ID=61042595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710598851.4A Active CN107589452B (zh) | 2017-07-21 | 2017-07-21 | 纵波和转换波的数据匹配方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107589452B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031261B (zh) * | 2018-06-01 | 2020-08-14 | 中国人民解放军战略支援部队信息工程大学 | 一种时差估计方法及装置 |
CN110618448B (zh) * | 2018-06-19 | 2021-02-09 | 中国石油化工股份有限公司 | 多波动态匹配方法及系统 |
CN112363223A (zh) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | 横波速度的获取方法及装置 |
CN112379439A (zh) * | 2020-11-02 | 2021-02-19 | 中国石油天然气集团有限公司 | 地震数据中纵波与横波匹配的方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445785A (zh) * | 2014-06-12 | 2016-03-30 | 中国石油化工股份有限公司 | 一种横波偏移速度建模方法 |
CN106226811A (zh) * | 2016-07-05 | 2016-12-14 | 中国石油天然气集团公司 | 一种基于属性的纵波数据和转换波数据匹配的方法及装置 |
CN106646611A (zh) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
-
2017
- 2017-07-21 CN CN201710598851.4A patent/CN107589452B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445785A (zh) * | 2014-06-12 | 2016-03-30 | 中国石油化工股份有限公司 | 一种横波偏移速度建模方法 |
CN106226811A (zh) * | 2016-07-05 | 2016-12-14 | 中国石油天然气集团公司 | 一种基于属性的纵波数据和转换波数据匹配的方法及装置 |
CN106646611A (zh) * | 2016-12-08 | 2017-05-10 | 中国石油天然气集团公司 | 一种纵波和转换波地震数据时间匹配的方法及装置 |
Non-Patent Citations (2)
Title |
---|
一种PS波与PP波反射时间匹配方法;李生杰等;《中山大学学报(自然科学版)》;20111130;第50卷(第6期);第43-47页 |
基于动态图像变形的PP与PS波层位直接匹配;蒋雪珍等;《地球物理学报》;20170331;第60卷(第3期);第1111-1114页 |
Also Published As
Publication number | Publication date |
---|---|
CN107589452A (zh) | 2018-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107589452B (zh) | 纵波和转换波的数据匹配方法和装置 | |
US9702993B2 (en) | Multi-parameter inversion through offset dependent elastic FWI | |
Huang et al. | A new spectral representation of earthquake data: Hilbert spectral analysis of station TCU129, Chi-Chi, Taiwan, 21 September 1999 | |
CN102466819B (zh) | 地震信号的频谱分析方法及装置 | |
CN108646293B (zh) | 基于黏声拟微分方程的黏声起伏地表正演模拟系统及方法 | |
CN112327358B (zh) | 一种粘滞性介质中声波地震数据正演模拟方法 | |
CN104808243B (zh) | 一种叠前地震贝叶斯反演方法和装置 | |
CN104570072A (zh) | 一种粘弹性介质中的球面pp波反射系数建模方法 | |
CN105425289A (zh) | 确定低频波阻抗的方法和装置 | |
CN104237938B (zh) | 利用动态时间调整进行纵波和转换波数据匹配的方法 | |
Sykulski et al. | Frequency-domain stochastic modeling of stationary bivariate or complex-valued signals | |
CN103777238A (zh) | 一种纯纵波各向异性波场模拟方法 | |
CN103064115A (zh) | 一种射线参数域纵波与转换波匹配方法 | |
CN109143340A (zh) | 一种基于常q模型的粘弹介质地震波模拟方法及系统 | |
CN107229075B (zh) | 深度域地震子波的确定方法和装置 | |
CN106569262A (zh) | 低频地震数据缺失下的背景速度模型重构方法 | |
Amundsen et al. | Elimination of temporal dispersion from the finite-difference solutions of wave equations in elastic and anelastic models | |
CN102866426A (zh) | 一种利用avo大角度道集分析岩体油气信息的方法 | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN110988990A (zh) | 一种高精度的地震属性反演方法 | |
Du et al. | Study on Optical Fiber Gas-Holdup Meter Signal Denoising Using Improved Threshold Wavelet Transform | |
CN111257930B (zh) | 一种黏弹各向异性双相介质区域变网格求解算子 | |
CN107561586B (zh) | 一种气泡压制的方法和装置 | |
CN106291675A (zh) | 一种基于基追踪技术的地震数据重构方法 | |
CN106291680A (zh) | 一种数据低频延拓方法 |
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 |