一种三维TTI介质射线追踪方法及系统

摘要

本发明公开了一种三维TTI介质射线追踪方法及系统，该方法包括基于声波近似方程推导出三维TTI介质下的程函方程；将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；基于解析后的程函方程建立三维TTI介质下的射线追踪方程；基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。

权利要求

1.一种三维TTI介质射线追踪方法，其特征在于，所述方法包括以下步骤：步骤一、基于声波近似方程推导出三维TTI介质下的程函方程，所述三维TTI介质下的程函方程为：其中，表示程函方程，τ表示射线旅行时， 表示表示射线在z轴上的位移量， 表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数；

步骤二、将所述三维TTI介质下的程函方程中的系数a0,a1,a2,a3,a4分别解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，其中，v表示动校正速度， η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角；

步骤三、基于解析后的程函方程建立三维TTI介质下的射线追踪方程；

步骤四、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。

2.如权利要求1所述的方法，其特征在于，根据如下表达式建立所述三维TTI介质下的射线追踪方程：其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度， 分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量， 分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量， 分别表示射线在x,y,z轴上的位移量。

3.如权利要求2所述的方法，其特征在于，在步骤四之前还包括对所述射线追踪方程进行初始化：设置初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度；

根据所述初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度获取相速度。

4.如权利要求3所述的方法，其特征在于，根据下式获取所述相速度V(θ′,θ,φ,φ′)：其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；

ε表示Thomsen各向异性参数。

5.如权利要求1至4中任意一项所述的方法，其特征在于，所述记录射线追踪到目标位置的射线旅行时，包括：设置追踪步长；

对于地下任一目标位置，记录射线追踪的起始位置到所述目标位置射线追踪使用的总的步长数目作为射线旅行时。

6.一种三维TTI介质射线追踪系统，其特征在于，所述系统包括：程函方程建模模块，其基于声波近似方程推导出三维TTI介质下的程函方程，所述三维TTI介质下的程函方程为：其中， 表示程函方程，τ表示射线旅行时， 表示表示射线在z轴上的位移量， 表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数；

程函方程解析模块，其将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，其中，v表示动校正速度， η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角；

射线追踪方程建模模块，其基于解析后的程函方程建立三维TTI介质下的射线追踪方程；

射线追踪模块，基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。

7.如权利要求6所述的系统，其特征在于，所述射线追踪方程建模模块基于解析后的程函方程建立三维TTI介质下的射线追踪方程为：其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度， 分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量， 分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量， 分别表示射线在x,y,z轴上的位移量。

说明书

一种三维TTI介质射线追踪方法及系统

技术领域

[0001] 本发明涉及地震资料处理领域，尤其涉及一种三维TTI介质射线追踪方法及系统。

背景技术

[0002] 在地震资料处理领域，对于横向各向同性(Transverse Isotropic，TI)介质，理想情况下，若介质对称轴是垂直的，则该横向各向同性介质为VTI(Transversely Isotropy With A Vertical Axis of Symmetry)介质，若介质对称轴是水平的，则该横向各向同性介质为HTI(Transversely Isotropy With A Horizontal Axis of Symmetry)介质。但是在实际情况中，TI介质的对称轴往往是倾斜的，此时称该横向各向同性介质称为TTI(Transversely Isotropy With A Tilt Axis of Symmetry)介质。

[0003] 目前的地震资料处理中，采用基于各向同性介质和VTI介质的理论假设前提条件对地下介质中地震波的传播进行简化表达，从而导致地震偏移成像数据处理中目的层的构造成像位置错误，甚至造成无论在深度上还是纵向上皆与测井数据不吻合的现象，给后期的圈闭规模计算、井位预测埋下隐患。

[0004] 因此，亟需一种算法稳健、精度较高的地震资料处理方法。

发明内容

[0005] 本发明所要解决的技术问题是提供一种算法稳健、精度较高的地震资料处理方法。

[0006] 为了解决上述技术问题，本发明提供了一种三维TTI介质射线追踪方法，该方法包括了以下步骤：步骤一、基于声波近似方程推导出三维TTI介质下的程函方程；步骤二、将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；步骤三、基于解析后的程函方程建立三维TTI介质下的射线追踪方程；步骤四、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。

[0007] 优选地，所述三维TTI介质下的程函方程 为：

[0008]

[0009] 其中，τ表示射线旅行时， 表示表示射线在z轴上的位移量， 表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数。

[0010] 优选地，将所述三维TTI介质下的程函方程中的系数a0,a1,a2,a3,a4分别解析为以下参数方程：

[0011]

[0012]

[0013]

