CN108139499A - Q-补偿的全波场反演 - Google Patents

Q-补偿的全波场反演 Download PDF

Info

Publication number
CN108139499A
CN108139499A CN201680057240.2A CN201680057240A CN108139499A CN 108139499 A CN108139499 A CN 108139499A CN 201680057240 A CN201680057240 A CN 201680057240A CN 108139499 A CN108139499 A CN 108139499A
Authority
CN
China
Prior art keywords
model
sound
wave field
full wave
viscous
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
CN201680057240.2A
Other languages
English (en)
Other versions
CN108139499B (zh
Inventor
孙宏川
E·G·韦尔德尔慕斯
J·刘
R·班赛尔
S·K·拉扎拉托斯
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.)
ExxonMobil Upstream Research Co
Original Assignee
Exxon Production Research Co
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 Exxon Production Research Co filed Critical Exxon Production Research Co
Publication of CN108139499A publication Critical patent/CN108139499A/zh
Application granted granted Critical
Publication of CN108139499B publication Critical patent/CN108139499B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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. 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/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/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. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/58—Media-related
    • G01V2210/584—Attenuation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

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)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

一种方法，包括：获得由声全波场反演过程生成的速度模型；利用计算机通过对地下区域的处理后的地震数据应用伪Q迁移来生成变量Q模型，其中所述速度模型被用作所述伪Q迁移中的引导约束；以及利用计算机生成最终地下速度模型，所述最终地下速度模型通过执行粘声全波场反演过程来恢复由地下区域中的气体异常引起的振幅衰减，其中所述变量Q模型在所述粘声全波场反演处理中是固定的。

Description

