[发明专利]一种粘声各向异性最小二乘逆时偏移成像方法

说明书

技术领域

本发明属于石油勘探领域，具体涉及一种粘声各向异性最小二乘逆时偏移成像方法。

背景技术

地球介质广泛存在着粘滞性和各向异性，而且这两种性质常常是同时存在的，介质的粘滞性影响地震波传播的能量，而介质的各向异性影响地震波传播的形态。因此，在进行地震波偏移成像方法研究时，需要同时考虑地下介质的粘滞性和各向异性，消除两者对地震波传播的影响，传统基于理想的弹性介质及其声学假设都是不成立的，需要对同时存在粘滞性和各向异性介质的成像方法进行研究。

传统的逆时偏移成像方法存在低频成像噪音、成像能量不均衡、假象严重及分辨率低等缺点，最小二乘逆时偏移方法是一种基于线性化反演理论的保幅成像方法，其基本思想是在宏观背景模型的基础上估计出一个最优化的扰动部分。很多模型及实际资料试验表明，它可以提高偏移剖面的分辨率，压制偏移假象，更重要的是，它可以给出真实地下反射界面的相对振幅。

发明内容

针对现有技术中存在的上述技术问题，本发明提出了一种粘声各向异性最小二乘逆时偏移成像方法，本发明能够对同时存在粘滞性和各向异性的复杂地下构造进行高精度成像，既能克服传统成像方法在同时处理同时存在粘滞性和各向异性介质的不足，又通过最小二乘方式消除了成像噪音，得到真振幅的成像剖面。

为了实现上述目的，本发明采用如下技术方案：

所述的粘声各向异性最小二乘逆时偏移成像方法，具体包括如下步骤：

步骤1：输入速度模型、各向异性Thomsen参数模型、品质因子Q模型及慢度扰动模型，并建立观测系统；

步骤2：采用有限差分方法进行粘声各向异性介质正演模拟：

正演模拟采用如下的波动方程：

$\left\{\begin{array}{c}\frac{{\∂}^{2}p}{\∂t^2}+\frac{{\mathrm{\τv}}_{px}}{2}\sqrt{-\frac{{\∂}^2}{\∂x^2}}\frac{\∂p}{\∂t}-v_{px}^2\frac{{\∂}^{2}p}{\∂x^2}-v_{pn}{v}_{pz}\frac{{\∂}^{2}q}{\∂z^2}=0\\ \frac{{\∂}^{2}q}{\∂t^2}+\frac{{\mathrm{\τv}}_{pz}}{2}\sqrt{-\frac{{\∂}^2}{\∂z^2}}\frac{\∂q}{\∂t}-v_{pn}{v}_{pz}\frac{{\∂}^{2}p}{\∂x^2}-v_{pz}^2\frac{{\∂}^{2}q}{\∂z^2}=0\end{array}\right.---\left(1\right)$