CN107783186A - 一种近地表速度建模方法及装置 - Google Patents

一种近地表速度建模方法及装置 Download PDF

Info

Publication number
CN107783186A
CN107783186A CN201710905940.9A CN201710905940A CN107783186A CN 107783186 A CN107783186 A CN 107783186A CN 201710905940 A CN201710905940 A CN 201710905940A CN 107783186 A CN107783186 A CN 107783186A
Authority
CN
China
Prior art keywords
earth
skew
shot point
correction value
speed
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
Application number
CN201710905940.9A
Other languages
English (en)
Other versions
CN107783186B (zh
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 CN201710905940.9A priority Critical patent/CN107783186B/zh
Publication of CN107783186A publication Critical patent/CN107783186A/zh
Application granted granted Critical
Publication of CN107783186B publication Critical patent/CN107783186B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01V1/282Application of seismic models, synthetic seismograms

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

本申请实施例提供了一种近地表速度建模方法及装置,该方法包括:获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;根据所述走时更新所述静校正量;根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;将所述偏移前CMP道集数据校正到所述偏移基准面上;根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。本申请实施例可提高建模精度和叠前深度偏移成像精度。

Description

一种近地表速度建模方法及装置
技术领域
本申请涉及地震资料处理技术领域,尤其是涉及一种近地表速度建模方法及装置。
背景技术
由于波场畸变等因素的影响,目前叠前深度偏移还无法实现真地表偏移,更多地还是在地表高程平滑面上建立速度模型进行偏移。而在建立叠前深度偏移速度模型时,一般不考虑近地表的影响。原因在于,其假定近地表的影响在偏前道集准备过程中通过基准面静校正和地表一致性剩余静校正方法已解决。这种假定减弱了地表高程平滑面上波场畸变的影响,有利于深度域速度模型的有效建立。但是,该方法建立的深度-速度模型会降低叠前深度偏移成像的精度,从而人为地降低了速度模型的精度,不利于深度偏移的准确成像。因此目前亟需一种可提高叠前深度偏移成像精度的近地表速度建模方案。
发明内容
本申请实施例的目的在于提供一种近地表速度建模方法及装置,以提高建模精度和叠前深度偏移成像精度。
为达到上述目的,一方面,本申请实施例提供了近地表速度建模方法,包括:
获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
根据所述走时更新所述静校正量;
根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
将所述偏移前CMP道集数据校正到所述偏移基准面上;
根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
优选的,所述获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量,包括:
基于初至层析反演静校正获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量。
优选的,所述确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时,包括:
根据公式ΔTs=∫dHs/Vw确定炮点从所述高速顶界面至所述地表高程平滑面所需走时ΔTs
根据公式ΔTG=∫dHG/Vw确定检波点从高速顶界面至所述地表高程平滑面所需走时ΔTG
其中,Hs为地表高程平滑面上的炮点至从所述高速顶界面的距离;Vw为初至反演所得风化层速度;HG为地表高程平滑面上的检波点至从所述高速顶界面的距离。
优选的,根据所述走时更新所述静校正量,包括:
根据公式sta_s'=sta_s-ΔT'S获得更新后炮点的静校正量;
根据公式sta_g'=sta_g-ΔT'G获得更新后检波点的静校正量;
其中,sta_s为炮点的静校正量;sta_s'为更新后炮点的静校正量;ΔT'S为炮点从地表经过风化层到高速层顶界面所需走时;sta_g为检波点的静校正量;sta_g'为更新后检波点的静校正量;ΔT'G为检波点从地表经过风化层到高速层顶界面所需走时。
优选的,在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时变化量为:
Δt2=-(ED-ES)S/Vr-(ED-ES)G/Vr
其中,Δt2为在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时;ED为统一基准面高程,ES为地表平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度。
优选的,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息达到预设阈值时,在所述高速层顶界面至所述偏移基准面间填充所述近地表反演速度模型,获得浅表层速度模型。
优选的,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息低于预设阈值时,根据所述高速层顶界面以下的速度信息,获得浅表层速度模型。
另一方面,本申请实施例还提供了一种近地表速度建模装置,包括:
信息获取模块,用于获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
偏移基准面确定模块,用于确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
走时获取模块,用于确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
静校正量更新模块,用于根据所述走时更新所述静校正量;
第一校正模块,用于根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
第二校正模块,用于将所述偏移前CMP道集数据校正到所述偏移基准面上;
模型获取模块,用于根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
优选的,所述获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量,包括:
基于初至层析反演静校正获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量。
优选的,所述确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时,包括:
根据公式ΔTs=∫dHs/Vw确定炮点从所述高速顶界面至所述地表高程平滑面所需走时ΔTs
根据公式ΔTG=∫dHG/Vw确定检波点从高速顶界面至所述地表高程平滑面所需走时ΔTG
其中,Hs为地表高程平滑面上的炮点至从所述高速顶界面的距离;Vw为初至反演所得风化层速度;HG为地表高程平滑面上的检波点至从所述高速顶界面的距离。
优选的,根据所述走时更新所述静校正量,包括:
根据公式sta_s'=sta_s-ΔT'S获得更新后炮点的静校正量;
根据公式sta_g'=sta_g-ΔT'G获得更新后检波点的静校正量;
其中,sta_s为炮点的静校正量;sta_s'为更新后炮点的静校正量;ΔT'S为炮点从地表经过风化层到高速层顶界面所需走时;sta_g为检波点的静校正量;sta_g'为更新后检波点的静校正量;ΔT'G为检波点从地表经过风化层到高速层顶界面所需走时。
优选的,在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数的走时变化量为:
Δt2=-(ED-ES)S/Vr-(ED-ES)G/Vr
其中,Δt2为在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时;ED为统一基准面高程,ES为地表平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度。
优选的,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息达到预设阈值时,在所述高速层顶界面至所述偏移基准面间填充所述近地表反演速度模型,获得浅表层速度模型。
优选的,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息低于预设阈值时,根据所述高速层顶界面以下的速度信息,获得浅表层速度模型。
再一方面,本申请实施例还提供了另一种近地表速度建模装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
根据所述走时更新所述静校正量;
根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
将所述偏移前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、确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时。
在本申请的一些实施方式中,所述确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时例如可以通过以下公式计算得到:
根据公式ΔTs=∫dHs/Vw确定炮点从所述高速顶界面至所述地表高程平滑面所需走时ΔTs
根据公式ΔTG=∫dHG/Vw确定检波点从高速顶界面至所述地表高程平滑面所需走时ΔTG
其中,Hs为地表高程平滑面上的炮点至从所述高速顶界面的距离;Vw为初至反演所得风化层速度;HG为地表高程平滑面上的检波点至从所述高速顶界面的距离。
S204、根据所述走时更新所述静校正量。
在本申请的一些实施方式中,所述根据所述走时更新所述静校正量例如可通过以下公式计算得到:
根据公式sta_s'=sta_s-ΔT'S获得更新后炮点的静校正量;
根据公式sta_g'=sta_g-ΔT'G获得更新后检波点的静校正量;
其中,sta_s为炮点的静校正量;sta_s'为更新后炮点的静校正量;ΔT'S为炮点从地表经过风化层到高速层顶界面所需走时;sta_g为检波点的静校正量;sta_g'为更新后检波点的静校正量;ΔT'G为检波点从地表经过风化层到高速层顶界面所需走时。
S205、根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据。其中,校正前后的走时变化量Δt1为Δt1=-sta_s'-sta_g'。其中,sta_s'为更新后炮点的静校正量,sta_g'为更新后检波点的静校正量。
S206、将所述偏移前CMP道集数据校正到所述偏移基准面上。
通过将所述偏移前CMP道集数校正到所述偏移基准面上,可以生成偏移基准面上的偏移前CMP道集数据,以便于后续使用。而在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时变化量为:
Δt2=-(ED-ES)S/Vr-(ED-ES)G/Vr
其中,Δt2为在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时;ED为统一基准面高程,ES为地表平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度。
S207、根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
在本申请的一些实施方式中,根据不同的情况,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型可以采用不同的方式。例如在本申请的一些实施方式中,当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息达到预设阈值时,在所述高速层顶界面至所述偏移基准面间填充所述近地表反演速度模型,获得浅表层速度模型。例如在本申请的另一些实施方式中,当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息低于预设阈值时,可根据所述高速层顶界面以下的速度信息,获得浅表层速度模型。也就是说当高速顶界面以下的速度信息无法用反射信息得到时,可以参考初至层析反演所得速度模型对高速顶界面以下的速度信息进行预测。这种预测包括速度区域变化规律和趋势的预测(定性)和速度值的预测(定量)。通过这种预测可得到反射信息无法得到的速度信息并将其填充到浅表层速度模型中。
本申请实施例将初至层析反演得到的近地表速度模型应用于叠前深度偏移浅表层速度建模过程中,从而有效避免了因浅层覆盖次数低,有效反射信息少等造成的浅表层偏移速度不准问题,因而提高了速度模型的精度,从而达到了进一步提高叠前深度偏移成像精度的目的。
图3和图4分别为库车盆地某三维地震资料在近地表反演速度模型填充前后的叠前深度偏移剖面。对比图3和图4可看出(尤其是黑色圆圈所在部分),图3速度精度低,偏移归位不准确,存在划弧现象。而在近地表反演速度模型填充后,浅层深度偏移成像精度明显改善,从而较好地解决了由于速度不准导致的偏移划弧问题。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参考图5所示,本申请实施例的近地表速度建模装置可以包括:
信息获取模块51,可以用于获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
偏移基准面确定模块52,可以用于确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
走时获取模块53,可以用于确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
静校正量更新模块54,可以用于根据所述走时更新所述静校正量;
第一校正模块55,可以用于根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
第二校正模块56,可以用于将所述偏移前CMP道集数据校正到所述偏移基准面上;
模型获取模块57,可以用于根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
结合图6所示,本申请另一实施例的近地表速度建模装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
根据所述走时更新所述静校正量;
根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
将所述偏移前CMP道集数据校正到所述偏移基准面上;
根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
同样,本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法和设备的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (15)

1.一种近地表速度建模方法,其特征在于,包括:
获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
根据所述走时更新所述静校正量;
根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
将所述偏移前CMP道集数据校正到所述偏移基准面上;
根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
2.如权利要求1所述的近地表速度建模方法,其特征在于,所述获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量,包括:
基于初至层析反演静校正获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量。
3.如权利要求1所述的近地表速度建模方法,其特征在于,所述确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时,包括:
根据公式ΔTs=∫dHs/Vw确定炮点从所述高速顶界面至所述地表高程平滑面所需走时ΔTs
根据公式ΔTG=∫dHG/Vw确定检波点从高速顶界面至所述地表高程平滑面所需走时ΔTG
其中,Hs为地表高程平滑面上的炮点至所述高速顶界面的距离;Vw为初至反演所得风化层速度;HG为地表高程平滑面上的检波点至所述高速顶界面的距离。
4.如权利要求1所述的近地表速度建模方法,其特征在于,根据所述走时更新所述静校正量,包括:
根据公式sta_s'=sta_s-ΔT′S获得更新后炮点的静校正量;
根据公式sta_g'=sta_g-ΔT′G获得更新后检波点的静校正量;
其中,sta_s为炮点的静校正量;sta_s'为更新后炮点的静校正量;ΔT′S为炮点从地表经过风化层到高速层顶界面所需走时;sta_g为检波点的静校正量;sta_g'为更新后检波点的静校正量;ΔT′G为检波点从地表经过风化层到高速层顶界面所需走时。
5.如权利要求1所述的近地表速度建模方法,其特征在于,在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时变化量为:
Δt2=-(ED-ES)S/Vr-(ED-ES)G/Vr
其中,Δt2为在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时;ED为统一基准面高程,ES为地表高程平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度。
6.如权利要求1所述的近地表速度建模方法,其特征在于,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息达到预设阈值时,在所述高速层顶界面至所述偏移基准面间填充所述近地表反演速度模型,获得浅表层速度模型。
7.如权利要求1所述的近地表速度建模方法,其特征在于,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息低于预设阈值时,根据所述高速层顶界面以下的速度信息,获得浅表层速度模型。
8.一种近地表速度建模装置,其特征在于,包括:
信息获取模块,用于获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
偏移基准面确定模块,用于确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
走时获取模块,用于确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
静校正量更新模块,用于根据所述走时更新所述静校正量;
第一校正模块,用于根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
第二校正模块,用于将所述偏移前CMP道集数据校正到所述偏移基准面上;
模型获取模块,用于根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
9.如权利要求8所述的近地表速度建模装置,其特征在于,所述获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量,包括:
基于初至层析反演静校正获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量。
10.如权利要求8所述的近地表速度建模装置,其特征在于,所述确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时,包括:
根据公式ΔTs=∫dHs/Vw确定炮点从所述高速顶界面至所述地表高程平滑面所需走时ΔTs
根据公式ΔTG=∫dHG/Vw确定检波点从高速顶界面至所述地表高程平滑面所需走时ΔTG
其中,Hs为地表高程平滑面上的炮点至从所述高速顶界面的距离;Vw为初至反演所得风化层速度;HG为地表高程平滑面上的检波点至从所述高速顶界面的距离。
11.如权利要求8所述的近地表速度建模装置,其特征在于,根据所述走时更新所述静校正量,包括:
根据公式sta_s'=sta_s-ΔT′S获得更新后炮点的静校正量;
根据公式sta_g'=sta_g-ΔT′G获得更新后检波点的静校正量;
其中,sta_s为炮点的静校正量;sta_s'为更新后炮点的静校正量;ΔT′S为炮点从地表经过风化层到高速层顶界面所需走时;sta_g为检波点的静校正量;sta_g'为更新后检波点的静校正量;ΔT′G为检波点从地表经过风化层到高速层顶界面所需走时。
12.如权利要求8所述的近地表速度建模装置,其特征在于,在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时变化量为:
Δt2=-(ED-ES)S/Vr-(ED-ES)G/Vr
其中,Δt2为在将所述偏移前CMP道集数据校正到所述偏移基准面上后,所述偏移前CMP道集数据的走时;ED为统一基准面高程,ES为地表平滑面高程;(ED-ES)S为炮点处ED和ES间的高程差;(ED-ES)G为检波点处ED和ES间的高程差;Vr为替换速度。
13.如权利要求8所述的近地表速度建模装置,其特征在于,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息达到预设阈值时,在所述高速层顶界面至所述偏移基准面间填充所述近地表反演速度模型,获得浅表层速度模型。
14.如权利要求8所述的近地表速度建模装置,其特征在于,所述根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型,包括:
当校正到所述偏移基准面上的偏移前CMP道集数据中所包含的高速层顶界面以下的反射信息低于预设阈值时,根据所述高速层顶界面以下的速度信息,获得浅表层速度模型。
15.一种装置近地表速度建模装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,其特征在于,所述计算机程序被所述处理器运行时执行如下步骤:
获取目标近地表的近地表反演速度模型及炮点、检波点的静校正量;
确定目标近地表的地表高程平滑面,并将所述地表高程平滑面确定为偏移基准面;
确定炮点、检波点从高速顶界面至所述地表高程平滑面所需走时;
根据所述走时更新所述静校正量;
根据更新后的静校正量校正CMP道集数据,获得偏移前CMP道集数据;
将所述偏移前CMP道集数据校正到所述偏移基准面上;
根据所述近地表反演速度模型及校正到所述偏移基准面上的偏移前CMP道集数据,获得浅表层速度模型。
CN201710905940.9A 2017-09-29 2017-09-29 一种近地表速度建模方法及装置 Active CN107783186B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710905940.9A CN107783186B (zh) 2017-09-29 2017-09-29 一种近地表速度建模方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710905940.9A CN107783186B (zh) 2017-09-29 2017-09-29 一种近地表速度建模方法及装置

Publications (2)

Publication Number Publication Date
CN107783186A true CN107783186A (zh) 2018-03-09
CN107783186B CN107783186B (zh) 2019-09-10

Family

ID=61433996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710905940.9A Active CN107783186B (zh) 2017-09-29 2017-09-29 一种近地表速度建模方法及装置

Country Status (1)

Country Link
CN (1) CN107783186B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109725354A (zh) * 2018-11-20 2019-05-07 中国石油天然气集团有限公司 各向异性速度建模方法及系统
CN109917454A (zh) * 2019-02-19 2019-06-21 中国石油天然气集团有限公司 基于双基准面的真地表叠前深度偏移成像方法及装置
CN110967752A (zh) * 2018-09-29 2020-04-07 中国石油化工股份有限公司 一种起伏地表偏移输入地震数据的静校正方法
CN111399037A (zh) * 2019-01-02 2020-07-10 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN111435171A (zh) * 2019-01-15 2020-07-21 中国石油集团东方地球物理勘探有限责任公司 海底节点地震数据的时钟漂移校正方法、装置及存储介质
CN112464145A (zh) * 2020-10-29 2021-03-09 中国石油天然气集团有限公司 复杂地表模型重构方法及装置
CN114076979A (zh) * 2020-08-17 2022-02-22 中国石油化工股份有限公司 近地表等效速度建模方法、装置、电子设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN104331594A (zh) * 2014-08-28 2015-02-04 中国石油天然气集团公司 一种叠前深度偏移速度建模中的结构模型生成方法及装置
CN104570102A (zh) * 2013-10-11 2015-04-29 中国石油化工股份有限公司 一种近地表速度模型与中深层速度模型的融合方法
CN106054252A (zh) * 2016-06-23 2016-10-26 中国石油天然气集团公司 一种叠前时间偏移的方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN104570102A (zh) * 2013-10-11 2015-04-29 中国石油化工股份有限公司 一种近地表速度模型与中深层速度模型的融合方法
CN104331594A (zh) * 2014-08-28 2015-02-04 中国石油天然气集团公司 一种叠前深度偏移速度建模中的结构模型生成方法及装置
CN106054252A (zh) * 2016-06-23 2016-10-26 中国石油天然气集团公司 一种叠前时间偏移的方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑鸿明 等: ""近似层替换的基准面静校正方法"", 《新疆地质》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110967752A (zh) * 2018-09-29 2020-04-07 中国石油化工股份有限公司 一种起伏地表偏移输入地震数据的静校正方法
CN110967752B (zh) * 2018-09-29 2022-05-03 中国石油化工股份有限公司 一种起伏地表偏移输入地震数据的静校正方法
CN109725354A (zh) * 2018-11-20 2019-05-07 中国石油天然气集团有限公司 各向异性速度建模方法及系统
CN109725354B (zh) * 2018-11-20 2020-07-10 中国石油天然气集团有限公司 各向异性速度建模方法及系统
CN111399037A (zh) * 2019-01-02 2020-07-10 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN111399037B (zh) * 2019-01-02 2022-12-02 中国石油天然气集团有限公司 高速顶界面提取的方法和装置
CN111435171A (zh) * 2019-01-15 2020-07-21 中国石油集团东方地球物理勘探有限责任公司 海底节点地震数据的时钟漂移校正方法、装置及存储介质
CN109917454A (zh) * 2019-02-19 2019-06-21 中国石油天然气集团有限公司 基于双基准面的真地表叠前深度偏移成像方法及装置
CN109917454B (zh) * 2019-02-19 2020-10-09 中国石油天然气集团有限公司 基于双基准面的真地表叠前深度偏移成像方法及装置
CN114076979A (zh) * 2020-08-17 2022-02-22 中国石油化工股份有限公司 近地表等效速度建模方法、装置、电子设备及介质
CN114076979B (zh) * 2020-08-17 2024-05-07 中国石油化工股份有限公司 近地表等效速度建模方法、装置、电子设备及介质
CN112464145A (zh) * 2020-10-29 2021-03-09 中国石油天然气集团有限公司 复杂地表模型重构方法及装置

Also Published As

Publication number Publication date
CN107783186B (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN107783186A (zh) 一种近地表速度建模方法及装置
CN107817525B (zh) 一种近地表常速填充方法及装置
CN105277978B (zh) 一种确定近地表速度模型的方法及装置
CN106990438A (zh) 复杂地表区域微测井点位及井深确定方法及装置
CN105301636B (zh) 速度模型的建立方法和装置
CN105301639A (zh) 基于vsp旅行时双加权层析反演速度场的方法及其装置
CN106896406A (zh) 基于阻抗域属性预测优质储层的方法和装置
US11125898B2 (en) Method and system for positioning seismic source in microseism monitoring
CN108828669B (zh) 一种二维相交测线静校正处理方法、装置及系统
CN108845350B (zh) 反演二维速度模型的方法及装置
CN106959467A (zh) 地震波速度反演方法和装置
CN107329174B (zh) 在地震剖面上显示测井信息的方法及装置
CN105093318B (zh) 一种自适应波动方程波场延拓静校正方法
CN106054252B (zh) 一种叠前时间偏移的方法及装置
CN104614762B (zh) 疏松砂岩气藏边界确定方法及装置
CN111413741A (zh) 一种砂岩型铀矿资源量计算方法和装置
CN106443829A (zh) 一种近地表模型构建方法及装置
CN111208558B (zh) 超深低幅度三维地质构造的建立方法及装置
CN111399037B (zh) 高速顶界面提取的方法和装置
CN107229074A (zh) 低勘探程度海相探区的变速成图方法和装置
CN107656314B (zh) 道集记录处理方法、装置及计算机存储介质
CN102749649B (zh) 一种优选风化层速度的确定方法
CN114185085B (zh) 推覆构造区的速度建模方法、装置、设备及可读存储介质
CN113495296B (zh) 层析静校正量的确定方法、装置、设备及可读存储介质
CN111435172A (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