Q-补偿的全波场反演
对相关申请的交叉引用
本申请要求于2015年10月2日提交的标题为Q-COMPENSATED FULL WAVEFIELDINVERSION的美国临时专利申请62/236,190的权益，该申请的全部内容通过引用被结合于此。
技术领域
本文描述的示例性实施例涉及地球物理勘探领域，并且更特别地涉及地球物理数据处理。具体而言，本文描述的实施例涉及通过补偿地下异常来改善地震数据的准确度的方法。
背景技术
本部分旨在介绍可以与本发明的示例性实施例相关联的现有技术的各个方面。这种讨论被认为有助于提供框架以便于更好地理解本发明的特定方面。因此，应当理解的是，本部分应该以此来阅读，而不必作为对现有技术的承认。
地震勘探的重要目标是准确地对通常称为反射体的地下结构进行成像。通过在地震勘测的执行期间获得原始地震数据有利于地震勘探。在地震勘测期间，地震能量通过例如受控爆炸在地面生成，并传递到地球。地震波由地下结构反射，并被许多称为地震检波器的传感器接收。由地震检波器接收到的地震数据被处理以便创建地下环境的准确映射。处理后的数据然后被检查，以识别可能含有油气(hydrocarbon)的地质构造。
以相对垂直方向传送到地球中的地震能量最有可能被反射器反射。这种能量提供关于地下结构有意义的信息。但是，地震能量可能由于在地下环境中常规发生的声阻抗中的异常而非期望地扩散。地震勘测期间地震能量的扩散可能导致地下特征在结果产生的地震数据中被错误地表示。
声阻抗是地震能量穿过地下环境的特定部分的容易程度的测量。本领域普通技术人员将认识到的是，声阻抗可以被定义为密度和地震速度的乘积。声阻抗通常由符号Z表示。
地震波在地下环境中行进时会因各种原因衰减。质量度量(有时称为质量因子)Q通常用于表示地下岩层(underground formation)的衰减特性。通常，Q与地震信号衰减成反比，并且可以从零值到无穷大。更具体而言，Q是无量纲质量因子，它是波的峰值能量与耗散能量的比值。随着波行进，由于球面发散和吸收，它们会随着距离和时间而失去能量。当恢复地震振幅以进行流体和岩性解释(诸如振幅与偏移量(AVO)分析)时，必须考虑到这种能量损失。具有相对高Q值的结构倾向于几乎没有衰减地传送地震波。倾向于更大程度衰减地震能量的结构具有较低的Q值。
与地下结构相关联的Q值用于在数学上改变地震数据值以更准确地表示地下环境中的结构。这个过程被本领域普通技术人员称为“Q迁移”。在Q迁移期间，表示地震能量行进通过具有相对低Q值的地下结构的地震数据值可以在频谱中被放大和扩大到比表示地震能量行进通过具有相对高Q值的地下结构的数据值更大的程度。改变与低Q值相关联的数据的振幅和相位考虑了当地震能量行进通过具有相对低Q值的结构时发生的较大信号衰减。
FWI是一种偏微分方程约束优化方法，它迭代地最小化测得的波场和计算出的波场之间的失配范数。地震FWI涉及多次迭代，并且单次迭代可以涉及以下计算：(1)正向方程的解，(2)伴随方程的解，以及(3)这些正向解和伴随解的卷积以产生成本函数的梯度。注意的是，对于二阶优化方法，诸如高斯牛顿法(Gauss-Newton)，还需要(4)扰动正向方程的解。例如，对于这种情况的更稳健的数学证明可以例如在美国专利公开2013/0238246中找到，其全部内容通过引用被结合于此。
用于在衰减声介质中模拟波的线性粘声波方程的常规一阶形式是：
具有压力p，速度v和存储器变量ml的适当的初始和边界条件。注意的是：
κ＝非松弛体块模量(limω→∞κ(ω)→κ)，
ρ＝质量密度，
v＝速度(3D空间中的v＝{vxvyvz}T)，
p＝压力,
ml＝机制l的存储器变量，
sp＝压力源，
sv＝速度源,
和
其中松弛参数τ∈l和τσl可以由方程(2)针对给定质量因子分布确定。
注意的是，连续的标量变量用斜体字符表示，并且向量和矩阵贯穿本文档用粗体非斜体字符表示。
其中
Q＝质量因子，
τ∈l＝SLS模型中机制l的应变松弛时间，
τσl＝SLS模型中机制l的应力松弛时间，
x＝空间坐标，
ω＝频率，
L＝SLS模型中使用的松弛机制的数量。
概念上，质量因子Q表示介质中储存能量与耗散能量的比率。确定应变和应力松弛时间以最佳地适合频带上的期望质量因子分布。
基于计算目标函数关于参数的梯度的全波场反演(FWI)方法通常通过使用伴随方法来高效地实现，伴随方法已被证明优于其它相关方法，诸如直接灵敏度分析、有限差异或复变量方法。
常规粘声系统的连续伴随(方程(1))是：
其中
以及
和分别是关于压力和速度的目标函数的导数。
用于地球物理的常用迭代反演方法是成本函数优化。成本函数优化涉及相对于模型θ的成本函数的值的迭代最小化或最大化。也称为目标函数的成本函数是模拟数据和观测数据之间不匹配程度的测量。模拟(模拟数据)通过首先利用适当的数值方法(诸如有限差分或有限元方法)将源信号在介质中的物理控制传播离散化，并且然后在计算机上使用当前的地球物理属性模型计算数值解来进行。
以下总结了FWI的局部成本函数优化过程：(1)选择起始模型；(2)计算搜索方向S(θ)；以及(3)搜索作为模型在搜索方向上的扰动的更新模型。
通过使用新的更新模型作为寻找另一个搜索方向的起始模型来迭代成本函数优化过程，然后将其用于扰动模型以便更好地解释观察到的数据。这个过程继续，一直到发现能够令人满意地解释观察到的数据的更新模型。常用的局部成本函数优化方法包括梯度搜索法、共轭梯度法、准牛顿法、高斯牛顿法和牛顿法。
声学近似中地震数据的局部成本函数优化是常见的地球物理反演任务，并且一般地说明其它类型的地球物理反演。当在声学近似中反演地震数据时，成本函数可以写为：
其中
θ＝N个参数的向量，(θ1,θ2,…θN)描述地下模型，
g＝集合索引，
wg＝集合g的源函数，它是空间坐标和时间的函数，对于点源，这是空间坐标的增量δ(delta)函数，
Ng＝集合的数量，
r＝集合内的接收器索引，
Nr＝集合中的接收器的数量，
t＝轨迹内的时间样本索引，
Nt＝时间样本的数量，
W＝范数函数(最小化函数，例如用于最小二乘函数(x)＝x2)，
ψcalc＝从模型θ的计算出的地震数据，
ψobs＝测得的地震数据(压力、应力、速度和/或加速度)。
集合(gathers)(来自共享共同几何形状的多个传感器的数据)可以是可以在地震正向建模程序的一次运行中被模拟的任何类型的集合(共同中点、共同源、共同偏移、共同接收器等)。通常，集合对应于地震炮(seismic shot)，但是地震炮可能比点源更为普遍。对于点源，集合索引g对应于各个点源的位置。这种广义的源数据ψobs可以或者在现场采集或者可以从使用点源采集到的数据进行合成。另一方面，计算出的数据ψcalc通常可以在正向建模时通过使用广义源函数直接计算。
FWI尝试更新离散模型θ，使得最小。这可以通过局部成本函数优化来完成，该优化如下更新给定模型θ(k)：
θ(i+1)＝θ(i)+γ(i)S(θ(i)), (5)
其中i是迭代次数，γ是模型更新的标量步长，并且S(θ)是搜索方向。对于最陡下降，这是关于模型参数取的失配函数的梯度的负值。在这种情况下，模型扰动或通过其更新模型的值通过将目标函数的梯度与步长γ相乘来计算，该步长γ必须重复计算。对于二阶优化技术，梯度由Hessian(目标函数相对于模型参数的二阶导数)进行缩放。的计算需要计算关于N个模型参数中的每一个的导数。在地球物理问题中，N通常非常大(超过一百万)，并且如果必须为每个单独的模型参数执行该计算，那么该计算可能非常耗时。幸运的是，伴随方法可以用于为所有模型参数一次高效地执行这种计算(Tarantola，1984)。
FWI通过地震数据的迭代反演生成用于叠前深度迁移和地质解释的高分辨率属性模型(Tarantola，1984；Pratt等人，1998)。随着计算机资源的不断增加和最近的技术进步，FWI能够处理大得多的数据集，并逐渐在3D实际数据应用中变得经济适用。但是，在常规的FWI中，被反演的数据通常被当作它们是在声地下介质(acoustic subsurface medium)中收集的，这与地球始终衰减的事实不一致。当介质中存在气体云时，控制衰减效应的质量因子(Q)在地震波传播中起重要作用，从而导致相位失真、振幅变小以及频率降低。因此，常规的声学FWI不能补偿Q效应，并且不能恢复气体异常下的振幅和带宽损失。
另一方面，粘声FWI使用波场传播中的中速和Q值。因此，当波前继续进行时，Q效应自然得到补偿。在一些情况下，如果浅层气体异常覆盖储层，那么严重的Q效应将屏蔽掉信号并且在声学FWI实现中引起跳周期问题。因此，粘声FWI算法和准确的Q模型是非常优选的。
但是，Q值并不容易确定。在许多方法中，基于射线的折射或反射Q层析成像已经被广泛研究。但是，在现场数据应用中，Q层析成像是繁琐的过程，并且反演很大程度上取决于如何将衰减信号与其未衰减的对等信号分开。近年来，已经提出诸如FWI的基于波的反演算法来反演Q值。理论上，这种基于波的方法更准确。但是，速度和Q反演可能以不同的步伐收敛，并且速度和Q梯度之间可能存在严重的能量泄漏，使得反演结果不可靠。
Zhou等人(2014)描述了如何使用声学FWI进行速度反演，以及然后如何使用FWI反演速度模型进行Q反演。但是，基于射线的Q层析成像非常耗时，并且它们没有进行真正的粘声波形反演。
Bai等人(2014)应用了粘声FWI进行速度反演，但是，他们也需要使用粘声FWI来反演Q模型。如业界普遍认为的那样，这种波形Q反演非常不稳定。速度反演中的误差可能容易泄漏到Q反演中，反之亦然。此外，这种方法很难做到面向目标的Q反演，并且Q的分辨率和幅度仍然是亟待解决的问题。
发明内容
一种方法，包括：获得由声全波场反演过程生成的速度模型；利用计算机通过对地下区域的处理后的地震数据应用伪Q迁移来生成变量Q模型，其中所述速度模型被用作所述伪Q迁移中的引导约束；以及利用计算机产生最终地下速度模型，所述最终地下速度模型通过执行粘声全波场反演过程来恢复由地下区域中的气体异常引起的振幅衰减，其中所述变量Q模型在所述粘声全波场反演过程中是固定的。
该方法还可以包括基于所述最终地下速度模型生成地下图像。
该方法还可以包括：生成处理后的地震数据，其中所述生成包括将基于声射线的叠前深度迁移应用于所述速度模型并且输出共同图像集合。
在该方法中，生成变量Q模型包括根据引导约束来平坦化所述共同图像集合。
在该方法中，所述引导约束定义所述速度模型内包含气体异常的区。
在该方法中，所述伪Q迁移仅应用于包含所述气体异常的区。
在该方法中，所述变量Q模型通过整个粘声全波场反演过程保持固定。
在该方法中，生成最终地下速度模型包括通过使粘声共同图像集合平坦化来应用伪Q迁移以构造另一个变量Q模型，并且由粘声全波场反演过程生成的速度模型被用作伪Q迁移中的引导约束。
该方法还可以包括迭代地重复以下操作直到达到预定的停止条件，所述操作包括：执行粘声全波场反演过程，然后从基于粘声射线的叠前深度迁移生成粘声共同图像集合，以及然后通过使粘声共同图像集合平坦化来生成另一个变量Q模型，其中由粘声全波场反演过程生成的速度模型被用作伪Q迁移中的引导约束。
该方法还可以包括进行地震勘测，其中至少一个源用于将声信号注入到地下并且至少一个接收器用于记录从地下特征反射的声信号。
该方法还可以包括使用最终地下速度模型和地下的图像来提取油气。
该方法还可以包括钻井以提取油气，其中所述井被布置在通过分析地下图像而确定的位置处。
在该方法中，所述引导约束由从声全波场反演过程反演的地质结构来引导。
附图说明
虽然本公开容许各种修改和替代形式，但是其具体示例实施例已经在附图中示出并且在本文中进行了详细描述。但是，应当理解的是，本文对具体示例实施例的描述不旨在将本公开限制于本文所公开的特定形式，而是相反，本公开要涵盖由所附权利要求限定的所有修改和等同物。还应当理解的是，附图不一定按比例绘制，而是将重点放在清楚地图示本发明的示例性实施例的原理上。此外，某些维度可能被夸大以帮助在视觉上传达这些原理。
图1是用于解释本发明的示例性实施例的操作的地下区域的图形表示。
图2图示了生成起动速度模型和起动Q模型的示例性方法。
图3图示了粘声FWI的示例性方法。
图4A图示了示例性速度模型。
图4B图示了与图4A中的气体云匹配的示例性Q模型。
图5A图示了通过声学正向建模生成的示例性共同偏移数据。
图5B图示了通过粘声正向建模生成的示例性共同偏移数据。
图6A图示了由声学FWI生成的示例性速度更新。
图6B图示了从具有固定的Q值的粘声FWI生成的示例性速度更新。
具体实施方式
本文描述了示例性实施例。但是，就以下描述特定于特定实施例而言，这仅仅是出于示例性目的，并且只是提供示例性实施例的描述。因此，本发明不限于下面描述的具体实施例，而是，它包括落入所附权利要求的真实精神和范围内的所有替代、修改和等同物。
本文描述的示例性实施例提供了全面的模型构建工作流程，其有效地补偿了FWI中的Q效应而不会遭受能量泄漏，并且能够生成具有很大改进的地下保真度的高分辨率属性曲线。
通过气体云层传播的地震波常常导致相位失真、振幅变小和频率降低。声波FWI不补偿这种Q效应，并且因此不能恢复气体异常下的振幅和带宽损失。本技术进步的非限制性实施例通过组合基于射线的Q模型构建与粘声FWI来补偿Q效应。在基于射线的Q模型构建中，伪Q迁移用于高效扫描所有可能的Q值以便获得最佳Q效应。与Q层析成像或全波场Q反演等其他Q估计方法相比，伪Q迁移可以在面向目标的局部气体区域内进行扫描，并且扫描过程既高效又稳定。
在粘声FWI中，Q效应在正向建模和伴随计算两者中都得以补偿。将Q补偿的FWI应用于复数合成数据揭露了气体区下面明显改善的结构，并且结果绝对可以有利于地质解释。只要Q效应被认为成问题，本技术进步通常就可以应用于任何现场数据。本技术进步最适用于其中地下存在强气体异常并且常规的声学模型构建工作流程不能恢复气体下面的潜在储层目标的振幅和带宽损失的数据集。
本技术进步的实施例提供了气体友好的模型构建工作流程，其可以包括：用于Q值的高效估计的基于射线的伪Q迁移；以及用于Q补偿速度更新的基于波的粘声FWI。
在粘声FWI中，Q值是固定的，并且只需要计算速度梯度。贯穿整个FWI速度反演固定Q值阻止了Q的能量泄漏转化为速度，这使得反演稳定并且反演的速度模型可靠。将Q补偿的FWI应用于具有气体云的合成和现场数据已表明气体下面改进的结构和连井，并且因此可以对地质解释产生积极影响。
图1是用于解释本技术进步的示例性实施例的操作的地下区域的图形表示。该图一般地由标号100表示。图100示出了水域101和沉积物区域102。本领域普通技术人员将认识到的是，水域101表现出通常由诸如999的大数字表示的非常高的Q值。因此，地震波以相对小的衰减行进通过水域101。沉积物区域102可以具有比水域101低得多的Q值。例如，沉积物区域102的Q值可以在大约150的范围内。
地下异常(诸如气帽等)可能具有非常低的Q值。在图100中，异常A 103表现出在大约20的范围内的Q值。类似地，异常B 104也表现出在大约20的范围内的Q值。异常A 103和异常B 104的低Q值导致与更深的地下结构对应的衰减的地震数据。作为说明，异常A 103负面影响在异常A 103下方延伸的异常衰减区域105中的地震数据的完整性。行进通过异常A103的任何地震能量在其返回到表面并被测量时将显著衰减。如果将异常A 103布置在油气的沉积物上方，那么可能隐藏可以识别油气的更深储层的存在的地震数据。这种现象在图100中由一系列反射器106、107、108和109示出。可能由显著衰减的地震数据表示的反射器106、107、108和109中的部分以虚线示出。与不可能由于异常A 103和异常B 104的存在而被显著衰减的反射器106、107、108和109的部分对应的数据在图1中以实线示出。根据本发明的示例性实施例的伪Q迁移的执行旨在恢复由于通过异常A 103和异常B 104而受到不利影响的地震能量的准确振幅、频率和相位数据。
图2图示了生成起始速度模型和起始Q模型的示例性方法，该起始速度模型和起始Q模型然后用作粘声FWI方法的输入(参见图3)。在步骤201中，生成初始速度模型V0 TOMO。初始速度模型可以通过基于射线的层析成像生成(Liu等人，2008)。但是，也可以使用其它技术。
在步骤203中，使用初始速度模型应用声学FWI。在步骤205中，从声学FWI生成高分辨率速度分布V1 FWI。任何FWI算法的关键可以被描述如下：使用起始地下物理属性模型，通过使用数值方案(例如，有限差分、有限元法等)求解波动方程来生成(即，建模或模拟)合成的地震数据。如本文所使用的术语速度模型或物理属性模型是指数字阵列，通常为3-D阵列，其中每个数字(其可以被称为模型参数)是单元格中的速度或另一个物理属性的值，其中出于计算目的，地下区域在概念上已被划分成离散单元格。将合成地震数据与现场地震数据进行比较，并使用两者之间的差计算误差或目标函数。使用目标函数，生成修改后的地下模型，用于模拟一组新的合成地震数据。这组新的合成地震数据与现场数据进行比较，以生成新的目标函数。重复该处理，直到目标函数被令人满意地最小化并且生成最终地下模型。使用全局或局部优化方法来最小化目标函数和更新地下模型。关于FWI的更多细节可以在Lee等人的美国专利公布2011/0194379中找到，其全部内容通过引用被结合于此。
在步骤207中，基于FWI反演速度模型V1 FWI，应用基于声射线的叠前深度迁移(PSDM)来计算其平坦度反映速度模型和Q模型的准确度的共同图像集合(common-image-gather)。几种叠前迁移方法可以用于执行步骤207，并且包括例如Kirchoff PSDM、单向波动方程迁移和反向时间迁移，其中的每一个都是本领域普通技术人员已知的。
在步骤209中，提供共同图像集合(CIG)作为伪Q迁移的输入。共同图像集合是共享一些共同几何属性(例如，共同偏移或共同角度)的地震轨迹的集合。
在步骤211中，从V1 FWI确定通常表示气体异常的低速度区，并将其提供为伪Q迁移的输入。在确定Q的后期阶段(参见步骤215)，这些低速度区被用作FWI引导约束。特别地，FWI引导约束描绘了可能的气体区域(目标气体区)，并且用在随后的伪Q迁移中以将伪Q迁移限制为仅可能包含气体的那些区域(如由V1 FWI中的低速度的区域所指示的)。引导约束由从声全波场反演过程反演的地质结构引导。
在步骤213中，获得Qo并将其提供作为伪Q迁移的输入。Qo可以是初始的齐次Q模型。
在步骤215中，将伪Q迁移应用于由FWI约束描绘的区域或区。伪Q迁移与Q层析成像不同。Q层析成像是非常繁琐的过程，其需要仔细准备不同的输入数据。此外，当共同图像角度的信噪比低时，Q层析成像是不稳定的。但是，伪Q迁移不仅高效，而且因为它与Q扫描相似，因此也是稳定的。下面是伪Q迁移理论的简要概述。在国际专利申请公开WO 2009/123790中可以找到伪Q迁移的更全面描述，该申请的全部内容通过引用被结合于此。
伪Q迁移的方面涉及构建Q积分表，该表可以用于恢复与迁移轨迹对应的数据的振幅、频率和相位值。Q积分数据是使用矩阵(相对于给定Q模型的Q积分的导数)和向量(Q模型的更新)的乘积来计算的。Q积分的导数可以部分地通过确定速度模型中的变化率来计算。在示例性实施例中，Q积分的导数有利地表示为稀疏矩阵。
导数矩阵独立于Q模型；因此，它可以被预先计算并存储，用于与随后开发的Q模型一起使用。然后计算每个特定图像点处以及针对其轨迹中的反射角度的Q积分的表。基本上，根据示例性实施例的伪Q迁移包括跟踪进(trace-in)和跟踪出(trace-out)操作，因此该计算实质上是1维处理操作。此外，伪Q迁移可以在目标定向中实现。只有预测到输入轨迹要受到低Q区(诸如异常A衰减区域105)影响(参见图1)时，输入轨迹才将为了振幅和相位恢复而被重新迁移。
作为示例，令c(x)为粘声速度场中的复数速度与频率变量ω的关系。在这些假设下，c(x)可以被表示如下：
其中c0是复数速度的声学部分，Q是表示衰减的质量因子，并且ω0是参考频率。
复数行进时间可以通过以下计算：
其中τ(x)是在声学介质c0中的行进时间，并且
方程(6)中的第一项包含迁移成像中的主要运动学信息，并且可以通过声学介质中的射线跟踪来计算。方程(6)中的第二项允许迁移补偿由于衰减而导致的振幅损失，并且方程(6)中的第三项允许迁移补偿由于分散而导致的相位失真。如方程(8)中所定义的，第二项和第三项都取决于T*，即沿着射线路径L的Q-1的积分。T*可以在与用于计算T的相同射线路径上计算。当Q被更新时，T保持不变并且T*的变化是：
而且，方程(9)可以被重写为矩阵形式：
ΔT*＝D*Δ(Q-1) (10)
在方程(10)中，D是关于Q-1(导数值)的T*的导数的矩阵。作为本文使用的术语的示例，D包含“基于速度模型的Q积分值的导数”，其中术语“Q积分值”由方程(8)表示，并且c0指速度模型。矩阵D可以被预先计算和存储，因为它不依赖于特定的Q模型。
伪Q迁移的示例性处理步骤可以陈述如下：(1)给出迁移的轨迹(共同图像集合或处理后的地震数据)，在迁移中使用的速度模型以及初始Q模型；(2)选择反射点并且估计那些点的反射角；(3)计算Q积分相对于Q模型的导数(d T*/dQ)并输出那些导数；以及(4)将Q积分的导数乘以初始Q模型以获得Q积分的表。
利用可用作初始Qo模型的简单齐次Q，图2中的过程将确定并输出优化变量Q模型Q1 Pseudo-Q(步骤217)，其通过最有效地减小图像集合曲率(即，平坦度)来判断。在伪Q优化的每次迭代中，射线路径矩阵是固定的，并且Q值被扰动。与任何基于反演的Q模型构建方法不同，伪Q迁移在合理的Q值范围内扫描可用的Q模型，并在FWI引导的目标区域内执行扫描(即，由FWI速度模型识别为是气体的区域)。因此，结果得到的Q1 Pseudo-Q没有局部最小值并且因此在粘声波场传播的后期阶段是优选的(参见图3)。
图2的共同特征是每个过程都可以是声学的。在图2中，整个过程有意保持与Q因子无关。这个声学过程有助于构建更可靠的起始Q模型，其从一开始就能够对声学模型进行改进。
获得图2中的方法的输出，并将其用作图3中的方法的第一次迭代的起始模型(步骤301和303)。在步骤305中，粘声FWI利用固定Q(固定并不意味着常量，因为Q1 Pseudo-Q是变量，而是固定是指在图3所示的方法的每次迭代i中使用相同的变量Q模型)来执行。如图3所示，图3使用了速度V1 FWI和Q模型Q1 Pseudo-Q，粘声FWI在其波场传播中自动包括Q，并能够有效地解释衰减效应。但是，Q模型在粘声FWI过程中固定，以避免Q反演与速度反演之间的能量泄漏。因此，只需要在粘声FWI过程中计算速度梯度。
在步骤307中，更新后的速度模型被输出用于在粘声PSDM中后续使用。
在步骤309中，获得Q模型Q1 Pseudo-Q用于在粘声PSDM中后续使用。
在步骤311中，使用来自步骤307的更新后的速度模型和来自步骤309的Q模型应用粘声PSDM。粘声Kirchoff迁移是生成粘声PSDM集合和堆叠的最广泛使用的方法。
在步骤313中，从粘声PSDM输出共同图像集合(CIG)。
类似于图2中的步骤211，步骤315包括使用FWI反演速度模型(来自步骤307)来确定气体异常(基于衰减速度的区域)，所述异常可用作伪Q反演的引导约束。
在步骤317中，获得Q模型Q1 Pseudo-Q用于在伪Q迁移中后续使用。
在步骤319，伪Q迁移稳定地扫描最佳Q值以进一步最小化图像集合曲率。通过扫描，该过程可以循环多个Q值并确定多个Q值中的哪一个是最优解。上文讨论了伪Q迁移，并且该讨论适用于这一步骤的执行。
图3中的工作流程迭代地更新速度和Q模型，并在可以由用户指定的预定迭代i＝N处停止。可替代地，可以使用其它停止标准。
注意的是，索引“i”是针对图3中的方法的大混合循环，而不是FWI步骤305的迭代次数。如本领域普通技术人员将理解的，FWI本身是迭代的过程。
图3的共同特征是每个过程可以是粘声的。在图3中，深度迁移和FWI两者都使用粘声引擎来保证波传播是一致的。
一旦完成对图3中整个过程的迭代(即，完成N次迭代)，图3的输出就是最终地下速度模型(或物理属性模型)和最终Q模型(步骤321和323)。最终物理属性地下模型可以用于生成用于解释地下的地下图像和/或管理油气勘探(步骤325)。将Q补偿的FWI应用于复数合成数据揭露了气体区下面明显改善的结构，并且结果绝对有利于地质解释。如本文所使用的，油气管理包括油气提取、油气生产、油气勘探、识别潜在的油气资源、识别井的位置、确定井注入和/或提取率、识别储层连接性、获取、处置和/放弃油气资源、回顾先前的油气管理决策、以及任何其它与油气有关的行为或活动。
以下描述本技术进步的应用的非限制性示例。
已经将图2和图3的方法应用于2-D合成拖缆数据。图4A图示了示例性速度模型。图4A是真正的地下模型，从中可以应用正向建模以生成合成数据。图4A中的速度模型根据深水油田地质情况进行了修改，在深水油田地质情况下，气体区位于水底附近并且使得常规深度成像和模型构建在下面的储层区域非常具有挑战性。
为了与低速度气体区401平行，如图4B所示，创建Q异常402。类似于图4A，图4B也在粘声建模中使用以计算合成数据。
图5A示出了使用图4A中的速度模型根据声学正向建模(共同偏移数据，其中偏移＝100m)生成的近偏移数据，而图5B示出了使用速度和Q异常两者根据粘声建模生成的类似数据。在圈出的区域501和502中，事件由于气体的吸收效应而衰减。与图5A相比，图5B中的唯一区别在于其中相位失真并且振幅衰减的气体区域。图5A和图5B之间的差异展示了Q衰减效应。由于Q，图5B明显地具有低的质量。
图6A示出了从声学FWI反演的速度更新(即，步骤203)，并且图6B示出了从具有固定Q值的粘声FWI反演的速度更新(即，步骤305)。在这两个FWI中，合成的粘声数据被用作观测数据。图6A和图6B展示了在FWI中包括Q或不包括Q对速度反演具有显著的影响。对于这种特定情况，在FWI中包括固定的Q明显改善了目标气体区域中的速度更新(比较601和602)。相同的工作流程也已应用于3维海洋油田数据并且具有固定Q值的FWI与声学FWI相比产生了明显更好的地质结构。
在所有实际应用中，本技术进步必须与根据本文公开编程的计算机结合使用。优选地，为了高效地执行FWI，计算机是本领域技术人员已知的高性能计算机(HPC)。这种高性能计算机通常涉及节点集群，每个节点具有多个CPU和允许并行计算的计算机存储器。可以使用任何交互式可视化程序和相关联硬件(诸如显示器和投影机)对模型进行可视化和编辑。系统的体系架构可以变化，并且可以由能够执行逻辑操作并根据本技术进步显示输出的任何数量的合适的硬件结构组成。本领域普通技术人员知道可从Cray或IBM获得的合适的超级计算机。
本技术可以容许各种修改和替代形式，并且上述示例仅仅作为示例示出。但是，本技术不旨在限于本文所公开的具体示例。实际上，本技术包括落入所附权利要求的精神和范围内的所有替代方案、修改和等同物。
参考文献
以下参考文献通过引用被整体结合于此：
Bai,J.,Yingst,D.,Bloor,R.,Leveille,J.,2014,Viscoacoustic waveforminversion of velocity structures in the time domain,Geophysics,79,R103-R119；
Liu,J.,2007,Method for performing pseudo-Q migration of seismic data,International Patent Publication WO2009/123790；
Liu,J.,Bear,L.,Krebs,J.,Montelli,R.,and Palacharla,G.,2008,Tomographic inversion by matrix transformation,Geophysics,73(5),VE35-VE38；
Pratt,R.G.,Shin,C.,and Hicks,G.J.,1998,Gauss-Newton and full Newtonmethods in frequency-space seismic waveform inversion:Geophysical JournalInternational,133,341–362；
Tarantola,A.,1984,Inversion of seismic reflection data in theacoustic approximation:Geophysics,49,1259–1266；and
Zhou,J.,Wu,X.,Teng,K.,Xie,Y.,Lefeuvre,F.,Anstey,I.,Sirgue,L.,2014,FWI-guided Q tomography for imaging in the presence of complex gas clouds,76th EAGE conference.

