CN112444861A - 速度模型更新方法、计算机存储介质及计算机系统 - Google Patents

速度模型更新方法、计算机存储介质及计算机系统 Download PDF

Info

Publication number
CN112444861A
CN112444861A CN201910798065.8A CN201910798065A CN112444861A CN 112444861 A CN112444861 A CN 112444861A CN 201910798065 A CN201910798065 A CN 201910798065A CN 112444861 A CN112444861 A CN 112444861A
Authority
CN
China
Prior art keywords
time difference
horizon
residual time
layer
velocity model
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
CN201910798065.8A
Other languages
English (en)
Other versions
CN112444861B (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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910798065.8A priority Critical patent/CN112444861B/zh
Publication of CN112444861A publication Critical patent/CN112444861A/zh
Application granted granted Critical
Publication of CN112444861B publication Critical patent/CN112444861B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • 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/282—Application of seismic models, synthetic seismograms
    • 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/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures

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

速度模型更新方法、计算机存储介质及计算机系统
技术领域
本发明属于地震资料处理领域，尤其涉及一种高精度的速度模型更新方法，还涉及存储该速度模型更新方法的计算机存储介质及计算机系统。
背景技术
速度问题是勘探地震学的核心问题，对于任何偏移方法，速度模型的正确性都是决定构造成像质量的关键因素，尤其是叠前深度偏移对速度模型的依赖性更强、反应更敏感、要求更高。建立叠前深度偏移速度模型一般包括构造建模、初始速度建模、目标线偏移成像、模型的优化迭代和更新、最终叠前深度偏移体偏移五个过程，其中速度模型的优化迭代基于共成像点道集拾取剩余时差，并利用网格层析反演技术更新速度模型，因而它是建立叠前深度偏移速度模型的关键和核心步骤。
在某些地区多次波在共成像点道集上与一次反射波混在一起，多次波压制方法无法将多次波压制干净，仅仅基于数据驱动拾取剩余时差，往往包含多次波的剩余时差，基于多次波的剩余时差所建立的矩阵方程，其反演出来的速度误差量往往是错误的。
发明内容
本发明要解决的技术问题是：常规的自动拾取剩余时差，基于数据驱动进行相干扫描，通过曲率、相干性、剩余深度差、拾取的深度范围等参数控制要拾取的同相轴。如果共成像点道集的多次波去除不干净，多次波被保留在数据中，过多的多次波落在自动拾取的范围内，那么多次波的剩余时差就会被拾取，基于多次波的剩余时差建立的矩阵方程反演出来的速度更新往往是错误的。
为解决上述技术问题，本发明提供了一种速度模型更新方法、计算机存储介质及计算机系统。
根据本发明的一个方面，提供了一种速度模型更新方法，其包括以下步骤：
确定地质层位；
根据所述地质层位、共成像点道集和待更新的速度模型对应的深度域偏移剖面，确定多次波发育的深度范围；
根据所述多次波发育的深度范围以及待更新的速度模型对应的深度域偏移剖面，获得与所述地质层位对应的深度域层位信息；
将所述深度域层位信息导入所述共成像点道集中，针对位于地质层位浅层的上层及地质层位深层的下层的地层分别自动拾取出第一剩余时差和第二剩余时差，针对位于地质层位浅层以下至地质层位深层以上的地层逐层拾取出第三剩余时差；
根据所述第一剩余时差、第二剩余时差和第三剩余时差，建立三维层析方程并进行网格层析反演；
判断网格层析反演后的速度模型是否合格，若是，则将网格层析反演后的速度模型作为更新后的速度模型输出，若否，则将网格层析反演后的作为所述待更新的速度模型进行更新。
优选的是，根据测井曲线和地质分层数据确定所述地质层位。
优选的是，在确定多次波发育的深度范围之前，所述方法还包括：对所述共成像点道集进行预处理。
优选的是，所述预处理包括对共成像点道集进行去噪处理和/或去除多次波处理。
优选的是，通过判断所述共成像点道集是否拉平来判断所述网格层析反演后的速度模型是否合格，其中，在共成像点道集拉平时所述网格层析反演后的速度模型合格，在共成像点道集未拉平时所述网格层析反演后的速度模型不合格。
优选的是，所述网格层析反演后的速度模型不合格时，通过图偏移方调整不合格的速度模型的深度域层位信息，并基于调整后的深度域层位信息对不合格的速度模型进行更新。
优选的是，通过数据驱动来自动拾取出所述第一剩余时差和所述第二剩余时差。
优选的是，所述逐层拾取出第三剩余时差包括：
在所述地质层位浅层以下至地质层位深层以上的每个层位上分别仅设置一个时窗；
沿层位逐一拾取每个时窗内的剩余时差，并将该剩余时差作为所述第三剩余时差。
根据本发明的另一个方面，提供了一种计算机存储介质，其中存储有计算机程序，当所述计算机程序被一个或多个处理器执行时执行上述速度模型更新方法。
根据本发明的再一个方面，提供了一种计算机系统，其包括处理器和存储有计算机程序的存储器，所述处理器在执行所述计算机程序时执行上述速度模型更新方法。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
应用本发明的速度模型更新方法，在多次波发育区域引入层位信息，强制性沿着层位追踪剩余时差，能够提高速度反演的精度。具体地，本发明属于避开多次波的混合型剩余时差拾取方法，针对没有多次波发育的范围基于数据驱动自动拾取，而针对有多次波发育的区域强制沿地质层位追踪剩余时差，通过混合型的层约束拾取方法，能够避免因多次波而引起的速度反演偏差，同时不会降低多次波没有发育区域速度反演精度，提高网格层析反演的准确性。
附图说明
通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
图1示出了本发明实施例一速度模型更新方法的流程图；
图2示出了本发明实施例二速度模型更新方法的流程图；
图3示出了剩余曲率与共成像点道集叠合图；
图4A示出了网格层析更新前的速度模型示意图；
图4B示出了网格层析更新后的速度模型示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
在现有技术中，常规的自动拾取剩余时差，基于数据驱动进行相干扫描，通过曲率、相干性、剩余深度差、拾取的深度范围等参数控制要拾取的同相轴。如果共成像点道集的多次波去除不干净，多次波被保留在数据中，过多的多次波落在自动拾取的范围内，那么多次波的剩余时差就会被拾取，基于多次波的剩余时差建立的矩阵方程反演出来的速度更新往往是错误的。
为了解决上述技术问题，本发明实施例提供了一种速度模型更新方法。
实施例一
图1示出了本发明实施例一速度模型更新方法的流程图。如图1所示，本发明实施例的速度模型更新方法主要包括步骤S101至步骤S108。
在步骤S101中，确定地质层位。
具体地，根据测井曲线和地质分层数据确定所述地质层位。
在步骤S102中，根据所述地质层位、共成像点道集和待更新的速度模型对应的深度域偏移剖面，确定多次波发育的深度范围。
在步骤S103中，根据所述多次波发育的深度范围以及待更新的速度模型对应的深度域偏移剖面，获得与所述地质层位对应的深度域层位信息。
在步骤S104中，将所述深度域层位信息导入所述共成像点道集中，针对位于地质层位浅层的上层及地质层位深层的下层的地层分别自动拾取出第一剩余时差和第二剩余时差，针对位于地质层位浅层以下至地质层位深层以上的地层逐层拾取出第三剩余时差。
具体地，针对位于地质层位浅层的上层及地质层位深层的下层的地层，按照传统方法进行自动拾取。即针对位于地质层位浅层的上层的地层自动拾取出第一剩余时差，针对位于地质层位深层的下层的地层自动拾取出第二剩余时差。在本发明一优选的实施例中，通过数据驱动来自动拾取出所述第一剩余时差和所述第二剩余时差。
针对位于地质层位浅层以下至地质层位深层以上的地层，逐层拾取出第三剩余时差，具体逐层拾取过程为：在所述地质层位浅层以下至地质层位深层以上的每个层位上分别仅设置一个时窗；沿层位逐一拾取每个时窗内的剩余时差，并将该剩余时差作为所述第三剩余时差。
值得注意的是，针对位于地质层位浅层以下至地质层位深层以上的每个层位，对每个层位设置且仅设置一个时窗，设置的时窗的规格不重要。针对此处的每个层位(即，位于地质层位浅层以下至地质层位深层以上的每个层位)，强制性沿该地质层位拾取剩余时差。这里，强制拾取指的是，此处的每个地质层位都要开时窗，且必须要经过每一个时窗拾取剩余时差。
在步骤S105中，根据所述第一剩余时差、第二剩余时差和第三剩余时差，建立三维层析方程并进行网格层析反演。
在步骤S106中，判断网格层析反演后的速度模型是否合格。
具体地，通过判断所述共成像点道集是否拉平来判断所述网格层析反演后的速度模型是否合格。在共成像点道集拉平时所述网格层析反演后的速度模型合格，在共成像点道集未拉平时所述网格层析反演后的速度模型不合格。
在步骤S107中，在判断出网格层析反演后的速度模型合格的情况下，将网格层析反演后的速度模型作为更新后的速度模型输出。
在步骤S108中，在判断出网格层析反演后的速度模型不合格的情况下，将网格层析反演后的作为所述待更新的速度模型进行更新。
具体地，在判断出网格层析反演后的速度模型不合格的情况下，首先通过图偏移方调整不合格的速度模型的深度域层位信息，然后基于调整后的深度域层位信息对不合格的速度模型进行更新，即将调整后的深度域层位信息带入步骤S104中，循环执行步骤S104至步骤S108，以对不合格的速度模型进行更新，直到更新后的速度模型合格时为止。
在本发明一优选的实施例中，在步骤S102之前，还包括对所述共成像点道集进行预处理。具体地，所述预处理包括对共成像点道集进行去噪处理和/或去除多次波处理。
应用本发明实施例的速度模型更新方法，在多次波发育区域引入层位信息，强制性沿着层位追踪剩余时差，能够提高速度反演的精度。具体地，本发明实施例提供了一种避开多次波的混合型剩余时差拾取方法，通过人工交互的方式，针对没有多次波发育的范围基于数据驱动自动拾取，而针对有多次波发育的区域强制沿地质层位追踪剩余时差，通过混合型的层约束拾取方法，能够避免因多次波而引起的速度反演偏差，同时不会降低多次波没有发育区域速度反演精度，提高网格层析反演的准确性。
实施例二
图2示出了本发明实施例二速度模型更新方法的流程图。如图2所示，本发明实施例的速度模型更新过程具体包括以下步骤。
步骤一、根据测井曲线信息和地质分层数据确定地质层位。
步骤二、对共成像点深度道集进行去噪、去除多次波等处理，提高共成像点深度道集的信噪比。
步骤三、结合步骤一中地质层位信息、步骤二中优化后的共成像点深度道集以及深度域偏移剖面识别残存多次波发育的深度范围。
步骤四、根据步骤三中获得多次波发育深度信息，在现有速度模型对应的深度偏移剖面上拾取几套地质层位，最浅的层位对应多次波发育的起始深度附近，最深的层位对应多次波发育的终止深度附近。
步骤五、将步骤四中解释的地质层位信息导入到共成像点深度道集中。
步骤六、对步骤五中的共成像点深度道集自动拾取剩余时差，其中地质层位浅层以上自动拾取，没有层位约束；地质层位浅层至地质层位深层沿每个层位开一个时窗，强制性沿地质层位拾取剩余时差；地质层位深层向下地层自动拾取，没有层位约束。
步骤七、基于步骤六中拾取的剩余时差，建立三维层析方程，进行网格层析反演，得到新一轮的速度模型。
步骤八、用步骤七中得到的速度模型进行目标线叠前深度偏移，判断共成像点道集是否拉平，道集拉平时停止速度模型迭代更新，不满足要求时应用图偏移技术计算与步骤七得到的新的速度模型相对应的新的深度域层位信息，然后迭代步骤五、六，基于新的层位追踪剩余时差。
下面通过实际数据来说明本发明实施例的效果。网格层析反演获得的速度模型避开了由于多次波剩余曲率而带来的速度偏差。本示例中，选取新疆地震资料，4.5km以下发育多次波，尤其是奥陶系内幕，层间多次波与一次反射波混在一起，为了保证大套地层速度的准确性，多次波发育区域强制性沿层位追踪剩余曲率。图3为剩余曲率与共成像点道集叠合图，图4A和图4B分别为网格层析反演更新前和更新后的速度模型。从图中可以看出，应用本发明实施例的速度模型更新方法，能够有效地避免由于多次波同相轴下拉而反演的低速异常。浅层拾取较密，网格层析反演精度高，从而保证了浅层速度的准确性。
应用本发明实施例的速度模型更新方法，在多次波发育区域引入层位信息，强制性沿着层位追踪剩余时差，能够提高速度反演的精度。具体地，本发明实施例提供了一种避开多次波的混合型剩余时差拾取方法，通过人工交互的方式，针对没有多次波发育的范围基于数据驱动自动拾取，而针对有多次波发育的区域强制沿地质层位追踪剩余时差，通过混合型的层约束拾取方法，能够避免因多次波而引起的速度反演偏差，同时不会降低多次波没有发育区域速度反演精度，提高网格层析反演的准确性。
实施例三
本实施例提供了一种计算机存储介质，其中存储有计算机程序，当所述计算机程序被一个或多个处理器执行时执行如实施例一或实施例二中所述的速度模型更新方法。
本实施例还提供了一种计算机系统，该计算机系统包括处理器和存储有计算机程序的存储器，所述处理器在执行所述计算机程序时执行如实施例一或实施例二所述的速度模型更新方法。
这里，实施例一涉及的速度模型更新方法包括以下步骤：
确定地质层位；
根据所述地质层位、共成像点道集和待更新的速度模型对应的深度域偏移剖面，确定多次波发育的深度范围；
根据所述多次波发育的深度范围以及待更新的速度模型对应的深度域偏移剖面，获得与所述地质层位对应的深度域层位信息；
将所述深度域层位信息导入所述共成像点道集中，针对位于地质层位浅层的上层及地质层位深层的下层的地层分别自动拾取出第一剩余时差和第二剩余时差，针对位于地质层位浅层以下至地质层位深层以上的地层逐层拾取出第三剩余时差；
根据所述第一剩余时差、第二剩余时差和第三剩余时差，建立三维层析方程并进行网格层析反演；
判断网格层析反演后的速度模型是否合格，若是，则将网格层析反演后的速度模型作为更新后的速度模型输出，若否，则将网格层析反演后的作为所述待更新的速度模型进行更新。
实施例二涉及的速度模型更新方法包括以下步骤：
步骤一、根据测井曲线信息和地质分层数据确定地质层位。
步骤二、对共成像点深度道集进行去噪、去除多次波等处理，提高共成像点深度道集的信噪比。
步骤三、结合步骤一中地质层位信息、步骤二中优化后的共成像点深度道集以及深度域偏移剖面识别残存多次波发育的深度范围。
步骤四、根据步骤三中获得多次波发育深度信息，在现有速度模型对应的深度偏移剖面上拾取几套地质层位，最浅的层位对应多次波发育的起始深度附近，最深的层位对应多次波发育的终止深度附近。
步骤五、将步骤四中解释的地质层位信息导入到共成像点深度道集中。
步骤六、对步骤五中的共成像点深度道集自动拾取剩余时差，其中地质层位浅层以上自动拾取，没有层位约束；地质层位浅层至地质层位深层沿每个层位开一个时窗，强制性沿地质层位拾取剩余时差；地质层位深层向下地层自动拾取，没有层位约束。
步骤七、基于步骤六中拾取的剩余时差，建立三维层析方程，进行网格层析反演，得到新一轮的速度模型。
步骤八、用步骤七中得到的速度模型进行目标线叠前深度偏移，判断共成像点道集是否拉平，道集拉平时停止速度模型迭代更新，不满足要求时应用图偏移技术计算与步骤七得到的新的速度模型相对应的新的深度域层位信息，然后迭代步骤五、六，基于新的层位追踪剩余时差。
应用本发明实施例，在多次波发育区域引入层位信息，强制性沿着层位追踪剩余时差，能够提高速度反演的精度。具体地，本发明实施例提供了一种避开多次波的混合型剩余时差拾取方法，通过人工交互的方式，针对没有多次波发育的范围基于数据驱动自动拾取，而针对有多次波发育的区域强制沿地质层位追踪剩余时差，通过混合型的层约束拾取方法，能够避免因多次波而引起的速度反演偏差，同时不会降低多次波没有发育区域速度反演精度，提高网格层析反演的准确性。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (10)

1.一种速度模型更新方法，其特征在于，包括以下步骤：
确定地质层位；
根据所述地质层位、共成像点道集和待更新的速度模型对应的深度域偏移剖面，确定多次波发育的深度范围；
根据所述多次波发育的深度范围以及待更新的速度模型对应的深度域偏移剖面，获得与所述地质层位对应的深度域层位信息；
将所述深度域层位信息导入所述共成像点道集中，针对位于地质层位浅层的上层及地质层位深层的下层的地层分别自动拾取出第一剩余时差和第二剩余时差，针对位于地质层位浅层以下至地质层位深层以上的地层逐层拾取出第三剩余时差；
根据所述第一剩余时差、第二剩余时差和第三剩余时差，建立三维层析方程并进行网格层析反演；
判断网格层析反演后的速度模型是否合格，若是，则将网格层析反演后的速度模型作为更新后的速度模型输出，若否，则将网格层析反演后的作为所述待更新的速度模型进行更新。
2.根据权利要求1所述的方法，其特征在于，根据测井曲线和地质分层数据确定所述地质层位。
3.根据权利要求1所述的方法，其特征在于，在确定多次波发育的深度范围之前，所述方法还包括：对所述共成像点道集进行预处理。
4.根据权利要求3所述的方法，其特征在于，所述预处理包括对共成像点道集进行去噪处理和/或去除多次波处理。
5.根据权利要求1所述的方法，其特征在于，通过判断所述共成像点道集是否拉平来判断所述网格层析反演后的速度模型是否合格，其中，在共成像点道集拉平时所述网格层析反演后的速度模型合格，在共成像点道集未拉平时所述网格层析反演后的速度模型不合格。
6.根据权利要求5所述的方法，其特征在于，所述网格层析反演后的速度模型不合格时，通过图偏移方调整不合格的速度模型的深度域层位信息，并基于调整后的深度域层位信息对不合格的速度模型进行更新。
7.根据权利要求1所述的方法，其特征在于，通过数据驱动来自动拾取出所述第一剩余时差和所述第二剩余时差。
8.根据权利要求1或7所述的方法，其特征在于，所述逐层拾取出第三剩余时差包括：
在所述地质层位浅层以下至地质层位深层以上的每个层位上分别仅设置一个时窗；
沿层位逐一拾取每个时窗内的剩余时差，并将该剩余时差作为所述第三剩余时差。
9.一种计算机存储介质，其特征在于，其中存储有计算机程序，当所述计算机程序被一个或多个处理器执行时，执行如权利要求1-8中任一项所述的速度模型更新方法。
10.一种计算机系统，其特征在于，包括处理器和存储有计算机程序的存储器，所述处理器在执行所述计算机程序时执行如权利要求1-8中任一项所述的速度模型更新方法。
CN201910798065.8A 2019-08-27 2019-08-27 速度模型更新方法、计算机存储介质及计算机系统 Active CN112444861B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910798065.8A CN112444861B (zh) 2019-08-27 2019-08-27 速度模型更新方法、计算机存储介质及计算机系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910798065.8A CN112444861B (zh) 2019-08-27 2019-08-27 速度模型更新方法、计算机存储介质及计算机系统

Publications (2)

Publication Number Publication Date
CN112444861A true CN112444861A (zh) 2021-03-05
CN112444861B CN112444861B (zh) 2024-06-18

Family

ID=74742104

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910798065.8A Active CN112444861B (zh) 2019-08-27 2019-08-27 速度模型更新方法、计算机存储介质及计算机系统

Country Status (1)

Country Link
CN (1) CN112444861B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020042678A1 (en) * 2000-08-07 2002-04-11 Dimitri Bevc Typing picks to horizons in migration velocity analysis
US20020049540A1 (en) * 2000-08-07 2002-04-25 Dimitri Bevc Velocity analysis using angle-domain common image gathers
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model
US20120241166A1 (en) * 2011-03-22 2012-09-27 Hongchuan Sun Residual Moveout Estimation Through Least Squares Inversion
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN102841375A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种复杂条件下基于角度域共成像点道集的层析速度反演方法
CN103675906A (zh) * 2012-09-19 2014-03-26 中国石油化工股份有限公司 用于提高复杂碎屑岩储层预测精度的多控储层预测方法
CN104570088A (zh) * 2013-10-29 2015-04-29 中国石油天然气股份有限公司 一种在共偏移距共方位角域进行多方位层析速度建模方法
CN104932021A (zh) * 2014-03-18 2015-09-23 中国石油化工股份有限公司 基于反向射线追踪的约束层析速度建模方法
CN105093281A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 一种反演框架下的地震多波建模方法
US20160109589A1 (en) * 2014-10-20 2016-04-21 Jonathan Liu Velocity tomography using property scans
US20160299243A1 (en) * 2014-10-24 2016-10-13 Landmark Graphics Corporation Structure tensor constrained tomographic velocity analysis
CN109116413A (zh) * 2018-07-30 2019-01-01 中国石油化工股份有限公司 成像域立体层析速度反演方法
CN109444959A (zh) * 2018-11-01 2019-03-08 科吉思石油技术咨询(北京)有限公司 全频高精度层速度场建立方法
CN109581521A (zh) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 Tti各向异性的局部层析方法及系统

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020042678A1 (en) * 2000-08-07 2002-04-11 Dimitri Bevc Typing picks to horizons in migration velocity analysis
US20020049540A1 (en) * 2000-08-07 2002-04-25 Dimitri Bevc Velocity analysis using angle-domain common image gathers
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model
US20120241166A1 (en) * 2011-03-22 2012-09-27 Hongchuan Sun Residual Moveout Estimation Through Least Squares Inversion
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN102841375A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种复杂条件下基于角度域共成像点道集的层析速度反演方法
CN103675906A (zh) * 2012-09-19 2014-03-26 中国石油化工股份有限公司 用于提高复杂碎屑岩储层预测精度的多控储层预测方法
CN104570088A (zh) * 2013-10-29 2015-04-29 中国石油天然气股份有限公司 一种在共偏移距共方位角域进行多方位层析速度建模方法
CN104932021A (zh) * 2014-03-18 2015-09-23 中国石油化工股份有限公司 基于反向射线追踪的约束层析速度建模方法
CN105093281A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 一种反演框架下的地震多波建模方法
US20160109589A1 (en) * 2014-10-20 2016-04-21 Jonathan Liu Velocity tomography using property scans
US20160299243A1 (en) * 2014-10-24 2016-10-13 Landmark Graphics Corporation Structure tensor constrained tomographic velocity analysis
CN109581521A (zh) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 Tti各向异性的局部层析方法及系统
CN109116413A (zh) * 2018-07-30 2019-01-01 中国石油化工股份有限公司 成像域立体层析速度反演方法
CN109444959A (zh) * 2018-11-01 2019-03-08 科吉思石油技术咨询(北京)有限公司 全频高精度层速度场建立方法

Also Published As

Publication number Publication date
CN112444861B (zh) 2024-06-18

Similar Documents

Publication Publication Date Title
US20100135115A1 (en) Multiple anisotropic parameter inversion for a tti earth model
CN102944896B (zh) 表层调查数据的模型法静校正方法
CN103765245A (zh) 混合式确定性－地质统计地球模型
US9366771B2 (en) System and method for migration velocity modeling
EP2631678A2 (en) Feature detection in seismic volumes
US11119233B2 (en) Method for estimating elastic parameters of subsoil
CN113945978A (zh) 一种深度域地震数据优化方法及装置
CN112444861A (zh) 速度模型更新方法、计算机存储介质及计算机系统
CN112415601B (zh) 表层品质因子q值的确定方法及装置
Zhang et al. Automatic first-arrival picking workflow by global path tracing
Santos et al. Fast estimation of common-reflection-surface parameters using local slopes
CN107193047A (zh) 非地表一致性时变动态剩余静校正方法和装置
CN113031072A (zh) 虚同相轴层间的多次波压制方法、装置及设备
CN113820745A (zh) 地震速度建模方法、装置、电子设备及介质
CN113219534B (zh) 一种叠前深度偏移速度质控方法、装置、介质及电子设备
CN113219533B (zh) 一种叠前时间偏移速度建模方法、装置、介质及电子设备
CN112505780B (zh) 地层深度数据的校正方法及装置
CN105549075B (zh) 求取浅层速度场的方法及装置
CN112415580B (zh) 消除速度模型突变界面的方法、叠前深度偏移处理方法
CN113267810B (zh) 地震勘探全深度速度建模方法及装置
CN113495303B (zh) 剥蚀量恢复的方法和装置
CN114185085B (zh) 推覆构造区的速度建模方法、装置、设备及可读存储介质
CN113466936A (zh) 一种断层阴影区crp道集获取方法、装置、设备及存储介质
CN111505711A (zh) 3d折射波波场偏移成像方法
CN113009577A (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