[发明专利]二维的横向各向异性介质HTI叠前深度偏移建模方法及装置

说明书

技术领域

本发明涉及地震勘探技术领域，具体地涉及石油地球物理勘探领域的一种新的叠前深度偏移技术，应用于HTI介质的偏移成像，特别涉及一种二维的横向各向异性介质HTI叠前深度偏移建模方法及装置。

背景技术

长期以来人们对于地下介质中的各向异性进行了大量的研究，但是直到Thomsen(1986)给出著名的VTI(transversely isotropic media with a vertical symmetry axis，带有垂直对称轴的横向各向同性介质)介质各向异性参数ε和δ之后才得以应用于地震勘探。在弱各性异性条件下，Thomsen(1986)推导了VTI介质中qP(准纵波)相速度表达式的近似公式：

V_P(θ)≈α₀(1+δsin²θcos²θ+εsin⁴θ)，(1)

这里α₀是qP波的垂直相速度；θ是相速度角；δ和ε是Thomsen引入的无量纲参数，它们是弹性刚度系数的线性组合。

通过使用扰动理论求解Christoffel方程，Mensch和Rasolofosaon(1997)将Thomsen的理论推广到任意对称轴的弱各向异性介质，同时引入了广义各向异性参数，作为Thomsen的VTI各向异性参数δ和ε的推广，前提仍然是在弱各向异性的条件下。Rasolofosaon(2003)证明对于任意各向异性介质，在弱各向异性条件下，给定任意方位λ上的二维剖面上的各向异性都可以用等效VTI介质来近似，也就是存在随方位变化的Thomsen参数δ(λ)和ε(λ)。Rasolofosaon(2000)对其给出的qP波相速度解析表达式进行了修改，具体如下：

V_P(θ，λ)＝α₀(1+δ(λ)sin²θcos²θ+ε(λ)sin⁴θ)+ΔE_ttriclinic(θ，λ)，(2)

其中：

δ(λ)＝δ_xcos²λ+2x_zcosλsinλ+δ_ysin²λ，(3)

ε(λ)＝δ_xcos⁴λ+δ_zcos²λsin²λ+2cosλsinλ(ε₁₆cos²λ+ε₂₆sin²λ)+ε_ysin⁴λ，(4)

ΔE_triclinic(θ，λ)＝2cosθsin³θ(ε₁₅cos³λ+ε₂₄sin³λ

+x_xcos²θsinλ+x_ysin²λcosλ)(5)

+2sinθcos³θ(ε₃₄sinλ+ε₃₅cosλ).

这里λ表示方位角，θ表示相位角，广义弱各向异性参数ε_x，ε_y，δ_x，δ_y和δ_z等是弹性刚度系数的线性组合，它们仍然是无量纲参数，具体的表达式为：