Claims (13)

1.一种方法，包括：
获得由声全波场反演过程生成的速度模型；
利用计算机通过对地下区域的处理后的地震数据应用伪Q迁移来生成变量Q模型，其中所述速度模型被用作所述伪Q迁移中的引导约束；以及
利用计算机产生最终地下速度模型，所述最终地下速度模型通过执行粘声全波场反演过程来恢复由地下区域中的气体异常引起的振幅衰减，其中所述变量Q模型在所述粘声全波场反演过程中是固定的。
2.如权利要求1所述的方法，还包括基于所述最终地下速度模型生成地下图像。
3.如任一前述权利要求所述的方法，还包括：生成处理后的地震数据，其中所述生成包括将基于声射线的叠前深度迁移应用于所述速度模型并且输出共同图像集合。
4.如权利要求3所述的方法，其中生成变量Q模型包括根据引导约束来平坦化所述共同图像集合。
5.如权利要求3所述的方法，其中所述引导约束定义所述速度模型内包含气体异常的区。
6.如权利要求5所述的方法，其中所述伪Q迁移仅应用于包含所述气体异常的区。
7.如任一前述权利要求所述的方法，其中所述变量Q模型通过整个粘声全波场反演过程保持固定。
8.如任一前述权利要求所述的方法，其中生成最终地下速度模型包括通过使粘声共同图像集合平坦化来应用伪Q迁移以构造另一个变量Q模型，并且由粘声全波场反演过程生成的速度模型被用作伪Q迁移中的引导约束。
9.如任一前述权利要求所述的方法，还包括迭代地重复以下操作直到达到预定的停止条件，所述操作包括：执行粘声全波场反演过程，然后从基于粘声射线的叠前深度迁移生成粘声共同图像集合，以及然后通过使粘声共同图像集合平坦化来生成另一个变量Q模型，其中由粘声全波场反演过程生成的速度模型被用作伪Q迁移中的引导约束。
10.如任一前述权利要求所述的方法，还包括进行地震勘测，其中至少一个源用于将声信号注入到地下，并且至少一个接收器用于记录从地下特征反射的声信号。
11.如权利要求2所述的方法，还包括使用最终地下速度模型和地下的图像来提取油气。
12.如权利要求11所述的方法，还包括钻井以提取油气，其中所述井被布置在通过分析地下图像而确定的位置处。
13.如任一前述权利要求所述的方法，其中所述引导约束由从声全波场反演过程反演的地质结构来引导。
CN201680057240.2A 2015-10-02 2016-08-30 Q-补偿的全波场反演 Expired - Fee Related CN108139499B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201562236190P 2015-10-02 2015-10-02
US62/236,190 2015-10-02
PCT/US2016/049392 WO2017058440A1 (en) 2015-10-02 2016-08-30 Q-compensated full wavefield inversion