[0014]

[0015]

[0016] 其中，v表示动校正速度， η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角。

[0017] 优选地，根据如下表达式建立所述三维TTI介质下的射线追踪方程：

[0018]

[0019]

[0020] 其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度， 分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量， 分别表示射线在x,y,z轴上的位移量。

[0021] 优选地，在步骤四之前还包括对所述射线追踪方程进行初始化：设置初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度；根据所述初始射线追踪的相角以及所述沿TTI介质对称轴方向传播的纵波地震波速度获取相速度。

[0022] 优选地，根据下式获取所述相速度V(θ′,θ,φ,φ′)：

[0023]

[0024] 其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为[0025]

[0026] 系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，

[0027] F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，[0028] D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；

[0029] ε表示Thomsen各向异性参数。

[0030] 优选地，所述记录射线追踪到目标位置的射线旅行时，包括：设置追踪步长；对于地下任一目标位置，记录射线追踪的起始位置到所述目标位置射线追踪使用的总的步长数目作为射线旅行时。

[0031] 本发明还提供了一种三维TTI介质射线追踪系统，所述系统包括：程函方程建模模块，其基于声波近似方程推导出三维TTI介质下的程函方程；程函方程解析模块，其将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式；射线追踪方程建模模块，其基于所述程函方程建立三维TTI介质下的射线追踪方程；射线追踪模块，其基于上述解析后的三维TTI介质下的程函方程以及射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的旅行时。

[0032] 优选地，所述程函方程建模模块基于所述程函方程建立三维TTI介质下的射线追踪方程 为：

[0033]

[0034] 其中，τ表示射线旅行时， 表示表示射线在z轴上的位移量， 表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数,如以下参数方程所示：

[0035]

[0036]

[0037]

[0038]

[0039]

[0040] 其中，v表示动校正速度， η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴的倾角，φ表示TTI介质对称轴方位角。

[0041] 优选地，所述射线追踪方程建模模块建立三维TTI介质下的射线追踪方程为：

[0042]

[0043]

[0044] 其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度， 分别表示射线追踪过程中x,y,z轴方向射线的角度信息， 表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量，分别表示射线追踪单位时间步长对应的射线参数在x,y,z三个方向的改变量， 分别表示射线在x,y,z轴上的位移量。

[0045] 与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：

[0046] 本发明实施例从声波近似方程出发，推导出三维TTI介质下的程函方程，并将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，进而基于解析后的程函方程建立三维TTI介质下的射线追踪方程并根据射线追踪方程进行射线追踪，该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。

[0047] 本发明的其他优点、目标，和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书，权利要求书，以及附图中所特别指出的结构来实现和获得。

附图说明

[0048] 附图用来提供对本申请的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本申请实施例的附图与本申请的实施例一起用于解释本申请的技术方案，但并不构成对本申请技术方案的限制。

[0049] 图1为根据本发明一实施例的三维TTI介质射线追踪方法的流程示意图；

[0050] 图2a-图2h为根据本发明一实施例的SEG管道模型示意图；

[0051] 图3a-图3f为根据本发明一实施例的盐丘模型示意图；

[0052] 图4为根据本发明另一实施例的三维TTI介质射线追踪系统的结构示意图。

具体实施方式

[0053] 以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。

[0054] 图1为根据本发明一实施例的三维TTI介质射线追踪方法的流程示意图，如图所示，包括以下步骤：

[0055] 步骤S110、基于声波近似方程推导出三维TTI介质下的程函方程。

[0056] 步骤S120、将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式。

[0057] 步骤S130、基于解析后的程函方程建立三维TTI介质下的射线追踪方程。

[0058] 步骤S140、基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的旅行时。

[0059] 具体的，在步骤S110中，基于声波近似方程推导出的TTI介质下的程函方程，如表达式(1)所示：

[0060]

[0061] 其中，表示程函方程，τ表示射线旅行时， 表示表示射线在z轴上的位移量，表示z轴方向上的射线参数的分量，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，a0,a1,a2,a3,a4均为方程的系数。

[0062] 接下来，在步骤S120中，将上述程函方程中的a0,a1,a2,a3,a4五个参数分别解析为下列参数方程：

[0063]

[0064]

[0065]

[0066]

[0067]

[0068] 其中，v表示动校正速度， η为非椭圆率，ε、δ表示Thomsen各向异性参数，θ表示TTI介质对称轴倾角，φ表示TTI介质对称轴方位角。

[0069] 进一步的，在上述程函方程和参数方程中，令 其中pz，py，px分别表示射线参数p在x,y,z轴方向的分量， 分别表示射线在x,y,z轴上的位移量。

