CN103238158A - 利用互相关目标函数进行的海洋拖缆数据同时源反演 - Google Patents

利用互相关目标函数进行的海洋拖缆数据同时源反演 Download PDF

Info

Publication number
CN103238158A
CN103238158A CN201180058068XA CN201180058068A CN103238158A CN 103238158 A CN103238158 A CN 103238158A CN 201180058068X A CN201180058068X A CN 201180058068XA CN 201180058068 A CN201180058068 A CN 201180058068A CN 103238158 A CN103238158 A CN 103238158A
Authority
CN
China
Prior art keywords
coding
data
source
collection
physical property
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
CN201180058068XA
Other languages
English (en)
Other versions
CN103238158B (zh
Inventor
P·S·如斯
J·R·克雷布斯
S·拉扎拉托斯
A·鲍姆斯泰因
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 CN103238158A publication Critical patent/CN103238158A/zh
Application granted granted Critical
Publication of CN103238158B publication Critical patent/CN103238158B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06G—ANALOGUE COMPUTERS
    • G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/48—Analogue computers for specific processes, systems or devices, e.g. simulators
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003—Seismic data acquisition in general, e.g. survey design
    • G01V1/005—Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
    • 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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/614—Synthetically generated data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

源（或接收器）编码（30）地球物理数据（80）道集的同时全波场反演方法，以确定地下区域的物理性质模型（20），尤其适于在数据获取（40）中不符合固定接收器几何条件的考察。反演包括优化互相关目标函数（100）。

Description

利用互相关目标函数进行的海洋拖缆数据同时源反演
相关申请的交叉参考
本申请要求2010年12月1日提交的发明名称为“Simultaneous Source Inversionfor Marine Streamer Data with Cross-Correlation Objective Function（利用互相关目标函数进行的海洋拖缆数据同时源反演）”的美国临时专利申请61/418,694和2011年7月20日提交的相同发明名称的美国临时专利申请61/509,904的权益，二者的全部内容均被引入本文作为参考。
发明领域
本发明总体涉及地球物理勘探领域，更具体地涉及地球物理数据处理。具体地，本发明是得自多种地球物理源如震源的数据反演方法，包括地球物理模拟——在进行一次模拟中计算来自多个同时活动的地球物理源的数据。
发明背景
即使应用现代的计算能力，地震全波场反演仍是计算昂贵的工作。但是，用此方法获得详细的地下表示的益处预期超过该阻碍。导致更快周转时间的算法和工作流程的发展是使此技术适于场大规模数据的关键步骤。地震全波形反演包括数据的正向和伴随模拟的数次迭代（重复，iteration）。因此，降低正向和伴随计算运行成本的技术将致使用户在合理的时间量内解决较大规模的问题。
地球物理反演[1，2]尝试找到最优解释观测数据并满足地质和地球物理限制的地下性质的模型。存在大量众所周知的地球物理反演方法。这些众所周知的方法属于两种类别即迭代反演和非迭代反演中的一种。以下是两种类别中每一种的普遍意义的定义：
非迭代反演——通过假设一些简单的背景模型和基于输入数据更新模型实现的反演。此方法不将更新模型用作另一反演步骤的输入。对于地震数据的情况，这些方法常被称为成象、偏移、绕射层析成象或博恩反演。
迭代反演——涉及重复改进地下性质模型从而建立理想地解释观测数据的模型的反演。如果反演收敛，则最终模型将更好地解释观测数据并将更紧密地接近实际的地下性质。迭代反演通常生成比非迭代反演更准确的模型，然而计算起来昂贵得多。
迭代反演总体上相对于非迭代反演是优选的，因为其生成更准确的地下参数模型。不幸地，迭代反演的计算费用如此昂贵，以致于将其应用于多种感兴趣的问题是不切实际的。这种高计算成本是源于所有反演技术均需要多个计算密集型模拟的事实。任何单独模拟的计算时间与所要反演的源数成比例，并且地球物理数据一般存在大量源，在此前述使用的术语源是指源设备的激活位置。迭代反演的问题加剧，因为必须计算的模拟数与反演的迭代数成比例，并且所需迭代数一般属于成百上千的级别。
用于地球物理学的最常用迭代反演方法是成本函数优化。成本函数优化包括成本函数S(M)的数值相对于模型M的迭代最小化或最大化，该成本函数S(M)是计算数据与观测数据之间错配度的度量（其有时也被称为目标函数），其中计算数据是通过计算机使用当前地球物理性质模型和源信号在由给定的地球物理性质模型表示的介质中的物理学支配传播（physics governing propagation）来模拟的。模拟计算可通过几种数值方法中任一种完成，该数值方法包括但不限于有限差分、有限元或射线追踪。模拟计算可在频率域或时间域中进行。
成本函数优化方法是局部或全局的[3]。全局方法简单地包括计算模型群体{M1，M2，M3,…}的成本函数S(M)和从大约最小化S(M)的群体选择具有一个或多个模型的组。如需进一步改进，则此新选定的模型组可被用作生成新模型群体的基础，该新模型群体可关于成本函数S(M)再次进行测试。对于全局方法，测试群体中的各模型可被认为是迭代，或在较高水平上各组所测群体可被认为是迭代。公知的全局反演方法包括蒙特卡罗（Monte Carlo）、模拟退火、遗传和进化算法。
不幸地，全局优化方法一般收敛极其缓慢，因此大多数地球物理反演是基于局部成本函数优化。算法1概括了局部成本函数优化。
算法1–进行局部成本函数优化的算法。.
此程序通过利用新更新的模型作为另一梯度搜索的起始模型进行迭代。该过程持续，直到发现理想地解释观测数据的更新模型。常用的局部成本函数反演方法包括梯度搜索、共轭梯度和牛顿法。
在声学近似中地震数据的局部成本函数优化是普遍的地球物理反演作业，并通常是其他类型的地球物理反演的示例。当在声学近似中反演地震数据时，成本函数可写作：
S ( M ) = Σ g = 1 N g Σ r = 1 N r Σ t = 1 N t W ( ψ calc ( M , r , t , w g ) - ψ obs ( r , t , w g ) )
（方程式1）
在此：
S=成本函数，
M=N个描述地下模型的参数（m1，m2，…mN）的向量，
g=道集（gather）指数，
wg=道集g的源函数，其是空间坐标和时间的函数，对于点源而言其是空间坐标的δ函数，
Ng=道集数，
r=道集中的接收器指数，
Nr=道集中的接收器数，
t=道（trace）中的时间样本指数，
Nt=时间样本数，
W=最小化标准函数（我们通常选择W（x）=x2，其是最小二乘方（L2）标准），
ψcalc=从模型M计算的地震压力数据，
ψobs=测量地震压力数据。
该道集可以是可在一轮地震正向建模程序中模拟的任何类型的道集。通常道集相应于地震炮点（shot），尽管炮点可比点源更普遍（general）。关于点源，道集指数g相应于各个点源的位置。关于平面波源，g将相应于不同的平面波传播方向。这种广义源（generalized source）数据ψobs可以在现场得到或可由利用点源得到的数据合成。另一方面，计算数据ψcalc通常可通过利用正向建模时广义源函数直接计算。关于多种类型的正向建模——包括有限差分建模，广义源所需的计算时间粗略地等同于点源所需的计算时间。
方程式（1）可简化为：
S ( M ) = Σ g = 1 N g W ( δ ( M , w g ) ) .    （方程式2）
其中，现暗指基于接收器和时间样本的总和，并且
δ(M,wg)=ψcalc(M,wg)-ψobs(wg) （方程式3）
反演试图更新模型M，以使S(M)是最小的。其可通过局部成本函数优化实现，该局部成本函数优化如下更新给定模型M(k)：
M ( k + 1 ) = M ( k ) - α ( k ) ▿ M S ( M )     （方程式4）
其中，k是迭代数，α是模型更新的标量大小，
Figure BDA00003285636000045
是关于模型参数取定的错配度函数的梯度。模型扰动或模型更新值通过目标函数梯度乘以步长α计算，其必须被反复计算。
基于方程式（2），可得到下列成本函数梯度方程式：
▿ M S ( M ) = Σ g = 1 N g ▿ M W ( δ ( M , g ) ) .     （方程式5）
因此，为计算成本函数的梯度，必须分别计算各道集对成本函数的贡献梯度，然后加和那些贡献。因此，计算
Figure BDA00003285636000046
所需的计算量（computational effort）是确定单个道集对梯度的贡献所需的计算量的Ng倍。对于地球物理学问题，Ng通常相应于地球物理源数，并且属于10,000至100,000的级别，大幅扩大了计算
Figure BDA00003285636000047
的成本。
注意，
Figure BDA00003285636000048
的计算需要计算关于N个模型参数中的每一个mi的导数W（δ）。由于对于地球物理学问题而言N通常很大（通常大于一百万），这种计算可能极其耗时——如果其必须对于每一个单独的模型参数进行。幸运地，可利用伴随（adjoint）法同时（at once）对所有模型参数有效地进行这种计算[1]。通过如下算法总结用于最小二乘方目标函数和网格模型参数化的伴随法：
Figure BDA00003285636000044
Figure BDA00003285636000051
算法2–利用伴随法计算网格模型的最小二乘方成本函数梯度的算法。
虽然利用伴随法计算梯度相对于其他方法有效，但是其仍十分昂贵。具体地，伴随法需要两种模拟，一种在时间上正向，一种在时间上反向，而对于地球物理学问题这些模拟通常是极其计算密集型的。而且，如上所述，这种伴随法计算必须分别对每一个测量数据道集进行，增加了Ng倍的计算成本。
所有反演种类的计算成本可通过自源组合反演数据而非单独地反演源而降低。这可被称为同时源反演。几种类型的源组合是已知的，包括：相干（coherently）加和近间距的源以生成有效源，该有效源生成一些所需形状（例如，平面波）的波前、加和宽间距的源、或完全或部分地叠加反演前的数据。
通过反演组合源实现的计算成本降低由于组合数据的反演通常产生不太准确的反演模型的事实而至少部分抵消。这种准确性的损失是缘于如下事实：当加和各个源时丢失信息，因此加和数据限制反演模型的强度不如未加和的数据。加和过程中的这种信息丢失通过在加和前编码各炮点记录被最小化。组合前进行编码显著保留同时源数据中的更多信息，因此更好地限制反演[4]。编码还允许组合近间距的源，因此允许给定计算区域的更多源组合。多种编码方案可用于这种技术，包括时移编码和随机相位编码。此背景部分的其余内容简要地评述多种已公开的地球物理同时源技术——编码和非编码的。
Van Manen[6]提出利用地震干涉测量方法加速正向模拟。地震干涉测量通过将源布置在目标区域的边界各处而工作。对这些源单独建模，并且记录需要格林函数（Green’s function）的所有位置的波场。任两个记录位置之间的格林函数然后可通过互相关在两个记录位置处得到的道并在所有边界源基础上加和而计算。如果待反演的数据在目标区域内具有很多的源和接收器（相对于在边界上具有一个或另一个），则这是计算所需格林函数十分有效的方法。但是，对于地震数据的情况，待反演数据的源和接收器均在目标区域中是很少的。因此，这种改进对于地震反演问题具有非常有限的应用性。
Berkhout[7]和Zhang[8]提出反演总体上可通过反演非编码同时源而改进，该非编码同时源在地下一些区域中被相干加和以生成一定程度上所需的波前。例如，点源数据可随时移而加和，该时移是源位置的线性函数，从而以相对于地表的一定特定角度生成下行平面波。这种技术可用于所有种类的反演。此方法的问题是，源道集的相干加和必定降低数据的信息量。因此，例如，加和生成平面波去除了地震数据中所有关于传播时间相对于源-接收器炮检距（source-receiver offset）的信息。该信息对于更新缓慢改变的背景速度模型是至关重要的，因此Berkhout法没有很好地进行限制。为克服这个问题，可反演数据的多个不同相干加和（例如，多个具有不同传播方向的平面波），但是然后由于反演成本与反演的不同加和数量成比例而失去效力。在本文中，这种相干加和的源被称为广义源。因此，广义源可以是点源或生成具有一定所需形状的波前的点源总和。
Van Riel[9]提出通过如下进行的反演：输入地震数据的非编码叠加或部分叠加（关于源-接收器炮检距），然后限定关于此叠加数据的成本函数，其将被优化。因此，此公开提出利用非编码的同时源改进基于成本函数的反演。与Berkhout[6]的同时源反演方法一样，通过此方法提出的叠加减少了待反演数据中的信息量，因此该反演比对于原始数据得到较不充分的限制。
Mora[10]提出反演数据，该数据是宽间距的源的总和。因此，此公开提出利用非编码的同时源模拟提高反演效力。加和宽间距的源具有这样的优势：比Berkhout提出的相干加和保留更多信息。但是，宽间距的源的加和意味着必须用于反演的孔隙（aperture）（反演的模型区域）必须增加以容纳所有宽间距的源。由于计算时间与此孔隙的面积成比例，Mora法没有产生如加和的源彼此接近时所可以实现的一样多的效力增加。
Ober[11]提出通过利用同时编码源加速地震偏移——非迭代反演的特殊情况。在测试不同编码方法后，Ober发现，所得偏移图像具有显著降低的信噪比，这是因为宽带编码函数必须仅大约正交的事实。因此，当加和大于16个炮点时，反演质量不理想。由于非迭代反演开始不是非常昂贵并且由于需要高信噪比反演，因此该技术在地球物理学工业中未被广泛实践。
Ikelle[12]提出通过同时模拟在不同时间间隔（在模拟中）激活的点源而快速正向模拟的方法。还讨论了将这些时移的同时源模拟数据解码返回单独的模拟的方法，该单独的模拟将已得自各个点源。然后这些解码数据可被用作任何常规反演程序的部分。Ikelle方法的问题是，提出的解码方法将生成单独的数据，其噪声水平与相邻源的数据之间的差异成比例。此噪声对于非横向恒定的地下模型将是显著的，例如包含倾斜反射器（dipping reflectors）的模型。此外，此噪声将与同时源数成比例增长。基于这些困难，如果用于反演非横向恒定的地下，Ikelle的同时源方法可导致不可接受水平的噪声。
Krebs等在PCT专利申请公开号WO2008/042081——本文将其被赋予的全部权限引入作为参考——中提出的源编码，是反演全波场数据非常成本有效的方法。（同时反演编码的道集的同样的方法将服务于接收器——通过源-接收器互易性或通过编码共源数据道集中的实际接收器位置）。对于固定接收器，单个有效源仅需进行正向和伴随计算；参见PCT专利申请公开号WO2009/117174，本文将其被赋予的全部权限引入作为参考。如果事实是针对一般的2D采集几何记录数百炮点和在3D考察的情况下记录数千炮点，此方法带来的计算保存（saving）相当大。在实践中，对于最常见的场数据采集几何，固定接收器假设不严格地有效。在海洋拖缆数据的情况下，每一个新炮点的源和接收器均发生移动。即使是在接收器位置固定的考察中，实践通常是：不是所有接收器都“监听”每一个炮点，并且正在监听的接收器可自炮点-至-炮点而不同。这也有悖于“固定接收器假设”。此外，由于逻辑问题，难以记录接近源的数据，并且这意味着一般缺少炮检距附近的数据。对于海洋和陆地考察也是这样的。这些因素均意味着，对于同时源道集，每一个接收器位置对于一些源炮点都将是漏测的数据。总之，在同时编码源反演中，对于给定的同时编码道集，需要每一个炮点的所有接收器位置处的数据，并且这可被称为同时编码源反演的固定接收器假设。在WO08/042081中，当固定接收器假设不理想时，一些公开的实施方式可比其它的好。因此，对同时编码源（和/或接收器）反演的直接应用进行调节或调整将是有利的，这在固定接收器假设被折衷时将增强其性能。本发明提供了如此实施的方式。Haber等[15]也描述了关于利用随机优化方法移动同时编码源反演中的接收器的问题的方法，并将其应用于直流电阻性问题。
发明概述
在一个实施方式中，本发明是利用同时源编码进行全波场反演以确定地下区域的物理性质模型的计算机执行方法，该全波场反演是来自不满足同时编码源反演固定接收器假设的考察的测量地球物理数据的全波场反演，该方法包括参考图2流程图描述的下列步骤——其中所有加和、模拟、计算和更新在计算机上进行：
（a）获得测量地球物理数据（40）的两个或更多个编码道集组，该编码道集在固定接收器假设无效的情况下获得，其中各道集与单个广义源或可选地与单个接收器相关，并且其中各道集用选自非同等编码函数（30）组的不同编码函数编码（60）；
（b）通过加和相应于单个接收器（60）或可选地单个源的各道集中的全部或选定的数据记录，加和组中的编码道集，并针对每个不同接收器或可选地每个不同源重复进行，生成同时编码道集（80）；或可选地在场考察（90）中获得编码数据，其中编码函数是测量源特征（50）；
（c）假设地下区域的物理性质模型（10），所述模型提供整个地下区域的位置处的至少一种物理性质的数值；
（d）利用假设的物理性质模型模拟数据的合成同时编码道集，其中模拟利用编码源特征（70），并且其中在单个模拟操作中模拟整个同时编码道集；
（e）计算度量测量数据同时编码道集与模拟同时编码道集之间互相关的目标函数（100）；
（f）通过优化目标函数（100）更新（110）物理性质模型；
（g）利用自前次迭代的更新的物理性质模型（20）作为步骤（c）的假设物理性质模型，再重复步骤（a）-（f）至少一次，生成进一步更新的物理性质模型；和
（h）在收敛（120）后，下载、显示或保存进一步更新的物理性质模型于计算机存储器（130）。
附图简述
由于专利法限制，一个或多个附图是彩色原图的黑白再现。彩色原图已提交于对应的美国申请中。具有彩色附图的此专利或专利申请公开副本可通过请求和支付必要的费用从美国专利商标局获得。
本发明及其优势将通过参考下文详细描述和附图而得到更好的理解，在附图中：
图1A-C是测试实例的数据，其中图1A显示应用中的“真实”速度模型，即，用于模拟合成考察数据的模型；图1B显示通过利用常规L2目标函数反演模拟考察数据获得的速度模型；和图1C显示通过利用本发明方法的互相关目标函数反演而获得的速度模型；
图2是显示本发明方法的一个示例性实施方式的基本步骤的流程图。
图3是显示间距20m的四个炮点的总和的模拟地震数据的截面；
图4是显示间距100m的四个炮点的总和的模拟地震数据的截面；和
图5A-5C显示模型研究的结果，显示本发明方法的实施方式利用互相关目标函数减少在迭代反演中影响收敛的局部极小化问题的优势。
本发明将结合实例实施方式进行描述。但是，在下文详述针对于本发明的具体实施方式或具体应用方面，这仅意为示例性，并不被解释为限制本发明的范围。相反，意为覆盖可包括在所附权利要求限定的本发明范围内的所有替代选择、修正和等同形式。
实例实施方式详述
本专利文本中提出的方法利用互相关目标函数反演编码同时源数据。经编码形成同时源数据的原始炮点道集具有漏测的近-和远-炮检距，这在海洋拖缆考察中是常见的。对于这种互相关目标函数如何用于非固定接收器数据的全波场反演（“FWI”）的直观理解可通过考虑成像问题而获得。在FWI中，反演一种或多种模式的地震波场，获得地球的物理性质模型。首先考虑同时源逆时偏移（“RTM”）的成像问题。
假设源于两个源的两个地下测量波场d1,d2，并且其通过如下编码，
Figure BDA00003285636000091
其中P是将数据投影到固定接收器几何所需的所有接收器上的投影算子。因此，如果所有炮点的算子P相同，则是固定接收器几何。编码源由如下给出，
S=α1S1+α2S2
将正向问题表示为Ld=S，其中L是偏微分算子，d是模拟响应，和S是编码源函数。如果将正向算子表示为L-1使得模拟数据d=L-1S并将反向传播算子表示为
Figure BDA00003285636000092
则通过互相关获得的图像由如下给出，
Figure BDA00003285636000093
其中P是使全部区域的计算数据映射至固定接收器组的投影算子。对于固定接收器情况，所有源的P相同，因为每个接收器监听所有源。
对于非固定接收器几何，假设Pk是第k个源到接收器子组上投影算子。则编码数据由如下给出，
Figure BDA00003285636000094
则图像由如下给出
如果采用几种编码实现方式，则第二项（second term）由于串音而消除并且其对图像的贡献减少，因而增强信噪比。成像与生成FWI梯度之间存在紧密的联系。如果考虑FWI的目标函数（通常被称为成本函数）为L2，即，最小二乘方，测量数据与模拟数据之间的范数（norm），然后L2目标函数的梯度，简单地是反向传播的数据剩余与模拟波场的互相关。对于固定接收器几何，这是良好的度量，因为特定接收器监听全部源。对于非固定接收器，这存在问题，因为编码测量数据相对于编码模拟数据是不足的。这种不足最终映射成模型假象。
但是，从接收器非固定的成像条件获得的认识可被认为是使用互相关目标函数进行FWI可有助于减轻此问题的线索。因此，替代L2，考虑互相关目标函数，其由如下给出：
Figure BDA00003285636000101
其中
Figure BDA00003285636000102
是非固定接收器的编码测量数据，暗示将接收器未监听源的道假设是零道。
Figure BDA00003285636000103
是来自同时源函数的模拟编码数据；这是固定接收器数据，因为所有源同时被爆破（fired），并且各接收器测量来自此同时源的信号。假设对于各炮点，假想地，我们能够分开考察中存在的接收器和漏测的接收器，从而可写成
Figure BDA00003285636000104
Figure BDA00003285636000105
贡献于
Figure BDA00003285636000106
的道即是测量数据中不可用的道，即，因此上述方程式中第二项的贡献相对于第一项很小。可将目标函数近似为，
Figure BDA00003285636000108
而且，如果不同的编码用于一些迭代，优选在每次迭代中，如Krebs等在PCT专利申请公开号WO2008/042081中教导，第二项最终减少。除每次迭代改变编码外，进一步的选择是在每次迭代应用编码组，该编码组被设计以获得目标函数的平均梯度，即，得自各编码的梯度在加和时给出平均梯度。因此，在每次迭代时，这种平均处理使串音减少——非常类似于同时源RTM成像的思想，在同时源RTM成像中，利用不同编码获得的图像叠加以消除串音噪声。梯度平均也在随机优化方法中被主张（Haber等，[15]）。
本发明方法的优选实施方式的总体目标是最大化——而非最小化——互相关目标函数。但是，可改写目标函数，例如，用其负数（negative）将其替代，在这样的方式下有效解决最小化问题。因此，本发明方法中此步骤的总体表述是优化目标函数。根据优化算法的选择，即，其解决最大化还是最小化算法，可选择目标函数，以使其保留对任何任意数据规模不敏感的互相关性的主要特征。互相关目标函数的优势是其使FWI对小波规模不敏感。有效地，波形的传播时间要素适于此目标函数。如果目标函数以相位是传播时间的函数的频率域撰写，则此可容易看到。.
合成实例
本发明方法的可行性通过利用图1A所示速度模型的合成实例得以证明。在该实例中，建立具有80m炮点间距的383个炮点，并且每20m布置接收器。为模拟非固定接收器几何，对于各炮点去除近和远炮检距，并采用正200m至3km的炮检距。如预期，近和远炮检距的漏测接收器数据有悖于固定接收器假设。那些炮点经编码形成测量编码数据。在非线性反演的每次迭代时，应用考察中全部连续炮点，并用不同的编码函将其编码8次，以形成8个同时源，从而计算平均梯度。取决于平均梯度中存在多少串音，可增加编码数量。此被执行用于L2目标函数和互相关目标函数。
利用带有时间窗的多分辨率法（multi-resolution method）实施反演。在多分辨率法中，数据通过频率带窗化，从而稳定化反演问题。一般，较低频率带的数据先反演，并且得自低频反演的模型用作下一个较高频率带的起始模型。该技术用于避免收敛于错误的解决方案——通常被称为局部极小化解决方案。图1B显示利用L2目标函数获得的模型。在模型的浅部和深部，假象均明确可见。图1C的模型是应用互相关目标函数的结果。模型清楚显示相对于图1B模型的显著改进，强化了互相关度量明显更适于处理非固定接收器几何同时源数据的观察结果。
应用互相关目标函数的反演在所增加的源记录间距过近时可遭遇局部极小化。进一步拉开炮点记录的间距解决了此问题，但使加入反演的信息量减少，因而降低最终反演模型的质量。解决此矛盾的方式是将所有可用的炮点记录分成几组，其中进一步拉开属于同组的炮点记录的间距。然后可在非线性反演的每次迭代中选择不同炮点组，直到所有组均已被使用。此过程可重复，直到收敛。这可进行的一种方式是利用使源特征与测量数据道集相乘的编码函数。所有不在某次迭代选定组中的数据记录的编码函数可被设定为等于零。事实上，可利用编码函数随机选择各组中的炮点，按照Boonyasiriwat and Schuster，2010表现的思想，这将导致炮点间距增加，其中，它们将这用于固定接收器传播和不同的目的。这可利用随机编码顺序[-1,0,+1]进行。这种增加炮点间距的技术可用于本发明方法，以减少非固定接收器的串音噪声。可设计其他实现较小炮点记录组的较大炮点间距的方法。
下文是应用上述策略的实例。
此方法的可行性通过合成实例得以证明。建立了383个具有20m炮点间距的炮点，并且每10m布置接收器。为模拟非固定接收器几何，各炮点采用正200m至3km的炮检距。那些炮点经编码形成测量编码数据。在非线性反演的每次迭代时，利用8个同时源计算平均梯度。图3显示4个近间距（在此为20m）炮点的总和，而图4显示进一步远离（在此为100m）的4个炮点的总和。基本思想是炮点间距增加使目标函数的串音要素减少；但是，过大间距导致数据空隙，其降低更新模型的能力。如上所述，这种情况可通过选择若干组形成编码炮点而减轻。各组具有明显隔离的炮点，但采用这些组中的几组填补了空间空隙，否则如果仅选择一组会保留该空间空隙。对于在此呈现的实例，每次迭代选择4组，并将2种实现方式用于各组，实现方式是用不同编码函数组进行同时源模拟。因此，每次迭代总计有8个编码同时炮点。为进行比较，每次迭代以8种实现方式用近间距编码炮点进行反演。图5A是此实例所用的真实模型。图5B显示在源近间距时同时源反演的结果。由于模型中存在串音，反演超过50次迭代也无法建立适当的更新。50次迭代的模型被作图于图5B。为进行比较，炮点明显隔离时50次迭代的模型被作图于图5C。反演过程中不存在收敛问题。图5C与5A之间比图5B与5A之间具有明显更好的一致性。
上述申请涉及本发明的具体实施方式，其目的是将本发明示例。但对于本领域技术人员明显的是，对本文所述实施方式的多种改动和变型是可以的。所有这种改动和变型均意为处于如所附权利要求限定的本发明范围内。本领域技术人员容易理解，在本发明的优选实施方式中，本发明方法的至少一些步骤在计算机中进行，即，本发明是计算机执行的。在这种情况下，所得更新物理性质模型可被下载、显示或保存于计算机存储器。
参考文献
1.Tarantola,A.,"Inversion of seismic reflection data in the acousticapproximation,"Geophysics49,1259-1266(1984).
2.Sirgue,L.,and Pratt G."Efficient waveform inversion and imaging:A strategyfor selecting temporal frequencies,"Geophysics69,231-248(2004).
3.Fallat,M.R.,Dosso,S.E.,“Geoacoustic inversion via local,global,and hybridalgorithms,”Journal of the Acoustical Society of America105,3219-3230(1999).
4.Hinkley,D.and Krebs,J.,“Gradient computation for simultaneous sourceinversion,”PCT Patent Application Publication No.WO2009/117174.
5.Krebs,J.R.,Anderson,J.A.,Neelamani,R.,Hinkley,D.,Jing,C.,Dickens,T.,Krohn,C.,Traynin,P.,“Iterative inversion of data from simultaneous geophysicalsources,”PCT Patent Application Publication No.WO2008/042081.
6.Van Manen,D.J.,Robertsson,J.O.A.,Curtis,A.,"Making wave by timereversal,"SEG International Exposition and75th Annual Meeting Expanded Abstracts,1763-1766(2005).
7.Berkhout,A.J.,"Areal shot record technology,"Journal of Seismic Exploration1,251-264(1992).
8.Zhang,Y.,Sun,J.,Notfors,C.,Gray,S.H.,Cherris,L.,Young,J.,"Delayed-shot3D depth migration,"Geophysics70,E21-E28(2005).
9.Van Riel,P.,and Hendrik,W.J.D.,“Method of estimating elastic andcompositional parameters from seismic and echo-acoustic data,”U.S.Patent No.6,876,928(2005).
10.Mora,P.,"Nonlinear two-dimensional elastic inversion of multi-offset seismicdata,"Geophysics52,1211-1228(1987).
11.Ober,C.C.,Romero,L.A.,Ghiglia,D.C.,"Method of Migrating SeismicRecords,"U.S.Patent No.6,021,094(2000).
12.Ikelle,L.T.,"Multi-shooting approach to seismic modeling and acquisition,"U.S.Patent No.6,327,537(2001).
13.Romero,L.A.,Ghiglia,D.C.,Ober,C.C.,Morton,S.A.,"Phase encoding ofshot records in prestack migration,"Geophysics65,426-436(2000).
14.Jing X.,Finn,C.J.,Dickens,T.A.,Willen,D.E.,"Encoding multiple shotgathers in prestack migration,"SEG International Exposition and70th Annual MeetingExpanded Abstracts,786-789(2000).
15.Haber,E.,Chung M.and Herrmann,“An effective method for parameterestimation with PDE constraints with multiple right hand sides,”Preprint-UBChttp://www.math.ubc.ca/～haber/pubs/PdeOptStochV5.pdf(2010).
16.Jerome R.Krebs,John E.Anderson,David Hinkley,Ramesh Neelamani,Sunwoong Lee,Anatoly Baumstein,and Martin-Daniel Lacasse,“Full-wavefieldseismic inversion using encoded sources,”Geophysics74-6,WCC177-WCC188(2009).
17.Boonyasiriwat,C.,and Schuster,G.,“3D Multisource full-waveform inversionusing dynamic random phase encoding,”SEG Expanded Abstracts29,1044-1049(2010).

Claims (27)

1.用于测量地球物理数据的同时编码源反演以确定地下区域的物理性质模型的计算机执行方法，所述测量地球物理数据在同时编码源反演的固定接收器假设无效的条件下获得，所述方法包括下列步骤，其中所有加和、模拟、计算和更新在经编程以执行它们的计算机中进行：
（a）获得具有所述测量地球物理数据的两个或更多个编码道集的组，所述测量地球物理数据在所述固定接收器假设无效的条件下获得，其中各道集与单个广义源相关，或可选地与单个接收器相关，并且其中各道集用选自非同等编码函数组的不同编码函数编码；
（b）通过加和相应于单个接收器或可选地相应于单个源的各道集中的所有或选定的数据记录，加和所述组中的编码道集，并对于每个不同的接收器或可选地对于每个不同的源进行反复，生成同时编码道集；
（c）假设所述地下区域的物理性质模型，所述模型提供整个所述地下区域中的位置处的至少一种物理性质的数值；
（d）利用所述假设物理性质模型模拟数据的合成同时编码道集，其中所述模拟应用编码源特征，并且其中整个同时编码道集在单个模拟操作中被模拟；
（e）计算度量所述测量数据的同时编码道集与所述模拟同时编码道集之间的互相关的目标函数；
（f）通过优化所述目标函数更新所述物理性质模型；
（g）利用来自前次迭代的所述更新物理性质模型作为步骤（c）所述的假设物理性质模型，再重复步骤（a）-（f）至少一次，生成进一步更新的物理性质模型；和
（h）下载、显示或保存所述进一步更新的物理性质模型于计算机存储器。
2.权利要求1所述的方法，其中所述目标函数
Figure FDA00003285635900012
由如下给出或数学上等同于：
Figure FDA00003285635900011
其中
Figure FDA00003285635900021
是测量数据的同时编码道集，而
Figure FDA00003285635900022
是包括所有接收器位置的模拟数据道的模拟同时编码道集，而无论在具体炮点中是活动还是非活动的，以及||...||表示选定的范数。
3.权利要求2所述的方法，进一步包括通过如下近似所述目标函数：
Figure FDA00003285635900023
其中
Figure FDA00003285635900024
是模拟同时编码道集，以使相应于在具体炮点中非活动的接收器的所述道在编码前被零替代。
4.权利要求1所述的方法，其中步骤（a）的道集所用的所述编码函数在至少一次所述迭代中被改变。
5.权利要求4所述的方法，其中两组或更多组编码函数被用于每次迭代，并被设计以获得所述目标函数的平均梯度。
6.权利要求1所述的方法，其中所有编码函数在每次迭代中被改变。
7.权利要求1所述的方法，进一步包括如步骤（a）所述获得至少一个另外的编码道集组，所述编码道集组具有所述测量地球物理数据的两个或更多个编码道集，和对于每个另外的组进行步骤（b），然后累积来自步骤（e）的所述相应计算目标函数，其中步骤（f）所述的所述物理性质模型更新通过最大化所述累积计算目标函数而确定。
8.权利要求1所述的方法，其中测量数据的所述编码道集通过如下被编码：利用对于所述道集选择的所述编码函数在时间上褶积来自所述道集的所有道。
9.权利要求1所述的方法，其中通过获得来自地球物理考察的数据的道集获得测量数据的所述两个或更多个编码道集，在所述地球物理考察中数据由多个同时运转、独特编码的源装置获得。
10.权利要求1所述的方法，其中所述测量地球物理数据来自所述地下区域的地震考察。
11.权利要求10所述的方法，其中广义震源是全部点源或全部平面波源。
12.权利要求1所述的方法，其中用于所述合成同时编码道集模拟的所述编码源特征采用与用于编码测量数据的所述同时编码道集相同的编码函数，或是通过利用与用于编码步骤（a）的相应测量道集相同的编码函数在时间上褶积测量或估测的源特征而作成的函数，其中所述测量地球物理数据包括每个源激活的测量或估测的源特征。
13.权利要求8所述的方法，其中所述编码函数是选自线性、随机相位、线性调频脉冲、改良线性调频脉冲、随机时移和频率不相关的相位编码的类型。
14.权利要求8所述的方法，其中所述编码函数对于一些源是一种类型，并且对于其他源是另一种类型。
15.权利要求1所述的方法，其中所述编码函数被优化以提高所述目标函数的质量。
16.权利要求1所述的方法，其中步骤（d）所述的模拟利用有限差分、有限元或有限体积模拟代码进行。
17.权利要求10所述的方法，其中所述物理性质模型是地震波速度、地震弹性参数、地震各向异性参数或地震滞弹性参数的模型。
18.权利要求1所述的方法，其中利用选自蒙特卡罗、模拟退火、遗传或进化算法的全局目标函数优化方法更新所述模型。
19.权利要求1所述的方法，其中利用选自梯度线搜索、共轭梯度或牛顿法的局部目标函数优化方法更新所述模型。
20.权利要求1所述的方法，其中最大化所述目标函数包括计算所述目标函数相对于所述物理性质模型的参数的梯度。
21.用于测量地球物理数据的同时编码源反演以确定地下区域的物理性质模型的计算机执行方法，所述测量地球物理数据在同时编码源反演的固定接收器假设无效的条件下获得，所述方法包括：
（a）从对所述固定接收器假设无效的所述地下区域的地球物理考察获得测量地球物理数据；
（b）利用编程计算机，通过迭代反演，反演所述测量数据，包括利用假设或更新的物理性质模型同时模拟代表多个考察源或可选地多个接收器的考察数据，其中所述模拟中的源或接收器特征被编码，生成地球物理数据的模拟同时编码道集，所述反演进一步包括，通过优化度量所述模拟同时编码道集与所述测量地球物理数据相应的同时编码道集之间的互相关的目标函数，确定用于下一次迭代的模型更新；和
（c）下载或显示更新的物理性质模型或将其保存于计算机存储器或数据存储器。
22.权利要求21所述的方法，其中一些或全部所述地球物理数据，以各数据组中源活动位置之间的间距与分组前的所述一些或全部所述地球物理数据相比增加的方式，被分成多个数据组，并且每次迭代的所述模拟同时编码道集相应于不同的数据组，并且所有组均在所述迭代的过程中被应用。
23.计算机程序产品，包括非临时性计算机可用介质，所述非临时性计算机可用介质中具有计算机可读程序代码，所述计算机可读程序代码适于被执行，以实施测量地球物理数据的同时编码源反演方法，从而确定地下区域的物理性质模型，所述测量地球物理数据在同时编码源反演的固定接收器假设无效的条件下获得，所述方法包括：
（a）输入具有所述测量地球物理数据的两个或更多个编码道集的组，其中各道集与单个广义源相关或可选地与单个接收器相关，并且其中各道集用选自非同等编码函数组的不同编码函数编码；
（b）通过加和相应于单个接收器或可选地相应于单个源的各道集中的所有数据记录，加和所述组中的所述编码道集，并对于每个不同的接收器或可选地对于每个不同的源进行反复，生成同时编码道集；
（c）输入所述地下区域的物理性质模型，所述模型提供整个所述地下区域中的位置处的至少一种物理性质的数值；
（d）利用假设物理性质模型模拟数据的合成同时编码道集，其中所述模拟应用编码源特征，并且其中整个同时编码道集在单个模拟操作中被模拟；
（e）计算度量所述测量数据的同时编码道集与所述模拟同时编码道集之间的互相关的目标函数；
（f）通过优化所述目标函数更新所述物理性质模型；和
（g）利用来自前次迭代的所述更新物理性质模型作为步骤（c）所述的假设物理性质模型，再重复步骤（a）-（f）至少一次，生成进一步更新的物理性质模型。
24.从地下区域生产烃的方法，包括：
（a）进行所述地下区域的地震考察，其中不满足同时编码源反演的固定接收器假设；
（b）获得所述地下区域的速度模型，所述速度模型通过包括如下的方法确定：
通过迭代反演，反演考察测量数据，所述迭代反演包括，利用假设或更新的速度模型同时模拟代表多个考察源或可选地多个接收器的考察数据，其中所述模拟中的源或接收器特征被编码，生成地球物理数据的模拟同时编码道集，所述反演进一步包括，通过优化度量所述模拟同时编码道集与所述测量数据相应的同时编码道集之间的互相关的目标函数，确定用于下一次迭代的速度模型更新；
（c）将井眼钻探到至少部分由利用（b）的更新速度模型完成的对所述地下区域结构的阐释而确定的所述地下区域中的层中；和
（d）从所述井眼生产烃。
25.权利要求1所述的方法，其中所述编码通过选定的编码函数将各道集和各源特征相乘而进行，并且其中在步骤（g）的每次迭代中，一个或多个编码函数被设定为等于零，以使不归零的所述数据记录相应于平均比编码前的所述测量地球物理数据相距更远的源激活位置，并且随着所述迭代进行所有数据记录均被应用，即，不归零。
26.权利要求25所述的方法，其中所述编码函数用于随机选择每次迭代中的所述炮点，其导致炮点间距增加。
27.权利要求1所述的方法，其中所述被反演的测量地球物理数据是全波场地震数据。
CN201180058068.XA 2010-12-01 2011-09-01 利用互相关目标函数进行的海洋拖缆数据同时源反演 Expired - Fee Related CN103238158B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US41869410P 2010-12-01 2010-12-01
US61/418,694 2010-12-01
US201161509904P 2011-07-20 2011-07-20
US61/509,904 2011-07-20
PCT/US2011/050209 WO2012074592A1 (en) 2010-12-01 2011-09-01 Simultaneous source inversion for marine streamer data with cross-correlation objective function

Publications (2)

Publication Number Publication Date
CN103238158A true CN103238158A (zh) 2013-08-07
CN103238158B CN103238158B (zh) 2016-08-17

Family

ID=46163017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180058068.XA Expired - Fee Related CN103238158B (zh) 2010-12-01 2011-09-01 利用互相关目标函数进行的海洋拖缆数据同时源反演

Country Status (11)

Country Link
US (1) US8688381B2 (zh)
EP (1) EP2646944A4 (zh)
KR (1) KR101797451B1 (zh)
CN (1) CN103238158B (zh)
AU (1) AU2011337143B2 (zh)
BR (1) BR112013008503A2 (zh)
CA (1) CA2815054C (zh)
MY (1) MY160148A (zh)
RU (1) RU2587498C2 (zh)
SG (1) SG189850A1 (zh)
WO (1) WO2012074592A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646615A (zh) * 2016-12-29 2017-05-10 中国石油天然气集团公司 一种面波频散曲线的数据处理方法及装置
CN113484914A (zh) * 2021-07-13 2021-10-08 中海石油(中国)有限公司 海上风浪一致性影响量板制作方法、系统、介质及设备

Families Citing this family (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112012009154A2 (pt) * 2009-10-23 2016-08-16 Exxonmobil Upstream Res Co método para melhorar um modelo geológico de uma região de subsuperfície, produto de programa de computador, e, método para controlar hidrocarbonetos em uma região de subsuperfície
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8892413B2 (en) 2011-03-30 2014-11-18 Exxonmobil Upstream Research Company Convergence rate of full wavefield inversion using spectral shaping
US9176930B2 (en) * 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
US9453928B2 (en) 2012-03-06 2016-09-27 Westerngeco L.L.C. Methods and computing systems for processing data
US20130242693A1 (en) * 2012-03-13 2013-09-19 Seoul National University R&Db Foundation Seismic imaging system using a reverse time migration algorithm
US10317548B2 (en) 2012-11-28 2019-06-11 Exxonmobil Upstream Research Company Reflection seismic data Q tomography
GB2510872A (en) * 2013-02-15 2014-08-20 Total Sa Method of modelling a subsurface volume
GB2510873A (en) 2013-02-15 2014-08-20 Total Sa Method of modelling a subsurface volume
US10360317B2 (en) * 2013-03-15 2019-07-23 Chevron U.S.A. Inc. Beam inversion by Monte Carlo back projection
SG11201508195PA (en) 2013-05-24 2015-12-30 Exxonmobil Upstream Res Co 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
US10871584B2 (en) 2013-06-17 2020-12-22 Westerngeco L.L.C. Seismic data processing
US9702998B2 (en) 2013-07-08 2017-07-11 Exxonmobil Upstream Research Company Full-wavefield inversion of primaries and multiples in marine environment
US20150019180A1 (en) * 2013-07-12 2015-01-15 Jonathan Liu Model Replacement in a Local Region by Inversion
WO2015026451A2 (en) 2013-08-23 2015-02-26 Exxonmobil Upstream Research Company Simultaneous sourcing during both seismic acquisition and seismic inversion
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
US10739478B2 (en) 2013-11-12 2020-08-11 Westerngeco L.L.C. Seismic data processing
US9910189B2 (en) 2014-04-09 2018-03-06 Exxonmobil Upstream Research Company Method for fast line search in frequency domain FWI
US10670752B2 (en) 2014-04-14 2020-06-02 Cgg Services Sas Method for iterative inversion of data from non-encoded composite sources
SG11201608175SA (en) * 2014-05-09 2016-11-29 Exxonmobil Upstream Res Co Efficient line search methods for multi-parameter full wavefield inversion
WO2015184549A1 (en) * 2014-06-03 2015-12-10 Mtt Innovation Incorporated Efficient, dynamic, high contrast lensing with applications to imaging, illumination and projection
US10185046B2 (en) 2014-06-09 2019-01-22 Exxonmobil Upstream Research Company Method for temporal dispersion correction for seismic simulation, RTM and FWI
AU2015280633B2 (en) 2014-06-17 2018-07-19 Exxonmobil Upstream Research Company Fast viscoacoustic and viscoelastic full-wavefield inversion
US10838092B2 (en) 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10422899B2 (en) * 2014-07-30 2019-09-24 Exxonmobil Upstream Research Company Harmonic encoding for FWI
US9921324B2 (en) 2014-08-13 2018-03-20 Chevron U.S.A. Inc. Systems and methods employing upward beam propagation for target-oriented seismic imaging
US20160061986A1 (en) * 2014-08-27 2016-03-03 Schlumberger Technology Corporation Formation Property Characteristic Determination Methods
US10386511B2 (en) 2014-10-03 2019-08-20 Exxonmobil Upstream Research Company Seismic survey design using full wavefield inversion
EP3210050A1 (en) 2014-10-20 2017-08-30 Exxonmobil Upstream Research Company Velocity tomography using property scans
US10359532B2 (en) 2014-12-10 2019-07-23 Schlumberger Technology Corporation Methods to characterize formation properties
US11163092B2 (en) 2014-12-18 2021-11-02 Exxonmobil Upstream Research Company Scalable scheduling of parallel iterative seismic jobs
US10520618B2 (en) 2015-02-04 2019-12-31 ExxohnMobil Upstream Research Company Poynting vector minimal reflection boundary conditions
US10317546B2 (en) 2015-02-13 2019-06-11 Exxonmobil Upstream Research Company Efficient and stable absorbing boundary condition in finite-difference calculations
MX2017007988A (es) 2015-02-17 2017-09-29 Exxonmobil Upstream Res Co Proceso de inversion de campo ondulatorio completo de multifase que genera un conjunto de datos libres de multiples.
CA2985738A1 (en) 2015-06-04 2016-12-08 Exxonmobil Upstream Research Company Method for generating multiple free seismic images
US10838093B2 (en) 2015-07-02 2020-11-17 Exxonmobil Upstream Research Company Krylov-space-based quasi-newton preconditioner for full-wavefield inversion
CN108139499B (zh) 2015-10-02 2020-02-14 埃克森美孚上游研究公司 Q-补偿的全波场反演
MX2018003495A (es) 2015-10-15 2018-06-06 Exxonmobil Upstream Res Co Apilados angulares de dominio de modelo de fwi con conservacion de amplitud.
CN105572742B (zh) * 2015-12-21 2018-08-10 中国石油天然气集团公司 一种确定海水深度的方法和装置
US10768324B2 (en) 2016-05-19 2020-09-08 Exxonmobil Upstream Research Company Method to predict pore pressure and seal integrity using full wavefield inversion
EP3458883A1 (en) * 2016-05-20 2019-03-27 Exxonmobil Research And Engineering Company Shape-based geophysical parameter inversion
WO2018013257A1 (en) * 2016-07-13 2018-01-18 Exxonmobil Upstream Research Company Joint full wavefield inversion of p-wave velocity and attenuation using an efficient first order optimization
US10871585B2 (en) * 2016-08-03 2020-12-22 Harris Corporation System for processing seismic data based upon linear optimization and related methods
EP3413092B1 (en) 2017-06-08 2022-06-01 Total Se Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product
AU2018297693A1 (en) * 2017-07-06 2019-12-12 Chevron U.S.A. Inc. System and method for full waveform inversion of seismic data
US11656377B2 (en) * 2018-03-30 2023-05-23 Cgg Services Sas Visco-acoustic full waveform inversion of velocity and Q
US11231516B2 (en) 2018-05-15 2022-01-25 Exxonmobil Upstream Research Company Direct migration of simultaneous-source survey data
US11372123B2 (en) 2019-10-07 2022-06-28 Exxonmobil Upstream Research Company Method for determining convergence in full wavefield inversion of 4D seismic data
CN112698389B (zh) * 2019-10-22 2024-02-20 中国石油化工股份有限公司 一种地震资料反演成像方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1947032A (zh) * 2004-03-17 2007-04-11 维斯特恩格科地震控股有限公司 海上地震测量方法和系统
US20070274155A1 (en) * 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
CN101233527A (zh) * 2005-07-27 2008-07-30 埃克森美孚上游研究公司 与从地下岩层抽取碳氢化合物相关的井建模
US20100018718A1 (en) * 2006-09-28 2010-01-28 Krebs Jerome R Iterative inversion of data from simultaneous geophysical sources
US20100088035A1 (en) * 2008-10-06 2010-04-08 Etgen John T Pseudo-analytical method for the solution of wave equations

Family Cites Families (142)

* 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
US4562540A (en) 1982-11-12 1985-12-31 Schlumberger Technology Corporation Diffraction tomography system and methods
FR2543306B1 (fr) 1983-03-23 1985-07-26 Elf Aquitaine Procede et dispositif pour l'optimisation des donnees sismiques
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
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
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
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
US5920828A (en) 1997-06-02 1999-07-06 Baker Hughes Incorporated Quality control seismic data processing system
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
US6549854B1 (en) * 1999-02-12 2003-04-15 Schlumberger Technology Corporation 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
DE69932932D1 (de) 1999-10-22 2006-10-05 Jason Geosystems B V Verfahren zur Bestimmung der elastischen Parameter und Felszusammensetzung von unterirdischen Formationen mit Hilfe von seismischen Daten
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
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
US6903999B2 (en) 2000-01-21 2005-06-07 Schlumberger Technology Corporation System and method for estimating seismic material properties
DE60112895D1 (de) 2000-01-21 2005-09-29 Schlumberger Holdings System und verfahren seismischer wellenfeldtrennung
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
AU1176802A (en) 2000-10-17 2002-04-29 Westerngeco Llc Method of using cascaded sweeps for source coding and harmonic cancellation
WO2002047011A1 (en) 2000-12-08 2002-06-13 Ortoleva Peter J 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
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
WO2004034088A2 (en) 2002-10-04 2004-04-22 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
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
WO2004095072A2 (en) 2003-03-27 2004-11-04 Exxonmobil Upstream Research Company Method to convert seismic traces into petrophysical property logs
US7072767B2 (en) 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
CA2520640C (en) 2003-04-01 2012-10-23 Exxonmobil Upstream Research Company Shaped high frequency vibratory source
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
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
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
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
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
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
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
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
CA2663662C (en) 2006-09-13 2016-07-05 Exxonmobil Upstream Research Company Rapid inversion of electromagnetic reconnaissance survey data
WO2008087505A2 (en) 2007-01-20 2008-07-24 Spectraseis Ag Time reverse reservoir localization
US7636275B2 (en) * 2007-02-06 2009-12-22 Conocophillips Company Direct time lapse inversion of seismic data
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
US20090164186A1 (en) 2007-12-20 2009-06-25 Bhp Billiton Innovation Pty Ltd. Method for determining improved estimates of properties of a model
EP2238474A4 (en) 2008-01-08 2018-06-20 Exxonmobil Upstream Research Company Spectral shaping inversion and migration of seismic data
US8577660B2 (en) 2008-01-23 2013-11-05 Schlumberger Technology Corporation Three-dimensional mechanical earth modeling
WO2009117174A1 (en) 2008-03-21 2009-09-24 Exxonmobil Upstream Research Company An efficient method for inversion of geophysical data
EP2105765A1 (en) 2008-03-28 2009-09-30 Schlumberger Holdings Limited Simultaneous inversion of induction data for dielectric permittivity and electric conductivity
AU2009229187C1 (en) 2008-03-28 2014-01-23 Exxonmobil Upstream Research Company Surface wave mitigation in spatially inhomogeneous media
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)
US20120095690A1 (en) * 2008-08-01 2012-04-19 Higginbotham Joseph H Methods and computer-readable medium to implement inversion of angle gathers for rock physics reflectivity attributes
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
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
US9052410B2 (en) 2009-02-12 2015-06-09 Conocophillips Company Multiple seismic signal inversion
US8352190B2 (en) 2009-02-20 2013-01-08 Exxonmobil Upstream Research Company Method for analyzing multiple geophysical data sets
US9075163B2 (en) 2009-04-17 2015-07-07 Westerngeco L.L.C. Interferometric seismic data processing

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1947032A (zh) * 2004-03-17 2007-04-11 维斯特恩格科地震控股有限公司 海上地震测量方法和系统
CN101233527A (zh) * 2005-07-27 2008-07-30 埃克森美孚上游研究公司 与从地下岩层抽取碳氢化合物相关的井建模
US20070274155A1 (en) * 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
US20100018718A1 (en) * 2006-09-28 2010-01-28 Krebs Jerome R Iterative inversion of data from simultaneous geophysical sources
US20100088035A1 (en) * 2008-10-06 2010-04-08 Etgen John T Pseudo-analytical method for the solution of wave equations

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王天禧: "石油地震勘探数据叠前速度反演稳定性初步分析", 《中国海上油气（地质）》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646615A (zh) * 2016-12-29 2017-05-10 中国石油天然气集团公司 一种面波频散曲线的数据处理方法及装置
CN106646615B (zh) * 2016-12-29 2018-12-25 中国石油天然气集团公司 一种面波频散曲线的数据处理方法及装置
CN113484914A (zh) * 2021-07-13 2021-10-08 中海石油(中国)有限公司 海上风浪一致性影响量板制作方法、系统、介质及设备
CN113484914B (zh) * 2021-07-13 2023-09-12 中海石油(中国)有限公司 海上风浪一致性影响量板制作方法、系统、介质及设备

Also Published As

Publication number Publication date
US8688381B2 (en) 2014-04-01
RU2587498C2 (ru) 2016-06-20
AU2011337143B2 (en) 2016-09-29
CA2815054A1 (en) 2012-06-07
CN103238158B (zh) 2016-08-17
EP2646944A1 (en) 2013-10-09
MY160148A (en) 2017-02-28
KR101797451B1 (ko) 2017-11-14
BR112013008503A2 (pt) 2016-08-16
SG189850A1 (en) 2013-06-28
US20120143506A1 (en) 2012-06-07
WO2012074592A1 (en) 2012-06-07
EP2646944A4 (en) 2017-02-22
AU2011337143A1 (en) 2013-06-20
KR20130121895A (ko) 2013-11-06
RU2013129772A (ru) 2015-01-10
CA2815054C (en) 2017-05-16

Similar Documents

Publication Publication Date Title
CN103238158B (zh) 利用互相关目标函数进行的海洋拖缆数据同时源反演
CN103119552B (zh) 同时源编码和源分离作为全波场反演的实际解决方案
Plessix Three-dimensional frequency-domain full-waveform inversion with an iterative solver
CN104204857B (zh) 正交源和接收器编码
Brossier et al. Which data residual norm for robust elastic frequency-domain full waveform inversion?
CN101952743B (zh) 地球物理数据反演的有效方法
Dou et al. Full-wavefield inversion of surface waves for mapping embedded low-velocity zones in permafrost
CN101681394B (zh) 来自并发地球物理源的数据的迭代反演
Xing et al. Modeling frequency‐independent Q viscoacoustic wave propagation in heterogeneous media
CN103703391A (zh) 使用频谱整形的全波场反演的收敛速度
MXPA06012781A (es) Inversion de forma de onda llena en pila previa en 3-d.
CN102892972A (zh) 地球物理数据的迭代反演方法中的伪迹减少
CN103119472B (zh) 利用同时和顺序源方法进行全波形反演的混合方法
Shen Early-arrival waveform inversion for near-surface velocity estimation
Aragao et al. Elastic full waveform inversion with probabilistic petrophysical clustering
van Leeuwen et al. Enabling affordable omnidirectional subsurface extended image volumes via probing
Almomin Tomographic full waveform inversion
Chen et al. 3‐DS Wave Imaging via Robust Neural Network Interpolation of 2‐D Profiles From Wave‐Equation Dispersion Inversion of Seismic Ambient Noise
Galuzzi Modelling and optimization techniques for acoustic Full Waveform Inversion in seismic exploration
Huang Modeling and inversion of seismic data using multiple scattering, renormalization and homotopy methods
Landa et al. Exploding reflectors revisited: 3D heterogeneous multiscale elastic media.
Hwang et al. Full-waveform inversion strategies using common-receiver gathers for ocean-bottom cable data
Liu Full-Waveform Inversion of Time-Lapse Seismic Data Using Physics-Based and Data-Driven Techniques
Izzatullah et al. Physics‐reliable frugal local uncertainty analysis for full waveform inversion
Kurzmann et al. Seismic applications of full waveform inversion

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20180901

CF01 Termination of patent right due to non-payment of annual fee