CN107817525A - 一种近地表常速填充方法及装置 - Google Patents
一种近地表常速填充方法及装置 Download PDFInfo
- Publication number
- CN107817525A CN107817525A CN201710905952.1A CN201710905952A CN107817525A CN 107817525 A CN107817525 A CN 107817525A CN 201710905952 A CN201710905952 A CN 201710905952A CN 107817525 A CN107817525 A CN 107817525A
- Authority
- CN
- China
- Prior art keywords
- layer top
- migration
- near surface
- speed
- speed layer
- 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
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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. analysis, for interpretation, for correction
Abstract
本申请实施例提供了一种近地表常速填充方法及装置,该方法包括:获取目标近地表的高速层顶界面信息及炮、检点静校正量;确定偏移基准面;根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;将所述统一基准面上的CMP道集数据校正到所述偏移基准面;确定所述目标近地表的叠前深度偏移初始速度模型;将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。本申请实施例可有利于提高叠前深度偏移成像精度。
Description
技术领域
本申请涉及地震资料处理技术领域,尤其是涉及一种近地表常速填充方法及装置。
背景技术
由于波场畸变等因素的影响,目前叠前深度偏移还无法实现真地表偏移,更多地还是在地表高程平滑面上建立速度模型进行偏移。而在建立叠前深度偏移速度模型时,一般不考虑近地表的影响。原因在于,其假定近地表的影响在偏前道集准备过程中通过基准面静校正和地表一致性剩余静校正方法已解决。这种假定减弱了地表高程平滑面上波场畸变的影响,有利于深度域速度模型的有效建立。但是,该方法建立的深度-速度模型会降低叠前深度偏移成像的精度,从而人为地降低了速度模型的精度,不利于深度偏移的准确成像。因此目前亟需一种可提高叠前深度偏移成像精度的近地表常速填充方案。
发明内容
本申请实施例的目的在于提供一种近地表常速填充方法及装置,以利于提高叠前深度偏移成像精度。
为达到上述目的,一方面,本申请实施例提供了一种近地表常速填充方法,包括:
获取目标近地表的高速层顶界面信息及炮、检点静校正量;
确定偏移基准面;
根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
确定所述目标近地表的叠前深度偏移初始速度模型;
将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
优选的,所述获取目标近地表的高速层顶界面信息及炮、检点静校正量,包括:
基于初至层析反演静校正获取目标近地表的高速层顶界面信息及炮、检点静校正量。
优选的,所述偏移基准面包括地表高程平滑面。
优选的,所述叠前深度偏移初始速度模型由所述目标近地表的叠前时间偏移速度转换得到。
优选的,在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时通过以下公式计算得到:
Tdata2smooth=Tdata-ΔTS-ΔTG+(ES-EV)S/Vr+(ES-EV)G/Vr;
其中,Tdata2smooth为在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时;Tdata为原始CMP道集数据的走时,ΔTS为炮点从地表经过风化层到高速层顶界面所需走时,ΔTG为检波点从地表经过风化层到高速层顶界面所需走时,ES为地表高程平滑面高程,EV为初至反演时得到的高速层顶界面高程;(ES-EV)S为炮点处ES和EV间的高程差;(ES-EV)G为检波点处ES和EV间的高程差;Vr为替换速度。
另一方面,本申请实施例还提供了一种近地表常速填充装置,包括:
信息及静校正量获取模块,用于获取目标近地表的高速层顶界面信息及炮、检点静校正量;
偏移基准面确定模块,用于确定偏移基准面;
第一校正模块,用于根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
第二校正模块,用于将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
初始模型确定模块,用于确定所述目标近地表的叠前深度偏移初始速度模型;
初始模型调整模块,用于将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
优选的,所述获取目标近地表的高速层顶界面信息及炮、检点静校正量,包括:
基于初至层析反演静校正获取目标近地表的高速层顶界面信息及炮、检点静校正量。
优选的,所述偏移基准面包括地表高程平滑面。
优选的,所述叠前深度偏移初始速度模型由所述目标近地表的叠前时间偏移速度转换得到。
优选的,在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时通过以下公式计算得到:
Tdata2smooth=Tdata-ΔTS-ΔTG+(ES-EV)S/Vr+(ES-EV)G/Vr;
其中,Tdata2smooth为在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时;Tdata为原始CMP道集数据的走时,ΔTS为炮点从地表经过风化层到高速层顶界面所需走时,ΔTG为检波点从地表经过风化层到高速层顶界面所需走时,ES为地表高程平滑面高程,EV为初至反演时得到的高速层顶界面高程;(ES-EV)S为炮点处ES和EV间的高程差;(ES-EV)G为检波点处ES和EV间的高程差;Vr为替换速度。
再一方面,本申请实施例还提供了另一种近地表常速填充装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的高速层顶界面信息及炮、检点静校正量;
确定偏移基准面;
根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
确定所述目标近地表的叠前深度偏移初始速度模型;
将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
由以上本申请实施例提供的技术方案可见,本申请实施例在建模时引入高速层顶界面作为一个控制层位,即将目标近地表位于高速层顶界面之上且位于偏移基准面之下的部分填充替换速度。这样就实现了表层速度模型与偏移前数据的匹配,而且表层速度建模问题也得到了有效解决,从而为后续浅层乃至中、深层速度建模奠定了坚实的基础,进而有利于进一步提高速度模型精度和偏移成像精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请一实施例的表层结构示意图;
图2为本申请一实施例中近地表常速填充方法的流程图;
图3为本申请一实施例中替换速度填充前叠前深度偏移剖面示意图;
图4为本申请一实施例中替换速度填充后叠前深度偏移剖面示意图;
图5为本申请一实施例中近地表常速填充装置的结构框图;
图6为本申请另一实施例中近地表常速填充装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
如图1所示,在以往表层速度建模准备偏移前数据时,一般先根据统一基准面静校正公式:
Tdata2datum=Tdata-ΔTS+(ED-EV)S/Vr-ΔTG+(ED-EV)G/Vr (1)
将地震数据校到统一基准面上,然后根据公式(2)将地震数据校正到偏移基准面上,再在此基础上应用反射信息迭代建立偏移速度场。
Tdata2smooth=Tdata2datum-(ED-ES)S/Vr-(ED-ES)G/Vr (2)
其中,Tdata2datum表示将地震数据校到统一基准面后的走时,Tdata表示原始地震数据走时,ΔTS为炮点从地表经过风化层到高速层顶界面所需走时,ΔTG为检波点从地表经过风化层到高速层顶界面所需走时,ED为统一基准面高程,Es为地表平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度(在陆上地震数据处理过程中,通常要将地震数据校正到一个统一基准面,而且这个统一基准面一般为水平面。将地震数据校正到基准面需要对炮点和检波点的高程差进行校正,这里需引入替换速度。当统一基准面高于地表时,相当于剥去基准面以上的地层;当基准面低于地表时,相当于在地表与基准面之间填充了地层;而剥去地层和填充地层的速度即为替换速度);Tdata2smooth为将地震数据校正到偏移基准面后的走时。此外,V1为高速层速度,Vw为初至反演所得风化层速度,S表示炮点,G表示检波点,CMP表示共中心点,CIP表示成像点,E为地表高程,Hs为地表平滑面到高速层顶界面间的距离,△T为地震波从地表到高速层顶界面所需旅行时,H为同一空间位置处不同地形面之间的高程差。
研究表明,现有技术的这种方法假定偏移基准面之下的速度可以由反射信息反演得到,但实际上由于浅层覆盖次数低,这种假定是不成立的,或者说反演出的速度精度较低。浅、表层速度不准会严重影响下伏地层速度反演工作,进而造成偏移速度不准确,偏移成像差。为解决这一问题本申请的发明人提出了本申请的技术方案。具体请详见以下部分的描述。
参考图2所示,本申请一实施例的近地表常速填充方法可以包括以下步骤:
S201、获取目标近地表的高速层顶界面信息及炮、检点静校正量。
在本申请的一些实施方式中,可基于初至层析反演静校正技术获取目标近地表的高速层顶界面信息及炮、检点静校正量。其中所述高速层顶界面信息例如可以包括高速层顶界面高程、高速层速度等。
S202、确定偏移基准面。
偏移基准面是人为定义的偏移参考面,在本申请的一些实施方式中,例如可以将地表高程平滑面作为偏移基准面。
S203、根据所述炮、检点静校正量将CMP道集数据校正到统一基准面。
在本申请的一些实施方式中,统一基准面是人为定义的参考面,它是地震剖面的起始零线,地震剖面的时间可以以这个基准面为参考。把数据校正到这个面上后,相当于激发点和检波点都位于这个基准面上。所述根据所述炮、检点静校正量将CMP道集数据校正到统一基准面,即根据所述炮、检点静校正量对CMP道集数据进行校正,从而将CMP道集数据正到统一基准面上。
S204、将所述统一基准面上的CMP道集数据校正到所述偏移基准面。
在本申请的一些实施方式中,通过将所述统一基准面上的CMP道集数据校正到所述偏移基准面,可以生成偏移基准面上的偏移前CMP道集数据,以便于后续使用。而在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时通过以下公式(3)计算得到:
Tdata2smooth=Tdata-ΔTS-ΔTG+(ES-EV)S/Vr+(ES-EV)G/Vr; (3)
其中,Tdata2smooth为在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时;Tdata为原始CMP道集数据的走时;ES为地表高程平滑面高程;(ES-EV)S为炮点处ES和EV间的高程差;(ES-EV)G为检波点处ES和EV间的高程差。
S205、确定所述目标近地表的叠前深度偏移初始速度模型。
在本申请的一些实施方式中,所述叠前深度偏移初始速度模型可由所述目标近地表的叠前时间偏移速度转换得到。
S206、将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
上由述公式(3)可以看出,公式(3)右侧前三项表示地震数据由地表经过风化层校正到高速层顶界面上,后两项表示通过替换速度将地震数据从高速层顶界面校正到偏移基准面上,意味着将风化层速度替换为替换速度,高速层顶界面以上为替换速度的各向同性介质,因此,基于该思路,本申请在建模时引入高速层顶界面作为一个控制层位,即将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。这样就实现了表层速度模型与偏移前数据的匹配,而且表层速度建模问题也得到了有效解决,从而为后续浅层乃至中、深层速度建模奠定了坚实的基础,进而有利于进一步提高速度模型精度和偏移成像精度。
图3和图4分别为库车盆地某三维地震资料在替换速度填充前后的偏移剖面。对比图3和图4可看出,在替换速度填充后,浅层深度偏移成像精度明显改善。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
结合图5所示,本申请一实施例的近地表常速填充装置可以包括:
信息及静校正量获取模块51,可以用于获取目标近地表的高速层顶界面信息及炮、检点静校正量;
偏移基准面确定模块52,可以用于确定偏移基准面;
第一校正模块53,可以用于根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
第二校正模块54,可以用于将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
初始模型确定模块55,可以用于确定所述目标近地表的叠前深度偏移初始速度模型;
初始模型调整模块56,可以用于将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
结合图6所示,本申请另一实施例的近地表常速填充装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的高速层顶界面信息及炮、检点静校正量;
确定偏移基准面;
根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
确定所述目标近地表的叠前深度偏移初始速度模型;
将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
同样,本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法和设备的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (11)
1.一种近地表常速填充方法,其特征在于,包括:
获取目标近地表的高速层顶界面信息及炮、检点静校正量;
确定偏移基准面;
根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
确定所述目标近地表的叠前深度偏移初始速度模型;
将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
2.如权利要求1所述的近地表常速填充方法,其特征在于,所述获取目标近地表的高速层顶界面信息及炮、检点静校正量,包括:
基于初至层析反演静校正获取目标近地表的高速层顶界面信息及炮、检点静校正量。
3.如权利要求1所述的近地表常速填充方法,其特征在于,所述偏移基准面包括地表高程平滑面。
4.如权利要求1所述的近地表常速填充方法,其特征在于,所述叠前深度偏移初始速度模型由所述目标近地表的叠前时间偏移速度转换得到。
5.如权利要求1所述的近地表常速填充方法,其特征在于,在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时通过以下公式计算得到:
Tdata2smooth=Tdata-ΔTS-ΔTG+(ES-EV)S/Vr+(ES-EV)G/Vr;
其中,Tdata2smooth为在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时;Tdata为原始CMP道集数据的走时,ΔTS为炮点从地表经过风化层到高速层顶界面所需走时,ΔTG为检波点从地表经过风化层到高速层顶界面所需走时,ES为地表高程平滑面高程,EV为初至反演时得到的高速层顶界面高程;(ES-EV)S为炮点处ES和EV间的高程差;(ES-EV)G为检波点处ES和EV间的高程差;Vr为替换速度。
6.一种近地表常速填充装置,其特征在于,包括:
信息及静校正量获取模块,用于获取目标近地表的高速层顶界面信息及炮、检点静校正量;
偏移基准面确定模块,用于确定偏移基准面;
第一校正模块,用于根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
第二校正模块,用于将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
初始模型确定模块,用于确定所述目标近地表的叠前深度偏移初始速度模型;
初始模型调整模块,用于将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
7.如权利要求6所述的近地表常速填充装置,其特征在于,所述获取目标近地表的高速层顶界面信息及炮、检点静校正量,包括:
基于初至层析反演静校正获取目标近地表的高速层顶界面信息及炮、检点静校正量。
8.如权利要求6所述的近地表常速填充装置,其特征在于,所述偏移基准面包括地表高程平滑面。
9.如权利要求6所述的近地表常速填充装置,其特征在于,所述叠前深度偏移初始速度模型由所述目标近地表的叠前时间偏移速度转换得到。
10.如权利要求6所述的近地表常速填充装置,其特征在于,在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时通过以下公式计算得到:
Tdata2smooth=Tdata-ΔTS-ΔTG+(ES-EV)S/Vr+(ES-EV)G/Vr;
其中,Tdata2smooth为在将所述统一基准面上的CMP道集数据校正到所述偏移基准面后,所述CMP道集数据的走时;Tdata为原始CMP道集数据的走时,ΔTS为炮点从地表经过风化层到高速层顶界面所需走时,ΔTG为检波点从地表经过风化层到高速层顶界面所需走时,ES为地表高程平滑面高程,EV为初至反演时得到的高速层顶界面高程;(ES-EV)S为炮点处ES和EV间的高程差;(ES-EV)G为检波点处ES和EV间的高程差;Vr为替换速度。
11.一种装置近地表常速填充装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,其特征在于,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的高速层顶界面信息及炮、检点静校正量;
确定偏移基准面;
根据所述炮、检点静校正量将CMP道集数据校正到统一基准面;
将所述统一基准面上的CMP道集数据校正到所述偏移基准面;
确定所述目标近地表的叠前深度偏移初始速度模型;
将所述高速层顶界面信息中的高速层顶界面加入所述叠前深度偏移初始速度模型,并将所述目标近地表位于所述高速层顶界面之上且位于所述偏移基准面之下的部分填充替换速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710905952.1A CN107817525B (zh) | 2017-09-29 | 2017-09-29 | 一种近地表常速填充方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710905952.1A CN107817525B (zh) | 2017-09-29 | 2017-09-29 | 一种近地表常速填充方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107817525A true CN107817525A (zh) | 2018-03-20 |
CN107817525B CN107817525B (zh) | 2019-09-10 |
Family
ID=61607036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710905952.1A Active CN107817525B (zh) | 2017-09-29 | 2017-09-29 | 一种近地表常速填充方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107817525B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109581495A (zh) * | 2018-10-24 | 2019-04-05 | 中国石油天然气集团有限公司 | 山地表层速度模型构建方法及系统 |
CN110967752A (zh) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | 一种起伏地表偏移输入地震数据的静校正方法 |
CN112630838A (zh) * | 2019-09-24 | 2021-04-09 | 中国石油化工股份有限公司 | 一种时间域浮动面转换深度域计算方法及计算系统 |
CN112764107A (zh) * | 2019-11-04 | 2021-05-07 | 中国石油天然气集团有限公司 | 地表高速覆盖地区的地震资料静校正方法及装置 |
CN112859167A (zh) * | 2019-11-28 | 2021-05-28 | 中国石油天然气股份有限公司 | 畸变地质体速度场的校正方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040015297A1 (en) * | 2002-04-17 | 2004-01-22 | Jean-Claude Dubois | Method for determining primary static corrections to be applied to seismic traces |
CN1580817A (zh) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | 分步静校正计算方法 |
US20050256648A1 (en) * | 2004-05-11 | 2005-11-17 | West Michael P | Velocity determination of the near-surface layers in the earth using exploration 2D or 3D seismic data |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
CN102736104A (zh) * | 2012-05-31 | 2012-10-17 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 近地表表层模型构建方法 |
CN103728666A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种确定静校正中间参考面替换速度的方法 |
-
2017
- 2017-09-29 CN CN201710905952.1A patent/CN107817525B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040015297A1 (en) * | 2002-04-17 | 2004-01-22 | Jean-Claude Dubois | Method for determining primary static corrections to be applied to seismic traces |
CN1580817A (zh) * | 2003-08-01 | 2005-02-16 | 中国石油天然气集团公司 | 分步静校正计算方法 |
US20050256648A1 (en) * | 2004-05-11 | 2005-11-17 | West Michael P | Velocity determination of the near-surface layers in the earth using exploration 2D or 3D seismic data |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
CN102736104A (zh) * | 2012-05-31 | 2012-10-17 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 近地表表层模型构建方法 |
CN103728666A (zh) * | 2012-10-16 | 2014-04-16 | 中国石油化工股份有限公司 | 一种确定静校正中间参考面替换速度的方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110967752A (zh) * | 2018-09-29 | 2020-04-07 | 中国石油化工股份有限公司 | 一种起伏地表偏移输入地震数据的静校正方法 |
CN110967752B (zh) * | 2018-09-29 | 2022-05-03 | 中国石油化工股份有限公司 | 一种起伏地表偏移输入地震数据的静校正方法 |
CN109581495A (zh) * | 2018-10-24 | 2019-04-05 | 中国石油天然气集团有限公司 | 山地表层速度模型构建方法及系统 |
CN112630838A (zh) * | 2019-09-24 | 2021-04-09 | 中国石油化工股份有限公司 | 一种时间域浮动面转换深度域计算方法及计算系统 |
CN112764107A (zh) * | 2019-11-04 | 2021-05-07 | 中国石油天然气集团有限公司 | 地表高速覆盖地区的地震资料静校正方法及装置 |
CN112859167A (zh) * | 2019-11-28 | 2021-05-28 | 中国石油天然气股份有限公司 | 畸变地质体速度场的校正方法及装置 |
CN112859167B (zh) * | 2019-11-28 | 2023-08-22 | 中国石油天然气股份有限公司 | 畸变地质体速度场的校正方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107817525B (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107817525B (zh) | 一种近地表常速填充方法及装置 | |
CN107783186B (zh) | 一种近地表速度建模方法及装置 | |
CN106990438B (zh) | 复杂地表区域微测井点位及井深确定方法及装置 | |
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN107765298B (zh) | 一种确定近地表速度模型的方法及装置 | |
CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN105301636B (zh) | 速度模型的建立方法和装置 | |
CN106959467B (zh) | 地震波速度反演方法和装置 | |
CN109917454A (zh) | 基于双基准面的真地表叠前深度偏移成像方法及装置 | |
CN104199103A (zh) | 转换波静校正方法及装置 | |
CN107329174B (zh) | 在地震剖面上显示测井信息的方法及装置 | |
Kanaujia et al. | 1D velocity structure and characteristics of contemporary local seismicity around the Tehri region, Garhwal Himalaya | |
CN108845350A (zh) | 反演二维速度模型的方法及装置 | |
CN108445532B (zh) | 一种深度域反偏移方法及装置 | |
CN108957539B (zh) | 层析偏移速度分析中的射线追踪方法及装置 | |
CN108828669B (zh) | 一种二维相交测线静校正处理方法、装置及系统 | |
US20200049844A1 (en) | Computer implemented method for improving a velocity model for seismic imaging | |
CN106443829A (zh) | 一种近地表模型构建方法及装置 | |
CN111399037B (zh) | 高速顶界面提取的方法和装置 | |
CN111208558B (zh) | 超深低幅度三维地质构造的建立方法及装置 | |
Chambers et al. | A practical implementation of wave front construction for 3-D isotropic media | |
CN106597536A (zh) | 一种地震成像的处理方法及装置 | |
CN113495296B (zh) | 层析静校正量的确定方法、装置、设备及可读存储介质 | |
CN107229074A (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 |