[0070] 对表达式(1)所示的程函方程进行化简，过程如下：

[0071] 在 中，将相应的a0，a1，a2，a3，a4五个参数方程中的 用 进行替换。

[0072] 接下来，将a0，a1，a2，a3，a4五个参数的表达式的各项用b系列参数进行化简，其具体方法如下：

[0073] 对参数a0表达式中各项用b系列参数表示为：

[0074]

[0075] b02＝v2(1+2η)cos2θ(cosφpx-sinφpy)2

[0076]

[0077] b04＝v2(1+2η)(sinφpx+cosφpy)2

[0078]

[0079] 则将a0化简为：a0＝b01-b02+b03-b04-b05+1，

[0080] 同理，对参数a1表达式中各项用b系列参数表示为：

[0081]

[0082] b12＝sin2θcosφpx(2vtsinφpy-1)

[0083] b13＝vt(px)2(cos2θcos2φ+sin2φ)

[0084] b14＝py(vtpy(cos2θsin2φ+cos2φ)+sin2θsinφ)

[0085] b15＝vt(2v2η(b12+b13+b14)-vt)+v2(1+2η)

[0086] 则将a1化简为：a1＝b11*b15。

[0087] 同理，对参数a2表达式中各项用b系列参数表示为：

[0088] b21＝-4sin2θ(3cos(2θ)+2)cos(2φ)((px)2-(py)2)

[0089] b22＝8sin2θ(3cos(2θ)+2)sin(2φ)pxpy

[0090] b23＝(2cos(2θ)+3cos(4θ)+3)((px)2+(py)2)

[0091] 则将a2化简为：

[0092] a3，a4暂可不用化简，即：

[0093]

[0094]

[0095] 在步骤S130中，基于解析后的程函方程建立三维TTI介质下的射线追踪方程的，如表达式(2)所示：

[0096]

[0097] 其中，dx表示单位时间步长射线追踪的路径长度在x轴方向的投影，ds表示单位时间步长射线所追踪的路径长度， 分别表示射线追踪过程中x,y,z轴方向射线的角度信息，表示程函方程，px,py,pz分别表示射线参数p在x,y,z轴上的分量， 分别表示射线追踪单位时间步长对应的射线参数在x,y,z轴三个方向的改变量， 分别表示射线在x,y,z轴上的位移量。

[0098] 在求解上述方程中的 时， 的求解步骤如下：

[0099] 因为在常规射线追踪系统中

[0100] 而 依 次 求 解 方 程 中即可得到 值，同理可求得

[0101] 的具体计算过程如下：

[0102]

[0103]

[0104] 进一步的，依次求解方程中 即：

[0105]

[0106]

[0107]

[0108]

[0109]

[0110] 同理， 具体计算过程如下：

[0111]

[0112]

[0113]

[0114]

[0115]

[0116]

[0117]

[0118] 具体计算过程如下：

[0119]

[0120]

[0121]

[0122] 在求解上述方程中的 时， 求解步骤如下：

[0123] 因为在常规射线追踪系统

[0124] 而

[0125] 依次求解方程中 即可得到 值，同理可求得

[0126] 的具体计算过程如下：

[0127]

[0128]

[0129]

[0130]

[0131]

[0132]

[0133]

[0134] 的具体计算过程如下：

[0135]

[0136]

[0137]

[0138]

[0139]

[0140]

[0141]

[0142]

[0143] 的具体计算过程如下：

[0144]

[0145]

[0146]

[0147]

[0148]

[0149]

[0150] 接下来，在得到三维TTI介质下的射线追踪方程之后，在步骤S140中，基于射线追踪方程进行射线追踪，具体包括基于射线追踪方程获取射线路径与射线参数，并在获取射线路径的同时记录射线追踪到目标位置的射线旅行时。

[0151] 在步骤S140之前，还包括对射线追踪方程进行初始化。

[0152] 具体的，设置初始射线追踪的相角以及TTI介质对称轴方向传播的纵波地震波速度vt，根据初始射线追踪的相角以及TTI介质对称轴方向传播的纵波地震波速度获取相速度，相速度如表达式(3)所示：

[0153]

[0154] 其中，θ为TTI介质对称轴倾角，θ′为TTI介质初始出射角，φ为TTI介质对称轴方位角，φ′为射线方位角，vt表示沿TTI介质对称轴方向传播的纵波地震波速度，相角γ表示为：

[0155]

[0156] 系数E＝sinθ′sinθcos(φ′-φ)+cosθ′cosθ，

