CN108139498A - 具有振幅保持的fwi模型域角度叠加 - Google Patents
具有振幅保持的fwi模型域角度叠加 Download PDFInfo
- Publication number
- CN108139498A CN108139498A CN201680059382.2A CN201680059382A CN108139498A CN 108139498 A CN108139498 A CN 108139498A CN 201680059382 A CN201680059382 A CN 201680059382A CN 108139498 A CN108139498 A CN 108139498A
- Authority
- CN
- China
- Prior art keywords
- subset
- computer
- reflectivity
- angle
- bandwidth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004321 preservation Methods 0.000 title description 6
- 238000000034 method Methods 0.000 claims abstract description 61
- 238000002310 reflectometry Methods 0.000 claims abstract description 38
- 230000008569 process Effects 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 10
- 229930195733 hydrocarbon Natural products 0.000 claims description 7
- 150000002430 hydrocarbons Chemical class 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000005553 drilling Methods 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 230000001052 transient effect Effects 0.000 claims description 2
- 230000000149 penetrating effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001881 scanning electron acoustic microscopy Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
-
- 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/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
-
- 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/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Mining & Mineral Resources (AREA)
- Business, Economics & Management (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- Entrepreneurship & Innovation (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Geophysics And Detection Of Objects (AREA)
- General Engineering & Computer Science (AREA)
Abstract
一种方法,包括:获得根据预定的地下反射角范围被分成子集的地震数据集;利用计算机分别对每个子集执行声波全波场反演过程以对密度进行反演并生成相应的密度模型;使用相应的密度模型,作为反射角的函数为每个子集生成声阻抗;以及使用计算机将每个子集的声阻抗变换成反射率剖面,其中所述变换包括通过它们相应的带宽对反射率剖面进行归一化。
Description
相关申请的交叉引用
本申请要求于2015年10月15日提交的标题为“FWI MODEL DOMAIN ANGLE STACKSWITH AMPLITUDE PRESERVATION”的美国临时专利申请62/241,780的权益,该申请的全部内容通过引用被结合于此。
技术领域
本文描述的示例性实施例涉及地球物理勘探领域,并且更特别地涉及地球物理数据处理。具体而言,本文描述的实施例涉及用于更高效地生成FWI模型域角度叠加(anglestack)的方法。
背景技术
本部分旨在介绍可以与本发明的示例性实施例相关联的现有技术的各个方面。这种讨论被认为有助于提供框架以便于更好地理解本发明的特定方面。因此,应当理解的是,本部分应该以此来阅读,而不必作为对现有技术的承认。
地震勘探的重要目标是准确地对通常称为反射体的地下结构进行成像。通过在地震勘测的执行期间获得原始地震数据有利于地震勘探。在地震勘测期间,地震能量通过例如受控***在地面生成,并传递到地球。地震波由地下结构反射,并被许多称为地震检波器的传感器接收。由地震检波器接收到的地震数据被处理以便创建地下环境的准确映射。处理后的数据然后被检查,以识别可能含有油气(hydrocarbon)的地质构造。
全波场反演(FWI)是用于估计地下属性(诸如速度或密度)的地球物理学方法。众所周知,与传统方法相比,它具有更高的分辨率和更准确的物理特性。FWI算法的基本组成部分可以被描述如下:使用起始地下物理属性模型,通过使用数值方案(例如,有限差分、有限元法等)求解波动方程来生成合成地震数据。将合成地震数据与现场地震数据进行比较,并使用两者之间的差计算目标函数的值。为了最小化目标函数,生成修改后的地下模型,其用于模拟一组新的合成地震数据。将这组新的合成地震数据与现场数据进行比较以重新计算目标函数的值。通过使用新的更新模型作为寻找另一个搜索方向的起始模型来迭代目标函数优化过程,然后将其用于干扰模型以便更好地解释观察到的数据。这个过程继续,一直到发现能够令人满意地解释观察到的数据的更新模型。可以使用全局或局部优化方法来最小化目标函数并更新地下模型。常用的局部目标函数优化方法包括但不限于梯度搜索法、共轭梯度法、准牛顿法、高斯牛顿法和牛顿法。常用的全局方法包括但不限于蒙特卡罗或网格搜索。
虽然FWI被期望提供地下属性,但是难以直接用FWI从地震数据中提取正确的粘弹性属性。由于FWI通过将数据与合成波形进行拟合来估计属性,因此它依赖于波动方程可以如何准确地解释实际物理特性,以及优化方法可以多好地将影响与不同的属性分开。当使用声波方程时,FWI可以基于数据集中的行进时间信息来生成P-波速度模型。但是,由于真实的地球是粘弹性的,因此振幅信息没有得到充分利用,并且声学模型无法解释所获取的数据中的所有振幅。如果FWI被期望提供可解释的产品,比如弹性阻抗,那么通常需要弹性模拟,但非常昂贵;一般而言,它是声学FWI计算的6到10倍。此外,由于有限的剪切波运动学信息以及采集中常常较差的信噪比,因此难以获得剪切波速度的初始模型。
使用弹性振幅信息的替代方法是形成角度叠加。可以对角度叠加执行振幅与角度(AVA)分析[5]以提取弹性属性。用Kirchhoff迁移生成的传统AVA叠加需要几何扩展校正来解决传播期间的振幅损失。但是,不能保证校正后的振幅将反映数据的真实振幅。此外,Kirchhoff迁移基于有利于平滑速度模型并可能在高对比度的介质中失败的射线追踪。角度计算是在一维(1D)假设下进行的,当地下结构复杂时,这种假设不够准确。基于逆时间迁移(RTM)的角度叠加[1,2]对于利用高分辨率速度模型更为先进。尽管如此,振幅保持仍然困难。Yu Zhang等人(2014)[3]报告了平衡图像振幅的最小二乘RTM;但是,它尚未被证明能够生成角度叠加。
发明内容
一种方法,包括:获得根据预定的地下反射角范围被分成子集的地震数据集;利用计算机分别对每个子集执行声波全波场反演过程以对密度进行反演并生成相应的密度模型;使用相应的密度模型,作为反射角的函数为每个子集生成声阻抗;以及使用计算机将每个子集的声阻抗变换成反射率剖面(reflectivity section),其中所述变换包括通过它们相应的带宽对反射率剖面进行归一化。
在该方法中,每个全波场反演过程从相同的速度模型开始。
在该方法中,每个全波场反演过程被独立地应用到子集。
在该方法中,所述获得包括通过使用包括反射体倾斜角和P-波速度的信息的数据屏蔽(data mask)将炮集(shot gather)划分为子集。
该方法还可以包括对于每个反射率剖面使用傅里叶变换、离散傅里叶变换或快速傅立叶变换以计算在相同位置处应用于全部反射率剖面的至少一个局部窗口内的平均频谱,并确定每个平均频谱的带宽。
在该方法中,所述确定带宽是基于10-dB点之间的距离。
在该方法中,所述确定带宽是基于具有最陡斜率的点之间的距离。
在该方法中,平均频谱在多个局部窗口内被计算,并且被平均。
该方法还可以包括确定多个角度的反射率值并通过内插构建角度与振幅的曲线。
该方法还可以包括使用反射率剖面管理油气生产。
在该方法中,管理油气生产包括在至少部分地由反射率剖面确定的位置处钻井。
一种利用指令编码的非瞬态计算机可读存储介质,当所述指令由计算机执行时使所述计算机实现一种方法,所述方法包括:获得根据预定的地下反射角范围被分成子集的地震数据集;利用计算机分别对每个子集执行声波全波场反演过程以对密度进行反演并生成相应的密度模型;使用相应的密度模型,作为反射角的函数为每个子集生成声阻抗;以及使用计算机将每个子集的声阻抗变换成反射率剖面,其中所述变换包括通过它们相应的带宽对反射率剖面进行归一化。
附图说明
虽然本公开容许各种修改和替代形式,但是其具体示例实施例已经在附图中示出并且在本文中进行了详细描述。但是,应当理解的是,本文对具体示例实施例的描述不旨在将本公开限制于本文所公开的特定形式,而是相反,本公开要涵盖由所附权利要求限定的所有修改和等同物。还应当理解的是,附图不一定按比例绘制,而是将重点放在清楚地图示本发明的示例性实施例的原理上。此外,某些维度可能被夸大以帮助在视觉上传达这些原理。
图1图示了用于生成FWI AVA叠加的示例性方法。
图2图示了由于散射角度引起的频谱扩展效应。
图3A图示了单炮集。
图3B图示了被静音成五个不同角度范围的单炮集。
图4A图示了具有声学FWI的角度叠加。
图4B图示了具有卷积的角度叠加。
图5图示了来自图4A和图4B中的每个角度叠加的垂直线,其被重叠以显示一致性。
具体实施方式
本文描述了示例性实施例。但是,就以下描述特定于特定实施例而言,这仅仅是出于示例性目的,并且只是提供示例性实施例的描述。因此,本发明不限于下面描述的具体实施例,而是,它包括落入所附权利要求的真实精神和范围内的所有替代、修改和等同物。
本文描述的示例性实施例提供了一种方法,该方法可以:1)对于复杂的地质情况是稳健的;2)保持振幅与角度信息;以及3)比弹性FWI更便宜。所提出的FWI模型域角度叠加可以通过对不同声学模型的不同角度范围的数据集进行反演来生成。通过数据拟合过程可以实现振幅保持,并且使用坡印廷(Poynting)向量角度计算更准确。坡印廷向量描述了在各向同性和各向异性介质中体波、界面波、导波和非均匀波的能量流。坡印廷向量自然地利用了高分辨率FWI速度模型。当作为FWI工作流程的组成部分实现时,目前的技术优势可以在不改变平台的情况下利用FWI产品角度叠加。更重要的是,示例性方法不限于建模引擎中的不完整物理特性。
本技术进步的示例性实施例使用FWI生成模型域振幅保持的角度叠加。有利地,目前的技术进步可以仅使用声学模拟,但是可以应用于所获取的地震数据的完整偏移。使用一个声学模型来拟合包含各种物理特性的所有数据是不可能的。但是,如果数据集通过反射角度分开,那么对于每个角度,都存在可以解释数据的声学模型。在所有模型组合的情况下,阻抗模型被形成为反射角的函数:I(θ)。声阻抗是地震能量穿过地下环境的特定部分的容易程度的量度。本领域的普通技术人员将认识到的是,声阻抗可以被定义为密度和地震速度的乘积。根据阻抗,反射率可以作为角度的函数被推导出来:R(θ),其正好是AVA的定义。
实际上,没有必要找到R(θ)的连续形式。相反,R(θ)可以在多个角度被确定,并且AVA曲线可以通过内插重建。
图1图示了用于生成FWI AVA叠加的示例性方法。在步骤101中,获得处理后的数据,其可以是从收集到的地震数据中生成的单炮集。这种处理后的数据是根据本领域普通技术人员已知的常规技术调整的地震数据。在步骤102中,炮集被划分成若干子集。每个子集都在相对小的反射角范围内。这可以通过使用精心设计的数据屏蔽来实现,这些数据屏蔽包括反射器倾角和P-波速度的信息。利用倾角和P-波速度,可以使用射线追踪方法来寻找在来自每个反射角的每个地下点的数据中的反射波的时间和偏移。数据屏蔽是隐藏原始数据的一部分使得只有期望的部分(即,期望角度范围)可用于后续处理的过程。附图中示出的角度范围仅仅是示例性的,并且可以使用其它角度范围。此外,角度范围的数量也是示例性的,因为炮集可以被划分成更多或更少的剖面。
在步骤103中,在步骤102中生成的每个数据子集上,应用声学FWI以独立地获得声阻抗。本领域的普通技术人员熟悉声学FWI,并且该过程的进一步细节被省略。所有的反演可以从相同的速度模型开始,并且假定速度模型建立已经完成并且足够准确,因此运动学在这个过程中不会被更新。模型更新仅用于解释数据振幅。在反演之后,用这些阻抗模拟的合成波形很好地拟合了真实数据,使得振幅信息被保持在模型域中。每个阻抗模型只能解释由数据屏蔽约束的一定范围的反射角的数据。在一个角度范围内,可以选择中间角度作为反射率的标称角度。在反演中形成梯度时,通过使用坡印廷向量[2],可以附加地保证这一点,因此梯度对标称角度的反射最为敏感。在有限差分模拟和梯度计算的过程中,使用坡印廷向量来分离波传播方向。数据分离可以基于射线理论进行。由于坡印廷向量基于波动理论,因此它可以帮助检查数据分离的准确性。
如本文使用的术语速度模型、密度模型或物理属性模型是指数字的数组,通常是3-D数组,其中可以被称为模型参数的每个数字是速度、密度或小区中其它物理属性的值,其中地下区域在概念上被划分为离散的小区以用于计算目的。
在反演之后,密度和速度两者都是已知的,并且步骤104可以根据声学FWI的结果来确定阻抗,因为阻抗是密度和速度的函数。
在声学FWI中,两个参数可以被反演以拟合数据振幅:P-波速度和密度。当使用L-2标准类型的目标函数时,通常选择P-波速度来拟合振幅和行进时间。但是,密度可能是用于反射率反演的更好的参数。密度具有比P-波速度简单得多的AVA响应。如在Aki-Richards方程中所描述的:
其中Δρ是密度扰动并且Δα是P-波速度扰动,密度具有不随角度变化的恒定AVA响应[6]。这表明,当密度扰动被反演以在某个角度θ拟合数据振幅时,扰动的值直接表示在该角度θ的反射率,而不管实际θ的值如何。相反,如果使用P-波速度扰动,那么为了获得R(θ),我们需要应用的校正。此外,方程[1]仅在扰动弱时是有效的,并且当扰动强时,校正不具有明确的形式。对于密度,恒定的AVA响应对于所有情况都是有效的。
在获得声阻抗之后,它们可以在步骤105中被成形或转换为反射率剖面(P-P反射率)。可以根据声阻抗相对于空间的导数(或者更一般地,垂直导数,其在一些情况下可能是时间)近似地确定反射率剖面。但是,还有一个更多的步骤来平衡跨不同角度的反射率谱,即,“拉伸”。虽然图1示出了在同一步骤中的“成形”和“拉伸”,但这些不一定被同时执行。
类似于迁移中的小波拉伸效应,从不同反射角的数据获得的反射率具有不同的分辨率。由于只有带限数据(ω0~ωf),其中ω0和ωf是在数据中存在的最小和最大频率,因此在每个角度θ,反射率谱仅在如图2所示的波数域中从到进行采样。波数域中的不同带宽导致空间域中的不同振幅。假定真实的反射率值R(θ1)和R(θ2)是相同的,那么反演的反射率之间的关系将是因此,为了在模型域中保持数据AVO,需要通过除以cosθ的因子来补偿频谱拉伸。实际上,使用单个反射角的数据是困难的。因此,获得补偿因子的替代方法是测量反射率的带宽。对于所有反射率剖面,可以使用傅里叶变换来计算在相同位置处应用到所有剖面的局部窗口内的平均频谱。只要位置在所有剖面中都是相同的,这就可以在多个位置处执行并求平均。傅立叶变换可以是优选的,但也可以使用其它变换,诸如FFT或DFT。带宽可以被定义为例如10-dB点之间的距离。但是,可以使用带宽的其它测量(即,3-dB点、半极大值全宽、最陡斜率点等),但是10-dB点之间的距离可能是优选的。然后,为了完成步骤105并补偿“拉伸”,每个反射率剖面通过其自己的带宽被归一化,使得对反射率振幅的频谱拉伸效应被校正。
图1中的方法(步骤106)的最终输出是不同角度的反射率叠加。
本技术进步应用于用从SEG SEAM阶段I模型提取的2-D切片生成的合成数据集。压力数据是用4公里最大偏移量的拖缆(streamer)模拟的。在水面上使用了吸收边界条件。因此,数据中不存在与自由表面相关的倍数。如图3A所示,通过与以下角度范围0-10;10-20;20-30;30-40;以及40-50度对应的虚线301、302、303、304和305将单个炮集划分成五个剖面。角度范围不一定需要重叠。基于这些线301-305设计五个数据屏蔽,每个屏蔽覆盖10度的反射角。因此,如图3B所示,生成了五个数据子集(参见步骤102)。从相同的速度模型开始,分别使用每个子集来执行声学FWI(参见步骤103)。在反演之后结果得到的声学阻抗(参见步骤104)被“成形”和“拉伸”(参见步骤105)成反射率剖面(参见步骤106),并且在图4A的面板中示出。最后面板中的间隙区域401是因为来自深部的大角度反射在采集偏移的外部。为了验证结果的质量,通过地震子波与反射率序列之间的卷积生成真实的反射率剖面,并且在图4B中示出。反射率序列使用Zoeppritz方程[4]和真实的合成模型进行计算。图4A和图4B中的AVA看起来是一致的。为了更仔细地研究,来自模型域叠加的五条垂直线402a、403a、404a、405a和406a以及来自卷积剖面402b、403b、404b、405b和406b的五条垂直线在图5中被重叠以示出一致性。很清楚,在所有角度和所有深度下,与卷积叠加相比,FWI模型域叠加具有非常相似的振幅。它表明工作流程在地质情况复杂的情况下是可靠且准确的。
最终反射率是可用于地下的解释和/或油气勘探的管理的地下图像的示例。如本文所使用的,油气管理包括油气提取、油气生产、油气勘探、识别潜在的油气资源、识别井的位置、确定井注入和/或提取率、识别储层连接性、获取、处置和/放弃油气资源、回顾先前的油气管理决策、以及其它与油气有关的行为或活动。
在所有实际应用中,本技术进步必须与根据本文公开内容编程的计算机结合使用。优选地,为了高效地执行FWI,计算机是如本领域技术人员已知的高性能计算机(HPC)。这种高性能计算机通常涉及节点集群,每个节点具有多个CPU和允许并行计算的计算机存储器。可以使用任何交互式可视化程序和相关联硬件(诸如显示器和投影机)对模型进行可视化和编辑。***的体系架构可以变化,并且可以由能够执行逻辑操作并根据本技术进步显示输出的任何数量的合适的硬件结构组成。本领域普通技术人员知道可从Cray或IBM获得的合适的超级计算机。
参考文献:
以下参考文献通过引用被整体地结合于此:
[1]Xu,S.,Y.Zhang和B.Tang,2011年,3D angle gathers from reverse timemigration:Geophysics,76:2,S77–S92.doi:10.1190/1.3536527;
[2]Thomas A.Dickens和Graham A.Winbow(2011)RTM angle gathers usingPoynting vectors.SEG Technical Program Expanded Abstracts 2011:pp.3109-3113;
[3]Yu Zhang,Lian Duan和Yi Xie(2013)A stable and practicalimplementation of least-squares reverse time migration.SEG Technical ProgramExpanded Abstracts 2013:pp.3716-3720;
[4]Encyclopedic Dictionary of Applied Geophysics,R.E.Sheriff,4thedition.,SEG,2002,p.400;
[5]Encyclopedic Dictionary of Applied Geophysics,R.E.Sheriff,4thedition.,SEG,2002,p.12;以及
[6]Aki,K和Richards,P(2002)Quantitative seismology,2nd edition,University Science Books,p.148.
Claims (12)
1.一种方法,包括:
获得根据预定的地下反射角范围被分成子集的地震数据集;
利用计算机分别对每个子集执行声波全波场反演过程以对密度进行反演并生成相应的密度模型;
使用相应的密度模型,作为反射角的函数为每个子集生成声阻抗;以及
使用计算机将每个子集的声阻抗变换成反射率剖面,其中所述变换包括通过它们相应的带宽对反射率剖面进行归一化。
2.如权利要求1所述的方法,其中每个全波场反演过程从相同的速度模型开始。
3.如任何一个前述权利要求所述的方法,其中每个全波场反演过程被独立地应用到子集。
4.如任何一个前述权利要求所述的方法,其中所述获得包括通过使用包括反射体倾斜角和P-波速度的信息的数据屏蔽将炮集划分为子集。
5.如任何一个前述权利要求所述的方法,还包括对于每个反射率剖面使用傅里叶变换、离散傅立叶变换或快速傅立叶变换以计算在相同位置处应用于全部反射率剖面的至少一个局部窗口内的平均频谱,并确定每个平均频谱的带宽。
6.如权利要求5所述的方法,其中确定带宽是基于10-dB点之间的距离。
7.如权利要求5所述的方法,其中确定带宽是基于具有最陡斜率的点之间的距离。
8.如权利要求5所述的方法,其中所述平均频谱在多个局部窗口内被计算并且被平均。
9.如任何一个前述权利要求所述的方法,还包括确定多个角度的反射率值并通过内插构建角度与振幅的曲线。
10.如任何一个前述权利要求所述的方法,还包括使用反射率剖面管理油气生产。
11.如任何一个前述权利要求所述的方法,其中管理油气生产包括在至少部分地由反射率剖面确定的位置处钻井。
12.一种利用指令编码的非瞬态计算机可读存储介质,当所述指令由计算机执行时使所述计算机实现一种方法,所述方法包括:
获得根据预定的地下反射角范围被分成子集的地震数据集;
利用计算机分别对每个子集执行声波全波场反演过程以对密度进行反演并生成相应的密度模型;
使用相应的密度模型,作为反射角的函数为每个子集生成声阻抗;以及
使用计算机将每个子集的声阻抗变换成反射率剖面,其中所述变换包括通过它们相应的带宽对反射率剖面进行归一化。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562241780P | 2015-10-15 | 2015-10-15 | |
US62/241,780 | 2015-10-15 | ||
PCT/US2016/049390 WO2017065889A1 (en) | 2015-10-15 | 2016-08-30 | Fwi model domain angle stacks with amplitude preservation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108139498A true CN108139498A (zh) | 2018-06-08 |
CN108139498B CN108139498B (zh) | 2019-12-03 |
Family
ID=56889232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680059382.2A Expired - Fee Related CN108139498B (zh) | 2015-10-15 | 2016-08-30 | 具有振幅保持的fwi模型域角度叠加 |
Country Status (9)
Country | Link |
---|---|
US (1) | US10520619B2 (zh) |
EP (1) | EP3362823B1 (zh) |
KR (1) | KR102021276B1 (zh) |
CN (1) | CN108139498B (zh) |
AU (1) | AU2016337084B2 (zh) |
BR (1) | BR112018004435A2 (zh) |
CA (1) | CA2998519A1 (zh) |
MX (1) | MX2018003495A (zh) |
WO (1) | WO2017065889A1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018031113A1 (en) * | 2016-08-12 | 2018-02-15 | Exxonmobil Upstream Research Company | Tomographically enhanced full wavefield inversion |
US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
CN109655891B (zh) * | 2017-10-11 | 2020-12-01 | 中国石油化工股份有限公司 | 克服全波形反演周波跳跃的方法及*** |
DE102018108544A1 (de) * | 2018-01-31 | 2019-08-01 | Bauhaus-Universität Weimar | Detektion und Lokalisierung von Bauwerksschäden mit vollständiger Wellenformumkehrung |
EP3749638B1 (en) | 2018-02-09 | 2023-04-05 | Vistagen Therapeutics, Inc. | Synthesis of 4-chlorokynurenines and intermediates |
US11340368B2 (en) | 2019-03-21 | 2022-05-24 | Saudi Arabian Oil Company | Generating a velocity model and a density model of a subsurface structure of a reservoir |
US11914101B2 (en) | 2020-01-31 | 2024-02-27 | ExxonMobil Technology and Engineering Company | Method for partitioning a search direction when using least squares reverse time migration |
US11320557B2 (en) | 2020-03-30 | 2022-05-03 | Saudi Arabian Oil Company | Post-stack time domain image with broadened spectrum |
US11193364B1 (en) * | 2020-06-03 | 2021-12-07 | Schlumberger Technology Corporation | Performance index using frequency or frequency-time domain |
NO346427B1 (en) * | 2021-02-03 | 2022-08-08 | Pre Stack Solutions Geo As | Method of seismic inversion using artificial kinematic constraints |
US11971513B2 (en) | 2021-05-21 | 2024-04-30 | Saudi Arabian Oil Company | System and method for forming a seismic velocity model and imaging a subterranean region |
WO2024150065A1 (en) | 2023-01-13 | 2024-07-18 | Cgg Services Sas | System and method for extracting amplitude-versus-angle or offset information from seismic data |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
US20040220743A1 (en) * | 2003-04-30 | 2004-11-04 | Conocophillips Company | Method for determining shear-wave velocity model for depth migration of mode-converted data |
US20040243313A1 (en) * | 2003-04-29 | 2004-12-02 | Karine Broto | Method for determining a seismic wave velocity model in a heterogeneous underground |
WO2014189679A1 (en) * | 2013-05-24 | 2014-11-27 | Exxonmobil Upstream Research Company | Multi-parameter inversion through offset dependent elastic fwi |
Family Cites Families (216)
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 |
US4545039A (en) | 1982-09-09 | 1985-10-01 | Western Geophysical Co. Of America | Methods for seismic exploration |
US4675851A (en) | 1982-09-09 | 1987-06-23 | Western Geophysical Co. | Method 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 | 内燃エンジンの作動状態制御方法 |
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 |
GB2322704B (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 |
US5924049A (en) | 1995-04-18 | 1999-07-13 | Western Atlas International, Inc. | Methods for acquiring and processing seismic data |
CA2188255C (en) | 1995-04-18 | 2003-03-25 | Craig J. Beasley | Method for providing uniform subsurface coverage in the presence of steep dips |
US5721710A (en) | 1995-09-29 | 1998-02-24 | Atlantic Richfield Company | High fidelity vibratory source seismic method with source separation |
US5719821A (en) | 1995-09-29 | 1998-02-17 | Atlantic Richfield Company | Method and apparatus for source separation of seismic vibratory signals |
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 |
US5715213A (en) | 1995-11-13 | 1998-02-03 | Mobil Oil Corporation | High fidelity vibratory source seismic method using a plurality of vibrator sources |
US5822269A (en) | 1995-11-13 | 1998-10-13 | Mobil Oil Corporation | Method for separation of a plurality of vibratory seismic energy source signals |
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 |
US6574564B2 (en) | 1998-10-01 | 2003-06-03 | Institut Francais Du Petrole | 3D prestack seismic data migration 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 |
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 |
EP1151326B1 (en) | 1999-02-12 | 2005-11-02 | Schlumberger Limited | 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 |
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 |
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 |
AU774883B2 (en) | 2000-01-21 | 2004-07-08 | Schlumberger Holdings Limited | 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 |
US6546339B2 (en) | 2000-08-07 | 2003-04-08 | 3D Geo Development, Inc. | Velocity analysis using angle-domain common image gathers |
AU2002211768B2 (en) | 2000-10-17 | 2005-08-04 | Westerngeco Seismic Holdings Limited | 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 |
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 |
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 |
US7072767B2 (en) | 2003-04-01 | 2006-07-04 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
WO2004095073A2 (en) | 2003-04-01 | 2004-11-04 | 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 |
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 |
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 |
EP1617309B1 (en) | 2004-07-15 | 2011-01-12 | Fujitsu Limited | Simulation technique with local grid refinement |
US7525873B1 (en) * | 2004-07-23 | 2009-04-28 | Bf Corporation North America Inc. | Seismic inversion of conditioned amplitude spectra |
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 |
EP1859301B1 (en) | 2005-02-22 | 2013-07-17 | 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 |
WO2006122146A2 (en) | 2005-05-10 | 2006-11-16 | William Marsh 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 |
US20090164756A1 (en) | 2005-10-18 | 2009-06-25 | Tor Dokken | Geological Response Data Imaging With Stream Processors |
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 |
BRPI0716853A2 (pt) | 2006-09-28 | 2013-10-01 | Exxonmobil Upstream Res Co | mÉtodos para determinar um modelo de propriedades fÍsicas para uma regiço de subsuperfÍcie, e para produzir hidrocarbonetes a partir de uma regiço de subsuperfÍcie |
EP2089738B1 (en) | 2006-12-07 | 2019-01-02 | Council of Scientific & Industrial Research | A method for computing an exact impulse response of a plane acoustic reflector at zero offset due to a point acoustic source |
WO2008087505A2 (en) | 2007-01-20 | 2008-07-24 | Spectraseis Ag | Time reverse reservoir localization |
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 |
WO2009088602A1 (en) | 2008-01-08 | 2009-07-16 | 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 |
ES2651923T3 (es) | 2008-03-21 | 2018-01-30 | Exxonmobil Upstream Research Company | Un método eficiente para la inversión de datos geofísicos |
EP2265975A4 (en) | 2008-03-28 | 2017-05-24 | Exxonmobil Upstream Research Company | Surface wave mitigation in spatially inhomogeneous media |
EP2105765A1 (en) | 2008-03-28 | 2009-09-30 | Schlumberger Holdings Limited | Simultaneous inversion of induction data for dielectric permittivity and electric conductivity |
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 |
CN102124374B (zh) | 2008-08-15 | 2013-07-17 | Bp北美公司 | 用于分离单独的同时震源的方法 |
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 |
WO2010085822A2 (en) | 2009-01-26 | 2010-07-29 | 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 |
WO2010095859A2 (ko) | 2009-02-17 | 2010-08-26 | Shin Changsoo | 지하구조 영상화 장치 및 방법 |
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 |
US7856528B1 (en) | 2009-08-11 | 2010-12-21 | Texas Memory Systems, Inc. | Method and apparatus for protecting data using variable size page stripes in a FLASH-based storage system |
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 |
CA2767757A1 (en) | 2009-09-09 | 2011-03-17 | 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 |
US8861308B2 (en) | 2009-12-07 | 2014-10-14 | Westerngeco L.L.C. | 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 |
WO2011091216A2 (en) | 2010-01-22 | 2011-07-28 | Schlumberger Canada Limited | Real-time formation anisotropy and dip evaluation using tri-axial induction measurements |
EP2529254A4 (en) | 2010-01-25 | 2017-08-02 | CGG Veritas Services (U.S.) Inc. | Methods and systems for estimating stress using seismic data |
US8265875B2 (en) | 2010-01-29 | 2012-09-11 | Westerngeco L.L.C. | Interpolation of periodic data |
AU2010344186B2 (en) | 2010-01-29 | 2016-04-14 | Exxonmobil Upstream Research Company | Temporary field storage of gas to optimize field development |
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 |
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 |
JP5812990B2 (ja) | 2010-06-15 | 2015-11-17 | 電気化学工業株式会社 | 透光性硬質基板積層体の製造方法 |
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 |
EP2606452A4 (en) | 2010-08-16 | 2017-08-16 | Exxonmobil Upstream Research Company | Reducing the dimensionality of the joint inversion problem |
US8243548B2 (en) * | 2010-08-27 | 2012-08-14 | Board Of Regents Of The University Of Texas System | Extracting SV shear data from P-wave seismic data |
US20120051176A1 (en) | 2010-08-31 | 2012-03-01 | Chevron U.S.A. Inc. | Reverse time migration back-scattering noise removal using decomposed wavefield directivity |
BR112013006423A2 (pt) | 2010-09-20 | 2016-07-26 | Chevron Usa Inc | sistema e método par agerar imagens de estruturas de sub-superfície |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
BR112013002842A2 (pt) | 2010-09-27 | 2016-06-07 | Exxonmobil Upstream Res Co | codificação de fonte e separação de fonte simultâneas como uma solução prática para inversão de campo de onda completa |
AU2011310635B2 (en) | 2010-09-28 | 2014-09-18 | Shell Internationale Research Maatschappij B.V. | Earth model estimation through an acoustic Full Waveform Inversion of seismic data |
CN103238158B (zh) | 2010-12-01 | 2016-08-17 | 埃克森美孚上游研究公司 | 利用互相关目标函数进行的海洋拖缆数据同时源反演 |
US9134442B2 (en) | 2010-12-16 | 2015-09-15 | Bp Corporation North America Inc. | Seismic acquisition using narrowband seismic sources |
US9702994B2 (en) | 2011-02-18 | 2017-07-11 | Westerngeco L.L.C. | Waveform inversion by multiple shot-encoding for non-fixed spread geometries |
AU2012233133B2 (en) | 2011-03-30 | 2014-11-20 | Exxonmobil Upstream Research Company | Convergence rate of full wavefield inversion using spectral shaping |
US20120275267A1 (en) | 2011-04-26 | 2012-11-01 | Ramesh Neelamani | Seismic Data Processing |
EP2707757B1 (en) | 2011-05-13 | 2016-08-10 | Saudi Arabian Oil Company | Couple time-distance dependent swept frequency source acquisition design and data de-noising |
US9075159B2 (en) | 2011-06-08 | 2015-07-07 | Chevron U.S.A., Inc. | System and method for seismic data inversion |
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 |
EP2751710B1 (en) | 2011-09-02 | 2017-08-02 | 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 |
SG11201404094RA (en) | 2012-03-08 | 2014-10-30 | Exxonmobil Upstream Res Co | Orthogonal source and receiver encoding |
US9541661B2 (en) | 2012-04-19 | 2017-01-10 | Cgg Services Sa | Device and method for deghosting variable depth streamer data |
US9435905B2 (en) | 2012-04-19 | 2016-09-06 | Cgg Services Sa | Premigration deghosting of seismic data with a bootstrap technique |
US20130311149A1 (en) | 2012-05-17 | 2013-11-21 | Yaxun Tang | Tomographically Enhanced Full Wavefield Inversion |
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 |
CA2943589C (en) * | 2014-03-28 | 2022-11-15 | Cgg Services Sa | Method and apparatus for characterizing azimuthal anisotropy using statistical moments |
US10591654B2 (en) * | 2015-05-21 | 2020-03-17 | Board Of Regents, The University Of Texas System | Ultra-sparse dielectric resonant polarizers and associated methods |
-
2016
- 2016-08-30 KR KR1020187013662A patent/KR102021276B1/ko active IP Right Grant
- 2016-08-30 CA CA2998519A patent/CA2998519A1/en not_active Abandoned
- 2016-08-30 EP EP16763157.1A patent/EP3362823B1/en active Active
- 2016-08-30 WO PCT/US2016/049390 patent/WO2017065889A1/en active Application Filing
- 2016-08-30 CN CN201680059382.2A patent/CN108139498B/zh not_active Expired - Fee Related
- 2016-08-30 AU AU2016337084A patent/AU2016337084B2/en not_active Ceased
- 2016-08-30 US US15/251,298 patent/US10520619B2/en active Active
- 2016-08-30 MX MX2018003495A patent/MX2018003495A/es unknown
- 2016-08-30 BR BR112018004435A patent/BR112018004435A2/pt not_active Application Discontinuation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
US20040243313A1 (en) * | 2003-04-29 | 2004-12-02 | Karine Broto | Method for determining a seismic wave velocity model in a heterogeneous underground |
US20040220743A1 (en) * | 2003-04-30 | 2004-11-04 | Conocophillips Company | Method for determining shear-wave velocity model for depth migration of mode-converted data |
WO2014189679A1 (en) * | 2013-05-24 | 2014-11-27 | Exxonmobil Upstream Research Company | Multi-parameter inversion through offset dependent elastic fwi |
Also Published As
Publication number | Publication date |
---|---|
US10520619B2 (en) | 2019-12-31 |
KR102021276B1 (ko) | 2019-09-16 |
AU2016337084A1 (en) | 2018-05-10 |
CN108139498B (zh) | 2019-12-03 |
EP3362823A1 (en) | 2018-08-22 |
US20170108602A1 (en) | 2017-04-20 |
MX2018003495A (es) | 2018-06-06 |
BR112018004435A2 (pt) | 2018-09-25 |
KR20180067650A (ko) | 2018-06-20 |
CA2998519A1 (en) | 2017-04-20 |
WO2017065889A1 (en) | 2017-04-20 |
EP3362823B1 (en) | 2019-10-09 |
AU2016337084B2 (en) | 2019-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108139498B (zh) | 具有振幅保持的fwi模型域角度叠加 | |
Martin et al. | Marmousi-2: An updated model for the investigation of AVO in structurally complex areas | |
Warner et al. | Anisotropic 3D full-waveform inversion | |
Stopin et al. | Multiparameter waveform inversion of a large wide-azimuth low-frequency land data set in Oman | |
US10317548B2 (en) | Reflection seismic data Q tomography | |
US20100131205A1 (en) | Method for identifying and analyzing faults/fractures using reflected and diffracted waves | |
CN106556861B (zh) | 一种基于全方位地震资料的方位avo反演方法 | |
CN107390270B (zh) | 一种基于弹性波逆时偏移ADCIGs的AVA分析方法 | |
Colombo et al. | Near-surface full-waveform inversion in a transmission surface-consistent scheme | |
Alam | Near-surface characterization using traveltime and full-waveform inversion with vertical and horizontal component seismic data | |
Kamei et al. | Application of waveform tomography to a crooked-line 2D land seismic data set | |
Duret et al. | Near-surface velocity modeling using a combined inversion of surface and refracted P-waves | |
Vigh et al. | Elastic full-waveform inversion using OBN data acquisition | |
Chen et al. | 3-D seismic diffraction separation and imaging using the local rank-reduction method | |
EA030770B1 (ru) | Система и способ адаптивной сейсмической оптики | |
CN109143398B (zh) | 一种自动网格层析深度域速度的建模方法 | |
Zhou et al. | Orthorhombic velocity model building and imaging of Luda field with WAZ OBC data | |
Gao et al. | Waveform tomography of two-dimensional three-component seismic data for HTI anisotropic media | |
Paul et al. | Interpretation Guided Seismic Data Processing of a Large 3D OBC Survey With Simultaneous Shooting: A Case Study from Offshore Abu Dhabi | |
Buia et al. | Depth imaging Coil data: Multi azimuthal tomography earth model building and depth imaging the full azimuth Tulip coil project | |
Chandola et al. | Seismic acquisition and imaging strategies for unlocking subsurface complexities in Malaysian basins | |
Huang et al. | Imaging Complex Subsurface Structures for Geothermal Exploration at Pirouette Mountain and Eleven-Mile Canyon in Nevada | |
Perrone et al. | Waveform tomography based on local image correlations | |
Glushchenko et al. | A broadband full azimuth land seismic case study Onshore Abu Dhabi | |
Ortega et al. | Processing of large offset data: experimental seismic line from Tenerife Field, Colombia |
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: 20191203 Termination date: 20200830 |