CN102239429A - Tti地球模型的多各向异性参数反演 - Google Patents

Tti地球模型的多各向异性参数反演 Download PDF

Info

Publication number
CN102239429A
CN102239429A CN200980148505XA CN200980148505A CN102239429A CN 102239429 A CN102239429 A CN 102239429A CN 200980148505X A CN200980148505X A CN 200980148505XA CN 200980148505 A CN200980148505 A CN 200980148505A CN 102239429 A CN102239429 A CN 102239429A
Authority
CN
China
Prior art keywords
data
model
tti
value
anisotropic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN200980148505XA
Other languages
English (en)
Other versions
CN102239429B (zh
Inventor
孙永和
王躍
徐佟
L·L·张
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Original Assignee
Chevron USA Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chevron USA Inc filed Critical Chevron USA Inc
Publication of CN102239429A publication Critical patent/CN102239429A/zh
Application granted granted Critical
Publication of CN102239429B publication Critical patent/CN102239429B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times

Landscapes

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

Abstract

一种用于确定倾斜横向各向同性(TTI)地球模型的各向异性模型参数的值的方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ和ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：为感兴趣地质体获取输入数据；确定输入数据与各向异性模型参数之间的理论关系；以及使用牵涉到下述过程的迭代或顺序组合的工作流，根据所述理论关系和所述输入数据，计算所述感兴趣地质体中的多个地下位置的每一个上的各向异性模型参数的值，所述过程包括输入数据预处理、传统层析反演、基于倾斜横向各向同性模型的三维层析反演、和使用倾斜横向各向同性模型的三维叠前深度偏移。

Description

TTI地球模型的多各向异性参数反演
技术领域
本发明涉及用于确定倾斜横向各向同性(Tilted TransverselyIsotropic，TTI)地球模型的各向异性模型参数的值的方法。
背景技术
为了分析地下地层的地质结构，勘探地球物理学家做了许多假设。其中之一是地下地层是各向同性的，而实际上，它基本上是各向异性的。这种错误假设可能导致地质结构的错误成像和解释。为了将地震处理技术扩展到各向异性介质，期望获得地质结构的各向异性的度量。
地震各向异性可以定义成地震速度与波传播方向的相关性。众所周知，倾斜轴横向各向同性地球模型或TTI地球模型可以用于在各向异性介质中模拟波的传播和获取地下地层的图像。描述TTI地球模型的物理参数包括(1)对称轴；(2)沿着对称轴的P-波(压缩)速度-Vp0；(3)为偏离对称轴的小角度-δ规定速度如何变化的参数；以及(4)在偏离对称轴的大角度-η上确定速度的参数(参见：Thomsen，″Weak Elastic Anisotropy″，Geophysics，vol.51，no.10，1986年10月，以及Alkhalifah和Tsvankin，″Velocity analysis for transverselyisotropic Media″，Geophysics，vol.60，1550-1566，1995)。
一些TTI地球模型还使用各向异性参数ε来描述各向异性介质中波的传播。参数ε满足如下关系η＝(ε-δ)/(1+2δ)。需要S(剪切)波速度来完整描述TTI地球模型，但在P-波处理中，S-波速度通常是使用与P-波速度的经验关系获得的。
通常，TTI地球模型是三向模型。模型中的每个点通过它的坐标和各向异性参数的值来描述。在某些情况下，如果各向异性介质的性质不是逐点变化的，则完全定义一个模型可能只需少量各向异性参数。但是，在大多数情况下，TTI地球模型需要各向异性参数的大量空间变化值来精确定义该模型。
TTI地球模型的各向异性参数可以从岩芯数据中直接测量。但是，钻井和取芯是非常昂贵的过程，直接测量只在很少的井位置上才有可能。对于3D成像，期望也使用横向扩展数据来确定TTI地球模型的各向异性参数。
发明内容
在本发明的一个方面中，提供了一种用于确定倾斜横向各向同性(TTI)地球模型的各向异性模型参数的值的计算机实现方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森(Thomsen)各向异性参数δ和ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：为感兴趣地质体获取输入数据；确定输入数据与各向异性参数之间的理论关系；以及使用牵涉到下述过程的迭代或顺序组合的工作流，根据所述理论关系和所述输入数据，计算所述感兴趣地质体中的多个地下位置的每一个上的各向异性模型参数的值，所述过程包括输入数据预处理、传统层析反演、基于倾斜横向各向同性模型的三维层析反演、和使用倾斜横向各向同性模型的三维叠前深度偏移(pre-stack depthmigration)。
在本发明的另一个方面中，提供了一种具有机器可执行指令的计算机产品，所述指令可被机器执行以执行确定TTI地球模型的各向异性参数的值的层析反演方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ和ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：确定输入数据与各向异性参数之间的关系，所述输入数据是为感兴趣地质体获取的；以及使用牵涉到下述过程的迭代或顺序组合的工作流，根据所述关系和所述输入数据，计算所述感兴趣地质体中的多个地下位置的每一个上的各向异性模型参数的值，所述过程包括输入数据预处理、三维层析反演、和三维TTI叠前深度偏移。
一种用于确定三维TTI地球模型的各向异性参数的值的计算机实现方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ，ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：(a)获取使共成像点道集(common-imaging-point gather)基本上平坦和使地震数据基本上与井数据相联系的初始偏移速度模型，所述初始偏移速度模型包括感兴趣地质体中的多个地下位置的每一个上的初始值Vp0，δ，ε(或η)；(b)将井位置上的检验炮数据，或VSP数据，或检验炮数据和VSP数据这两者输入三维层析反演中以确定井位置附近的Vp0的更新值，所述Vp0的值通过校正量ΔVp0来更新，其中Vp0＝Vp0+ΔVp0；(c)使用相对变化量Δδ＝(VP0)/Vp0确定δ(initial)的改进增量Δδ；(d)在考虑到地质一致性和正则化的情况下将相对变化量Δδ从近井位置外推到多个地下位置的每一个上的整个三维TTI地球模型，以确定δ的更新值，其中δ＝δ+Δδ；(e)使用多个地下位置的每一个上的外推三维Δδ来确定Vp0的更新值＝Vp0(1-Δδ)并获取三维扩展更新增量ΔVp0＝-ΔδVp0；(f)将利用通过Vp0和δ的更新值定义的改进偏移速度模型获取的共成像点道集中的从近到中偏距/角度剩余时差信息(near-to-mid-offset/angle residual moveout information)输入TTI层析反演过程中，以进一步提供多个地下位置的每一个上的δ的更新值；以及(g)将利用改进模型获取的共成像点道集中的从近到远偏距/角度剩余时差信息输入TTI层析反演过程中，以提供多个地下位置的每一个上的η的更新值。
当通过参照形成本说明书的一部分、相同标号在各个图形中表示相应部件的附图考虑如下描述和所附权利要求书，本发明的这些和其它目的、特征、和特性，以及操作方法、结构的相关元件的功能、部件的组合、和制造经济性将变得更加显而易见。但是，不言而喻，这些附图只用于例示和描述的目的，而无意作为限制本发明的定义。正如用在说明书和权利要求书中的那样，除非上下文另有明确指明，单数形式“一个”、“一种”、和“该”也包括复数指示物。
附图说明
图1示出了根据本发明的一个实施例确定倾斜横向各向同性(TTI)地球模型的各向异性模型参数的值的方法；
图2示出了根据本发明的一个实施例确定各向异性参数的值的工作流；以及
图3示出了可以根据本发明的一个实施例用在图2的工作流中的各种情形。
具体实施方式
本发明的实施例提供了一种用于精确确定感兴趣地质体中的多个位置的每一个上的倾斜横向各向同性(TTI)地球模型的各向异性参数的值的方法。在一个实施例中，使用井位置上的直接测量值校准估计的各向异性参数。本发明的实施例还涉及具有机器可执行指令的计算机产品，该指令可被机器执行以执行用于确定TTI地球模型的各向异性参数的值的方法。
在一个实施例中，用于精确确定TTI地球模型的各向异性参数的值的方法包括三个主要方面：(1)使用输入数据确定各向异性参数；(2)输入数据与各向异性参数之间的关系；以及(3)用于将输入数据转换成各向异性参数的值的过程。被用于将输入数据转换成各向异性参数的值的过程被称为反演过程(inversion process)。
现在参照图1，该图示出了根据本发明的一个实施例确定TTI地球模型的各向异性参数的值的方法10。借助于图1的方法，为TTI地球模型涵盖的感兴趣地质体中的每个坐标(x，y，z)确定各向异性参数(Vp0，δ，η和/或ε)的值。该方法从为感兴趣地质体获取输入数据的过程15开始。输入数据可以包括，例如，表面地震数据、垂直地震剖面数据、VSP(垂直地震剖面)数据、检验炮数据、测井数据、解释数据、区域趋势、先验数据、或前面数据的任何组合。解释数据包括层位位置(horizon position)和形状的挑选。
记录来自地球的反射的表面地震数据具有很大的覆盖范围，并且可容易地在大多数区域中获得。在不同角度上反射的波对各向异性参数敏感。因此，表面地震数据可以用作用于各向异性参数估计的输入数据。可测量数据包括不同偏距(offset)上的行进时间和不同偏距/角度上偏移道集的深度误差。记录来自不同方向的波的直接到达的VSP和检验炮数据也对各向异性参数敏感。测量值是不同钻孔位置上的行进时间。在现有技术中众所周知，VSP数据是通过将波接收器放在钻孔下面而将波发射器放在表面附近而获取的。检验炮数据是通过将波发射器和接收器两者都放在钻孔下面获取的。
在获得输入数据之后，该方法10转到过程20，其中确定输入数据与各向异性参数之间的理论关系。射线理论可以用于推导记录的行进时间(输入数据)与TTI地球模型的各向异性参数之间的理论关系。具体地说，首先，通过实现描述在如下文献中的TTI指定技术可以获取输入数据与各向异性参数之间的理论关系，以指定各向异性地球模型速度：Thomsen，″Weak elastic anisotropy，″Geophysics，vol.51，pgs.1954-66(1986)；和/或Alkahalifah等人，″Velocity analysis fortransversely isotropic media，″Geophysics，vol.60，pgs.1550-1566(1995)。然后，通过描述在如下文献中的一种或多种技术进行各向异性射线追踪：Cerveny，″Seismic rays and ray intensities in inhomo-geneous anisotropic media，″Geophysical Journal，vol.29，pgs.1-13(1972)和/或Gajewski等人，″Vertical seismic profile synthetics bydynamic ray tracing in laterally varying layered anisotropicstructures，″Journal of Geophysics Res.，vol.95，pgs.11301-11315，(1990)。有关使用射线理论推导输入数据与各向异性参数之间的理论关系的其它信息可以从2008年3月24日提交的发明名称为“偏移地震数据的系统和方法”的美国专利申请第12/079,170号中找到。
在确定了输入数据与各向异性参数之间的理论关系之后，该方法转到过程25，其中根据理论关系和输入数据计算感兴趣地质体中的多个地下位置的每一个上的各向异性模型参数的值。在一个实施例中，各向异性参数的计算使用牵涉到下述过程的迭代或顺序组合的工作流，所述过程包括输入数据预处理、传统层析反演或基于倾斜横向各向同性模型的三维层析反演、和使用倾斜横向各向同性模型的三维叠前深度偏移。并且，所述过程的迭代或顺序组合可以包括解释性挑选(interpretive picking)。
过程25的各向异性参数的估计是使用所记录数据和理论关系求解模型参数的计算过程。这种过程叫做层析反演。该反演可以使用各种工作流不同地执行。例如，可以顺序地(即，每次一个地)估计各向异性参数，或优选地，可以同时估计各种参数。也可以使用一种类型的数据估计一个或多个参数。在另一个实施例中，所有可用输入数据被用于同时估计所有各向异性参数。并且，各向异性模型参数的值可以使用具有输入数据的全集或子集和输出数据的全集或子集的反演的各种组合来确定。
过程25的层析反演是递归的和迭代的，可以使用三维倾斜横向各向同性射线追踪来模拟三维波传播。输入数据预处理可以包括在执行三维TTI层析反演之前对输入进行归一化。输入数据预处理还可以包括在执行三维TTI层析反演之前，使用3D TTI叠前深度偏移算法来偏移地震数据，将偏移地震数据按照它们的地下位置和它们的偏移偏距/角度(migration offset/angle)分类成共成像点(CIP)道集，并在共成像点道集域中选择和量化剩余时差。TTI叠前深度偏移算法在现有技术中是众所周知的。
各向异性参数的值通过迭代地选择Vp0，δ和η的值来计算。具体地说，计算各向异性参数的值直到(a)感兴趣地质体中的地下结构的地震图像位置与它们在钻井中识别的空间位置相联系，和(b)共成像点道集中的剩余时差在模型中的每个点上最小。地下结构包括岩石边界。
现在参照图2，这个图形示出了根据本发明的一个实施例确定倾斜横向各向同性(TTI)地球模型的各向异性参数的值的工作流100。该迭代工作流可以从测量三维地震反射数据的过程105开始。然后，该工作流转到将地震反射数据用于为感兴趣地质体确定初始TTI地球模型的过程110。初始TTI地球模型可以使用各种手段来确定。例如，初始TTI地球模型可以从速度信息和对测量的数据进行的各种处理，以及有关感兴趣地质体所在的区域的局部或全局知识中确定。而且，初始TTI地球模型构建还包括估计地下反射体结构倾角和获取感兴趣地质体中的TTI对称轴。在初始TTI地球模型中，在感兴趣地质体中的多个位置(x，y，z)的每一个上提供每个参数Vp0，δ和η的值。
在确定了初始TTI地球模型之后，该工作流转到框115和过程120，其中根据射线追踪，将井数据(框115)用于执行三维(3D)TTI层析(过程120)，以更新井附近的TTI地球模型的Vp0的值。井数据可以包括代表有关井中波速的信息的VSP数据或检验炮数据。井附近感兴趣地质体中的Vp0的更新值对应于Vp0(initial)+ΔVp0。借助于这种运算，在初始TTI地球模型中更新井附近Vp0的值。具体地说，在图2的实施例中使用VSP数据和检验炮数据将地震数据与井数据联系在一起。
在更新了井附近Vp0的值之后，该工作流100转到过程125，其中按照如下变换Δδ＝ΔVp0/Vp0以及δ(updated)＝δ(initial)+Δδ，使用Vp0的更新值更新井附近δ的初始值。过程110-125的变换使得能够通过提供井附近更精确的TTI模型来细化初始3D TTI地球模型。
一旦利用Vp0和δ的更新值在井附近更新了初始TTI地球模型，工作流100就转到过程130，其中填充模型中的其余点，即，井附近的位置以外的点。在图2的实施例中，使用三维外推过程和井附近的Vp0和δ的更新值填充模型中的其余点。三维外推过程通过首先外推模型的其余位置中的δ或Δδ，然后使用变换ΔVp0＝-Δδ＊Vp0获取Vp0的三维更新值来执行。应该理解，可以使用各种类型的算法来外推井位置以外的δ或Δδ的值，并获取Vp0和δ的扩展更新值。外推过程的结果(参见框135)提供了具有井联系(在模型的每个点上具有Vp0和δ的更新值)的三维TTI地球模型。框135的三维TTI地球模型包括模型的每个点(x，y，z)上的如下参数：Vp0(initial)+ΔVp0、δ(initial)+Δδ和η(initial)。
在过程140，将在框135上确定的模型用于执行地震反射数据的三维(3D)叠前TTI深度成像。借助于这个过程，通过模型处理/偏移已经为感兴趣地质体记录的地震反射数据，以提供感兴趣地质体地下新图像。应该理解，地震数据的偏移通常使得能够从偏移地震数据形成感兴趣地质体的图像，该图像比利用各向同性算法从偏移数据形成的图像更精确地描绘存在于感兴趣地震体中的地质特征。过程140的地震反射数据的三维叠前TTI深度成像提供了共成像点道集(框145)。在本领域中众所周知，共成像点道集对应于终止在感兴趣地质体中的相同图像位置上的偏移地震数据。
一旦通过3D叠前TTI深度成像分析识别了共成像点道集，工作流100接着转到执行剩余时差分析的过程150和/或155。在本领域中众所周知，剩余时差分析是速度模型细化中的基本步骤。该分析一般通过对内联线和参照线位置的网格使用相似性来执行。作为深度和偏距的函数，为每个共成像点道集计算一个相似性板。在过程150中为从近到中偏距执行对共成像点道集的剩余时差分析，而在过程155中为从近到远偏距执行对共成像点道集的剩余时差分析。近和远偏距分别对应于源与接收器之间的小距离和大距离。
应该理解，图2的过程150，155旨在涵盖多个时差分析。也就是说，在本发明的实施例中可以只为从近到中偏距(过程150)或只为从近到远偏距(过程155)执行图2的剩余时差分析。可替代地，在本发明的一个实施例中，如图2所示，可以同时为从近到中偏距和从近到远偏距执行剩余时差分析。
剩余时差分析为在过程135中定义的3D TTI模型的每个点提供了更新波速Vp0。在为模型的每个点确定了更新波速Vp0之后，在框165a和/或165b中确定是否达到收敛，即，过程150和/或过程155的剩余时差分析的结果是否在预定容限之内。如果达到收敛，则在框166中结束工作流100。
如果未达到收敛，则工作流100转到框170和/或框180，其中确定是否应该为模型的每个点更新δ和/或η的值。如果在框170询问的结果是肯定的，工作流就转到过程185，其中使用对从近到中偏距的剩余时差分析的结果执行三维(3D)TTI层析，以便为模型的每个点更新δ的值。只更新δ的值是非常有益的，因为这样的更新不会修改模型中波速Vp0的值。其结果是，在小倾角设置中不会在性质上破坏事先利用在框135中构建的模型获得的井数据与地震数据之间的联系。那样，可以为TTI模型的值实现更快速和更健壮的转换。
可替代地或附加地，确定是否应该为模型的每个点更新η的值。如果在框180上询问的结果是肯定的，工作流就转到过程195，其中使用对从远到中偏距的剩余时差分析的结果执行三维(3D)TTI层析，以便为模型的每个点更新η的值。
如果不更新δ和/或η的值，则确定是否应该更新波速Vp0的值(框175)。如果询问的结果是否定的，则在框166中结束工作流。如果询问的结果是肯定的，则工作流转到过程190，其中执行3D TTI层析以更新感兴趣地质体中速度Vp0的值。
在框197上，过程185，190和195的3D TTI层析分析的结果提供了在模型的每个点(x，y，z)上具有δ，Vp0和η的更新值的新三维(3D)TTI模型。然后，对框197的新模型中δ，Vp0和η的更新值确定是否达到收敛(框198)。可以使用各种判据来确定是否达到收敛。例如，确定获得的值是否在预定阈值之下。如果询问的结果是肯定的，则在过程199中结束图2的工作流。如果询问的结果是否定的，则工作流100返回到框115，其中根据射线追踪，将井数据用于执行三维(3D)TTI层析(过程120)，以更新井附近的在框197中获得的TTI地球模型的Vp0的值。因此，将在框197中获得的TTI地球模型用作新初始模型用于随后迭代。在一个实施例中，这个新初始TTI地球模型比用在第一次迭代中的初始模型更精细，因为它使共成像点道集基本上平坦和使地震数据基本上与井数据相联系。在确定了在井附近具有Vp0的更新值的修正模型之后，工作流100转到确定δ，和/或Vp0和/或η的更新值的过程125-198。然后迭代该工作流直到δ，和/或Vp0和/或η的值都达到收敛(即，δ，和/或Vp0和/或η的值在两次随后迭代之间基本上没有变化)。
应该理解，图2的工作流旨在包含优化三维(3D)TTI模型的若干情形。现在参照图3，这个图形示出了每次迭代时可以应用的各种情形。在第一种情形下，δ是在模型中的每个位置上更新的唯一参数。第一种情形是优选的，因为在小倾角设置中δ的更新不改变Vp0的值，其结果是，在优化δ期间不会破坏井数据与地震数据之间的定量联系。在第二种情形下，更新参数δ和Vp0。在第三种情形下，更新参数δ和η。在第四种情形下，更新速度Vp0。在第五种情形下，更新Vp0和η。在第六种情形下，更新η，以及在第七种情形下，更新δ，η和Vp0。
应该理解，在给定迭代中更新的一个或多个参数可能与在随后迭代中更新的一个或多个参数不同。换句话说，在每次迭代时优化的各向异性参数(δ，η和Vp0)可能不同。因此，在一个实施例中，设想可以使用各种情形来优化三维TTI模型。并且，设想在按照第一种情形优化了第一参数，例如δ之后，可以沿着图2的工作流优化第二参数，例如η，然后第三参数，例如Vp0。因此，可以同时或顺序优化参数δ，η和Vp0。并且，应该理解，每次迭代时更新的参数的选择在很大程度上取决于包括深度成像和层析的整个模型构建过程的成熟度和用于构成3D TTI模型的井数据和地震数据的类型和质量。
尽管为了例示的目的，根据当前认为是最实用优选实施例的实施例对本发明作了详细描述，但应该明白，这样的细节仅仅为了该目的，本发明不局限于所公开的实施例，而是相反，打算涵盖在所附权利要求书的精神和范围之内的所有修改和等同安排。
应该理解，确定倾斜横向各向同性(TTI)地球模型的各向异性模型参数的值所涉及的不同动作可以按照机器可执行指令或代码来执行。这些机器可执行指令可以嵌入数据存储媒体中。可以将处理器配置成执行这些指令。
包括可执行代码的牵涉到编程的计算机系统的软件功能可以用于实现上述模型。软件代码可由通用计算机执行。在运行时，可以将代码和可能相关数据记录存储在通用计算机平台内。但是，在其它时候，可以将软件存储在其它地方和/或进行传输以便装入适当通用计算机系统中。因此，上面讨论的实施例牵涉到具有至少一种机器可读媒体携带的代码的一个或多个模块形式的一种或多种软件或计算机产品。计算机系统的处理器对这样代码的执行使计算机平台能够基本上以在本文所讨论和例示的实施例中执行的方式实现各种功能。
正如本文所使用的那样，像计算机或机器“可读媒体”那样的术语指的是参与将指令提供给处理器加以执行的任何媒体。计算机或机器“可读媒体”可以广义地称为“计算机产品”。这样的媒体可以呈现许多形式，包括但不限于这些非易失性媒体、易失性媒体、和传输媒体。非易失性媒体包括，例如，像如上所讨论运行的任何计算机中的任何存储设备那样的光盘或磁盘。易失性媒体包括像计算机系统的主存储器那样的动态存储器。物理传输媒体包括同轴电缆、铜线和光纤，包括包含计算机系统内的总线的电线。载波传输媒体可以呈现像在射频(RF)和红外(IR)数据通信期间生成的那样的电或电磁信号，或声或光波的形式。因此，计算机可读媒体的常见形式包括，例如，软盘、柔性盘、硬盘、磁带、任何其它磁性媒体、CD-ROM、DVD、任何其它光学媒体、像打孔卡、纸带、带有穿孔图案的任何其它物体媒体那样的不常使用媒体、RAM、PROM、EPROM、FLASH-EPROM、任何其它存储芯片或盒、传输数据或指令的载波、传输这样载波的电缆或链路、或计算机可以从中读取或发送编程代码和/或数据的任何其它媒体。在将一条或多条指令的一个或多个序列传送到处理器加以执行时可能牵涉到这些计算机可读媒体形式的许多形式。
应理解，本发明构想可以尽可能地将任何实施例的一个或多个特征与任何其它实施例的一个或多个特征组合。

Claims (16)

1.一种确定倾斜横向各向同性(TTI)地球模型的各向异性模型参数的值的计算机实现方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ和ε(或η＝(ε-δ)/(1+2δ))，所述方法包含：
为感兴趣地质体获取输入数据；
确定输入数据与各向异性模型参数之间的理论关系；以及
使用牵涉到下述过程的迭代或顺序组合的工作流，根据所述理论关系和所述输入数据，计算所述感兴趣地质体中的多个地下位置的每一个上的各向异性模型参数的值，从而提供存在于所述感兴趣地质体中的地质特征的精确模型，所述过程包括输入数据预处理、传统层析反演、基于倾斜横向各向同性模型的三维层析反演、和使用倾斜横向各向同性模型的三维叠前深度偏移。
2.如权利要求1所述的方法，其中，所述过程的迭代或顺序组合包括解释性挑选。
3.如权利要求1所述的方法，其中，所述层析反演是递归的和迭代的，并且使用三维倾斜横向各向同性射线追踪来模拟三维波传播。
4.如权利要求1所述的方法，其中，所述输入数据包括表面地震数据、垂直地震剖面数据、检验炮数据、VSP数据、测井数据、解释数据、区域趋势、先验数据、或前面数据的任何组合。
5.如权利要求4所述的方法，其中，所述解释数据包括层位位置和形状的挑选。
6.如权利要求1所述的方法，进一步包含使用输入数据创建初始TTI模型，所述初始TTI模型是通过估计地下反射体结构倾角和获取感兴趣地质体中的TTI对称轴创建的。
7.如权利要求1所述的方法，进一步包含在执行反演之前归一化输入数据。
8.如权利要求1所述的方法，其中，所述输入数据包括地震数据，所述方法进一步包含，在执行三维TTI层析反演之前，使用3D TTI叠前深度偏移算法偏移地震数据，将偏移地震数据按照它们的地下位置和它们的偏移偏距/角度分类成共成像点(CIP)道集，并在共成像点道集域中选择和量化剩余时差。
9.如权利要求1所述的方法，其中，使用具有输入数据的全集或子集和输出数据的全集或子集的反演的各种组合，使用三维TTI层析反演同时地或者按顺序分别地计算多个地下位置的每一个上的各向异性模型参数的值。
10.如权利要求9所述的方法，其中，所述输出数据包括Vp0，δ和η，和/或ε。
11.如权利要求1所述的方法，进一步包含使用具有Vp0，δ和η的TTI模型执行TTI叠前深度偏移过程。
12.如权利要求1所述的方法，其中，所述计算包括通过迭代地选择Vp0，δ和η的值来计算各向异性参数的值，直到(a)感兴趣地质体中的地下结构的地震图像位置与它们在钻井中识别的空间位置相联系，和(b)共成像点道集中的剩余时差在模型中的每个点上最小。
13.如权利要求12所述的方法，其中，地下结构包括岩石边界。
14.一种确定三维TTI地球模型的各向异性模型参数的值的计算机实现方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ，ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：
a.获取使共成像点道集基本上平坦和使地震数据基本上与井数据相联系的初始TTI地球模型，初始偏移速度模型包括感兴趣地质体中的多个地下位置的每一个上的初始值Vp0，δ，ε(或η)；
b.将井位置上的检验炮数据、或VSP数据、或检验炮数据和VSP数据这两者输入三维层析反演中以确定井位置附近的Vp0的更新值，所述Vp0的值通过校正量ΔVp0来更新，其中Vp0＝Vp0+ΔVp0；
c.使用相对变化量Δδ＝(Vp0)/Vp0确定δ(initial)的改进增量Δδ；
d.在考虑到地质一致性和正则化的情况下将相对变化量Δδ从近井位置外推到多个地下位置的每一个上的整个三维TTI地球模型，以确定δ的更新值，其中δ＝δ+Δδ；
e.使用多个地下位置的每一个上的外推三维Δδ来确定Vp0的更新值＝Vp0(1-Δδ)并获取三维扩展更新增量ΔVp0＝-ΔδVp0；
f.将利用通过Vp0和δ的更新值定义的改进偏移速度模型获取的共成像点道集中的从近到中偏距/角度剩余时差信息输入TTI层析反演过程中，以进一步提供多个地下位置的每一个上的δ的更新值；以及
g.将利用改进模型获取的共成像点道集中的从近到远偏距/角度剩余时差信息输入TTI层析反演过程中，以提供多个地下位置的每一个上的η的更新值。
15.如权利要求14所述的方法，进一步包含迭代地重复(a)，(b)，(c)，(d)，(e)，(f)和(g)，其中在迭代结束时确定的TTI模型的改进参数对应于下一次迭代时TTI模型的初始参数。
16.一种具有机器可执行指令的计算机产品，所述指令可被机器执行以执行用于确定TTI地球模型的各向异性参数的值的层析反演方法，所述各向异性参数包括沿着倾斜对称轴的P-波速度(Vp0)、和作为偏离对称轴的波传播角度的函数的代表波速度变化的汤姆森各向异性参数δ和ε(或η＝(ε-δ)/(1+2δ))，所述方法包括：
确定输入数据与各向异性参数之间的关系，所述输入数据是为感兴趣地质体获取的；以及
使用牵涉到下述过程的迭代或顺序组合的工作流，根据所述关系和所述输入数据，计算所述感兴趣地质体中的多个地下位置的每一个上的各向异性参数的值，所述过程包括输入数据预处理、三维层析反演、和三维TTI叠前深度偏移。
CN200980148505XA 2008-12-03 2009-10-22 Tti地球模型的多各向异性参数反演 Expired - Fee Related CN102239429B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US12/327,671 2008-12-03
US12/327,671 US20100135115A1 (en) 2008-12-03 2008-12-03 Multiple anisotropic parameter inversion for a tti earth model
PCT/US2009/061612 WO2010065209A2 (en) 2008-12-03 2009-10-22 Multiple anisotropic parameter inversion for a tti earth model

Publications (2)

Publication Number Publication Date
CN102239429A true CN102239429A (zh) 2011-11-09
CN102239429B CN102239429B (zh) 2013-11-20

Family

ID=42222700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200980148505XA Expired - Fee Related CN102239429B (zh) 2008-12-03 2009-10-22 Tti地球模型的多各向异性参数反演

Country Status (9)

Country Link
US (2) US20100135115A1 (zh)
EP (1) EP2370924A2 (zh)
CN (1) CN102239429B (zh)
AU (1) AU2009322868A1 (zh)
BR (1) BRPI0922702A2 (zh)
CA (1) CA2745076A1 (zh)
EA (1) EA021972B1 (zh)
SG (1) SG171913A1 (zh)
WO (1) WO2010065209A2 (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104076391A (zh) * 2014-04-16 2014-10-01 孙赞东 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法
CN105493100A (zh) * 2013-08-29 2016-04-13 兰德马克绘图国际公司 静态地球模型校准方法和系统
CN106709031A (zh) * 2016-12-29 2017-05-24 天津瀚海星云数字科技有限公司 高精度3D GIS Earch的实现方法及系统
CN107966729A (zh) * 2016-10-19 2018-04-27 中国石油化工股份有限公司 一种三维tti介质射线追踪方法及系统
CN108732623A (zh) * 2017-04-17 2018-11-02 中国石油化工股份有限公司 Tti介质层析反演成像方法及计算机可读存储介质
CN109581499A (zh) * 2018-11-08 2019-04-05 成都捷科思石油天然气技术发展有限公司 一种利用各向异性叠前深度偏移生成构造图的方法
CN109581521A (zh) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 Tti各向异性的局部层析方法及系统
CN109655889A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 一种各向异性参数联合反演方法及系统
CN113009580A (zh) * 2019-12-20 2021-06-22 中国石油天然气集团有限公司 变偏移距vsp初至反演vti各向异性参数方法和装置

Families Citing this family (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2011252076B2 (en) * 2010-05-12 2014-04-10 Shell Internationale Research Maatschappij B.V. Seismic P-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis
EP2591381A4 (en) * 2010-07-08 2015-12-02 Geco Technology Bv METHODS AND DEVICES FOR TRANSFORMING COLLECTED DATA TO IMPROVE VISUALIZATION CAPACITY
US20120099396A1 (en) * 2010-10-22 2012-04-26 Chevron U.S.A. Inc. System and method for characterization with non-unique solutions of anisotropic velocities
CN102141633B (zh) * 2010-12-10 2012-08-29 中国科学院地质与地球物理研究所 各向异性三维叠前时间偏移方法
US8861309B2 (en) * 2011-01-31 2014-10-14 Chevron U.S.A. Inc. Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context
FR2976088B1 (fr) * 2011-05-30 2014-03-07 Cggveritas Services Sa Procede de tomographie non lineaire pour un axe de symetrie principal d'un modele de vitesse anisotrope et dispositif
US9207342B2 (en) * 2012-03-09 2015-12-08 Chevron U.S.A. Inc. Correction of shear log for elastic anisotropy
CN102721977B (zh) * 2012-05-31 2014-08-06 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法
CN102830433B (zh) * 2012-08-16 2014-10-01 中国石油化工股份有限公司 基于频率域用偶极横波测井资料计算各向异性方位角的方法
RU2596182C1 (ru) * 2012-08-17 2016-08-27 Лэндмарк Графикс Корпорейшн Система и способ для отображения сейсмических данных
CA2906362A1 (en) * 2013-04-02 2014-10-09 Halliburton Energy Services, Inc. Anisotropy analysis using direct and reflected arrivals in seismic survey data
US20140301165A1 (en) * 2013-04-03 2014-10-09 Westerngeco L.L.C. Seismic data processing using joint tomography
CN104216011B (zh) * 2013-06-05 2017-08-04 上海青凤致远地球物理地质勘探科技有限公司 一种TTI介质稳定的qP波逆时偏移方法
EP2997401A4 (en) * 2013-06-18 2017-02-22 Halliburton Energy Services, Inc. Methods and systems for seismic data analysis using a tilted transversely isotropic (tti) model
CN103352691B (zh) * 2013-07-05 2015-11-11 天津大学 一种正交偶极子声波测井接收声系装置
US10529144B2 (en) 2013-08-23 2020-01-07 Landmark Graphics Corporation Local updating of 3D geocellular model
CN104570119B (zh) * 2013-10-29 2017-11-17 中国石油化工股份有限公司 一种三维垂直地震剖面反射波拉伸校正方法
US20160187512A1 (en) * 2014-01-14 2016-06-30 Cgg Services Sa Detecting and estimating anisotropy errors using full waveform inversion and ray based tomography
US10267937B2 (en) 2014-04-17 2019-04-23 Saudi Arabian Oil Company Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data
US9562983B2 (en) * 2014-04-17 2017-02-07 Saudi Arabian Oil Company Generating subterranean imaging data based on vertical seismic profile data
WO2015164917A1 (en) * 2014-05-01 2015-11-05 Curtin University Of Technology A method of estimating a volume fraction of coal in a layered geological structure
DE112015006388T5 (de) * 2015-05-01 2017-12-14 Halliburton Energy Services, Inc. Schätzung anisotroper Parameter aus Walkaway-VSP-Daten mittels Differenzialevolution
US9442207B1 (en) * 2015-06-12 2016-09-13 Chevron U.S.A. Inc. System and method for computing residual moveout from seismic images
WO2018080460A1 (en) * 2016-10-25 2018-05-03 Landmark Graphics Corporation Estimating interval anisotropy parameter for pre-stack depth migration using a least-squares method
CN106896408B (zh) * 2017-03-23 2018-10-16 中国石油天然气股份有限公司 一种角度域叠前时间偏移方法
US10901103B2 (en) 2018-03-20 2021-01-26 Chevron U.S.A. Inc. Determining anisotropy for a build section of a wellbore
CN108710148B (zh) * 2018-05-29 2019-05-24 中国科学院地质与地球物理研究所 三维倾角域稳相叠前深度偏移方法和装置
CN111596346B (zh) * 2019-02-20 2023-04-25 中国石油天然气集团有限公司 弹性波速度反演方法和装置
CN110161562B (zh) * 2019-04-29 2020-09-18 中南大学 倾斜横向各向同性介质中的裂缝弱度反演方法及系统
CN112444861B (zh) * 2019-08-27 2024-06-18 中国石油化工股份有限公司 速度模型更新方法、计算机存储介质及计算机系统
EP3859405A1 (en) * 2020-01-28 2021-08-04 Repsol Exploración, S.A. Computer-implemented method for determining a velocity image of a domain of the subsurface structural geology in an oil and gas reservoir
CN113296153B (zh) * 2020-02-24 2023-04-25 中国石油天然气集团有限公司 轴对称介质各向异性参数确定方法及装置
CN111624647B (zh) * 2020-06-05 2022-06-24 中油奥博(成都)科技有限公司 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置
CN113970784B (zh) * 2020-07-22 2024-09-06 中国石油化工股份有限公司 一种用于更新地震各向异性参数的方法及系统
CN112130202B (zh) * 2020-08-14 2023-06-30 中国石油天然气集团有限公司 一种正交各向异性速度反演的方法及系统
CN114185092B (zh) * 2020-09-15 2024-05-07 中国石油化工股份有限公司 储层水平缝发育程度评价方法、装置、电子设备及介质
CN112305589A (zh) * 2020-09-22 2021-02-02 中国石油天然气集团有限公司 各向异性介质深度域成像的方法及装置
CN112540407B (zh) * 2020-12-01 2023-04-25 中国煤炭地质总局地球物理勘探研究院 一种叠前深度偏移各向异性场建立方法
CN114594515B (zh) * 2020-12-07 2024-03-29 中国石油化工股份有限公司 基于缓变各向异性的井控速度反演方法
CN113221228B (zh) * 2021-06-04 2022-09-16 中国电建集团成都勘测设计研究院有限公司 水电站地下洞群围岩力学参数反演方法
CN113568045B (zh) * 2021-07-26 2022-11-11 中南大学 一种tti介质裂缝倾角的反演方法
CN114325832B (zh) * 2021-12-23 2024-10-11 太原理工大学 一种裂缝参数和弹性参数同步反演方法及系统
CN117471538A (zh) * 2022-07-22 2024-01-30 大庆油田有限责任公司 多井联控VTI介质各向异性参数δ场迭代更新方法
CN115951405B (zh) * 2023-03-08 2023-05-09 中国地震局地球物理研究所 一种基于稀疏化卷积和自适应矩估计的速度变化成像方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
US6920084B2 (en) * 2003-01-15 2005-07-19 Western Geco, L.L.C. Offset transformation to zero dip that preserves angle of incidence
US7072767B2 (en) * 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
US20080273421A1 (en) * 2007-03-05 2008-11-06 Zvi Koren Model-based time-preserving tomography

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6714480B2 (en) * 2002-03-06 2004-03-30 Schlumberger Technology Corporation Determination of anisotropic moduli of earth formations
US7924652B2 (en) * 2007-06-01 2011-04-12 Baker Hughes Incorporated Method for determining seismic anisotropy
CA2710607A1 (en) * 2008-02-28 2009-09-03 Exxonmobil Upstream Research Company Rock physics model for simulating seismic response in layered fractured rocks
US8750074B2 (en) * 2008-11-25 2014-06-10 Baker Hughes Incorporated Anisotropic parameter determination
US8798967B2 (en) * 2011-03-30 2014-08-05 Chevron U.S.A. Inc. System and method for computations utilizing optimized earth model representations

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6920084B2 (en) * 2003-01-15 2005-07-19 Western Geco, L.L.C. Offset transformation to zero dip that preserves angle of incidence
US7072767B2 (en) * 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
US20080273421A1 (en) * 2007-03-05 2008-11-06 Zvi Koren Model-based time-preserving tomography

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LAXMIDHAR BEHERA ET AL.: "Migration velocity analysis for tilted transversely isotropic media", 《GEOPHYSICAL PROSPECTING》 *
SYLVESTRE CHARLES ET AL: "To TTI or not to TTI? SemiAutomated 3D Tomographic Velocity Analysis in the Canadian Foothills: A Case History", 《2006 CSPG-CSEG-CWLS CONVENTION》 *
YAPING ZHU ET AL.: "Physical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media", 《GEOPHYSICS》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105493100A (zh) * 2013-08-29 2016-04-13 兰德马克绘图国际公司 静态地球模型校准方法和系统
CN104076391A (zh) * 2014-04-16 2014-10-01 孙赞东 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法
CN104076391B (zh) * 2014-04-16 2015-12-02 孙学凯 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法
CN107966729A (zh) * 2016-10-19 2018-04-27 中国石油化工股份有限公司 一种三维tti介质射线追踪方法及系统
CN106709031A (zh) * 2016-12-29 2017-05-24 天津瀚海星云数字科技有限公司 高精度3D GIS Earch的实现方法及系统
CN108732623A (zh) * 2017-04-17 2018-11-02 中国石油化工股份有限公司 Tti介质层析反演成像方法及计算机可读存储介质
CN109581521A (zh) * 2017-09-28 2019-04-05 中国石油化工股份有限公司 Tti各向异性的局部层析方法及系统
CN109655889A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 一种各向异性参数联合反演方法及系统
CN109655889B (zh) * 2017-10-11 2021-04-02 中国石油化工股份有限公司 一种各向异性参数联合反演方法及系统
CN109581499A (zh) * 2018-11-08 2019-04-05 成都捷科思石油天然气技术发展有限公司 一种利用各向异性叠前深度偏移生成构造图的方法
CN113009580A (zh) * 2019-12-20 2021-06-22 中国石油天然气集团有限公司 变偏移距vsp初至反演vti各向异性参数方法和装置

Also Published As

Publication number Publication date
US8666668B2 (en) 2014-03-04
US20120185171A1 (en) 2012-07-19
CA2745076A1 (en) 2010-06-10
US20100135115A1 (en) 2010-06-03
BRPI0922702A2 (pt) 2016-01-05
EA021972B1 (ru) 2015-10-30
WO2010065209A3 (en) 2010-08-12
SG171913A1 (en) 2011-07-28
CN102239429B (zh) 2013-11-20
AU2009322868A1 (en) 2010-06-10
EP2370924A2 (en) 2011-10-05
WO2010065209A2 (en) 2010-06-10
EA201170747A1 (ru) 2011-12-30

Similar Documents

Publication Publication Date Title
CN102239429B (zh) Tti地球模型的多各向异性参数反演
Taillandier et al. First-arrival traveltime tomography based on the adjoint-state method
EP2846175B1 (en) Seismic survey analysis
Hansen et al. Linear inverse Gaussian theory and geostatistics
US8694299B2 (en) Artifact reduction in iterative inversion of geophysical data
CN101999086B (zh) 用于确定地震数据质量的方法
US20130223187A1 (en) Geological Structure Contour Modeling and Imaging
US7768870B2 (en) Method for adjusting a seismic wave velocity model according to information recorded in wells
WO2017035104A1 (en) Velocity model seismic static correction
MX2011007084A (es) Adquisicion de datos y migracion antes del apilamiento con base en analisis de visibilidad sismica.
US20110098996A1 (en) Sifting Models of a Subsurface Structure
Wang et al. Multiparameter TTI tomography of P-wave reflection and VSP data
KR101931092B1 (ko) 동시 및 순차 소스 방법을 이용하는 풀 파형 반전을 위한 하이브리드 방법
Nguyen et al. Reconstructing disturbance zones ahead of the tunnel face by elastic waveform inversion supported by a parametric level-set representation
US9103933B2 (en) Estimating a property by assimilating prior information and survey data
Pan et al. Random objective waveform inversion of surface waves
US6868039B2 (en) Method for determining a seismic wave velocity model in a heterogeneous underground
CN108508481B (zh) 一种纵波转换波地震数据时间匹配的方法、装置及系统
GB2503640A (en) Quality Assurance in a Full Waveform Inversion Process
US20230099919A1 (en) System and method for stochastic full waveform inversion
Zhou First-break vertical seismic profiling tomography for Vinton Salt Dome
CN102792187A (zh) 基于零化子的波反演
Duchesne et al. Improving seismic velocity estimation for 2D poststack time migration of regional seismic data using kriging with an external drift
CN104345337B (zh) 一种用于地震反演的时控储层参数建模方法
Jardin et al. Depth seismic imaging using reflection and first arrival traveltime tomography: Application to a deep profile across the Northern Emirates Foothills

Legal Events

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

Granted publication date: 20131120

Termination date: 20151022

EXPY Termination of patent right or utility model