基于各向异性全变分约束的弹性波直接包络反演方法技术

本发明专利技术涉及一种基于各向异性全变分约束的弹性波直接包络反演方法，是通过各向异性全变分约束与弹性波直接包络反演获得弹性强散射介质的高精度纵、横波速度结构。首先，对弹性波场进行波场模式分解分别获得纵、横波场，计算包络场分别获得正传纵、横波包络场；然后，根据梯度表达式，可得到弹性波直接包络反演的纵、横波速度梯度，并可更新纵、横波速度模型；之后，对当前纵、横波速度模型施加各向异性全变分约束可得到约束后的纵、横波速度更新模型；最后，将各向异性全变分约束弹性波直接包络反演结果作为初始模型，进行各向异性全变分约束弹性波全波形反演可得到弹性强散射介质的高精度纵、横波速度模型。横波速度模型。横波速度模型。

Elastic wave direct envelope inversion method based on anisotropic total variational constraints

【技术实现步骤摘要】
基于各向异性全变分约束的弹性波直接包络反演方法


[0001]本专利技术涉及一种在弹性波直接包络反演中施加各向异性全变分约束以提升强散射介质纵、横波速度结构反演效果的方法。

技术介绍

[0002]强散射介质通常能够形成油气资源的良好封存介质，因此，强散射介质的高精度参数建模方法已经在勘探地球物理领域引起高度重视。全波形反演方法是目前勘探地震领域参数建模精度最高的方法，其能充分利用地震波场的运动学和动力学信息。近年来，人们希望利用全波形反演方法对实际勘探中遇到的强散射介质(如盐丘等)进行高精度参数建模，以提高强散射体及其屏蔽区的成像质量。但是，常规全波形反演方法基于波恩近似这一弱散射近似，且实际地震数据往往缺失有效低频信息，这使得直接利用常规全波形反演方法进行强散射介质的参数建模存在一定难度。近年来，国内外学者从不同角度对全波形反演的方法理论框架进行了改进，以适应强散射介质的参数反演，这些研究绝大多数均是基于声波介质，主要包括Laplace
‑
Fourier域波形反演方法、标准设置方法、全变分约束方法、深度学习方法、以及直接包络反演方法等。关于弹性强散射介质的多参数建模研究仍处于起步探索阶段，现有方法对于弹性强散射体内部以及下方速度建模效果均不十分理想。总之，当前缺乏一种能够在低频数据缺失且无先验信息情况下对强散射介质进行高精度纵、横波速度建模的方法。
[0003]直接包络反演方法能够在地震数据缺失低频且无先验信息情况下，对强散射介质进行大尺度参数建模，且其具有相对较高的计算效率。直接包络反演方法与全波形反演方法的最大区别是其定义了直接包络敏感核函数(区别于常规包络反演的波形敏感核函数)，能够将低频包络数据扰动直接映射为强散射介质的大尺度参数扰动。但是，当前的弹性波直接包络反演方法在对弹性强散射介质进行纵、横波速度建模时，难以获得理想的强散射体内部速度信息，从而影响强散射体下边界及其下方的速度建模效果。

技术实现思路