Publications (2)

Publication Number Publication Date
CN108139499A true CN108139499A (zh) 2018-06-08
CN108139499B CN108139499B (zh) 2020-02-14

Family

ID=56853935

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680057240.2A Expired - Fee Related CN108139499B (zh) 2015-10-02 2016-08-30 Q-补偿的全波场反演

Country Status (10)

Country Link
US (1) US10310113B2 (zh)
EP (1) EP3356862A1 (zh)
KR (1) KR102020759B1 (zh)
CN (1) CN108139499B (zh)
AU (1) AU2016331881B8 (zh)
BR (1) BR112018003117A2 (zh)
CA (1) CA2998522A1 (zh)
MX (1) MX2018002797A (zh)
RU (1) RU2693495C1 (zh)
WO (1) WO2017058440A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703318A (zh) * 2019-09-24 2020-01-17 自然资源部第一海洋研究所 叠前地震数据的直接反演方法

Families Citing this family (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112018070577A2 (pt) 2016-04-07 2019-02-12 Bp Exploration Operating Company Limited detecção de localizações de ingresso de areia de fundo de poço
AU2017246521B2 (en) 2016-04-07 2023-02-02 Bp Exploration Operating Company Limited Detecting downhole sand ingress locations
US10459096B2 (en) 2016-07-13 2019-10-29 Exxonmobil Upstream Research Company Joint full wavefield inversion of P-wave velocity and attenuation using an efficient first order optimization
US10698126B2 (en) * 2016-08-12 2020-06-30 Exxonmobil Upstream Research Company Tomographically enhanced full wavefield inversion
EP3608503B1 (en) 2017-03-31 2022-05-04 BP Exploration Operating Company Limited Well and overburden monitoring using distributed acoustic sensors
US10295685B2 (en) 2017-04-06 2019-05-21 Saudi Arabian Oil Company Generating common image gather using wave-field separation
US11016212B2 (en) 2017-04-11 2021-05-25 Saudi Arabian Oil Company Compressing seismic wavefields in three-dimensional reverse time migration
EA202090528A1 (ru) 2017-08-23 2020-07-10 Бп Эксплорейшн Оперейтинг Компани Лимитед Обнаружение мест скважинных пескопроявлений
US10996372B2 (en) * 2017-08-25 2021-05-04 Exxonmobil Upstream Research Company Geophysical inversion with convolutional neural networks
WO2019072899A2 (en) 2017-10-11 2019-04-18 Bp Exploration Operating Company Limited EVENT DETECTION USING FREQUENCY DOMAIN ACOUSTIC CHARACTERISTICS
DE102018108544A1 (de) * 2018-01-31 2019-08-01 Bauhaus-Universität Weimar Detektion und Lokalisierung von Bauwerksschäden mit vollständiger Wellenformumkehrung
US11656377B2 (en) * 2018-03-30 2023-05-23 Cgg Services Sas Visco-acoustic full waveform inversion of velocity and Q
US11275190B2 (en) 2018-05-16 2022-03-15 Saudi Arabian Oil Company Generating diffraction images based on wave equations
US11340366B2 (en) * 2018-10-17 2022-05-24 Exxonmobil Upstream Research Company Accurate velocity model estimation and imaging in the presence of localized attenuation (Q) anomalies
CA3120493A1 (en) 2018-11-29 2020-06-04 Bp Exploration Operating Company Limited Das data processing to identify fluid inflow locations and fluid type
GB201820331D0 (en) 2018-12-13 2019-01-30 Bp Exploration Operating Co Ltd Distributed acoustic sensing autocalibration
CN112255674A (zh) * 2019-07-22 2021-01-22 中国石油天然气集团有限公司 时空变表层q场建立方法及系统
US11681043B2 (en) * 2019-09-03 2023-06-20 Saudi Arabian Oil Company Diffraction imaging using pseudo dip-angle gather
EP4045766A1 (en) 2019-10-17 2022-08-24 Lytt Limited Fluid inflow characterization using hybrid das/dts measurements
WO2021073740A1 (en) 2019-10-17 2021-04-22 Lytt Limited Inflow detection using dts features
WO2021093974A1 (en) 2019-11-15 2021-05-20 Lytt Limited Systems and methods for draw down improvements across wellbores
US11313988B2 (en) 2019-12-13 2022-04-26 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using seismic diffraction imaging
US11402529B2 (en) 2020-01-09 2022-08-02 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging
US11947062B2 (en) 2020-02-13 2024-04-02 ExxonMobil Technology and Engineering Company Velocity tomography using time lags of wave equation migration
US11467303B2 (en) 2020-03-09 2022-10-11 Saudi Arabian Oil Company Identifying geologic features in a subterranean formation using a post-stack seismic diffraction imaging condition
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum
US11662493B2 (en) 2020-04-08 2023-05-30 ExxonMobil Technology and Engineering Comany Methodology for enhancing properties of geophysical data with deep learning networks
CN113589379B (zh) * 2020-04-30 2023-10-20 中国石油化工股份有限公司 叠前道集保真性振幅补偿优化方法
CN111736220B (zh) * 2020-05-13 2023-02-10 中国石油天然气集团有限公司 一种逆时偏移成像方法及装置
US11656378B2 (en) 2020-06-08 2023-05-23 Saudi Arabian Oil Company Seismic imaging by visco-acoustic reverse time migration
EP4165284A1 (en) 2020-06-11 2023-04-19 Lytt Limited Systems and methods for subterranean fluid flow characterization
EP4168647A1 (en) 2020-06-18 2023-04-26 Lytt Limited Event model training using in situ data
CN112305590B (zh) * 2020-09-23 2023-12-22 中国石油天然气集团有限公司 粘声介质叠前时间偏移计算方法及装置
CN112415601A (zh) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 表层品质因子q值的确定方法及装置
WO2022245555A1 (en) 2021-05-21 2022-11-24 Exxonmobil Upstream Research Company Method and system for rock physics constrained data integration of elastic fwi properties and seismic stacks
CN113219536B (zh) * 2021-06-25 2022-03-01 成都理工大学 依赖频率的纵横波衰减参数的叠前地震反演方法
CN113866827B (zh) * 2021-09-29 2023-05-09 中国石油大学(华东) 一种解释性速度建模地震成像方法、系统、介质和设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1516814A (zh) * 2001-06-16 2004-07-28 ά˹�ض���Ƶ���ع����޹�˾ 处理数据的方法
CN101334483A (zh) * 2008-06-13 2008-12-31 徐基祥 一种在地震数据处理中衰减瑞雷波散射噪声的方法
WO2009123790A1 (en) * 2008-03-31 2009-10-08 Exxonmobil Upstream Research Company Method for performing pseudo-q migration of seismic data
US20120155217A1 (en) * 2010-12-16 2012-06-21 Bp Corporation North America Inc. Seismic acquisition using narrowband seismic sources

Family Cites Families (229)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3812457A (en) 1969-11-17 1974-05-21 Shell Oil Co Seismic exploration method
US3864667A (en) 1970-09-11 1975-02-04 Continental Oil Co Apparatus for surface wave parameter determination
US3984805A (en) 1973-10-18 1976-10-05 Daniel Silverman Parallel operation of seismic vibrators without phase control
US4168485A (en) 1974-08-12 1979-09-18 Continental Oil Company Simultaneous use of pseudo-random control signals in vibrational exploration methods
US4675851A (en) 1982-09-09 1987-06-23 Western Geophysical Co. Method for seismic exploration
US4545039A (en) 1982-09-09 1985-10-01 Western Geophysical Co. Of America Methods for seismic exploration
US4575830A (en) 1982-10-15 1986-03-11 Schlumberger Technology Corporation Indirect shearwave determination
US4594662A (en) 1982-11-12 1986-06-10 Schlumberger Technology Corporation Diffraction tomography systems and methods with fixed detector arrays
JPS59189278A (ja) 1983-03-23 1984-10-26 橋本電機工業株式会社 ウイケツト型平板乾燥機
JPS606032A (ja) 1983-06-22 1985-01-12 Honda Motor Co Ltd 内燃エンジンの作動状態制御方法
US4694404A (en) 1984-01-12 1987-09-15 Key Bank N.A. High-speed image generation of complex solid objects using octree encoding
US4924390A (en) 1985-03-04 1990-05-08 Conoco, Inc. Method for determination of earth stratum elastic parameters using seismic energy
US4715020A (en) 1986-10-29 1987-12-22 Western Atlas International, Inc. Simultaneous performance of multiple seismic vibratory surveys
FR2589587B1 (fr) 1985-10-30 1988-02-05 Inst Francais Du Petrole Procede de prospection sismique marine utilisant un signal vibratoire code et dispositif pour sa mise en oeuvre
US4707812A (en) 1985-12-09 1987-11-17 Atlantic Richfield Company Method of suppressing vibration seismic signal correlation noise
US4823326A (en) 1986-07-21 1989-04-18 The Standard Oil Company Seismic data acquisition technique having superposed signals
US4686654A (en) 1986-07-31 1987-08-11 Western Geophysical Company Of America Method for generating orthogonal sweep signals
US4766574A (en) 1987-03-31 1988-08-23 Amoco Corporation Method for depth imaging multicomponent seismic data
US4953657A (en) 1987-11-30 1990-09-04 Halliburton Geophysical Services, Inc. Time delay source coding
US4969129A (en) 1989-09-20 1990-11-06 Texaco Inc. Coding seismic sources
US4982374A (en) 1989-10-23 1991-01-01 Halliburton Geophysical Services, Inc. Method of source coding and harmonic cancellation for vibrational geophysical survey sources
GB9011836D0 (en) 1990-05-25 1990-07-18 Mason Iain M Seismic surveying
US5469062A (en) 1994-03-11 1995-11-21 Baker Hughes, Inc. Multiple depths and frequencies for simultaneous inversion of electromagnetic borehole measurements
GB2293010B (en) 1994-07-07 1998-12-09 Geco As Method of processing seismic data
US5583825A (en) 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
DE69625978T2 (de) 1995-04-18 2003-11-13 Western Atlas Int Inc Gleichmässige unteroberflächendeckung bei anwesenheit von steilen neigungen
US5924049A (en) 1995-04-18 1999-07-13 Western Atlas International, Inc. Methods for acquiring and processing seismic data
US5719821A (en) 1995-09-29 1998-02-17 Atlantic Richfield Company Method and apparatus for source separation of seismic vibratory signals
US5721710A (en) 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation
US5822269A (en) 1995-11-13 1998-10-13 Mobil Oil Corporation Method for separation of a plurality of vibratory seismic energy source signals
US5715213A (en) 1995-11-13 1998-02-03 Mobil Oil Corporation High fidelity vibratory source seismic method using a plurality of vibrator sources
US5790473A (en) 1995-11-13 1998-08-04 Mobil Oil Corporation High fidelity vibratory source seismic method for use in vertical seismic profile data gathering with a plurality of vibratory seismic energy sources
US5838634A (en) 1996-04-04 1998-11-17 Exxon Production Research Company Method of generating 3-D geologic models incorporating geologic and geophysical constraints
US5798982A (en) 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
GB9612471D0 (en) 1996-06-14 1996-08-14 Geco As Method and apparatus for multiple seismic vibratory surveys
US5878372A (en) 1997-03-04 1999-03-02 Western Atlas International, Inc. Method for simultaneous inversion processing of well log data using a plurality of earth models
US6014342A (en) 1997-03-21 2000-01-11 Tomo Seis, Inc. Method of evaluating a subsurface region using gather sensitive data discrimination
US5999489A (en) 1997-03-21 1999-12-07 Tomoseis Inc. High vertical resolution crosswell seismic imaging
US5920838A (en) 1997-06-02 1999-07-06 Carnegie Mellon University Reading and pronunciation tutor
FR2765692B1 (fr) 1997-07-04 1999-09-10 Inst Francais Du Petrole Methode pour modeliser en 3d l'impedance d'un milieu heterogene
GB2329043B (en) 1997-09-05 2000-04-26 Geco As Method of determining the response caused by model alterations in seismic simulations
US5999488A (en) 1998-04-27 1999-12-07 Phillips Petroleum Company Method and apparatus for migration by finite differences
US6219621B1 (en) 1998-06-30 2001-04-17 Exxonmobil Upstream Research Co. Sparse hyperbolic inversion of seismic data
US6388947B1 (en) 1998-09-14 2002-05-14 Tomoseis, Inc. Multi-crosswell profile 3D imaging and method
FR2784195B1 (fr) 1998-10-01 2000-11-17 Inst Francais Du Petrole Methode pour realiser en 3d avant sommation, une migration de donnees sismiques
US6574564B2 (en) 1998-10-01 2003-06-03 Institut Francais Du Petrole 3D prestack seismic data migration method
US6225803B1 (en) 1998-10-29 2001-05-01 Baker Hughes Incorporated NMR log processing using wavelet filter and iterative inversion
US6021094A (en) 1998-12-03 2000-02-01 Sandia Corporation Method of migrating seismic records
US6754588B2 (en) 1999-01-29 2004-06-22 Platte River Associates, Inc. Method of predicting three-dimensional stratigraphy using inverse optimization techniques
AU3229900A (en) 1999-02-12 2000-08-29 Prange, Michael Uncertainty constrained subsurface modeling
US6058073A (en) 1999-03-30 2000-05-02 Atlantic Richfield Company Elastic impedance estimation for inversion of far offset seismic sections
FR2792419B1 (fr) 1999-04-16 2001-09-07 Inst Francais Du Petrole Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol
GB9927395D0 (en) 1999-05-19 2000-01-19 Schlumberger Holdings Improved seismic data acquisition method
US6327537B1 (en) 1999-07-19 2001-12-04 Luc T. Ikelle Multi-shooting approach to seismic modeling and acquisition
FR2798197B1 (fr) 1999-09-02 2001-10-05 Inst Francais Du Petrole Methode pour former un modele d'une formation geologique, contraint par des donnees dynamiques et statiques
EP2296013B1 (en) 1999-10-22 2016-03-30 CGG Services (NL) B.V. Method of estimating elastic and compositional parameters from seismic and echo-acoustic data
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2800473B1 (fr) 1999-10-29 2001-11-30 Inst Francais Du Petrole Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques
DE19954866A1 (de) 1999-11-15 2001-05-31 Infineon Technologies Ag Verfahren zur Behandlung einer durch Epitaxie hergestellten Oberfläche eines SiC-Halbleiterkörpers und danach hergestellten Schottkykontakt
CN1188710C (zh) 2000-01-21 2005-02-09 施鲁博格控股有限公司 估算地震介质特性的系统和方法
CN1188711C (zh) 2000-01-21 2005-02-09 施鲁博格控股有限公司 用于地震波场分离的系统和方法
US6826486B1 (en) 2000-02-11 2004-11-30 Schlumberger Technology Corporation Methods and apparatus for predicting pore and fracture pressures of a subsurface formation
FR2805051B1 (fr) 2000-02-14 2002-12-06 Geophysique Cie Gle Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques
GB2359363B (en) 2000-02-15 2002-04-03 Geco Prakla Processing simultaneous vibratory seismic data
US6687659B1 (en) 2000-03-24 2004-02-03 Conocophillips Company Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications
US6317695B1 (en) 2000-03-30 2001-11-13 Nutec Sciences, Inc. Seismic data processing method
US6546339B2 (en) 2000-08-07 2003-04-08 3D Geo Development, Inc. Velocity analysis using angle-domain common image gathers
US6950786B1 (en) 2000-10-10 2005-09-27 Schlumberger Technology Corporation Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot
CA2426160A1 (en) 2000-10-17 2002-04-25 David Lee Nyland Method of using cascaded sweeps for source coding and harmonic cancellation
AU2002239619A1 (en) 2000-12-08 2002-06-18 Peter J. Ortoleva Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories
FR2818753B1 (fr) 2000-12-21 2003-03-21 Inst Francais Du Petrole Methode et dispositif de prospection sismique par emission simultanee de signaux sismisques obtenus en codant un signal par des sequences pseudo aleatoires
FR2821677B1 (fr) 2001-03-05 2004-04-30 Geophysique Cie Gle Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees
US6751558B2 (en) 2001-03-13 2004-06-15 Conoco Inc. Method and process for prediction of subsurface fluid and rock pressures in the earth
US6927698B2 (en) 2001-08-27 2005-08-09 Larry G. Stolarczyk Shuttle-in receiver for radio-imaging underground geologic structures
US6545944B2 (en) 2001-05-30 2003-04-08 Westerngeco L.L.C. Method for acquiring and processing of data from two or more simultaneously fired sources
US6882958B2 (en) 2001-06-28 2005-04-19 National Instruments Corporation System and method for curve fitting using randomized techniques
GB2379013B (en) 2001-08-07 2005-04-20 Abb Offshore Systems Ltd Microseismic signal processing
US6593746B2 (en) 2001-08-27 2003-07-15 Larry G. Stolarczyk Method and system for radio-imaging underground geologic structures
US7672824B2 (en) 2001-12-10 2010-03-02 Westerngeco L.L.C. Method for shallow water flow detection
US7069149B2 (en) 2001-12-14 2006-06-27 Chevron U.S.A. Inc. Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume
US7330799B2 (en) 2001-12-21 2008-02-12 Société de commercialisation des produits de la recherche appliquée-Socpra Sciences et Génie s.e.c. Method and algorithm for using surface waves
US6842701B2 (en) 2002-02-25 2005-01-11 Westerngeco L.L.C. Method of noise removal for cascaded sweep data
GB2387226C (en) 2002-04-06 2008-05-12 Westerngeco Ltd A method of seismic surveying
FR2839368B1 (fr) 2002-05-06 2004-10-01 Total Fina Elf S A Methode de decimation de traces sismiques pilotee par le trajet sismique
US6832159B2 (en) 2002-07-11 2004-12-14 Schlumberger Technology Corporation Intelligent diagnosis of environmental influence on well logs with model-based inversion
US6906981B2 (en) 2002-07-17 2005-06-14 Pgs Americas, Inc. Method and system for acquiring marine seismic data using multiple seismic sources
FR2843202B1 (fr) 2002-08-05 2004-09-10 Inst Francais Du Petrole Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration
US6832155B2 (en) 2002-09-23 2004-12-14 Itt Manufacturing Enterprises, Inc. Methods and apparatus for determining phase ambiguities in ranging and navigation systems
AU2003279870A1 (en) 2002-10-04 2004-05-04 Paradigm Geophysical Corporation Method and system for limited frequency seismic imaging
GB2396448B (en) 2002-12-21 2005-03-02 Schlumberger Holdings System and method for representing and processing and modeling subterranean surfaces
US7027927B2 (en) 2002-12-23 2006-04-11 Schlumberger Technology Corporation Methods for determining formation and borehole parameters using fresnel volume tomography
US6735527B1 (en) 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
US6999880B2 (en) 2003-03-18 2006-02-14 The Regents Of The University Of California Source-independent full waveform inversion of seismic data
US7184367B2 (en) 2003-03-27 2007-02-27 Exxonmobil Upstream Research Company Method to convert seismic traces into petrophysical property logs
WO2004095073A2 (en) 2003-04-01 2004-11-04 Exxonmobil Upstream Research Company Shaped high frequency vibratory source
US7072767B2 (en) 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
NO322089B1 (no) 2003-04-09 2006-08-14 Norsar V Daglig Leder Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder
GB2400438B (en) 2003-04-11 2005-06-01 Westerngeco Ltd Determination of waveguide parameters
US6970397B2 (en) 2003-07-09 2005-11-29 Gas Technology Institute Determination of fluid properties of earth formations using stochastic inversion
US6882938B2 (en) 2003-07-30 2005-04-19 Pgs Americas, Inc. Method for separating seismic signals from two or more distinct sources
GB2405473B (en) 2003-08-23 2005-10-05 Westerngeco Ltd Multiple attenuation method
US6944546B2 (en) 2003-10-01 2005-09-13 Halliburton Energy Services, Inc. Method and apparatus for inversion processing of well logging data in a selected pattern space
US6901333B2 (en) 2003-10-27 2005-05-31 Fugro N.V. Method and device for the generation and application of anisotropic elastic parameters
US7046581B2 (en) 2003-12-01 2006-05-16 Shell Oil Company Well-to-well tomography
US7716251B2 (en) 2003-12-02 2010-05-11 International Business Machines Corporation Systems and method for indexing, searching and retrieving semantic objects
US20050128874A1 (en) 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US7359283B2 (en) 2004-03-03 2008-04-15 Pgs Americas, Inc. System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers
US7791980B2 (en) 2004-05-21 2010-09-07 Westerngeco L.L.C. Interpolation and extrapolation method for seismic recordings
FR2872584B1 (fr) 2004-06-30 2006-08-11 Inst Francais Du Petrole Methode pour simuler le depot sedimentaire dans un bassin respectant les epaisseurs des sequences sedimentaires
MX2007000170A (es) 2004-07-07 2007-03-30 Exxonmobil Upstream Res Co Aplicaciones de red bayesiana para geologia y geografia.
EP1617309B1 (en) 2004-07-15 2011-01-12 Fujitsu Limited Simulation technique with local grid refinement
US7646924B2 (en) 2004-08-09 2010-01-12 David Leigh Donoho Method and apparatus for compressed sensing
US7480206B2 (en) 2004-09-13 2009-01-20 Chevron U.S.A. Inc. Methods for earth modeling and seismic imaging using interactive and selective updating
FR2876458B1 (fr) 2004-10-08 2007-01-19 Geophysique Cie Gle Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples
US7088639B2 (en) * 2004-11-17 2006-08-08 Rdsp I L.P. Method for determining formation quality factor from well log data and its application to seismic reservoir characterization
GB2422433B (en) 2004-12-21 2008-03-19 Sondex Wireline Ltd Method and apparatus for determining the permeability of earth formations
US7373251B2 (en) 2004-12-22 2008-05-13 Marathon Oil Company Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data
US7230879B2 (en) 2005-02-12 2007-06-12 Chevron U.S.A. Inc. Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects
WO2006090374A2 (en) 2005-02-22 2006-08-31 Paradigm Geophysical Ltd. Multiple suppression in angle domain time and depth migration
US7840625B2 (en) 2005-04-07 2010-11-23 California Institute Of Technology Methods for performing fast discrete curvelet transforms of data
US7271747B2 (en) 2005-05-10 2007-09-18 Rice University Method and apparatus for distributed compressed sensing
US7405997B2 (en) 2005-08-11 2008-07-29 Conocophillips Company Method of accounting for wavelet stretch in seismic data
US7397728B2 (en) * 2005-08-26 2008-07-08 Westerngeco L.L.C. Method for processing a record of seismic traces
RU2440604C2 (ru) 2005-10-18 2012-01-20 Синвент Ас Визуализация данных отклика геологической среды с использованием потоковых процессоров
US7373252B2 (en) 2005-11-04 2008-05-13 Western Geco L.L.C. 3D pre-stack full waveform inversion
FR2895091B1 (fr) 2005-12-21 2008-02-22 Inst Francais Du Petrole Methode pour mettre a jour un modele geologique par des donnees sismiques
GB2436626B (en) 2006-03-28 2008-08-06 Westerngeco Seismic Holdings Method of evaluating the interaction between a wavefield and a solid body
US7620534B2 (en) 2006-04-28 2009-11-17 Saudi Aramco Sound enabling computerized system for real time reservoir model calibration using field surveillance data
US20070274155A1 (en) 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
US7725266B2 (en) 2006-05-31 2010-05-25 Bp Corporation North America Inc. System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling
US7599798B2 (en) 2006-09-11 2009-10-06 Westerngeco L.L.C. Migrating composite seismic response data to produce a representation of a seismic volume
US8121823B2 (en) 2006-09-28 2012-02-21 Exxonmobil Upstream Research Company Iterative inversion of data from simultaneous geophysical sources
MX2009001681A (es) 2006-12-07 2009-02-26 Council Scient Ind Res Un metodo para computar una respuesta de impulso exacta de un reflector acustico plano a offset cero debido a una fuente acustica puntual.
ATE543109T1 (de) 2007-01-20 2012-02-15 Spectraseis Ag Zeitumkehr-reservoir-lokalisierung
US8248886B2 (en) 2007-04-10 2012-08-21 Exxonmobil Upstream Research Company Separation and noise removal for multiple vibratory source seismic data
US7715986B2 (en) 2007-05-22 2010-05-11 Chevron U.S.A. Inc. Method for identifying and removing multiples for imaging with beams
US7974824B2 (en) 2007-06-29 2011-07-05 Westerngeco L. L. C. Seismic inversion of data containing surface-related multiples
JP2009063942A (ja) 2007-09-10 2009-03-26 Sumitomo Electric Ind Ltd 遠赤外線カメラ用レンズ、レンズユニット及び撮像装置
US20090070042A1 (en) 2007-09-11 2009-03-12 Richard Birchwood Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state
US20090083006A1 (en) 2007-09-20 2009-03-26 Randall Mackie Methods and apparatus for three-dimensional inversion of electromagnetic data
WO2009067041A1 (en) 2007-11-19 2009-05-28 Steklov Mathematical Institute Ras Method and system for evaluating the characteristic properties of two contacting media and of the interface between them based on mixed surface waves propagating along the interface
US7732381B2 (en) 2007-11-30 2010-06-08 Schlumberger Technology Corporation Conductive cement formulation and application for use in wells
US20090164186A1 (en) 2007-12-20 2009-06-25 Bhp Billiton Innovation Pty Ltd. Method for determining improved estimates of properties of a model
CN101910871A (zh) 2008-01-08 2010-12-08 埃克森美孚上游研究公司 地震数据的频谱整形反演和偏移
US8577660B2 (en) 2008-01-23 2013-11-05 Schlumberger Technology Corporation Three-dimensional mechanical earth modeling
EA017177B1 (ru) 2008-03-21 2012-10-30 Эксонмобил Апстрим Рисерч Компани Эффективный способ инверсии геофизических данных
EP2105765A1 (en) 2008-03-28 2009-09-30 Schlumberger Holdings Limited Simultaneous inversion of induction data for dielectric permittivity and electric conductivity
EP2269094A4 (en) 2008-03-28 2017-06-14 Exxonmobil Upstream Research Company Characterizing spatial variablility of surface waves in seismic processing
US8275592B2 (en) 2008-04-07 2012-09-25 Westerngeco L.L.C. Joint inversion of time domain controlled source electromagnetic (TD-CSEM) data and further data
US8494777B2 (en) 2008-04-09 2013-07-23 Schlumberger Technology Corporation Continuous microseismic mapping for real-time 3D event detection and location
US8345510B2 (en) 2008-06-02 2013-01-01 Pgs Geophysical As Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s)
US20110182141A1 (en) 2008-08-14 2011-07-28 Schlumberger Technology Corporation Method and system for monitoring a logging tool position in a borehole
WO2010019957A1 (en) 2008-08-15 2010-02-18 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US8559270B2 (en) 2008-08-15 2013-10-15 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US20100054082A1 (en) 2008-08-29 2010-03-04 Acceleware Corp. Reverse-time depth migration with reduced memory requirements
US8296069B2 (en) 2008-10-06 2012-10-23 Bp Corporation North America Inc. Pseudo-analytical method for the solution of wave equations
US7616523B1 (en) 2008-10-22 2009-11-10 Pgs Geophysical As Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth
US9213119B2 (en) 2008-10-29 2015-12-15 Conocophillips Company Marine seismic acquisition
US20100118651A1 (en) 2008-11-10 2010-05-13 Chevron U.S.A. Inc. Method for generation of images related to a subsurface region of interest
US20100142316A1 (en) 2008-12-07 2010-06-10 Henk Keers Using waveform inversion to determine properties of a subsurface medium
US8095345B2 (en) 2009-01-20 2012-01-10 Chevron U.S.A. Inc Stochastic inversion of geophysical data for estimating earth model parameters
US8369184B2 (en) 2009-01-26 2013-02-05 Shotspotter, Inc. Systems and methods with improved three-dimensional source location processing including constraint of location solutions to a two-dimensional plane
US9052410B2 (en) 2009-02-12 2015-06-09 Conocophillips Company Multiple seismic signal inversion
US20110299361A1 (en) 2009-02-17 2011-12-08 Changsoo Shin Apparatus and method for imaging subsurface structure
US9110191B2 (en) 2009-03-30 2015-08-18 Westerngeco L.L.C. Multiple attenuation for ocean-bottom seismic data
US8547794B2 (en) 2009-04-16 2013-10-01 Baker Hughes Incorporated Extending the coverage of VSP/CDP imaging by using first-order downgoing multiples
US9075163B2 (en) 2009-04-17 2015-07-07 Westerngeco L.L.C. Interferometric seismic data processing
GB2483202B (en) 2009-07-06 2013-10-16 Exxonmobil Upstream Res Co Method for seismic interpretation using seismic texture attributes
US8176284B2 (en) 2009-08-11 2012-05-08 Texas Memory Systems, Inc. FLASH-based memory system with variable length page stripes including data protection information
US20110044127A1 (en) 2009-08-19 2011-02-24 Clement Kostov Removing free-surface effects from seismic data acquired in a towed survey
US8923093B2 (en) 2009-08-25 2014-12-30 Westerngeco L.L.C. Determining the quality of a seismic inversion
US20110131020A1 (en) * 2009-09-09 2011-06-02 Conocophillips Company Dip guided full waveform inversion
US8406081B2 (en) 2009-09-25 2013-03-26 Landmark Graphics Corporation Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers
GB2486121B (en) 2009-10-01 2014-08-13 Halliburton Energy Serv Inc Apparatus and methods of locating downhole anomalies
US9244181B2 (en) 2009-10-19 2016-01-26 Westerngeco L.L.C. Full-waveform inversion in the traveltime domain
WO2011071812A2 (en) 2009-12-07 2011-06-16 Geco Technology B.V. Simultaneous joint inversion of surface wave and refraction data
FR2955396B1 (fr) 2010-01-15 2013-03-01 Cggveritas Services Sa Dispositif de traitement de donnees sismiques marines
US9482775B2 (en) 2010-01-22 2016-11-01 Schlumberger Technology Corporation Real-time formation anisotropy and dip evaluation using tri-axial induction measurements
CA2787693A1 (en) 2010-01-25 2011-07-28 CGGVeritas Services (U.S.) Inc. Methods and systems for estimating stress using seismic data
US20130098608A1 (en) 2010-01-29 2013-04-25 Robert Barnum Temporary field storage of gas to optimize field development
US8265875B2 (en) 2010-01-29 2012-09-11 Westerngeco L.L.C. Interpolation of periodic data
US8537638B2 (en) 2010-02-10 2013-09-17 Exxonmobil Upstream Research Company Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration
US8612156B2 (en) 2010-03-05 2013-12-17 Vialogy Llc Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization
CA2791694C (en) 2010-03-12 2022-07-12 CGGVeritas Services (U.S.) Inc. Methods and systems for performing azimuthal simultaneous elastic inversion
US8680865B2 (en) 2010-03-19 2014-03-25 Schlumberger Technology Corporation Single well reservoir imaging apparatus and methods
US20110235464A1 (en) 2010-03-24 2011-09-29 John Brittan Method of imaging the earth's subsurface during marine seismic data acquisition
US8223587B2 (en) 2010-03-29 2012-07-17 Exxonmobil Upstream Research Company Full wavefield inversion using time varying filters
US9176244B2 (en) 2010-03-31 2015-11-03 Schlumberger Technology Corporation Data set inversion using source-receiver compression
KR101167715B1 (ko) 2010-04-30 2012-07-20 서울대학교산학협력단 복소 구배 최소자승법에의한 파형 역산을 이용한 지하 구조의 영상화 장치 및 방법
US8576663B2 (en) 2010-04-30 2013-11-05 Schlumberger Technology Corporation Multicomponent seismic inversion of VSP data
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8756042B2 (en) 2010-05-19 2014-06-17 Exxonmobile Upstream Research Company Method and system for checkpointing during simulations
US9061485B2 (en) 2010-06-15 2015-06-23 Denki Kagaku Kogyo Kabushiki Kaisha Method of manufacturing translucent rigid substrate laminate
US20110320180A1 (en) 2010-06-29 2011-12-29 Al-Saleh Saleh M Migration Velocity Analysis of Seismic Data Using Common Image Cube and Green's Functions
US8612188B2 (en) 2010-07-12 2013-12-17 The University Of Manchester Wave modelling
US9195783B2 (en) 2010-08-16 2015-11-24 Exxonmobil Upstream Research Company Reducing the dimensionality of the joint inversion problem
US20120051176A1 (en) 2010-08-31 2012-03-01 Chevron U.S.A. Inc. Reverse time migration back-scattering noise removal using decomposed wavefield directivity
CN103119471A (zh) 2010-09-20 2013-05-22 雪佛龙美国公司 用于生成地下构造的图像的系统和方法
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
KR101908278B1 (ko) 2010-09-27 2018-10-17 엑손모빌 업스트림 리서치 캄파니 풀 파동장 반전을 위한 실제적인 해결책으로서의 동시 소스 인코딩 및 소스 분리
US20130311151A1 (en) 2010-09-28 2013-11-21 René-Edouard André Michel Plessix Earth model estimation through an acoustic full waveform inversion of seismic data
RU2587498C2 (ru) 2010-12-01 2016-06-20 Эксонмобил Апстрим Рисерч Компани Инверсия одновременных источников для данных сейсмоприемной косы с взаимнокорреляционной целевой функцией
AU2012212520B2 (en) 2011-01-31 2014-09-11 Chevron U.S.A. Inc. Extracting geologic information from multiple offset stacks and/or angle stacks
US9702994B2 (en) 2011-02-18 2017-07-11 Westerngeco L.L.C. Waveform inversion by multiple shot-encoding for non-fixed spread geometries
EP2691795A4 (en) 2011-03-30 2015-12-09 CONVERGENCE SPEED OF COMPLETE WAVELENGTH INVERSION USING SPECTRAL SHAPING
US20120275267A1 (en) 2011-04-26 2012-11-01 Ramesh Neelamani Seismic Data Processing
US9513389B2 (en) 2011-05-13 2016-12-06 Saudi Arabian Oil Company Coupled time-distance dependent swept frequency source acquisition design and data de-noising
US20120316790A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase extrapolation
US20120316844A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase unwrapping
US20120316791A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion by non-linear model update
US9075159B2 (en) 2011-06-08 2015-07-07 Chevron U.S.A., Inc. System and method for seismic data inversion
CA2839277C (en) 2011-09-02 2018-02-27 Exxonmobil Upstream Research Company Using projection onto convex sets to constrain full-wavefield inversion
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
RU2612896C2 (ru) 2012-03-08 2017-03-13 Эксонмобил Апстрим Рисерч Компани Ортогональное кодирование источника и приемника
US9435905B2 (en) 2012-04-19 2016-09-06 Cgg Services Sa Premigration deghosting of seismic data with a bootstrap technique
US9541661B2 (en) 2012-04-19 2017-01-10 Cgg Services Sa Device and method for deghosting variable depth streamer data
US20130311149A1 (en) 2012-05-17 2013-11-21 Yaxun Tang Tomographically Enhanced Full Wavefield Inversion
US9442204B2 (en) * 2012-08-06 2016-09-13 Exxonmobil Upstream Research Company Seismic inversion for formation properties and attenuation effects
US9952340B2 (en) 2013-03-15 2018-04-24 General Electric Company Context based geo-seismic object identification
US9702993B2 (en) 2013-05-24 2017-07-11 Exxonmobil Upstream Research Company Multi-parameter inversion through offset dependent elastic FWI
US10459117B2 (en) 2013-06-03 2019-10-29 Exxonmobil Upstream Research Company Extended subspace method for cross-talk mitigation in multi-parameter inversion
US20140372043A1 (en) 2013-06-17 2014-12-18 Wenyi Hu Full Waveform Inversion Using Perfectly Reflectionless Subgridding
WO2015100544A1 (zh) 2013-12-30 2015-07-09 中国石油天然气股份有限公司 基于零偏垂直地震剖面数据估计品质因子的方法和装置
US20150253445A1 (en) 2014-03-06 2015-09-10 Westerngeco L.L.C. Visualization of seismic attributes
CA2947410A1 (en) 2014-06-17 2015-12-30 Exxonmobil Upstream Research Company Fast viscoacoustic and viscoelastic full-wavefield inversion
WO2016001750A2 (en) * 2014-06-30 2016-01-07 Cgg Services Sa Seismic data processing using matching filter based cost function optimization
US9829592B2 (en) * 2014-12-16 2017-11-28 Pgs Geophysical As Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method
US10670751B2 (en) * 2015-03-27 2020-06-02 Cgg Services Sas Full waveform inversion method for seismic data processing using preserved amplitude reverse time migration
EP3076205B1 (en) * 2015-03-31 2023-06-14 CGG Services SAS Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1516814A (zh) * 2001-06-16 2004-07-28 ά˹�ض���Ƶ���ع����޹�˾ 处理数据的方法
WO2009123790A1 (en) * 2008-03-31 2009-10-08 Exxonmobil Upstream Research Company Method for performing pseudo-q migration of seismic data
CN101334483A (zh) * 2008-06-13 2008-12-31 徐基祥 一种在地震数据处理中衰减瑞雷波散射噪声的方法
US20120155217A1 (en) * 2010-12-16 2012-06-21 Bp Corporation North America Inc. Seismic acquisition using narrowband seismic sources

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703318A (zh) * 2019-09-24 2020-01-17 自然资源部第一海洋研究所 叠前地震数据的直接反演方法
CN110703318B (zh) * 2019-09-24 2021-06-11 自然资源部第一海洋研究所 叠前地震数据的直接反演方法

