一种三维地震各向异性介质逆时偏移成像方法及装置

摘要

一种三维地震各向异性介质逆时偏移成像方法及装置，应用于反射波地震数据处理技术领域，该方法包括：确定需要进行逆时偏移成像的所有炮(S11)；针对每一炮，执行如下处理步骤(S12)：在该炮对应的炮点位置放置子波，并应用耦合二阶偏微分方程对该炮进行炮点波场模拟(S121)；通过对该炮对应的炮数据应用所述耦合二阶偏微分方程，对该炮进行检波点波场模拟(S122)；应用互相关成像条件对所述炮点波场模拟的结果及所述检波点波场模拟的结果进行成像，得到该炮的单炮逆时偏移结果(S123)；对所有炮都执行上述处理步骤后，将所有炮的单炮逆时偏移结果叠加起来形成逆时偏移成像剖面(S13)。本发明采用了稳定的耦合二阶偏微分方程实现波动方程的差分求解，能够解决TTI介质对称轴倾角突变引起的计算不稳定问题，并最终解决速度急剧变化的三维复杂构造成像问题。

主权项

一种三维地震各向异性介质逆时偏移成像方法，其特征在于，包括：确定需要进行逆时偏移成像的所有炮；针对每一炮，执行如下处理步骤：在该炮对应的炮点位置放置子波，并应用耦合二阶偏微分方程对该炮进行炮点波场模拟；通过对该炮对应的炮数据应用所述耦合二阶偏微分方程，对该炮进行检波点波场模拟；应用互相关成像条件对所述炮点波场模拟的结果及所述检波点波场模拟的结果进行成像，得到该炮的单炮逆时偏移结果；对所有炮都执行上述处理步骤后，将所有炮的单炮逆时偏移结果叠加起来形成逆时偏移成像剖面；其中，所述耦合二阶偏微分方程为：$\left\{\begin{array}{c}\frac{1}{v_{\mathrm{pz}}^2}\frac{{\partial }^{2}p}{{\partial t}^2}=(1+2ϵ)H_{2}p+H_{1}q+\frac{ϵ-\delta }{\sigma }{H}_1(p-q)\\ \frac{1}{v_{\mathrm{pz}}^2}\frac{{\partial }^{2}q}{\partial t^2}=(1+2\delta )H_{2}p+H_{1}q-\frac{ϵ-\delta }{\sigma }{H}_2(p-q)\\ H_1={\sin }^2\theta {\cos }^2\phi \frac{{\partial }^2}{\partial x^2}+{\sin }^2\theta {\sin }^2\phi \frac{{\partial }^2}{\partial y^2}+{\cos }^2\theta \frac{{\partial }^2}{\partial z^2}+\\ {\sin }^2\theta \sin 2\phi \frac{{\partial }^2}{\partial x\partial y}+\sin 2\theta \sin \phi \frac{{\partial }^2}{\partial y\partial z}+\sin 2\theta \cos \phi \frac{{\partial }^2}{\partial x\partial z}\\ H_2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial z^2}-H_1\end{array}\right.$所述方程中，δ、ε分别为成像空间对应的Thomson各向异性参数；θ、φ分别为成像空间对应的对称轴倾角参数和对称轴方位角参数。