[0004]本专利技术的目的就是针对上述现有技术的不足，提供一种新的强散射介质弹性波多参数反演方法，进而解决油气资源地震勘探中强散射储层纵、横波速度高精度建模问题，为复杂强散射储层成像提供高精度的纵、横波速度模型。
[0005]本专利技术的思想是：充分利用各向异性全变分约束能够锐化模型边界和增强层内均匀性的优势，将其引入弹性波直接包络反演过程中。在弹性波直接包络反演每次迭代中，对迭代更新后的纵、横波速度模型分别施加各向异性全变分约束，使得强散射体内部趋于均匀且边界更为明显。最终结合各向异性全变分约束弹性波全波形反演技术获得强散射介质的高精度纵、横波速度模型，克服现有技术的不足。
[0006]本专利技术的目的是通过以下技术方案实现的：
[0007]首先，准备好预处理后的弹性波多分量观测地震数据；利用观测地震数据估计震
源子波，给出纵、横波速度的初始模型(无需先验信息)；在初始速度模型上进行正演模拟，得到模拟地震数据与正传地震波场；对正传地震波场进行波场模式分解，分别得到纵、横波正传波场，求取相应的包络场，分别得到纵、横波正传包络场；对观测包络数据和模拟包络数据作差得到包络伴随源；在初始纵、横波速度模型上进行伴随包络场计算与模式分解，得到纵、横波伴随包络场；利用正传包络场与伴随包络场分别计算纵、横波速度梯度；求取步长，并进行纵、横波速度模型的更新；对更新后的纵、横波速度模型分别施加各向异性全变分约束，获得约束后的纵、横波速度模型；进行迭代更新，直到满足停止条件，得到强散射介质的大尺度纵、横波速度结构；以强散射介质的大尺度纵、横波速度模型为初始模型，进行各向异性全变分约束弹性波全波形反演，得到强散射介质的高精度纵、横波速度建模结果。
[0008]本专利技术所述的基于各向异性全变分约束的弹性波直接包络反演方法是通过MATLAB平台实现的；
[0009]本专利技术所述的基于各向异性全变分约束的弹性波直接包络反演方法，包括以下步骤：
[0010]a、安装MATLAB软件平台；
[0011]b、对数据进行静校正、去噪预处理，得到高质量的弹性波多分量观测地震数据；
[0012]c、对地震数据进行子波估计，提取每一炮数据的震源子波；
[0013]d、通过背景速度分析，得到大致的背景纵、横波速度范围，生成背景速度模型，其中不含有强散射体的任何先验信息，分别作为反演的初始纵波速度模型v
p0
和初始横波速度模型v
s0
；
[0014]e、在初始纵、横波速度模型上计算弹性波多分量模拟地震数据，对模拟数据取包络得到模拟包络数据上角标i表示i方向的分量，对于二维情况，指x(水平)和z(垂直)方向；对观测地震数据取包络得到观测包络数据通过公式(1)计算弹性波直接包络反演的目标函数σ
EDEI
：
[0015][0016]式中，求和符号下角标sr表示对所有震源和检波点进行积分，t表示时间，T表示总的记录时间长度；
[0017]f、在初始模型上计算模拟地震波场，对模拟地震波场进行波场模式分解得到正传纵波波场和正传横波波场，分别取包络得到纵波正传包络场和横波正传包络场
[0018]g、计算模拟包络数据与观测包络数据的差，得到伴随源，伴随源反传得到伴随包络场；对伴随包络场进行波场模式分解得到纵波伴随包络场和横波伴随包络场
[0019]h、纵波正传包络场与纵波伴随包络场进行零延迟互相关得到弹性波直接包络反演的纵波速度梯度，如公式(2)所示：
[0020][0021]式中，v
p
表示纵波速度，ρ表示密度。弹性波直接包络反演的横波速度梯度可由公式(3)计算：
[0022][0023]式中，v
s
表示横波速度，
·
表示点乘，μ表示剪切模量；
[0024]i、选择合适的步长，采用最速下降法对纵、横波速度模型进行更新。假设当前迭代次数为m，当前迭代更新得到的纵波速度和横波速度模型分别用和表示，对二者施加各向异性全变分约束的过程等价于求解公式(4)和(5)所示的最优化问题：
[0025][0026][0027]其中，J1和J2分别表示对纵、横波速度模型施加各向异性全变分约束的目标函数，α1和α2分别表示纵、横波速度的更新步长，和分别为第m次迭代的纵、横波速度的更新量，λ1和λ2分别为对纵、横波速度施加各向异性全变分约束的权系数，||
·
||表示二范数，||
·
||
ATV
表示各向异性全变分范数。对速度模型v计算各向异性全变分范数的具体表达式如公式(6)所示：
[0028][0029]其中，nz和nx分别表示模型垂向和横向的网格点数。求解公式(4)和(5)所示优化问题，得到各向异性全变分约束后的纵、横波速度模型作为当前迭代更新后的速度模型；
[0030]j、在更新后的模型上，进行迭代停止条件判断；如果不满足停止条件，将更新后的纵、横波速度模型作为初始模型本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于各向异性全变分约束的弹性波直接包络反演方法，其特征在于，是在地震数据缺失低频信息且没有模型先验信息的情况下，通过对弹性波直接包络反演过程施加各向异性全变分约束以获得弹性强散射介质的高质量大尺度纵、横波速度结构；将获得的高质量大尺度纵、横波速度结构作为各向异性全变分约束弹性波全波形反演的初始模型，获得弹性强散射介质的高精度纵、横波速度结构；包括以下步骤：a、安装MATLAB软件平台；b、对数据进行静校正、去噪预处理，得到高质量的弹性波多分量观测地震数据；c、对地震数据进行子波估计，提取每一炮数据的震源子波；d、通过背景速度分析，得到大致的背景纵、横波速度范围，生成背景速度模型，其中不含有强散射体的任何先验信息，分别作为反演的初始纵波速度模型v
p0
和初始横波速度模型v
s0
；e、在初始纵、横波速度模型上计算弹性波多分量模拟地震数据，对模拟数据取包络得到模拟包络数据上角标i表示i方向的分量，对于二维情况，指x(水平)和z(垂直)方向；对观测地震数据取包络得到观测包络数据通过公式(1)计算弹性波直接包络反演的目标函数σ
EDEI
：式中，求和符号下角标sr表示对所有震源和检波点进行积分，t表示时间，T表示总的记录时间长度；f、在初始模型上计算模拟地震波场，对模拟地震波场进行波场模式分解得到正传纵波波场和正传横波波场，分别取包络得到纵波正传包络场和横波正传包络场g、计算模拟包络数据与观测包络数据的差，得到伴随源，伴随源反传得到伴随包络场；对伴随包络场进行波场模式分解得到纵波伴随包络场和横波伴随包络场h、纵波正传包络场与纵波伴随包络场进行零延迟互相关得到弹性波直接包络反演的纵波速度梯度，如公式(2)所示：式中，v

【专利技术属性】
技术研发人员：张盼，韩立国，周奕秀，
申请(专利权)人：吉林大学，
类型：发明
国别省市：