[0157] F＝(sinθ′cosθcos(φ′-φ)-cosθ′sinθ)2+sin2θ′sin2(φ′-φ)，[0158] D＝[(1+2ε)F-E2]2+4(1+2δ)E2F；

[0159] ε表示Thomsen各向异性参数。

[0160] 进一步地，射线旅行时的记录方法如下：设置固定追踪步长，对于地下任一目标位置，记录射线追踪起始位置到目标位置射线追踪总共使用的步长数目来记录旅行时。

[0161] 本发明实施例从声波近似方程出发，推导出三维TTI介质下的程函方程，并将所述三维TTI介质下的程函方程解析为表征沿TTI介质对称轴方向传播的纵波地震波速度、Thomsen各向异性参数、TTI介质对称轴倾角及TTI介质对称轴方位角的形式，进而基于解析后的程函方程建立三维TTI介质下的射线追踪方程并根据射线追踪方程进行射线追踪，该方法能够获得精确的旅行时，利用射线追踪的结果对地震资料进行处理，能够有效解决目前地震资料处理方法精度低的问题。

[0162] 在一个实际的示例中，将本发明实施例提供的三维TTI介质射线追踪方法应用于TTI高斯束偏移成像，具体为通过加窗局部倾斜叠加并计算由震源和接收点出射的各向异性高斯束。基于三维TTI介质射线追踪方法获取射线参数、射线路径和射线旅行时，并根据射线旅行时获得高斯束的走时。根据射线参数、射线路径和高斯束旅行时计算动力学参量Q(s)、P(s)。根据动力学参量Q(s)、P(s)以及高斯束函数计算各条射线周围设定范围内不同频率的波场值。累加上述不同频率的波场值，并根据累加的结果进行成像。

[0163] 图2a-图2h为根据本发明一实施例的SEG管道模型示意图，通过SEG管道模型可以看到利用本发明实施例的射线追踪方法进行高斯束偏移成像的效果。如图所示，一方面各向异性地层下面的地层分界面被拉平，与理论参数模型所呈现图像一致，而各向同性介质射线追踪应用于高斯束偏移成像，各向异性地层下面的地层分界面未拉平。另一方面，通过成像道集对比，各向同性成像道集道集不平，而TTI成像道集道集拉平。高斯束偏移成像与理论参数模型所呈现图像一致，成像精度高。

[0164] 在另一个实际的示例中，利用本发明实施例提供的三维TTI介质射线追踪方法进行TTI参数建模，并将该TTI参数模型用于偏移成像。

[0165] 具体为，利用初始参数模型进行TTI偏移，提供成像道集与偏移剖面，如果成像道集同相轴拉平，则输出初始模型参数及成像道集和偏移剖面数据，如果成像道集同相轴未被拉平，则在偏移剖面上进行属性参数拾取，在成像道集中进行剩余深度拾取，利用拾取得到的属性参数在初始参数场上进行TTI射线追踪，提供射线旅行时及路径信息，与剩余深度一起作为矩阵建立模块的输入数据构建矩阵，最终通过求解矩阵得到更新后的参数场数据。

[0166] 进一步的，在进行深度拾取时，对于不同的反演参数需要的道集偏移距和角度范围是有区别的。具体的，如果单独反演vt，仅仅需要角度小于30°的道集或者偏移距深度比值小于1范围内的数据。如果单独反演ε，则需要角度小于45°的道集或者偏移距深度比值大1范围内的数据。如果同时反演三个各向异性参数vt,ε,δ，则需要对所有数据进行剩余深度拾取。对于不同的反演参数，矩阵求解过程中需要求解的对象不同。

[0167] 将上述成像道集参数和偏移剖面参数用于各向异性深度偏移成像。

[0168] 图3a-图3f为根据本发明一实施例的盐丘模型示意图，通过盐丘模型可以看到利用本发明实施例的射线追踪方法进行TTI参数建模偏移成像的效果。图3a-图3f依次为初始偏移剖面、参数建模之后的偏移剖面、理论值对应的偏移剖面、初始成像道集、参数建模之后的成像道集、理论值对应的成像道集，如图所示，得到的图像与理论模型精度相差不大，成像精度高。

[0169] 图4为根据本发明另一实施例的三维TTI介质射线追踪系统的结构示意图，如图所示，该系统包括，

[0170] 程函方程建模模块41，其执行步骤S110中的操作，程函方程解析模块42，其执行步骤S120中的操作，射线追踪方程建模模块43，其执行步骤S130中的操作，射线追踪模块44，其执行步骤S140中的操作，此处不再赘述。

[0171] 虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。