CN106461802A - 用于多参数全波场反演的有效的线性搜索方法 - Google Patents
用于多参数全波场反演的有效的线性搜索方法 Download PDFInfo
- Publication number
- CN106461802A CN106461802A CN201580022137.XA CN201580022137A CN106461802A CN 106461802 A CN106461802 A CN 106461802A CN 201580022137 A CN201580022137 A CN 201580022137A CN 106461802 A CN106461802 A CN 106461802A
- Authority
- CN
- China
- Prior art keywords
- search
- model
- linear search
- physical characteristics
- linear
- 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
- 238000000034 method Methods 0.000 title claims abstract description 88
- 230000005284 excitation Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000000704 physical effect Effects 0.000 abstract description 2
- 239000013598 vector Substances 0.000 description 10
- 238000012986 modification Methods 0.000 description 8
- 230000004048 modification Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 7
- 238000013508 migration Methods 0.000 description 6
- 230000005012 migration Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000007430 reference method Methods 0.000 description 1
- 230000017105 transposition Effects 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/32—Transforming one recording into another or 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/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- 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/364—Seismic filtering
- G01V1/368—Inverse filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/42—Waveform, i.e. using raw or pre-filtered trace data
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Analysis (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
用于同时对多个类别的物理特性参数(例如,速度和各向异性)的全波场地震数据进行反演的方法,该方法通过进行以下操作:针对每个类别的参数计算目标函数的梯度(即,搜索方向);然后将(优选地,穷尽的)第一遍独立线性搜索应用于每个参数类别,以沿着每个参数类别的搜索方向获得对应的步长;然后尚未更新模型,使用步长定义在所有参数类别的梯度之间的相对缩放。接下来,重新结合每个经缩放的搜索方向,以形成新的搜索方向,并且沿着新的搜索方向实行新的第二遍线性搜索,并且用所得的步长同时更新所有的参数。作为前述供选择的两遍实施例的供选择的方案,可以在每个第一遍线性搜索之后更新模型,并且不实行第二遍线性搜索。
Description
相关申请的交叉引用
本申请要求2014年5月9日提交的标题为“用于多参数全波场反演的有效的线性搜索方法(EFFICIENT LINE SEARCH METHODS FOR MULTI-PARAMETER FULL WAVEFIELDINVERSION)”的美国临时专利申请61/990,860的权益,该专利申请的全部内容以引用方式并入本文。
技术领域
一般来说,本公开涉及用于碳氢化合物的地球物理勘探的领域,并且更具体地,涉及地震数据处理。具体地,本公开涉及用于在地震数据的多参数全波场反演(“FWI”)中实行有效的线性搜索以推断地下物理特性模型的方法。在碳氢化合物的探测或碳氢化合物的生产中这样的模型可以是有用的。
背景技术
全波场反演是非线性反演技术,该非线性反演技术通过使所模拟的地震波场和所观测的地震波场之间的失配最小化,恢复地球模型。由于与FWI相关联的高计算成本,常规的实施方式利用局部优化技术来估计最优模型参数。广泛使用的局部优化技术是基于梯度的一阶方法(例如,最速下降或非线性共轭梯度),该方法仅利用目标函数的梯度信息定义搜索方向。虽然唯梯度一阶方法是相对有效的(该方法要求仅计算目标函数的梯度),但是其收敛性一般是缓慢的。通过使用二阶方法,能够显著提高FWI的收敛性。因为二阶方法利用目标函数的梯度和曲率信息两者来确定模型参数空间中的最优搜索方向,所以实现该提高的收敛性。(搜索方向单位矢量s通过m更新的=m+αs与模型更新过程有关,其中α(标量)是步长。)
一阶方法和二阶方法之间的主要差异是二阶方法用逆矩阵海赛函数(inverseHessian)(例如,高斯牛顿/牛顿方法)或用所投影的海赛函数的逆矩阵(例如,子空间方法)预处理梯度。海赛函数是目标函数关于模型参数的二阶偏导数的矩阵。一般来说,二阶方法具有吸引力,不仅是因为它们相对快速的收敛速率,还因为在多参数反演的情况下平衡不同的参数类别的梯度和为具有不同的数据灵敏度的参数类别(例如,速度、各向异性、衰减等)提供有意义的更新的能力。在二阶方法中,如果同时对此参数类别进行反演,则使用海赛函数的参数类别的最优缩放在多参数反演中是关键的。然而,因为计算海赛函数的逆矩阵非常昂贵,所以这是在实践中广泛采用二阶方法是主要的障碍。二阶方法的另一个缺点是如果目标函数不是二次的或凸性的(convex)(例如,其中初始模型远不是真实模型),则海赛函数或其近似不能准确地预测目标函数的形状。从此,不能适当地缩放不同参数类别的梯度,从而得到次优搜索方向。
发明内容
在一个实施例中,本发明是用于对地震数据进行迭代反演以同时推断地下的至少两个物理特性的模型的方法,所述方法包括:
(a)对于每个物理特性,针对物理特性的参数,计算目标函数的梯度,所述目标函数测量地震数据中的全部或一部分和对应模型所模拟的地震数据之间的不匹配;
(b)对于每个物理特性,根据梯度计算模型空间中的搜索方向;
(c)针对至少两个物理特性在搜索方向之间或之中进行交替线性搜索,以确定沿着搜索方向中的每个搜索方向的最优步长;以及
(d)使用最优步长更新模型。
附图说明
通过参考以下详细描述和附图可以更好地理解本发明及其优点,在附图中:
图1例示在本发明的交替的一遍线性搜索实施例中的最优搜索方向;
图2例示可以由本发明的交替的一遍线性搜索实施例得到的次优搜索方向;
图3例示使用本发明的交替的两遍线性搜索实施例的最优搜索方向;
图4A至图4F示出,对于本发明的方法的测试示例,用于两个参数(速度和各向异性)的初始模型和最终反演模型,并且使用初始模型和最终模型比较在迁移之后的数据;
图5是示出在本发明的方法的交替的一遍线性搜索中的基本步骤的流程图;
图6是示出在本发明的方法的交替的两遍线性搜索中的基本步骤的流程图;以及
图7是示出在本文中所公开的扩展的交替的两遍线性搜索中的基本步骤的流程图,其中级联反演方法与交替的两遍线性搜索结合。
由于对附图中颜色使用的专利规则限制,所以图4A至图4F中的一些是原始彩色生成的附图的黑白再现。
将结合示例实施例描述本发明。然而,在一定程度上,以下详细描述针对本发明的特定实施例或特定用途,这旨在仅是说明性的,并且不被解释为限制本发明的范围。相反,旨在覆盖可以包含在如由随附权利要求所限定的本发明的范围内的所有的供选择的方案、修改和等价物。
具体实施方式
通过使用可以被称为交替的一遍/两遍线性搜索方法(该方法不要求来自海赛函数矩阵的明示信息,而是通过连续的线性搜索近似二阶信息),本发明弥合一阶和二阶最优方法之间的差距。将示出本发明的方法能够适当地缩放具有不同的数据灵敏度的参数类别的梯度,并且能够同时为多个参数类别提供有意义的更新。在实践中,因为本发明的方法没有假设目标函数是二次的,并且如果能够有效地实施每个线性搜索,则本发明的方法还能够显著地更廉价,所以与基于海赛函数的二阶方法相比,本发明的方法能够更加强健。
虽然为了简单起见,使用两个参数类别描述了理论,但是本发明适用于任何数量的参数类别的同时反演,并且该方法扩展到超过两个参数类别是简单明了的。对于其中对两个参数类别进行反演的情况,模型能够表示为包含两个不同子模型的矢量,即,m=(m1m2)T,其中m1和m2分别是第一模型参数类别和第二模型参数类别,其中T代表转置。当前迭代s的搜索方向是这两个参数类别的搜索方向的级联,并且能够被写为:
其中,s1和s2分别是第一参数类别和第二参数类别的搜索方向。基于一阶的方法通常沿着方向s更新模型m。该方法的一个主要问题是如果模型参数类别是具有非常不同的单位和对FWI目标函数的非常不同的灵敏度的物理量(例如,速度、各向异性、衰减等),则从FWI梯度导出的得到的搜索方向将通常具有非常不同的量值,即,s1的量值可以与s2的量值显著不同。这常常导致FWI选取收敛路径,FWI沿着该收敛路径主要更新对目标函数更加灵敏的参数类别,同时使对目标函数不太灵敏的参数类别保持基本不更新。因此,FWI可以朝次优解收敛。下面是所描述的本发明的基于交替的线性搜索的方法的两个供选择的实施例,所述实施例解决了该问题,其目标是为所有的参数类别同时提供最优更新,同时不会由于不同的单位或数据灵敏度而偏向某些参数类别。
方法I:交替的一遍线性搜索
为了开始,定义两个基矢量,两个基矢量中的每个基矢量包含具体模型参数类别的搜索方向:
其中,0指代包含零的矢量。然后,原始搜索方向能够被写为上面的两个基矢量的和,上面的两个基矢量彼此正交(即,):
在该第一方法中,用图5中所示的基本步骤以交替的方式更新两个模型参数类别。当前模型52被用于模拟预测的数据,并且其与所测量的数据51结合,以计算目标函数,并且然后在步骤53计算关于第一模型参数的目标函数的梯度,以在步骤54生成搜索方向在步骤55,沿着方向实行线性搜索,并且更新模型。由于仅对于第一参数类别来说,是非零的,所以将仅更新参数类别m1。在已经更新模型之后,然后,沿着方向实行线性搜索,并且再次更新模型(步骤56)。在该情况下,由于仅对于第二参数类别来说,是非零的,所以将仅更新参数m2,如图1所示。在该方法中,一个FWI迭代是指以该交替的方式更新两个参数类别的过程,并且在步骤57,针对另一个迭代重复该过程,除非根据一些预定收敛准则另一个迭代已经收敛或已经达到另一个停止点。因为以独立的方式更新m1和m2,所以克服了在上面所描述的搜索方向(或梯度)中的缩放失衡。然后,所有的参数类别以交替的方式接收显著的更新。该方法可以被称为交替的一遍线性搜索。该方法的一个潜在的缺陷是连续的线性搜索的次序能够影响收敛路径,并且如果目标函数具有相对复杂的空间,则这可以导致之字形收敛路径(图2)。
更详细地,线性搜索涉及在高度多维模型参数空间中沿着搜索方向针对各种不同的“步”长的模型模拟的地震数据。选择使模型模拟的数据和所测量的数据之间的不匹配最小化的步长。在上面所描述的方法I中,在模型中更新第一物理特性之后,然后所更新的模型被用于模拟用于线性搜索的数据,以更新第二物理特性。尽管第二线性搜索仅更新第二物理特性,但是根据第一线性搜索更新第一物理特性将影响用于第二线性搜索的模型模拟的数据,并且因而将影响对第二物理特性的所得的更新。
为了进一步解释图1,由在目标函数轮廓的中心中的星号表示解,因为与模型参数m2相比,模型参数m1对目标函数更加灵敏,所以目标函数轮廓是椭圆而不是圆。图1示出在这种情形下通过使用本发明的方法的交替的一遍线性搜索实施例进行两个独立步骤能够在单个迭代过程周期中获得该解。图1还示出传统的线性搜索s由不同的灵敏度偏离右边方向,并且将要求更多的迭代以获得该解。图2示出其中交替的一遍线性搜索可能是有问题的,即,当m1和m2耦合时,如由是倾斜的目标函数轮廓所示的,意味着这两个参数类别之间的强相关性。在此类状况下,交替的一遍线性搜索可能要求额外的迭代以找到最优搜索方向。图3示出本发明的交替的两遍线性搜索实施例(接下来所解释的)如何能够解决该问题。
方法II:交替的两遍线性搜索
通过对方法I的修改,能够减轻上面所识别的方法I的缺陷中的一些。如在图6中所指示的,所观测的数据61和当前地下模型62被用于针对每个参数类别计算目标函数的梯度(步骤63),得到每个参数类别的搜索方向(64)。在步骤65,针对每个参数类别,使用在步骤64针对此类参数类别所计算的搜索方向,实行第一遍独立线性搜索。保存所计算的步长,但尚不实行模型更新。在步骤66,该修改的方法在其方面中的一个方面上类似于所谓的子空间方法(Kennett等人,1998),但不是以交替的方式更新所有模型类别,该修改的方法首先使用两个缩放因子α和β,即,步长,重新结合两个基矢量和以形成如下的新的搜索方向:
一旦获得该新的结合的搜索方向,就实行第二遍线性搜索(利用新的搜索方向),以同时更新两个参数类别(步骤67)。该方法和已知的子空间方法(Kennett等人,1998)之间的主要差异是如何估计缩放因子α和β。在子空间方法中,通过对以下所投影的海赛函数矩阵求逆矩阵,确定缩放因子:
其中,H是海赛函数矩阵。代替使用海赛函数,方法II的优选的实施例通过线性搜索估计最优缩放因子(步骤65)。在本公开的方法II的该两个参数示例中,如下首先实行两个独立的线性搜索:
(i)首先,确定最优步长α,使得该最优步长α沿着由定义的搜索方向,使目标函数最小化,但实际上不更新模型;
(ii)然后,确定最优步长β,使得最优步长β沿着由定义的搜索方向,使目标函数最小化,实际上也不更新模型。
注意,该方法不用所估计的步长更新模型参数的事实表示该两遍线性搜索方法和先前所描述的一遍线性搜索方法之间的主要差异中的一个。一旦确定了缩放因子,则使用等式4形成新的搜索方向。然后,使用新的搜索方向s新的实行第二遍线性搜索,以更新两个参数类别,即,我们尝试找到步距λ,使得沿着方向s新的,使目标函数最小化。该方法可以被称为交替的两遍线性搜索方法:第一遍线性搜索确定不同的搜索分量之中的相对缩放,并且所估计的缩放有效地使原始搜索方向旋转,以获得新的搜索方向;然后,第二遍线性搜索使用新的搜索方向更新模型参数。由于通过独立的线性搜索确定缩放因子,所以然后能够使用所估计的缩放因子,校正原始搜索方向(或梯度)的缩放。主要通过沿着基矢量每个参数类别有多灵敏(即,沿着s1,m1的灵敏度和沿着s2,m2的灵敏度),确定所结合的搜索方向的每个分量的相对量值。
实践考虑和扩展
通过穷尽的线性搜索确定缩放
在本文中所公开的交替的两遍线性搜索方法中,一个目标是第一遍线性搜索为每个参数类别找到适当的缩放因子(步长)。为了适当地这样做,当实行第一遍线性搜索时,可以实行穷尽的线性搜索,以找到拐点。采取两个参数类别的反演,例如,意味着通过在扰动强度的范围内沿着搜索方向使级联模型m扰动,迭代地扫描目标函数值(实质上,这意味着沿着s1扫描m1,因为的第二分量是零)。当目标函数值大于先前的目标函数值时,终止该线性搜索。返回的α被用作缩放因子,以为第一参数类别的搜索方向进行加权。接下来描述了例示穷尽的线性搜索的算法:
do while
αi←αi+Δα
计算
如果则停止且返回α=αi
i←i+1
end
其中,Δα是用户提供的增量值,并且J是目标函数值。类似地,能够通过沿着搜索方向使m扰动,找到第二参数类别的最优缩放因子β。
扩展的方法-合并数据集的级联分解
还能够通过并入在发明人Ayeni等人的母案专利申请“一种用于通过波场分量的级联反演估计多个地下参数的方法(A Method for Estimating Multiple SubsurfaceParameters by Cascaded Inversion of Wavefield Components)”中所描述的级联反演方法来扩展本发明的方法,该母案专利申请以允许其的所有司法权以引用方式并入本文。该扩展允许使用对各个参数类别最灵敏的数据的部分,容易地调节所估计的梯度(并且因此调节搜索方向)。在图7的流程图中概括了本发明的方法的扩展的实施例。本文中所公开的该扩展的方法与Ayeni等人的文章中所描述的不同,至少因为:在本发明中,在每次迭代中同时更新所有的参数类别,而在由Ayeni等人教导的反演中,在每次迭代时仅更新一个参数类别。
多个线性搜索的成本考虑
与常规的基于梯度的一阶方法相比较,交替的一/两遍方法要求实行更多的线性搜索,这可能是高成本的。为了减轻该问题,可以仅使用震源激发的(一个或更多个)小的子集,可以优选地从整个调查随机选择震源激发的(一个或更多个)小的子集,以实行线性搜索。这些随机选择的激发能够被用于第一遍和/或第二遍线性搜索。
示例
在测试示例中,本文中所公开的交替的两遍线性搜索方法被应用于3D场数据示例的各向异性VTI(垂直横向各向同性)反演。在该示例中,同时对两个参数类别进行反演:Thompson的各向异性参数η和正常时差(Normal Moveout)速度VNMO。在反演期间,其它Thompson各向异性参数δ被设置为等于零,并在整个反演中,其它Thompson各向异性参数δ是固定的。图4A-图4F将模型与从初始模型和反演模型导出的迁移集合进行比较。图4A和图4D示出从初始(4A)模型和最终(4D)模型导出的基尔霍夫深度迁移集合。图4B和图4E分别示出初始模型和最终η模型,而图4C和图4F示出初始模型和反演VNMO模型。注意,相对于使用指示所导出的模型的准确性的开始模型迁移的图4A中所示的集合,使用从交替的两遍线性搜索方法(图4D中所示的迁移的集合)导出的VNMO和η的FWI反演模型迁移的所迁移的集合是扁平的。
上述描述针对本发明的特定的实施例,是出于例示其的目的。然而,对于本领域中的技术人员将明显的是,对本文中所描述的实施例进行很多修改和变化是可能的。所有此类修改和变化旨在在如由随附权利要求书限定的本发明的范围内。如对于在该技术领域中工作的人员来说将明显的是,使用根据本文中的本公开编程的计算机实行本发明的方法的所有实践应用。
参考文献
Kennett、B.L.N.、M.S.Sambridge和P.R.Williamson,“用于多个参数类别的大型逆矩阵问题的子空间方法(Subspace methods for large inverse problems withmultiple parameter classes)”,地理物理学杂志94,237-247(1988)。
Claims (11)
1.一种用于对地震数据进行迭代反演以同时推断地下的至少两个物理特性的模型的方法,所述方法包括:
(a)对于每个物理特性,针对所述物理特性的参数计算目标函数的梯度,所述目标函数测量所述地震数据中的全部或部分和对应模型模拟的地震数据之间的不匹配;
(b)对于每个物理特性,根据所述梯度计算模型空间中的搜索方向;
(c)针对所述至少两个物理特性,在所述搜索方向之间或之中进行交替线性搜索,以确定沿着所述搜索方向中的每个搜索方向的最优步长;以及
(d)使用所述最优步长更新所述模型。
2.根据权利要求1所述的方法,其中(c)-(d)包括:
(i)针对第一物理特性实行第一线性搜索,并且使用来自所述第一线性搜索的所述最优步长更新所述模型,生成第一更新的模型;
(ii)然后,使用所述第一更新的模型针对第二物理特性实行第二线性搜索,并且然后使用来自所述第二线性搜索的所述最优步长再次更新所述模型,生成第二更新的模型;以及
(iii)针对正被建模的任何额外的物理特性,如在(ii)中的进一步更新所述模型,为反演的当前迭代生成最终更新的模型。
3.根据权利要求2所述的方法,其中所述第一线性搜索和所述第二线性搜索以及任何额外的线性搜索是穷尽的线性搜索,其被实行以根据所述步长在所述目标函数中找到拐点。
4.根据权利要求1所述的方法,其中(c)-(d)包括:
(i)使用初始模型针对第一物理特性实行第一线性搜索,并且保存来自所述第一线性搜索的所述最优步长,但尚不更新所述模型;
(ii)使用所述初始模型针对第二物理特性实行第二线性搜索,并且保存来自所述第二线性搜索的所述最优步长,但尚不更新所述模型;
(iii)针对正被建模的任何额外的物理特性,如在(ii)中的使用所述初始模型实行线性搜索并保存所述最优步长;以及
(iv)使用来自(i)-(iii)的所保存的步长计算新的搜索方向,然后沿着所述新的搜索方向实行第二遍线性搜索,确定第二遍最优步长,并且使用所述第二遍最优步长为所述反演的当前迭代生成更新的模型。
5.根据权利要求4所述的方法,其中在(a)中,所述地震数据中的一个或更多个部分被选择和被用于针对每个物理特性计算所述梯度,基于所述物理特性对所述数据的灵敏度,根据所选择的准则确定所述一个或更多个部分。
6.根据权利要求5所述的方法,进一步包括,在所述迭代反演的下一个迭代中,针对所述至少两个物理特性中的至少一个物理特性,扩展所述地震数据的所选择的一个或更多个部分,以包含所述地震数据的更多部分。
7.根据权利要求4所述的方法,其中在(iii)中的所述第一线性搜索和所述第二线性搜索以及任何额外的线性搜索是穷尽的线性搜索,其被实行以根据所述步长在所述目标函数中找到拐点。
8.根据权利要求1所述的方法,其中所述地震数据中的一些震源激发被丢弃并且不在所述方法中使用,以便减少计算成本和时间。
9.根据权利要求1所述的方法,其中数据反演是全波场反演。
10.根据权利要求1所述的方法,其中所述至少两个物理特性包括P波速、S波速、正常时差速度、一个或更多个各向异性参数和衰减中的至少两个。
11.根据权利要求1所述的方法,其中计算所述搜索方向而不用所述目标函数的海赛函数预处理所述梯度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201461990860P | 2014-05-09 | 2014-05-09 | |
US61/990,860 | 2014-05-09 | ||
PCT/US2015/022786 WO2015171215A1 (en) | 2014-05-09 | 2015-03-26 | Efficient line search methods for multi-parameter full wavefield inversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106461802A true CN106461802A (zh) | 2017-02-22 |
CN106461802B CN106461802B (zh) | 2019-05-31 |
Family
ID=52824600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580022137.XA Expired - Fee Related CN106461802B (zh) | 2014-05-09 | 2015-03-26 | 用于多参数全波场反演的有效的线性搜索方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US9977142B2 (zh) |
EP (1) | EP3140675A1 (zh) |
KR (1) | KR101915451B1 (zh) |
CN (1) | CN106461802B (zh) |
AU (1) | AU2015256626B2 (zh) |
CA (1) | CA2947847C (zh) |
MX (1) | MX2016013366A (zh) |
RU (1) | RU2649214C1 (zh) |
SG (1) | SG11201608175SA (zh) |
WO (1) | WO2015171215A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150323689A1 (en) * | 2014-05-09 | 2015-11-12 | Yaxun Tang | Efficient Line Search Methods for Multi-Parameter Full Wavefield Inversion |
CN107886571A (zh) * | 2017-11-03 | 2018-04-06 | 中原工学院 | 一种利用计算机多维空间的数学建模方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN110361742B (zh) * | 2019-06-21 | 2021-03-26 | 中国人民解放军国防科技大学 | 一种基于星载三频毫米波雷达的云雨微物理参数反演方法 |
CN113376257B (zh) * | 2021-06-07 | 2023-03-17 | 北京理工大学 | 一种基于共轭梯度优化方法反演的多参数测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102906599A (zh) * | 2010-02-10 | 2013-01-30 | 埃克森美孚上游研究公司 | 用于全波场反演和逆时偏移中的地下参数估计的方法 |
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
CN103270430A (zh) * | 2011-06-08 | 2013-08-28 | 雪佛龙美国公司 | 用于地震数据反演的***和方法 |
CN103713315A (zh) * | 2012-09-28 | 2014-04-09 | 中国石油化工股份有限公司 | 一种地震各向异性参数全波形反演方法及装置 |
Family Cites Families (215)
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 |
US5719821A (en) | 1995-09-29 | 1998-02-17 | Atlantic Richfield Company | Method and apparatus for source separation of seismic vibratory signals |
US5721710A (en) | 1995-09-29 | 1998-02-24 | Atlantic Richfield Company | High fidelity vibratory source seismic method with source separation |
US5822269A (en) | 1995-11-13 | 1998-10-13 | Mobil Oil Corporation | Method for separation of a plurality of vibratory seismic energy source signals |
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 |
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 |
US5999489A (en) | 1997-03-21 | 1999-12-07 | Tomoseis Inc. | High vertical resolution crosswell seismic imaging |
US6014342A (en) | 1997-03-21 | 2000-01-11 | Tomo Seis, Inc. | Method of evaluating a subsurface region using gather sensitive data discrimination |
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 |
US6021379A (en) * | 1997-07-29 | 2000-02-01 | Exxon Production Research Company | Method for reconstructing seismic wavefields |
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 |
US6480790B1 (en) | 1999-10-29 | 2002-11-12 | Exxonmobil Upstream Research Company | Process for constructing three-dimensional geologic models having adjustable geologic interfaces |
FR2800473B1 (fr) | 1999-10-29 | 2001-11-30 | Inst Francais Du Petrole | Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques |
DE19954866A1 (de) | 1999-11-15 | 2001-05-31 | Infineon Technologies Ag | Verfahren zur Behandlung einer durch Epitaxie hergestellten Oberfläche eines SiC-Halbleiterkörpers und danach hergestellten Schottkykontakt |
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 |
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 |
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 |
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 |
EP2105765A1 (en) | 2008-03-28 | 2009-09-30 | Schlumberger Holdings Limited | Simultaneous inversion of induction data for dielectric permittivity and electric conductivity |
EP2265975A4 (en) | 2008-03-28 | 2017-05-24 | 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) |
US20110182141A1 (en) | 2008-08-14 | 2011-07-28 | Schlumberger Technology Corporation | Method and system for monitoring a logging tool position in a borehole |
US8559270B2 (en) | 2008-08-15 | 2013-10-15 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
CN102124374B (zh) | 2008-08-15 | 2013-07-17 | Bp北美公司 | 用于分离单独的同时震源的方法 |
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 |
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 |
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 |
US8576663B2 (en) | 2010-04-30 | 2013-11-05 | Schlumberger Technology Corporation | Multicomponent seismic inversion of VSP data |
KR101167715B1 (ko) * | 2010-04-30 | 2012-07-20 | 서울대학교산학협력단 | 복소 구배 최소자승법에의한 파형 역산을 이용한 지하 구조의 영상화 장치 및 방법 |
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 |
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 |
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 |
US8437998B2 (en) * | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
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 |
US20120316790A1 (en) * | 2011-06-08 | 2012-12-13 | Chevron U.S.A. Inc. | System and method for data inversion with phase extrapolation |
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 |
US20120316844A1 (en) * | 2011-06-08 | 2012-12-13 | Chevron U.S.A. Inc. | System and method for data inversion with phase unwrapping |
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 |
US9348050B2 (en) * | 2012-05-23 | 2016-05-24 | Exxonmobil Upstream Research Company | Near-surface noise prediction and removal for data recorded with simultaneous seismic sources |
US20140200816A1 (en) * | 2013-01-14 | 2014-07-17 | Westerngeco L.L.C. | Seismic data processing |
US9702993B2 (en) | 2013-05-24 | 2017-07-11 | Exxonmobil Upstream Research Company | Multi-parameter inversion through offset dependent elastic FWI |
US10459117B2 (en) * | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
US20140372043A1 (en) * | 2013-06-17 | 2014-12-18 | Wenyi Hu | Full Waveform Inversion Using Perfectly Reflectionless Subgridding |
WO2015073487A1 (en) * | 2013-11-12 | 2015-05-21 | Westerngeco Llc | Seismic data processing |
US20150185345A1 (en) * | 2013-12-31 | 2015-07-02 | Chevron U.S.A. Inc. | System and method for seismic imaging of a complex subsurface |
WO2015171215A1 (en) * | 2014-05-09 | 2015-11-12 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10401516B2 (en) * | 2014-10-24 | 2019-09-03 | Westerngeco L.L.C. | Travel-time objective function for full waveform inversion |
-
2015
- 2015-03-26 WO PCT/US2015/022786 patent/WO2015171215A1/en active Application Filing
- 2015-03-26 SG SG11201608175SA patent/SG11201608175SA/en unknown
- 2015-03-26 AU AU2015256626A patent/AU2015256626B2/en not_active Ceased
- 2015-03-26 US US14/670,039 patent/US9977142B2/en active Active
- 2015-03-26 EP EP15715936.9A patent/EP3140675A1/en not_active Withdrawn
- 2015-03-26 CA CA2947847A patent/CA2947847C/en not_active Expired - Fee Related
- 2015-03-26 KR KR1020167034474A patent/KR101915451B1/ko active IP Right Grant
- 2015-03-26 RU RU2016148177A patent/RU2649214C1/ru not_active IP Right Cessation
- 2015-03-26 MX MX2016013366A patent/MX2016013366A/es unknown
- 2015-03-26 CN CN201580022137.XA patent/CN106461802B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102906599A (zh) * | 2010-02-10 | 2013-01-30 | 埃克森美孚上游研究公司 | 用于全波场反演和逆时偏移中的地下参数估计的方法 |
CN103270430A (zh) * | 2011-06-08 | 2013-08-28 | 雪佛龙美国公司 | 用于地震数据反演的***和方法 |
CN103713315A (zh) * | 2012-09-28 | 2014-04-09 | 中国石油化工股份有限公司 | 一种地震各向异性参数全波形反演方法及装置 |
CN103091711A (zh) * | 2013-01-24 | 2013-05-08 | 中国石油天然气集团公司 | 全波形反演方法及装置 |
Non-Patent Citations (2)
Title |
---|
LIN YOUZUO 等: "Spatially-Variant Tikhonov Regularization for Double-Difference Waveform Inversion", 《TENTH ANNUAL CONFERENCE ON CARBON CAPTURE AND SEQUESTRATION》 * |
高凤霞 等: "几种优化方法在频率域全波形反演中的应用效果及对比分析研究", 《地球物理学进展》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150323689A1 (en) * | 2014-05-09 | 2015-11-12 | Yaxun Tang | Efficient Line Search Methods for Multi-Parameter Full Wavefield Inversion |
US9977142B2 (en) * | 2014-05-09 | 2018-05-22 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
CN107886571A (zh) * | 2017-11-03 | 2018-04-06 | 中原工学院 | 一种利用计算机多维空间的数学建模方法 |
Also Published As
Publication number | Publication date |
---|---|
EP3140675A1 (en) | 2017-03-15 |
RU2649214C1 (ru) | 2018-03-30 |
MX2016013366A (es) | 2017-01-26 |
AU2015256626B2 (en) | 2017-10-12 |
WO2015171215A1 (en) | 2015-11-12 |
CN106461802B (zh) | 2019-05-31 |
CA2947847C (en) | 2018-08-14 |
US9977142B2 (en) | 2018-05-22 |
US20150323689A1 (en) | 2015-11-12 |
KR101915451B1 (ko) | 2018-11-06 |
CA2947847A1 (en) | 2015-11-12 |
AU2015256626A1 (en) | 2016-11-17 |
KR20160149308A (ko) | 2016-12-27 |
SG11201608175SA (en) | 2016-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109709603A (zh) | 地震层位识别与追踪方法、*** | |
US8571842B2 (en) | Method of determining parameter from sparse measurement data | |
US9494711B2 (en) | Adaptive weighting of geophysical data types in joint inversion | |
CN106461802A (zh) | 用于多参数全波场反演的有效的线性搜索方法 | |
CN108645994A (zh) | 一种基于多点地质统计学的地质随机反演方法及装置 | |
CN108139498A (zh) | 具有振幅保持的fwi模型域角度叠加 | |
CN106324674A (zh) | 一种页岩气toc叠前地震反演预测方法 | |
Santiago et al. | A methodology for the characterization of flow conductivity through the identification of communities in samples of fractured rocks | |
Sajeva et al. | Estimation of velocity macro-models using stochastic full-waveform inversion | |
Roberts et al. | Joint stochastic constraint of a large data set from a salt dome | |
CN108351432A (zh) | 数据驱动的聚焦反演 | |
Ashena et al. | A Novel 2.5 D Deep Network Inversion of Gravity Anomalies to Estimate Basement Topography | |
Yang et al. | Enhancing data-driven seismic inversion using physics-guided spatiotemporal data augmentation | |
Masmoudi et al. | Multi-parameters scanning in HTI media | |
Zhdanov et al. | 3D potential field migration for rapid imaging of gravity gradiometry data-A case study from Broken Hill, Australia, with comparison to 3D regularized inversion | |
Niri* et al. | A multi-objective optimization method for creating reservoir models that simultaneously match seismic and geologic data | |
Jiang et al. | Imaging complex near-surface land area with joint traveltime and waveform inversion | |
Zou et al. | Machine learning based fluids and lithofacies prediction based on the integration of well logging data and seismic inversion | |
Perrone et al. | Waveform tomography based on local image correlations | |
Alfayez et al. | Semi-automated prestack seismic inversion workflow using temporal convolutional networks | |
Washbourne et al. | Crosswell seismic imaging in three dimensions | |
Phillips et al. | Machine Learning Applications and Examples for Natural Resources | |
Sun et al. | Enabling Uncertainty Quantification in a standard Full Waveform Inversion method using Normalizing Flows | |
Bijani et al. | 2-D tomography of first-arrivals using the genetic algorithm with elitism | |
Wei et al. | Training an Inverse Operator on Shingling Seismic Data based on cGANs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20190531 Termination date: 20210326 |