Also Published As

Publication number Publication date
US10310113B2 (en) 2019-06-04
RU2693495C1 (ru) 2019-07-03
US20170097428A1 (en) 2017-04-06
AU2016331881B8 (en) 2019-07-18
AU2016331881A1 (en) 2018-04-26
KR102020759B1 (ko) 2019-09-11
MX2018002797A (es) 2018-04-13
EP3356862A1 (en) 2018-08-08
AU2016331881B2 (en) 2019-07-11
CN108139499B (zh) 2020-02-14
BR112018003117A2 (pt) 2018-09-25
KR20180059539A (ko) 2018-06-04
CA2998522A1 (en) 2017-04-06
WO2017058440A1 (en) 2017-04-06
AU2016331881A8 (en) 2019-07-18

Similar Documents

Publication Publication Date Title
CN108139499B (zh) Q-补偿的全波场反演
CA3122509C (en) Machine learning-augmented geophysical inversion
US10459096B2 (en) Joint full wavefield inversion of P-wave velocity and attenuation using an efficient first order optimization
US9829592B2 (en) Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method
CA2964893C (en) Structure tensor constrained tomographic velocity analysis
EP3803469A1 (en) Method for generating an image of a subsurface of an area of interest from seismic data
US10310117B2 (en) Efficient seismic attribute gather generation with data synthesis and expectation method
Thiel et al. Comparison of acoustic and elastic full‐waveform inversion of 2D towed‐streamer data in the presence of salt
EP3243089A1 (en) Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume over a time period
CN115877449A (zh) 用于在勘测区域内获得地下堆叠图像的计算机实现方法
WO2005019866A2 (en) Method and apparatus for kinematically linking multiple seismic domains
EA030770B1 (ru) Система и способ адаптивной сейсмической оптики
Jaimes-Osorio et al. Amplitude variation with offset inversion using acoustic-elastic local solver
WO2015094458A1 (en) A method of correcting velocity for complex surface topography
US11255993B2 (en) Variable aperture estimation using bottom-up ray tracing
US11231514B2 (en) Method for attenuation compensation utilizing non-stationary matching filters
Kim et al. Locating microseismic events using both head and direct waves travelling in vertical transverse isotropic media
US20240159930A1 (en) Method and apparatus for implementing full waveform inversion using angle gathers
Aftab et al. The effect of the number, dimensions and location of anomalies on the crosswell Traveltime tomography results
Zuniga et al. Velocity Inversion of PSP-waves on Marine Seismic Data using OBN Technology
Dzulkefli et al. Full Wavefield Redatuming: Accurate Velocity Modelling for Imaging Beneath Complex Overburden
Brittan et al. Dense multi-offset reflection tomography
EA044564B1 (ru) Построение модели скорости
McAdow et al. High Resolution Full Wavefield Inversion: Impact on Imaging and Reservoir Characterisation
Re et al. Letting the Data Drive the Initial Model Building for Refraction Tomography

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200214

Termination date: 20200830